SEOBNRv2_opt_dp1.h

Go to the documentation of this file.

REAL8 tmp53=sigmaKerr->data[2]*sigmaKerr->data[2];
REAL8 tmp54=sqrt(tmp53);
REAL8 tmp83=coeffs->KK*eta;
REAL8 tmp145=-1.+tmp83;
REAL8 tmp56=1/x->data[0];
REAL8 tmp60=x->data[0]*x->data[0];
REAL8 tmp153=(1.0/(x->data[0]*x->data[0]));
REAL8 tmp5=1.*e3z;
REAL8 tmp48=0.+tmp5;
REAL8 tmp81=tmp53+tmp60;
REAL8 tmp82=tmp81*tmp81;
REAL8 tmp149=(1.0/(tmp145*tmp145));
REAL8 tmp152=1.*tmp149;
REAL8 tmp156=1.*tmp153*tmp53;
REAL8 tmp157=1/tmp145;
REAL8 tmp158=2.*tmp157*tmp56;
REAL8 tmp159=tmp152+tmp156+tmp158;
REAL8 tmp160=coeffs->k0*eta;
REAL8 tmp161=1./(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp164=1.*coeffs->k4*tmp161;
REAL8 tmp165=1./(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp166=1.*coeffs->k3*tmp165;
REAL8 tmp167=1.*coeffs->k2*tmp153;
REAL8 tmp168=1.*coeffs->k1*tmp56;
REAL8 tmp169=pow(x->data[0],-5.);
REAL8 tmp170=1.*tmp56;
REAL8 tmp190=log(tmp170);
REAL8 tmp191=coeffs->k5l*tmp190;
REAL8 tmp192=coeffs->k5+tmp191;
REAL8 tmp193=1.*tmp169*tmp192;
REAL8 tmp196=1.+tmp164+tmp166+tmp167+tmp168+tmp193;
REAL8 tmp198=log(tmp196);
REAL8 tmp201=eta*tmp198;
REAL8 tmp202=1.+tmp160+tmp201;
REAL8 tmp204=-(tmp159*tmp202*tmp53*tmp60);
REAL8 tmp205=tmp204+tmp82;
REAL8 tmp206=1/tmp205;
REAL8 tmp211=0.+tmp60;
REAL8 tmp213=1/tmp211;
REAL8 tmp217=-3.*eta;
REAL8 tmp208=0.+p->data[1];
REAL8 tmp209=tmp208*tmp48;
REAL8 tmp210=0.+tmp209;
REAL8 tmp57=0.;
REAL8 tmp58=2.*tmp54*x->data[0];
REAL8 tmp59=tmp57+tmp58;
REAL8 tmp246=(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp252=eta*eta;
REAL8 tmp230=0.+p->data[0];
REAL8 tmp231=1.*tmp230;
REAL8 tmp232=0.+tmp231;
REAL8 tmp233=tmp232*tmp232;
REAL8 tmp234=26.+tmp217;
REAL8 tmp235=2.*eta*tmp165*tmp234;
REAL8 tmp236=6.*eta*tmp153;
REAL8 tmp237=1.+tmp235+tmp236;
REAL8 tmp238=log(tmp237);
REAL8 tmp239=1.+tmp238;
REAL8 tmp223=1.*tmp48;
REAL8 tmp224=0.+tmp223;
REAL8 tmp225=0.+p->data[2];
REAL8 tmp226=tmp224*tmp225;
REAL8 tmp227=0.+tmp226;
REAL8 tmp228=tmp227*tmp227;
REAL8 tmp229=1.*tmp213*tmp228*tmp60;
REAL8 tmp240=1.*tmp159*tmp202*tmp213*tmp233*tmp239*tmp60;
REAL8 tmp241=tmp210*tmp210;
REAL8 tmp242=1.*tmp206*tmp211*tmp241*tmp60;
REAL8 tmp259=pow(x->data[0],6.);
REAL8 tmp261=(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp262=0.+tmp229+tmp240+tmp242;
REAL8 tmp281=sqrt(tmp211);
REAL8 tmp247=6.*sigmaKerr->data[2]*tmp206*tmp210*tmp211*tmp246*tmp48;
REAL8 tmp248=8.*sigmaStar->data[2]*tmp206*tmp210*tmp211*tmp246*tmp48;
REAL8 tmp249=tmp247+tmp248;
REAL8 tmp250=0.08333333333333333*eta*tmp249*tmp56;
REAL8 tmp251=-109.*eta;
REAL8 tmp253=51.*tmp252;
REAL8 tmp254=tmp251+tmp253;
REAL8 tmp255=8.*tmp206*tmp210*tmp211*tmp246*tmp254*tmp48;
REAL8 tmp256=-6.*eta;
REAL8 tmp257=39.*tmp252;
REAL8 tmp258=tmp256+tmp257;
REAL8 tmp260=-12.*tmp159*tmp202*tmp206*tmp210*tmp233*tmp239*tmp258*tmp259*tmp48;
REAL8 tmp263=-180.*eta*tmp206*tmp210*tmp211*tmp261*tmp262*tmp48;
REAL8 tmp264=tmp255+tmp260+tmp263;
REAL8 tmp265=0.006944444444444444*sigmaKerr->data[2]*tmp153*tmp264;
REAL8 tmp266=103.*eta;
REAL8 tmp267=-60.*tmp252;
REAL8 tmp268=tmp266+tmp267;
REAL8 tmp269=-4.*tmp206*tmp210*tmp211*tmp246*tmp268*tmp48;
REAL8 tmp270=-16.*eta;
REAL8 tmp271=21.*tmp252;
REAL8 tmp272=tmp270+tmp271;
REAL8 tmp273=-12.*tmp159*tmp202*tmp206*tmp210*tmp233*tmp239*tmp259*tmp272*tmp48;
REAL8 tmp274=3.*eta;
REAL8 tmp275=23.+tmp274;
REAL8 tmp276=-4.*eta*tmp206*tmp210*tmp211*tmp261*tmp262*tmp275*tmp48;
REAL8 tmp277=tmp269+tmp273+tmp276;
REAL8 tmp278=0.013888888888888888*sigmaStar->data[2]*tmp153*tmp277;
REAL8 tmp279=tmp250+tmp265+tmp278;
REAL8 tmp304=(tmp232*tmp232*tmp232*tmp232);
REAL8 tmp305=tmp159*tmp159;
REAL8 tmp306=(1.0/(tmp211*tmp211));
REAL8 tmp307=tmp239*tmp239;
REAL8 tmp308=tmp202*tmp202;
REAL8 tmp311=tmp262*tmp262;
REAL8 tmp282=tmp159*tmp202*tmp60;
REAL8 tmp283=sqrt(tmp282);
REAL8 tmp284=-tmp283;
REAL8 tmp285=tmp159*tmp202*tmp206*tmp211*tmp60;
REAL8 tmp286=sqrt(tmp285);
REAL8 tmp287=1.*tmp281*tmp286;
REAL8 tmp288=tmp284+tmp287;
REAL8 tmp289=1.+tmp229+tmp240+tmp242;
REAL8 tmp292=1.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp165;
REAL8 tmp293=-8.*sigmaKerr->data[2];
REAL8 tmp294=14.*sigmaStar->data[2];
REAL8 tmp295=-36.*sigmaKerr->data[2]*tmp159*tmp202*tmp213*tmp233*tmp239*tmp246;
REAL8 tmp296=-30.*sigmaStar->data[2]*tmp159*tmp202*tmp213*tmp233*tmp239*tmp246;
REAL8 tmp297=3.*sigmaKerr->data[2]*tmp262*x->data[0];
REAL8 tmp298=4.*sigmaStar->data[2]*tmp262*x->data[0];
REAL8 tmp299=tmp293+tmp294+tmp295+tmp296+tmp297+tmp298;
REAL8 tmp300=0.08333333333333333*eta*tmp299*tmp56;
REAL8 tmp301=27.*eta;
REAL8 tmp302=-353.+tmp301;
REAL8 tmp303=-2.*eta*tmp302;
REAL8 tmp309=360.*tmp252*tmp259*tmp304*tmp305*tmp306*tmp307*tmp308;
REAL8 tmp310=-2.*tmp262*tmp268*x->data[0];
REAL8 tmp312=-(eta*tmp275*tmp311*tmp60);
REAL8 tmp313=-47.*eta;
REAL8 tmp314=54.*tmp252;
REAL8 tmp315=tmp262*tmp272*x->data[0];
REAL8 tmp316=tmp313+tmp314+tmp315;
REAL8 tmp317=-6.*tmp159*tmp202*tmp213*tmp233*tmp239*tmp246*tmp316;
REAL8 tmp318=tmp303+tmp309+tmp310+tmp312+tmp317;
REAL8 tmp319=0.013888888888888888*sigmaStar->data[2]*tmp153*tmp318;
REAL8 tmp320=8.+tmp274;
REAL8 tmp321=-112.*eta*tmp320;
REAL8 tmp322=810.*tmp252*tmp259*tmp304*tmp305*tmp306*tmp307*tmp308;
REAL8 tmp323=4.*tmp254*tmp262*x->data[0];
REAL8 tmp324=-45.*eta*tmp311*tmp60;
REAL8 tmp325=16.*eta;
REAL8 tmp326=147.*tmp252;
REAL8 tmp327=tmp258*tmp262*x->data[0];
REAL8 tmp328=tmp325+tmp326+tmp327;
REAL8 tmp329=-6.*tmp159*tmp202*tmp213*tmp233*tmp239*tmp246*tmp328;
REAL8 tmp330=tmp321+tmp322+tmp323+tmp324+tmp329;
REAL8 tmp331=0.006944444444444444*sigmaKerr->data[2]*tmp153*tmp330;
REAL8 tmp332=0.+sigmaStar->data[2]+tmp292+tmp300+tmp319+tmp331;
REAL8 tmp290=1./sqrt(tmp289);
REAL8 tmp337=e3z*tmp332;
REAL8 tmp338=0.+tmp337;
REAL8 tmp335=(1.0/(tmp205*tmp205));
REAL8 tmp377=tmp211*tmp211;
REAL8 tmp374=sqrt(tmp289);
REAL8 tmp379=tmp224*tmp332;
REAL8 tmp380=0.+tmp379;
REAL8 tmp342=1./sqrt(tmp282);
REAL8 tmp343=tmp159*tmp202*tmp239*tmp60;
REAL8 tmp344=sqrt(tmp343);
REAL8 tmp345=1./sqrt(tmp285);
REAL8 tmp346=2.*tmp54;
REAL8 tmp347=0.;
REAL8 tmp348=0.+tmp346+tmp347;
REAL8 tmp349=tmp205*tmp348;
REAL8 tmp350=-4.*tmp81*x->data[0];
REAL8 tmp351=1/tmp196;
REAL8 tmp352=2.*coeffs->k2;
REAL8 tmp353=3.*coeffs->k3;
REAL8 tmp354=4.*coeffs->k4;
REAL8 tmp355=5.*tmp192*tmp56;
REAL8 tmp356=tmp354+tmp355;
REAL8 tmp357=1.*tmp356*tmp56;
REAL8 tmp358=tmp353+tmp357;
REAL8 tmp359=1.*tmp358*tmp56;
REAL8 tmp360=tmp352+tmp359;
REAL8 tmp361=1.*tmp360*tmp56;
REAL8 tmp362=coeffs->k1+tmp361;
REAL8 tmp363=-(eta*tmp159*tmp351*tmp362);
REAL8 tmp364=1.*tmp157;
REAL8 tmp365=1.*tmp53*tmp56;
REAL8 tmp366=tmp364+tmp365;
REAL8 tmp367=-2.*tmp202*tmp366;
REAL8 tmp368=2.*tmp159*tmp202*x->data[0];
REAL8 tmp369=tmp363+tmp367+tmp368;
REAL8 tmp370=1.*tmp369*tmp53;
REAL8 tmp371=tmp350+tmp370;
REAL8 tmp372=tmp371*tmp59;
REAL8 tmp373=tmp349+tmp372;
REAL8 tmp375=1.+tmp374;
REAL8 tmp383=1.+tmp229+tmp240+tmp242+tmp374;
REAL8 tmp393=1./sqrt(tmp211);
REAL8 tmp394=1./(tmp205*tmp205*tmp205);
REAL8 tmp336=(1.0/sqrt(tmp289*tmp289*tmp289));
REAL8 tmp376=1/tmp375;
REAL8 tmp397=1.*tmp159*tmp202*tmp206*tmp241*tmp261*tmp377*tmp380;
REAL8 tmp398=-(tmp159*tmp202*tmp233*tmp239*tmp380*tmp60);
REAL8 tmp399=tmp211*tmp380*tmp383;
REAL8 tmp400=0.+tmp398+tmp399;
REAL8 tmp401=1.*tmp159*tmp202*tmp400*tmp60;
REAL8 tmp402=0.+tmp397+tmp401;
REAL8 tmp212=1/tmp159;
REAL8 tmp214=1/tmp202;
REAL8 tmp334=sqrt(tmp211*tmp211*tmp211);
REAL8 tmp406=tmp239*tmp261*tmp305*tmp308;
REAL8 tmp407=1./sqrt(tmp406);
REAL8 tmp408=-2.*tmp159*tmp202*tmp60;
REAL8 tmp409=tmp344*tmp369;
REAL8 tmp410=tmp408+tmp409;
REAL8 tmp414=1.*tmp213*x->data[0];
REAL8 tmp415=-4.*tmp159*tmp202*tmp246;
REAL8 tmp416=tmp369*tmp81;
REAL8 tmp417=tmp415+tmp416;
REAL8 tmp418=0.5*tmp153*tmp206*tmp212*tmp214*tmp417*tmp81;
REAL8 tmp419=tmp414+tmp418;
REAL8 tmp420=2.*tmp374;
REAL8 tmp421=1.+tmp420;
REAL8 tmp385=2.*tmp206*tmp210*tmp211*tmp48*tmp60;
REAL8 tmp386=1.*tmp206*tmp210*tmp211*tmp290*tmp48*tmp60;
REAL8 tmp387=tmp385+tmp386;
REAL8 tmp215=tmp153*tmp205*tmp212*tmp213*tmp214;
REAL8 tmp216=1./sqrt(tmp215);
REAL8 tmp218=4.+tmp217;
REAL8 tmp219=1.*p->data[0];
REAL8 tmp220=0.+tmp219;
REAL8 tmp221=(tmp220*tmp220*tmp220*tmp220);
REAL8 tmp222=2.*eta*tmp153*tmp218*tmp221;
REAL8 tmp243=1.+tmp222+tmp229+tmp240+tmp242;
REAL8 tmp341=(1.0/sqrt(tmp211*tmp211*tmp211));
REAL8 tmp405=1/tmp383;
REAL8 tmp430=-0.5*tmp210*tmp281*tmp286*tmp375*tmp380*tmp407*tmp410*x->data[0];
REAL8 tmp431=1.*tmp210*tmp281*tmp286*tmp380*tmp419*tmp421*x->data[0];
REAL8 tmp432=0.+tmp431;
REAL8 tmp433=tmp283*tmp432;
REAL8 tmp434=tmp430+tmp433;
REAL8 tmp435=tmp344*tmp434;
REAL8 tmp436=0.+tmp435;
REAL8 dHdp1=(1.*eta*(-(tmp165*tmp279*tmp332)-tmp241*tmp246*tmp288*tmp334*tmp335*tmp336*tmp338*tmp48+1.*e3z*tmp206*tmp279*tmp59-0.5*tmp210*tmp336*tmp342*tmp344*tmp345*tmp373*tmp376*tmp393*tmp394*tmp402*tmp48*tmp60+0.5*tmp290*tmp335*tmp341*tmp342*tmp344*tmp345*tmp373*tmp376*(1.*tmp159*tmp202*tmp206*tmp224*tmp241*tmp261*tmp279*tmp377+2.*tmp159*tmp202*tmp206*tmp210*tmp261*tmp377*tmp380*tmp48+1.*tmp159*tmp202*tmp60*(tmp211*tmp224*tmp279*tmp383+tmp211*tmp380*tmp387-tmp159*tmp202*tmp224*tmp233*tmp239*tmp279*tmp60))+1.*e3z*tmp206*tmp210*tmp279*tmp281*tmp288*tmp290*x->data[0]+1.*tmp206*tmp281*tmp288*tmp290*tmp338*tmp48*x->data[0]+1.*tmp206*tmp48*tmp59*x->data[0]+1.*tmp153*tmp212*tmp213*tmp214*tmp286*tmp344*tmp405*(-0.5*tmp206*tmp241*tmp246*tmp286*tmp290*tmp334*tmp380*tmp407*tmp410*tmp48-0.5*tmp210*tmp224*tmp279*tmp281*tmp286*tmp375*tmp407*tmp410*x->data[0]-0.5*tmp281*tmp286*tmp375*tmp380*tmp407*tmp410*tmp48*x->data[0]+tmp283*(2.*tmp206*tmp241*tmp246*tmp286*tmp290*tmp334*tmp380*tmp419*tmp48+1.*tmp210*tmp224*tmp279*tmp281*tmp286*tmp419*tmp421*x->data[0]+1.*tmp281*tmp286*tmp380*tmp419*tmp421*tmp48*x->data[0]))-(0.5*tmp210*tmp342*tmp344*tmp345*tmp373*tmp393*tmp394*tmp402*tmp48*tmp60*(1.0/(tmp375*tmp375)))/tmp289-tmp153*tmp212*tmp213*tmp214*tmp286*tmp387*tmp436*(1.0/(tmp383*tmp383))+(1.*tmp206*tmp210*tmp211*tmp216*tmp48*tmp60)/sqrt(tmp243)))/sqrt(1.+2.*eta*(-1.+0.5*tmp290*tmp335*tmp341*tmp342*tmp344*tmp345*tmp373*tmp376*tmp402+1.*tmp153*tmp212*tmp213*tmp214*tmp286*tmp405*tmp436+1.*tmp206*tmp338*tmp59+1.*tmp206*tmp210*tmp281*tmp288*tmp290*tmp338*x->data[0]+1.*tmp206*tmp210*tmp59*x->data[0]+1.*coeffs->dheffSSv2*eta*tmp161*(s1Vec->data[2]*s1Vec->data[2]+s2Vec->data[2]*s2Vec->data[2])-0.5*tmp165*(0.+tmp332*tmp332)+tmp216*sqrt(tmp243)));