SEOBNRv2_opt_dp2tortoise.h

Go to the documentation of this file.

REAL8 tmp1=1.*e3z;
REAL8 tmp3=0.+tmp1;
REAL8 tmp4=1.*tmp3;
REAL8 tmp5=0.+tmp4;
REAL8 tmp9=x->data[0]*x->data[0];
REAL8 tmp12=coeffs->KK*eta;
REAL8 tmp13=-1.+tmp12;
REAL8 tmp17=(1.0/(x->data[0]*x->data[0]));
REAL8 tmp10=0.+tmp9;
REAL8 tmp11=1/tmp10;
REAL8 tmp20=1/x->data[0];
REAL8 tmp16=sigmaKerr->data[2]*sigmaKerr->data[2];
REAL8 tmp14=(1.0/(tmp13*tmp13));
REAL8 tmp15=1.*tmp14;
REAL8 tmp18=1.*tmp16*tmp17;
REAL8 tmp19=1/tmp13;
REAL8 tmp21=2.*tmp19*tmp20;
REAL8 tmp22=tmp15+tmp18+tmp21;
REAL8 tmp24=coeffs->k0*eta;
REAL8 tmp25=1./(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp26=1.*coeffs->k4*tmp25;
REAL8 tmp27=1./(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp28=1.*coeffs->k3*tmp27;
REAL8 tmp29=1.*coeffs->k2*tmp17;
REAL8 tmp30=1.*coeffs->k1*tmp20;
REAL8 tmp31=pow(x->data[0],-5.);
REAL8 tmp32=1.*tmp20;
REAL8 tmp33=log(tmp32);
REAL8 tmp34=coeffs->k5l*tmp33;
REAL8 tmp35=coeffs->k5+tmp34;
REAL8 tmp36=1.*tmp31*tmp35;
REAL8 tmp37=1.+tmp26+tmp28+tmp29+tmp30+tmp36;
REAL8 tmp38=log(tmp37);
REAL8 tmp39=eta*tmp38;
REAL8 tmp40=1.+tmp24+tmp39;
REAL8 tmp6=0.+p->data[2];
REAL8 tmp7=tmp5*tmp6;
REAL8 tmp8=0.+tmp7;
REAL8 tmp42=tmp16+tmp9;
REAL8 tmp43=tmp42*tmp42;
REAL8 tmp44=-(tmp16*tmp22*tmp40*tmp9);
REAL8 tmp45=tmp43+tmp44;
REAL8 tmp48=-3.*eta;
REAL8 tmp50=1.*p->data[0];
REAL8 tmp51=0.+tmp50;
REAL8 tmp62=26.+tmp48;
REAL8 tmp63=2.*eta*tmp27*tmp62;
REAL8 tmp64=6.*eta*tmp17;
REAL8 tmp65=1.+tmp63+tmp64;
REAL8 tmp66=log(tmp65);
REAL8 tmp67=1.+tmp66;
REAL8 tmp68=tmp22*tmp22;
REAL8 tmp69=(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp70=tmp40*tmp40;
REAL8 tmp71=tmp67*tmp68*tmp69*tmp70;
REAL8 tmp86=sqrt(tmp16);
REAL8 tmp60=1/tmp45;
REAL8 tmp90=(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp96=eta*eta;
REAL8 tmp72=1./sqrt(tmp71);
REAL8 tmp73=1/tmp42;
REAL8 tmp74=sqrt(tmp71);
REAL8 tmp75=tmp73*tmp74;
REAL8 tmp76=-1.+tmp75;
REAL8 tmp77=-(tmp42*tmp51*tmp72*tmp76);
REAL8 tmp78=p->data[0]+tmp77;
REAL8 tmp79=1.*tmp78;
REAL8 tmp80=0.+tmp79;
REAL8 tmp81=tmp80*tmp80;
REAL8 tmp54=tmp8*tmp8;
REAL8 tmp55=1.*tmp11*tmp54*tmp9;
REAL8 tmp56=0.+p->data[1];
REAL8 tmp57=tmp3*tmp56;
REAL8 tmp58=0.+tmp57;
REAL8 tmp59=tmp58*tmp58;
REAL8 tmp61=1.*tmp10*tmp59*tmp60*tmp9;
REAL8 tmp82=1.*tmp11*tmp22*tmp40*tmp67*tmp81*tmp9;
REAL8 tmp103=pow(x->data[0],6.);
REAL8 tmp104=(1.0/(tmp10*tmp10));
REAL8 tmp106=0.+tmp55+tmp61+tmp82;
REAL8 tmp125=sqrt(tmp10);
REAL8 tmp91=6.*sigmaKerr->data[2]*tmp11*tmp5*tmp8*tmp90;
REAL8 tmp92=8.*sigmaStar->data[2]*tmp11*tmp5*tmp8*tmp90;
REAL8 tmp93=tmp91+tmp92;
REAL8 tmp94=0.08333333333333333*eta*tmp20*tmp93;
REAL8 tmp95=-109.*eta;
REAL8 tmp97=51.*tmp96;
REAL8 tmp98=tmp95+tmp97;
REAL8 tmp99=8.*tmp11*tmp5*tmp8*tmp90*tmp98;
REAL8 tmp100=-6.*eta;
REAL8 tmp101=39.*tmp96;
REAL8 tmp102=tmp100+tmp101;
REAL8 tmp105=-12.*tmp102*tmp103*tmp104*tmp22*tmp40*tmp5*tmp67*tmp8*tmp81;
REAL8 tmp107=-180.*eta*tmp106*tmp11*tmp5*tmp69*tmp8;
REAL8 tmp108=tmp105+tmp107+tmp99;
REAL8 tmp109=0.006944444444444444*sigmaKerr->data[2]*tmp108*tmp17;
REAL8 tmp110=103.*eta;
REAL8 tmp111=-60.*tmp96;
REAL8 tmp112=tmp110+tmp111;
REAL8 tmp113=-4.*tmp11*tmp112*tmp5*tmp8*tmp90;
REAL8 tmp114=-16.*eta;
REAL8 tmp115=21.*tmp96;
REAL8 tmp116=tmp114+tmp115;
REAL8 tmp117=-12.*tmp103*tmp104*tmp116*tmp22*tmp40*tmp5*tmp67*tmp8*tmp81;
REAL8 tmp118=3.*eta;
REAL8 tmp119=23.+tmp118;
REAL8 tmp120=-4.*eta*tmp106*tmp11*tmp119*tmp5*tmp69*tmp8;
REAL8 tmp121=tmp113+tmp117+tmp120;
REAL8 tmp122=0.013888888888888888*sigmaStar->data[2]*tmp121*tmp17;
REAL8 tmp123=tmp109+tmp122+tmp94;
REAL8 tmp148=tmp67*tmp67;
REAL8 tmp149=(tmp80*tmp80*tmp80*tmp80);
REAL8 tmp152=tmp106*tmp106;
REAL8 tmp126=tmp22*tmp40*tmp9;
REAL8 tmp127=sqrt(tmp126);
REAL8 tmp128=-tmp127;
REAL8 tmp129=tmp10*tmp22*tmp40*tmp60*tmp9;
REAL8 tmp130=sqrt(tmp129);
REAL8 tmp131=1.*tmp125*tmp130;
REAL8 tmp132=tmp128+tmp131;
REAL8 tmp133=1.+tmp55+tmp61+tmp82;
REAL8 tmp136=1.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp27;
REAL8 tmp137=-8.*sigmaKerr->data[2];
REAL8 tmp138=14.*sigmaStar->data[2];
REAL8 tmp139=-36.*sigmaKerr->data[2]*tmp11*tmp22*tmp40*tmp67*tmp81*tmp90;
REAL8 tmp140=-30.*sigmaStar->data[2]*tmp11*tmp22*tmp40*tmp67*tmp81*tmp90;
REAL8 tmp141=3.*sigmaKerr->data[2]*tmp106*x->data[0];
REAL8 tmp142=4.*sigmaStar->data[2]*tmp106*x->data[0];
REAL8 tmp143=tmp137+tmp138+tmp139+tmp140+tmp141+tmp142;
REAL8 tmp144=0.08333333333333333*eta*tmp143*tmp20;
REAL8 tmp145=27.*eta;
REAL8 tmp146=-353.+tmp145;
REAL8 tmp147=-2.*eta*tmp146;
REAL8 tmp150=360.*tmp103*tmp104*tmp148*tmp149*tmp68*tmp70*tmp96;
REAL8 tmp151=-2.*tmp106*tmp112*x->data[0];
REAL8 tmp153=-(eta*tmp119*tmp152*tmp9);
REAL8 tmp154=-47.*eta;
REAL8 tmp155=54.*tmp96;
REAL8 tmp156=tmp106*tmp116*x->data[0];
REAL8 tmp157=tmp154+tmp155+tmp156;
REAL8 tmp158=-6.*tmp11*tmp157*tmp22*tmp40*tmp67*tmp81*tmp90;
REAL8 tmp159=tmp147+tmp150+tmp151+tmp153+tmp158;
REAL8 tmp160=0.013888888888888888*sigmaStar->data[2]*tmp159*tmp17;
REAL8 tmp161=8.+tmp118;
REAL8 tmp162=-112.*eta*tmp161;
REAL8 tmp163=810.*tmp103*tmp104*tmp148*tmp149*tmp68*tmp70*tmp96;
REAL8 tmp164=4.*tmp106*tmp98*x->data[0];
REAL8 tmp165=-45.*eta*tmp152*tmp9;
REAL8 tmp166=16.*eta;
REAL8 tmp167=147.*tmp96;
REAL8 tmp168=tmp102*tmp106*x->data[0];
REAL8 tmp169=tmp166+tmp167+tmp168;
REAL8 tmp170=-6.*tmp11*tmp169*tmp22*tmp40*tmp67*tmp81*tmp90;
REAL8 tmp171=tmp162+tmp163+tmp164+tmp165+tmp170;
REAL8 tmp172=0.006944444444444444*sigmaKerr->data[2]*tmp17*tmp171;
REAL8 tmp173=0.+sigmaStar->data[2]+tmp136+tmp144+tmp160+tmp172;
REAL8 tmp23=1/tmp22;
REAL8 tmp41=1/tmp40;
REAL8 tmp180=tmp22*tmp40*tmp67*tmp9;
REAL8 tmp181=sqrt(tmp180);
REAL8 tmp182=sqrt(tmp133);
REAL8 tmp175=1./sqrt(tmp10);
REAL8 tmp185=-2.*tmp22*tmp40*tmp9;
REAL8 tmp186=1/tmp37;
REAL8 tmp187=2.*coeffs->k2;
REAL8 tmp188=3.*coeffs->k3;
REAL8 tmp189=4.*coeffs->k4;
REAL8 tmp190=5.*tmp20*tmp35;
REAL8 tmp191=tmp189+tmp190;
REAL8 tmp192=1.*tmp191*tmp20;
REAL8 tmp193=tmp188+tmp192;
REAL8 tmp194=1.*tmp193*tmp20;
REAL8 tmp195=tmp187+tmp194;
REAL8 tmp196=1.*tmp195*tmp20;
REAL8 tmp197=coeffs->k1+tmp196;
REAL8 tmp198=-(eta*tmp186*tmp197*tmp22);
REAL8 tmp199=1.*tmp19;
REAL8 tmp200=1.*tmp16*tmp20;
REAL8 tmp201=tmp199+tmp200;
REAL8 tmp202=-2.*tmp201*tmp40;
REAL8 tmp203=2.*tmp22*tmp40*x->data[0];
REAL8 tmp204=tmp198+tmp202+tmp203;
REAL8 tmp205=tmp181*tmp204;
REAL8 tmp206=tmp185+tmp205;
REAL8 tmp134=1./sqrt(tmp133);
REAL8 tmp212=1.*tmp11*x->data[0];
REAL8 tmp213=-4.*tmp22*tmp40*tmp90;
REAL8 tmp214=tmp204*tmp42;
REAL8 tmp215=tmp213+tmp214;
REAL8 tmp216=0.5*tmp17*tmp215*tmp23*tmp41*tmp42*tmp60;
REAL8 tmp217=tmp212+tmp216;
REAL8 tmp209=tmp173*tmp5;
REAL8 tmp210=0.+tmp209;
REAL8 tmp87=0.;
REAL8 tmp88=2.*tmp86*x->data[0];
REAL8 tmp89=tmp87+tmp88;
REAL8 tmp207=1.+tmp182;
REAL8 tmp183=1.+tmp182+tmp55+tmp61+tmp82;
REAL8 tmp227=1./sqrt(tmp126);
REAL8 tmp228=1./sqrt(tmp129);
REAL8 tmp229=(1.0/(tmp45*tmp45));
REAL8 tmp230=2.*tmp86;
REAL8 tmp231=0.;
REAL8 tmp232=0.+tmp230+tmp231;
REAL8 tmp233=tmp232*tmp45;
REAL8 tmp234=-4.*tmp42*x->data[0];
REAL8 tmp235=1.*tmp16*tmp204;
REAL8 tmp236=tmp234+tmp235;
REAL8 tmp237=tmp236*tmp89;
REAL8 tmp238=tmp233+tmp237;
REAL8 tmp240=tmp10*tmp10;
REAL8 tmp252=pow(tmp10,-2.5);
REAL8 tmp176=(1.0/sqrt(tmp133*tmp133*tmp133));
REAL8 tmp239=1/tmp207;
REAL8 tmp255=1.*tmp210*tmp22*tmp240*tmp40*tmp59*tmp60*tmp69;
REAL8 tmp256=-(tmp210*tmp22*tmp40*tmp67*tmp81*tmp9);
REAL8 tmp257=tmp10*tmp183*tmp210;
REAL8 tmp258=0.+tmp256+tmp257;
REAL8 tmp259=1.*tmp22*tmp258*tmp40*tmp9;
REAL8 tmp260=0.+tmp255+tmp259;
REAL8 tmp244=2.*tmp11*tmp5*tmp8*tmp9;
REAL8 tmp245=1.*tmp11*tmp134*tmp5*tmp8*tmp9;
REAL8 tmp246=tmp244+tmp245;
REAL8 tmp218=2.*tmp182;
REAL8 tmp219=1.+tmp218;
REAL8 tmp46=tmp11*tmp17*tmp23*tmp41*tmp45;
REAL8 tmp47=1./sqrt(tmp46);
REAL8 tmp49=4.+tmp48;
REAL8 tmp52=(tmp51*tmp51*tmp51*tmp51);
REAL8 tmp53=2.*eta*tmp17*tmp49*tmp52;
REAL8 tmp83=1.+tmp53+tmp55+tmp61+tmp82;
REAL8 tmp177=e3z*tmp173;
REAL8 tmp178=0.+tmp177;
REAL8 tmp226=(1.0/sqrt(tmp10*tmp10*tmp10));
REAL8 tmp184=1/tmp183;
REAL8 tmp264=-0.5*tmp125*tmp130*tmp206*tmp207*tmp210*tmp58*tmp72*x->data[0];
REAL8 tmp265=1.*tmp125*tmp130*tmp210*tmp217*tmp219*tmp58*x->data[0];
REAL8 tmp266=0.+tmp265;
REAL8 tmp267=tmp127*tmp266;
REAL8 tmp268=tmp264+tmp267;
REAL8 tmp269=tmp181*tmp268;
REAL8 tmp270=0.+tmp269;
REAL8 dHdp2=(1.*eta*(-(tmp123*tmp173*tmp27)+1.*e3z*tmp123*tmp60*tmp89-0.5*tmp176*tmp181*tmp227*tmp228*tmp229*tmp238*tmp239*tmp252*tmp260*tmp5*tmp8*tmp9+0.5*tmp134*tmp181*tmp226*tmp227*tmp228*tmp229*tmp238*tmp239*(1.*tmp123*tmp22*tmp240*tmp40*tmp5*tmp59*tmp60*tmp69+1.*tmp22*tmp40*tmp9*(tmp10*tmp210*tmp246+tmp10*tmp123*tmp183*tmp5-tmp123*tmp22*tmp40*tmp5*tmp67*tmp81*tmp9))-tmp132*tmp175*tmp176*tmp178*tmp5*tmp58*tmp60*tmp8*tmp90+1.*e3z*tmp123*tmp125*tmp132*tmp134*tmp58*tmp60*x->data[0]+1.*tmp11*tmp130*tmp17*tmp181*tmp184*tmp23*tmp41*(-0.5*tmp130*tmp134*tmp175*tmp206*tmp210*tmp5*tmp58*tmp72*tmp8*tmp90-0.5*tmp123*tmp125*tmp130*tmp206*tmp207*tmp5*tmp58*tmp72*x->data[0]+tmp127*(2.*tmp130*tmp134*tmp175*tmp210*tmp217*tmp5*tmp58*tmp8*tmp90+1.*tmp123*tmp125*tmp130*tmp217*tmp219*tmp5*tmp58*x->data[0]))-tmp11*tmp130*tmp17*tmp23*tmp246*tmp270*tmp41*(1.0/(tmp183*tmp183))-(0.5*tmp181*tmp227*tmp228*tmp229*tmp238*tmp252*tmp260*tmp5*tmp8*tmp9*(1.0/(tmp207*tmp207)))/tmp133+(1.*tmp11*tmp47*tmp5*tmp8*tmp9)/sqrt(tmp83)))/sqrt(1.+2.*eta*(-1.+0.5*tmp134*tmp181*tmp226*tmp227*tmp228*tmp229*tmp238*tmp239*tmp260+1.*tmp11*tmp130*tmp17*tmp184*tmp23*tmp270*tmp41+1.*tmp178*tmp60*tmp89+1.*tmp125*tmp132*tmp134*tmp178*tmp58*tmp60*x->data[0]+1.*tmp58*tmp60*tmp89*x->data[0]+1.*coeffs->dheffSSv2*eta*tmp25*(s1Vec->data[2]*s1Vec->data[2]+s2Vec->data[2]*s2Vec->data[2])-0.5*tmp27*(0.+tmp173*tmp173)+tmp47*sqrt(tmp83)));