SEOBNRv2_opt_.h

Go to the documentation of this file.

REAL8 tmp11=sigmaKerr->data[2]*sigmaKerr->data[2];
REAL8 tmp12=sqrt(tmp11);
REAL8 tmp20=coeffs->KK*eta;
REAL8 tmp21=-1.+tmp20;
REAL8 tmp13=1/x->data[0];
REAL8 tmp17=x->data[0]*x->data[0];
REAL8 tmp4=1./(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp24=(1.0/(x->data[0]*x->data[0]));
REAL8 tmp22=(1.0/(tmp21*tmp21));
REAL8 tmp23=1.*tmp22;
REAL8 tmp25=1.*tmp11*tmp24;
REAL8 tmp26=1/tmp21;
REAL8 tmp27=2.*tmp13*tmp26;
REAL8 tmp28=tmp23+tmp25+tmp27;
REAL8 tmp29=coeffs->k0*eta;
REAL8 tmp30=1.*coeffs->k4*tmp4;
REAL8 tmp31=1./(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp32=1.*coeffs->k3*tmp31;
REAL8 tmp33=1.*coeffs->k2*tmp24;
REAL8 tmp34=1.*coeffs->k1*tmp13;
REAL8 tmp35=pow(x->data[0],-5.);
REAL8 tmp36=1.*tmp13;
REAL8 tmp37=log(tmp36);
REAL8 tmp38=coeffs->k5l*tmp37;
REAL8 tmp39=coeffs->k5+tmp38;
REAL8 tmp40=1.*tmp35*tmp39;
REAL8 tmp41=1.+tmp30+tmp32+tmp33+tmp34+tmp40;
REAL8 tmp42=log(tmp41);
REAL8 tmp43=eta*tmp42;
REAL8 tmp44=1.+tmp29+tmp43;
REAL8 tmp18=tmp11+tmp17;
REAL8 tmp19=tmp18*tmp18;
REAL8 tmp45=-(tmp11*tmp17*tmp28*tmp44);
REAL8 tmp46=tmp19+tmp45;
REAL8 tmp6=1.*e3z;
REAL8 tmp7=0.+tmp6;
REAL8 tmp50=0.+tmp17;
REAL8 tmp51=1/tmp50;
REAL8 tmp55=-3.*eta;
REAL8 tmp8=0.+p->data[1];
REAL8 tmp9=tmp7*tmp8;
REAL8 tmp10=0.+tmp9;
REAL8 tmp47=1/tmp46;
REAL8 tmp14=0.;
REAL8 tmp15=2.*tmp12*x->data[0];
REAL8 tmp16=tmp14+tmp15;
REAL8 tmp68=0.+p->data[0];
REAL8 tmp69=1.*tmp68;
REAL8 tmp70=0.+tmp69;
REAL8 tmp71=tmp70*tmp70;
REAL8 tmp72=26.+tmp55;
REAL8 tmp73=2.*eta*tmp31*tmp72;
REAL8 tmp74=6.*eta*tmp24;
REAL8 tmp75=1.+tmp73+tmp74;
REAL8 tmp76=log(tmp75);
REAL8 tmp77=1.+tmp76;
REAL8 tmp87=(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp61=1.*tmp7;
REAL8 tmp62=0.+tmp61;
REAL8 tmp63=0.+p->data[2];
REAL8 tmp64=tmp62*tmp63;
REAL8 tmp65=0.+tmp64;
REAL8 tmp66=tmp65*tmp65;
REAL8 tmp67=1.*tmp17*tmp51*tmp66;
REAL8 tmp78=1.*tmp17*tmp28*tmp44*tmp51*tmp71*tmp77;
REAL8 tmp79=tmp10*tmp10;
REAL8 tmp80=1.*tmp17*tmp47*tmp50*tmp79;
REAL8 tmp90=0.+tmp67+tmp78+tmp80;
REAL8 tmp98=eta*eta;
REAL8 tmp110=3.*eta;
REAL8 tmp99=(tmp70*tmp70*tmp70*tmp70);
REAL8 tmp100=tmp28*tmp28;
REAL8 tmp101=pow(x->data[0],6.);
REAL8 tmp102=(1.0/(tmp50*tmp50));
REAL8 tmp103=tmp77*tmp77;
REAL8 tmp104=tmp44*tmp44;
REAL8 tmp112=tmp90*tmp90;
REAL8 tmp146=sqrt(tmp50);
REAL8 tmp84=1.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp31;
REAL8 tmp85=-8.*sigmaKerr->data[2];
REAL8 tmp86=14.*sigmaStar->data[2];
REAL8 tmp88=-36.*sigmaKerr->data[2]*tmp28*tmp44*tmp51*tmp71*tmp77*tmp87;
REAL8 tmp89=-30.*sigmaStar->data[2]*tmp28*tmp44*tmp51*tmp71*tmp77*tmp87;
REAL8 tmp91=3.*sigmaKerr->data[2]*tmp90*x->data[0];
REAL8 tmp92=4.*sigmaStar->data[2]*tmp90*x->data[0];
REAL8 tmp93=tmp85+tmp86+tmp88+tmp89+tmp91+tmp92;
REAL8 tmp94=0.08333333333333333*eta*tmp13*tmp93;
REAL8 tmp95=27.*eta;
REAL8 tmp96=-353.+tmp95;
REAL8 tmp97=-2.*eta*tmp96;
REAL8 tmp105=360.*tmp100*tmp101*tmp102*tmp103*tmp104*tmp98*tmp99;
REAL8 tmp106=103.*eta;
REAL8 tmp107=-60.*tmp98;
REAL8 tmp108=tmp106+tmp107;
REAL8 tmp109=-2.*tmp108*tmp90*x->data[0];
REAL8 tmp111=23.+tmp110;
REAL8 tmp113=-(eta*tmp111*tmp112*tmp17);
REAL8 tmp114=-47.*eta;
REAL8 tmp115=54.*tmp98;
REAL8 tmp116=-16.*eta;
REAL8 tmp117=21.*tmp98;
REAL8 tmp118=tmp116+tmp117;
REAL8 tmp119=tmp118*tmp90*x->data[0];
REAL8 tmp120=tmp114+tmp115+tmp119;
REAL8 tmp121=-6.*tmp120*tmp28*tmp44*tmp51*tmp71*tmp77*tmp87;
REAL8 tmp122=tmp105+tmp109+tmp113+tmp121+tmp97;
REAL8 tmp123=0.013888888888888888*sigmaStar->data[2]*tmp122*tmp24;
REAL8 tmp124=8.+tmp110;
REAL8 tmp125=-112.*eta*tmp124;
REAL8 tmp126=810.*tmp100*tmp101*tmp102*tmp103*tmp104*tmp98*tmp99;
REAL8 tmp127=-109.*eta;
REAL8 tmp128=51.*tmp98;
REAL8 tmp129=tmp127+tmp128;
REAL8 tmp130=4.*tmp129*tmp90*x->data[0];
REAL8 tmp131=-45.*eta*tmp112*tmp17;
REAL8 tmp132=16.*eta;
REAL8 tmp133=147.*tmp98;
REAL8 tmp134=-6.*eta;
REAL8 tmp135=39.*tmp98;
REAL8 tmp136=tmp134+tmp135;
REAL8 tmp137=tmp136*tmp90*x->data[0];
REAL8 tmp138=tmp132+tmp133+tmp137;
REAL8 tmp139=-6.*tmp138*tmp28*tmp44*tmp51*tmp71*tmp77*tmp87;
REAL8 tmp140=tmp125+tmp126+tmp130+tmp131+tmp139;
REAL8 tmp141=0.006944444444444444*sigmaKerr->data[2]*tmp140*tmp24;
REAL8 tmp142=0.+sigmaStar->data[2]+tmp123+tmp141+tmp84+tmp94;
REAL8 tmp143=e3z*tmp142;
REAL8 tmp144=0.+tmp143;
REAL8 tmp147=tmp17*tmp28*tmp44;
REAL8 tmp150=tmp17*tmp28*tmp44*tmp47*tmp50;
REAL8 tmp154=1.+tmp67+tmp78+tmp80;
REAL8 tmp155=1./sqrt(tmp154);
REAL8 tmp199=tmp142*tmp62;
REAL8 tmp200=0.+tmp199;
REAL8 tmp194=sqrt(tmp154);
REAL8 tmp49=1/tmp28;
REAL8 tmp52=1/tmp44;
REAL8 tmp151=sqrt(tmp150);
REAL8 tmp203=1.+tmp194+tmp67+tmp78+tmp80;
REAL8 tmp162=tmp17*tmp28*tmp44*tmp77;
REAL8 tmp163=sqrt(tmp162);
REAL8 tmp197=(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp171=1/tmp41;
REAL8 tmp172=2.*coeffs->k2;
REAL8 tmp173=3.*coeffs->k3;
REAL8 tmp174=4.*coeffs->k4;
REAL8 tmp175=5.*tmp13*tmp39;
REAL8 tmp176=tmp174+tmp175;
REAL8 tmp177=1.*tmp13*tmp176;
REAL8 tmp178=tmp173+tmp177;
REAL8 tmp179=1.*tmp13*tmp178;
REAL8 tmp180=tmp172+tmp179;
REAL8 tmp181=1.*tmp13*tmp180;
REAL8 tmp182=coeffs->k1+tmp181;
REAL8 tmp183=-(eta*tmp171*tmp182*tmp28);
REAL8 tmp184=1.*tmp26;
REAL8 tmp185=1.*tmp11*tmp13;
REAL8 tmp186=tmp184+tmp185;
REAL8 tmp187=-2.*tmp186*tmp44;
REAL8 tmp188=2.*tmp28*tmp44*x->data[0];
REAL8 tmp189=tmp183+tmp187+tmp188;
REAL8 tmp195=1.+tmp194;
REAL8 tmp148=sqrt(tmp147);
REAL8 Hreal=sqrt(1.+2.*eta*(-1.+1.*tmp144*tmp16*tmp47+1.*tmp10*tmp144*tmp146*(-tmp148+1.*tmp146*tmp151)*tmp155*tmp47*x->data[0]+1.*tmp10*tmp16*tmp47*x->data[0]+1.*coeffs->dheffSSv2*eta*tmp4*(s1Vec->data[2]*s1Vec->data[2]+s2Vec->data[2]*s2Vec->data[2])-0.5*tmp31*(0.+tmp142*tmp142)+(0.5*tmp155*tmp163*((0.+2.*tmp12)*tmp46+tmp16*(1.*tmp11*tmp189-4.*tmp18*x->data[0]))*(0.+1.*tmp17*tmp28*tmp44*(0.+tmp200*tmp203*tmp50-tmp17*tmp200*tmp28*tmp44*tmp71*tmp77)+1.*tmp197*tmp200*tmp28*tmp44*tmp47*tmp79*(tmp50*tmp50))*(1.0/(tmp46*tmp46))*(1.0/sqrt(tmp50*tmp50*tmp50)))/(tmp195*sqrt(tmp147)*sqrt(tmp150))+(1.*tmp151*tmp24*tmp49*tmp51*tmp52*(0.+tmp163*(tmp148*(0.+1.*tmp10*tmp146*tmp151*(1.+2.*tmp194)*tmp200*x->data[0]*(0.5*tmp18*tmp24*tmp47*tmp49*tmp52*(tmp18*tmp189-4.*tmp28*tmp44*tmp87)+1.*tmp51*x->data[0]))-(0.5*tmp10*tmp146*tmp151*tmp195*tmp200*(tmp163*tmp189-2.*tmp17*tmp28*tmp44)*x->data[0])/sqrt(tmp100*tmp104*tmp197*tmp77))))/tmp203+(1.*sqrt(1.+tmp67+tmp78+tmp80+2.*eta*tmp24*(4.+tmp55)*pow(0.+1.*p->data[0],4.)))/sqrt(tmp24*tmp46*tmp49*tmp51*tmp52)));