SEOBNRv2_opt_pdata1.h

Go to the documentation of this file.

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