SEOBNRv2_opt_tortoise.h

Go to the documentation of this file.

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