SEOBNRv2_opt_dx0_sph.h

Go to the documentation of this file.

REAL8 tmp9=1/x->data[0];
REAL8 tmp14=sigmaKerr->data[2]*sigmaKerr->data[2];
REAL8 tmp15=sqrt(tmp14);
REAL8 tmp22=coeffs->KK*eta;
REAL8 tmp23=-1.+tmp22;
REAL8 tmp19=x->data[0]*x->data[0];
REAL8 tmp5=pow(x->data[0],-5.);
REAL8 tmp26=(1.0/(x->data[0]*x->data[0]));
REAL8 tmp31=pow(x->data[0],-6.);
REAL8 tmp34=1./(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp36=1./(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp39=1.*tmp9;
REAL8 tmp40=log(tmp39);
REAL8 tmp41=coeffs->k5l*tmp40;
REAL8 tmp42=coeffs->k5+tmp41;
REAL8 tmp24=(1.0/(tmp23*tmp23));
REAL8 tmp25=1.*tmp24;
REAL8 tmp27=1.*tmp14*tmp26;
REAL8 tmp28=1/tmp23;
REAL8 tmp29=2.*tmp28*tmp9;
REAL8 tmp30=tmp25+tmp27+tmp29;
REAL8 tmp45=1.*coeffs->k4*tmp34;
REAL8 tmp46=1.*coeffs->k3*tmp36;
REAL8 tmp47=1.*coeffs->k2*tmp26;
REAL8 tmp48=1.*coeffs->k1*tmp9;
REAL8 tmp49=1.*tmp42*tmp5;
REAL8 tmp50=1.+tmp45+tmp46+tmp47+tmp48+tmp49;
REAL8 tmp53=coeffs->k0*eta;
REAL8 tmp54=log(tmp50);
REAL8 tmp55=eta*tmp54;
REAL8 tmp56=1.+tmp53+tmp55;
REAL8 tmp20=tmp14+tmp19;
REAL8 tmp7=1.*e3z;
REAL8 tmp8=0.+tmp7;
REAL8 tmp10=p->data[2]*tmp9;
REAL8 tmp11=0.+tmp10;
REAL8 tmp12=tmp11*tmp8;
REAL8 tmp13=0.+tmp12;
REAL8 tmp63=tmp20*tmp20;
REAL8 tmp64=-(tmp14*tmp19*tmp30*tmp56);
REAL8 tmp65=tmp63+tmp64;
REAL8 tmp16=0.;
REAL8 tmp17=2.*tmp15*x->data[0];
REAL8 tmp18=tmp16+tmp17;
REAL8 tmp71=1/tmp65;
REAL8 tmp76=0.+tmp19;
REAL8 tmp87=1.*tmp8;
REAL8 tmp88=0.+tmp87;
REAL8 tmp89=-(p->data[1]*tmp9);
REAL8 tmp90=0.+tmp89;
REAL8 tmp91=tmp88*tmp90;
REAL8 tmp92=0.+tmp91;
REAL8 tmp77=1/tmp76;
REAL8 tmp93=tmp92*tmp92;
REAL8 tmp81=-3.*eta;
REAL8 tmp32=-(coeffs->k5l*tmp31);
REAL8 tmp33=-4.*coeffs->k4*tmp5;
REAL8 tmp35=-3.*coeffs->k3*tmp34;
REAL8 tmp37=-2.*coeffs->k2*tmp36;
REAL8 tmp38=-(coeffs->k1*tmp26);
REAL8 tmp43=-5.*tmp31*tmp42;
REAL8 tmp44=tmp32+tmp33+tmp35+tmp37+tmp38+tmp43;
REAL8 tmp51=1/tmp50;
REAL8 tmp99=0.+p->data[0];
REAL8 tmp100=1.*tmp99;
REAL8 tmp101=0.+tmp100;
REAL8 tmp102=tmp101*tmp101;
REAL8 tmp103=26.+tmp81;
REAL8 tmp104=2.*eta*tmp103*tmp36;
REAL8 tmp105=6.*eta*tmp26;
REAL8 tmp106=1.+tmp104+tmp105;
REAL8 tmp94=(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp95=(1.0/(tmp76*tmp76));
REAL8 tmp107=log(tmp106);
REAL8 tmp108=1.+tmp107;
REAL8 tmp58=-2.*tmp14*tmp36;
REAL8 tmp59=-2.*tmp26*tmp28;
REAL8 tmp60=tmp58+tmp59;
REAL8 tmp21=4.*tmp20*x->data[0];
REAL8 tmp52=-(eta*tmp14*tmp19*tmp30*tmp44*tmp51);
REAL8 tmp57=-2.*tmp14*tmp30*tmp56*x->data[0];
REAL8 tmp61=-(tmp14*tmp19*tmp56*tmp60);
REAL8 tmp62=tmp21+tmp52+tmp57+tmp61;
REAL8 tmp66=(1.0/(tmp65*tmp65));
REAL8 tmp118=tmp13*tmp13;
REAL8 tmp82=4.+tmp81;
REAL8 tmp83=1.*p->data[0];
REAL8 tmp84=0.+tmp83;
REAL8 tmp85=(tmp84*tmp84*tmp84*tmp84);
REAL8 tmp75=1/tmp30;
REAL8 tmp78=1/tmp56;
REAL8 tmp79=tmp26*tmp65*tmp75*tmp77*tmp78;
REAL8 tmp124=2.*eta*tmp26*tmp82*tmp85;
REAL8 tmp125=1.*tmp19*tmp77*tmp93;
REAL8 tmp126=1.*tmp102*tmp108*tmp19*tmp30*tmp56*tmp77;
REAL8 tmp127=1.*tmp118*tmp19*tmp71*tmp76;
REAL8 tmp128=1.+tmp124+tmp125+tmp126+tmp127;
REAL8 tmp110=-6.*eta*tmp103*tmp34;
REAL8 tmp111=-12.*eta*tmp36;
REAL8 tmp112=tmp110+tmp111;
REAL8 tmp113=1/tmp106;
REAL8 tmp147=(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp96=-2.*tmp93*tmp94*tmp95;
REAL8 tmp97=2.*p->data[1]*tmp77*tmp88*tmp92;
REAL8 tmp98=2.*tmp77*tmp93*x->data[0];
REAL8 tmp109=1.*eta*tmp102*tmp108*tmp19*tmp30*tmp44*tmp51*tmp77;
REAL8 tmp114=1.*tmp102*tmp112*tmp113*tmp19*tmp30*tmp56*tmp77;
REAL8 tmp115=-2.*tmp102*tmp108*tmp30*tmp56*tmp94*tmp95;
REAL8 tmp116=2.*tmp102*tmp108*tmp30*tmp56*tmp77*x->data[0];
REAL8 tmp117=1.*tmp102*tmp108*tmp19*tmp56*tmp60*tmp77;
REAL8 tmp119=-(tmp118*tmp19*tmp62*tmp66*tmp76);
REAL8 tmp120=2.*tmp118*tmp71*tmp94;
REAL8 tmp121=-2.*p->data[2]*tmp13*tmp71*tmp76*tmp8;
REAL8 tmp122=2.*tmp118*tmp71*tmp76*x->data[0];
REAL8 tmp154=tmp109+tmp114+tmp115+tmp116+tmp117+tmp119+tmp120+tmp121+tmp122+tmp96+tmp97+tmp98;
REAL8 tmp157=0.+tmp125+tmp126+tmp127;
REAL8 tmp171=(tmp101*tmp101*tmp101*tmp101);
REAL8 tmp172=tmp30*tmp30;
REAL8 tmp173=pow(x->data[0],6.);
REAL8 tmp176=eta*eta;
REAL8 tmp174=tmp108*tmp108;
REAL8 tmp177=tmp56*tmp56;
REAL8 tmp185=103.*eta;
REAL8 tmp186=-60.*tmp176;
REAL8 tmp187=tmp185+tmp186;
REAL8 tmp190=3.*eta;
REAL8 tmp191=23.+tmp190;
REAL8 tmp195=-16.*eta;
REAL8 tmp196=21.*tmp176;
REAL8 tmp197=tmp195+tmp196;
REAL8 tmp202=-47.*eta;
REAL8 tmp203=54.*tmp176;
REAL8 tmp204=tmp157*tmp197*x->data[0];
REAL8 tmp205=tmp202+tmp203+tmp204;
REAL8 tmp193=tmp157*tmp157;
REAL8 tmp170=(eta*eta*eta);
REAL8 tmp179=pow(x->data[0],7.);
REAL8 tmp180=1./(tmp76*tmp76*tmp76);
REAL8 tmp182=pow(x->data[0],5.);
REAL8 tmp227=-109.*eta;
REAL8 tmp228=51.*tmp176;
REAL8 tmp229=tmp227+tmp228;
REAL8 tmp234=-6.*eta;
REAL8 tmp235=39.*tmp176;
REAL8 tmp236=tmp234+tmp235;
REAL8 tmp241=16.*eta;
REAL8 tmp242=147.*tmp176;
REAL8 tmp243=tmp157*tmp236*x->data[0];
REAL8 tmp244=tmp241+tmp242+tmp243;
REAL8 tmp262=sqrt(tmp76);
REAL8 tmp142=-3.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp34;
REAL8 tmp143=-36.*eta*sigmaKerr->data[2]*tmp102*tmp108*tmp30*tmp44*tmp51*tmp77*tmp94;
REAL8 tmp144=-30.*eta*sigmaStar->data[2]*tmp102*tmp108*tmp30*tmp44*tmp51*tmp77*tmp94;
REAL8 tmp145=-36.*sigmaKerr->data[2]*tmp102*tmp112*tmp113*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp146=-30.*sigmaStar->data[2]*tmp102*tmp112*tmp113*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp148=72.*sigmaKerr->data[2]*tmp102*tmp108*tmp147*tmp30*tmp56*tmp95;
REAL8 tmp149=60.*sigmaStar->data[2]*tmp102*tmp108*tmp147*tmp30*tmp56*tmp95;
REAL8 tmp150=-108.*sigmaKerr->data[2]*tmp102*tmp108*tmp19*tmp30*tmp56*tmp77;
REAL8 tmp151=-90.*sigmaStar->data[2]*tmp102*tmp108*tmp19*tmp30*tmp56*tmp77;
REAL8 tmp152=-36.*sigmaKerr->data[2]*tmp102*tmp108*tmp56*tmp60*tmp77*tmp94;
REAL8 tmp153=-30.*sigmaStar->data[2]*tmp102*tmp108*tmp56*tmp60*tmp77*tmp94;
REAL8 tmp155=3.*sigmaKerr->data[2]*tmp154*x->data[0];
REAL8 tmp156=4.*sigmaStar->data[2]*tmp154*x->data[0];
REAL8 tmp158=3.*sigmaKerr->data[2]*tmp157;
REAL8 tmp159=4.*sigmaStar->data[2]*tmp157;
REAL8 tmp160=tmp143+tmp144+tmp145+tmp146+tmp148+tmp149+tmp150+tmp151+tmp152+tmp153+tmp155+tmp156+tmp158+tmp159;
REAL8 tmp161=0.08333333333333333*eta*tmp160*tmp9;
REAL8 tmp162=-8.*sigmaKerr->data[2];
REAL8 tmp163=14.*sigmaStar->data[2];
REAL8 tmp164=-36.*sigmaKerr->data[2]*tmp102*tmp108*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp165=-30.*sigmaStar->data[2]*tmp102*tmp108*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp166=3.*sigmaKerr->data[2]*tmp157*x->data[0];
REAL8 tmp167=4.*sigmaStar->data[2]*tmp157*x->data[0];
REAL8 tmp168=tmp162+tmp163+tmp164+tmp165+tmp166+tmp167;
REAL8 tmp169=-0.08333333333333333*eta*tmp168*tmp26;
REAL8 tmp175=720.*tmp170*tmp171*tmp172*tmp173*tmp174*tmp44*tmp51*tmp56*tmp95;
REAL8 tmp178=720.*tmp108*tmp112*tmp113*tmp171*tmp172*tmp173*tmp176*tmp177*tmp95;
REAL8 tmp181=-1440.*tmp171*tmp172*tmp174*tmp176*tmp177*tmp179*tmp180;
REAL8 tmp183=2160.*tmp171*tmp172*tmp174*tmp176*tmp177*tmp182*tmp95;
REAL8 tmp184=720.*tmp171*tmp173*tmp174*tmp176*tmp177*tmp30*tmp60*tmp95;
REAL8 tmp188=-2.*tmp154*tmp187*x->data[0];
REAL8 tmp189=-2.*tmp157*tmp187;
REAL8 tmp192=-2.*eta*tmp154*tmp157*tmp19*tmp191;
REAL8 tmp194=-2.*eta*tmp191*tmp193*x->data[0];
REAL8 tmp198=tmp154*tmp197*x->data[0];
REAL8 tmp199=tmp157*tmp197;
REAL8 tmp200=tmp198+tmp199;
REAL8 tmp201=-6.*tmp102*tmp108*tmp200*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp206=-6.*eta*tmp102*tmp108*tmp205*tmp30*tmp44*tmp51*tmp77*tmp94;
REAL8 tmp207=-6.*tmp102*tmp112*tmp113*tmp205*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp208=12.*tmp102*tmp108*tmp147*tmp205*tmp30*tmp56*tmp95;
REAL8 tmp209=-18.*tmp102*tmp108*tmp19*tmp205*tmp30*tmp56*tmp77;
REAL8 tmp210=-6.*tmp102*tmp108*tmp205*tmp56*tmp60*tmp77*tmp94;
REAL8 tmp211=tmp175+tmp178+tmp181+tmp183+tmp184+tmp188+tmp189+tmp192+tmp194+tmp201+tmp206+tmp207+tmp208+tmp209+tmp210;
REAL8 tmp212=0.013888888888888888*sigmaStar->data[2]*tmp211*tmp26;
REAL8 tmp213=27.*eta;
REAL8 tmp214=-353.+tmp213;
REAL8 tmp215=-2.*eta*tmp214;
REAL8 tmp216=360.*tmp171*tmp172*tmp173*tmp174*tmp176*tmp177*tmp95;
REAL8 tmp217=-2.*tmp157*tmp187*x->data[0];
REAL8 tmp218=-(eta*tmp19*tmp191*tmp193);
REAL8 tmp219=-6.*tmp102*tmp108*tmp205*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp220=tmp215+tmp216+tmp217+tmp218+tmp219;
REAL8 tmp221=-0.027777777777777776*sigmaStar->data[2]*tmp220*tmp36;
REAL8 tmp222=1620.*tmp170*tmp171*tmp172*tmp173*tmp174*tmp44*tmp51*tmp56*tmp95;
REAL8 tmp223=1620.*tmp108*tmp112*tmp113*tmp171*tmp172*tmp173*tmp176*tmp177*tmp95;
REAL8 tmp224=-3240.*tmp171*tmp172*tmp174*tmp176*tmp177*tmp179*tmp180;
REAL8 tmp225=4860.*tmp171*tmp172*tmp174*tmp176*tmp177*tmp182*tmp95;
REAL8 tmp226=1620.*tmp171*tmp173*tmp174*tmp176*tmp177*tmp30*tmp60*tmp95;
REAL8 tmp230=4.*tmp154*tmp229*x->data[0];
REAL8 tmp231=4.*tmp157*tmp229;
REAL8 tmp232=-90.*eta*tmp154*tmp157*tmp19;
REAL8 tmp233=-90.*eta*tmp193*x->data[0];
REAL8 tmp237=tmp154*tmp236*x->data[0];
REAL8 tmp238=tmp157*tmp236;
REAL8 tmp239=tmp237+tmp238;
REAL8 tmp240=-6.*tmp102*tmp108*tmp239*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp245=-6.*eta*tmp102*tmp108*tmp244*tmp30*tmp44*tmp51*tmp77*tmp94;
REAL8 tmp246=-6.*tmp102*tmp112*tmp113*tmp244*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp247=12.*tmp102*tmp108*tmp147*tmp244*tmp30*tmp56*tmp95;
REAL8 tmp248=-18.*tmp102*tmp108*tmp19*tmp244*tmp30*tmp56*tmp77;
REAL8 tmp249=-6.*tmp102*tmp108*tmp244*tmp56*tmp60*tmp77*tmp94;
REAL8 tmp250=tmp222+tmp223+tmp224+tmp225+tmp226+tmp230+tmp231+tmp232+tmp233+tmp240+tmp245+tmp246+tmp247+tmp248+tmp249;
REAL8 tmp251=0.006944444444444444*sigmaKerr->data[2]*tmp250*tmp26;
REAL8 tmp252=8.+tmp190;
REAL8 tmp253=-112.*eta*tmp252;
REAL8 tmp254=810.*tmp171*tmp172*tmp173*tmp174*tmp176*tmp177*tmp95;
REAL8 tmp255=4.*tmp157*tmp229*x->data[0];
REAL8 tmp256=-45.*eta*tmp19*tmp193;
REAL8 tmp257=-6.*tmp102*tmp108*tmp244*tmp30*tmp56*tmp77*tmp94;
REAL8 tmp258=tmp253+tmp254+tmp255+tmp256+tmp257;
REAL8 tmp259=-0.013888888888888888*sigmaKerr->data[2]*tmp258*tmp36;
REAL8 tmp260=tmp142+tmp161+tmp169+tmp212+tmp221+tmp251+tmp259;
REAL8 tmp273=1.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp36;
REAL8 tmp274=0.08333333333333333*eta*tmp168*tmp9;
REAL8 tmp275=0.013888888888888888*sigmaStar->data[2]*tmp220*tmp26;
REAL8 tmp276=0.006944444444444444*sigmaKerr->data[2]*tmp258*tmp26;
REAL8 tmp277=0.+sigmaStar->data[2]+tmp273+tmp274+tmp275+tmp276;
REAL8 tmp68=2.*tmp15;
REAL8 tmp69=0.;
REAL8 tmp70=tmp68+tmp69;
REAL8 tmp279=e3z*tmp277;
REAL8 tmp280=0.+tmp279;
REAL8 tmp263=tmp19*tmp30*tmp56;
REAL8 tmp264=sqrt(tmp263);
REAL8 tmp265=-tmp264;
REAL8 tmp266=tmp19*tmp30*tmp56*tmp71*tmp76;
REAL8 tmp267=sqrt(tmp266);
REAL8 tmp268=1.*tmp262*tmp267;
REAL8 tmp269=tmp265+tmp268;
REAL8 tmp270=1.+tmp125+tmp126+tmp127;
REAL8 tmp271=1./sqrt(tmp270);
REAL8 tmp286=1./sqrt(tmp76);
REAL8 tmp290=1./sqrt(tmp263);
REAL8 tmp297=1./sqrt(tmp266);
REAL8 tmp292=2.*tmp30*tmp56*x->data[0];
REAL8 tmp340=tmp76*tmp76;
REAL8 tmp342=tmp277*tmp88;
REAL8 tmp343=0.+tmp342;
REAL8 tmp337=sqrt(tmp270);
REAL8 tmp351=1.+tmp125+tmp126+tmp127+tmp337;
REAL8 tmp363=-(tmp102*tmp108*tmp19*tmp30*tmp343*tmp56);
REAL8 tmp364=tmp343*tmp351*tmp76;
REAL8 tmp365=0.+tmp363+tmp364;
REAL8 tmp310=(1.0/sqrt(tmp76*tmp76*tmp76));
REAL8 tmp311=tmp108*tmp19*tmp30*tmp56;
REAL8 tmp312=sqrt(tmp311);
REAL8 tmp313=0.+tmp68+tmp69;
REAL8 tmp314=tmp313*tmp65;
REAL8 tmp315=-4.*tmp20*x->data[0];
REAL8 tmp316=2.*coeffs->k2;
REAL8 tmp317=3.*coeffs->k3;
REAL8 tmp318=4.*coeffs->k4;
REAL8 tmp319=5.*tmp42*tmp9;
REAL8 tmp320=tmp318+tmp319;
REAL8 tmp321=1.*tmp320*tmp9;
REAL8 tmp322=tmp317+tmp321;
REAL8 tmp323=1.*tmp322*tmp9;
REAL8 tmp324=tmp316+tmp323;
REAL8 tmp325=1.*tmp324*tmp9;
REAL8 tmp326=coeffs->k1+tmp325;
REAL8 tmp327=-(eta*tmp30*tmp326*tmp51);
REAL8 tmp328=1.*tmp28;
REAL8 tmp329=1.*tmp14*tmp9;
REAL8 tmp330=tmp328+tmp329;
REAL8 tmp331=-2.*tmp330*tmp56;
REAL8 tmp332=tmp292+tmp327+tmp331;
REAL8 tmp333=1.*tmp14*tmp332;
REAL8 tmp334=tmp315+tmp333;
REAL8 tmp335=tmp18*tmp334;
REAL8 tmp336=tmp314+tmp335;
REAL8 tmp338=1.+tmp337;
REAL8 tmp339=1/tmp338;
REAL8 tmp373=1.*tmp118*tmp147*tmp30*tmp340*tmp343*tmp56*tmp71;
REAL8 tmp374=1.*tmp19*tmp30*tmp365*tmp56;
REAL8 tmp375=0.+tmp373+tmp374;
REAL8 tmp283=(1.0/sqrt(tmp270*tmp270*tmp270));
REAL8 tmp291=eta*tmp19*tmp30*tmp44*tmp51;
REAL8 tmp293=tmp19*tmp56*tmp60;
REAL8 tmp294=tmp291+tmp292+tmp293;
REAL8 tmp298=-(tmp19*tmp30*tmp56*tmp62*tmp66*tmp76);
REAL8 tmp299=eta*tmp19*tmp30*tmp44*tmp51*tmp71*tmp76;
REAL8 tmp300=2.*tmp30*tmp56*tmp71*tmp94;
REAL8 tmp301=2.*tmp30*tmp56*tmp71*tmp76*x->data[0];
REAL8 tmp302=tmp19*tmp56*tmp60*tmp71*tmp76;
REAL8 tmp303=tmp298+tmp299+tmp300+tmp301+tmp302;
REAL8 tmp357=0.5*tmp154*tmp271;
REAL8 tmp358=tmp109+tmp114+tmp115+tmp116+tmp117+tmp119+tmp120+tmp121+tmp122+tmp357+tmp96+tmp97+tmp98;
REAL8 tmp134=(1.0/(tmp56*tmp56));
REAL8 tmp426=tmp108*tmp147*tmp172*tmp177;
REAL8 tmp427=1./sqrt(tmp426);
REAL8 tmp428=-2.*tmp19*tmp30*tmp56;
REAL8 tmp429=tmp312*tmp332;
REAL8 tmp430=tmp428+tmp429;
REAL8 tmp431=-0.5*tmp13*tmp262*tmp267*tmp338*tmp343*tmp427*tmp430*x->data[0];
REAL8 tmp432=1.*tmp77*x->data[0];
REAL8 tmp433=-4.*tmp30*tmp56*tmp94;
REAL8 tmp434=tmp20*tmp332;
REAL8 tmp435=tmp433+tmp434;
REAL8 tmp436=0.5*tmp20*tmp26*tmp435*tmp71*tmp75*tmp78;
REAL8 tmp437=tmp432+tmp436;
REAL8 tmp438=2.*tmp337;
REAL8 tmp439=1.+tmp438;
REAL8 tmp440=1.*tmp13*tmp262*tmp267*tmp343*tmp437*tmp439*x->data[0];
REAL8 tmp441=0.+tmp440;
REAL8 tmp442=tmp264*tmp441;
REAL8 tmp443=tmp431+tmp442;
REAL8 tmp444=tmp312*tmp443;
REAL8 tmp445=0.+tmp444;
REAL8 tmp447=1/tmp351;
REAL8 tmp138=(1.0/(tmp30*tmp30));
REAL8 tmp416=1./sqrt(tmp311);
REAL8 tmp417=eta*tmp108*tmp19*tmp30*tmp44*tmp51;
REAL8 tmp418=tmp112*tmp113*tmp19*tmp30*tmp56;
REAL8 tmp419=2.*tmp108*tmp30*tmp56*x->data[0];
REAL8 tmp420=tmp108*tmp19*tmp56*tmp60;
REAL8 tmp421=tmp417+tmp418+tmp419+tmp420;
REAL8 tmp381=-2.*eta*tmp330*tmp44*tmp51;
REAL8 tmp382=2.*eta*tmp30*tmp44*tmp51*x->data[0];
REAL8 tmp383=-5.*coeffs->k5l*tmp26;
REAL8 tmp384=-5.*tmp26*tmp42;
REAL8 tmp385=tmp383+tmp384;
REAL8 tmp386=1.*tmp385*tmp9;
REAL8 tmp387=-(tmp26*tmp320);
REAL8 tmp388=tmp386+tmp387;
REAL8 tmp389=1.*tmp388*tmp9;
REAL8 tmp390=-(tmp26*tmp322);
REAL8 tmp391=tmp389+tmp390;
REAL8 tmp392=1.*tmp391*tmp9;
REAL8 tmp393=-(tmp26*tmp324);
REAL8 tmp394=tmp392+tmp393;
REAL8 tmp395=-(eta*tmp30*tmp394*tmp51);
REAL8 tmp396=(1.0/(tmp50*tmp50));
REAL8 tmp397=eta*tmp30*tmp326*tmp396*tmp44;
REAL8 tmp398=-(eta*tmp326*tmp51*tmp60);
REAL8 tmp399=2.*tmp30*tmp56;
REAL8 tmp400=2.*tmp14*tmp26*tmp56;
REAL8 tmp401=2.*tmp56*tmp60*x->data[0];
REAL8 tmp402=tmp381+tmp382+tmp395+tmp397+tmp398+tmp399+tmp400+tmp401;
REAL8 tmp1=s1Vec->data[2]*s1Vec->data[2];
REAL8 tmp3=s2Vec->data[2]*s2Vec->data[2];
REAL8 tmp4=tmp1+tmp3;
REAL8 tmp80=1./sqrt(tmp79);
REAL8 tmp132=sqrt(tmp128);
REAL8 tmp307=tmp277*tmp277;
REAL8 tmp308=0.+tmp307;
REAL8 dHdx0=(1.*eta*(1.5*tmp308*tmp34-tmp260*tmp277*tmp36-4.*coeffs->dheffSSv2*eta*tmp4*tmp5-0.25*tmp154*tmp283*tmp290*tmp297*tmp310*tmp312*tmp336*tmp339*tmp375*tmp66+0.25*tmp271*tmp290*tmp297*tmp310*tmp336*tmp339*tmp375*tmp416*tmp421*tmp66-tmp18*tmp280*tmp62*tmp66+0.5*tmp271*tmp290*tmp297*tmp310*tmp312*tmp339*tmp375*tmp66*(0.+tmp18*(-8.*tmp19-4.*tmp20+1.*tmp14*tmp402)+tmp313*tmp62+tmp334*tmp70)+1.*tmp13*tmp18*tmp71+1.*e3z*tmp18*tmp260*tmp71+1.*tmp13*tmp262*tmp269*tmp271*tmp280*tmp71+1.*tmp13*tmp19*tmp269*tmp271*tmp280*tmp286*tmp71+1.*tmp280*tmp70*tmp71-eta*tmp134*tmp26*tmp267*tmp44*tmp445*tmp447*tmp51*tmp75*tmp77-tmp138*tmp26*tmp267*tmp445*tmp447*tmp60*tmp77*tmp78+0.5*tmp26*tmp297*tmp303*tmp445*tmp447*tmp75*tmp77*tmp78-2.*tmp267*tmp36*tmp445*tmp447*tmp75*tmp77*tmp78-p->data[2]*tmp18*tmp71*tmp8*tmp9-p->data[2]*tmp262*tmp269*tmp271*tmp280*tmp71*tmp8*tmp9-2.*tmp267*tmp445*tmp447*tmp75*tmp78*tmp9*tmp95-tmp13*tmp18*tmp62*tmp66*x->data[0]-tmp13*tmp262*tmp269*tmp271*tmp280*tmp62*tmp66*x->data[0]+1.*e3z*tmp13*tmp260*tmp262*tmp269*tmp271*tmp71*x->data[0]-0.5*tmp13*tmp154*tmp262*tmp269*tmp280*tmp283*tmp71*x->data[0]+1.*tmp13*tmp70*tmp71*x->data[0]+1.*tmp13*tmp262*tmp271*tmp280*tmp71*x->data[0]*(-0.5*tmp290*tmp294+0.5*tmp262*tmp297*tmp303+1.*tmp267*tmp286*x->data[0])+0.5*tmp271*tmp290*tmp297*tmp310*tmp312*tmp336*tmp339*tmp66*(1.*eta*tmp19*tmp30*tmp365*tmp44*tmp51+1.*tmp19*tmp365*tmp56*tmp60-tmp118*tmp147*tmp30*tmp340*tmp343*tmp56*tmp62*tmp66+1.*eta*tmp118*tmp147*tmp30*tmp340*tmp343*tmp44*tmp51*tmp71+1.*tmp118*tmp147*tmp340*tmp343*tmp56*tmp60*tmp71+4.*tmp118*tmp182*tmp30*tmp343*tmp56*tmp71*tmp76-2.*p->data[2]*tmp13*tmp19*tmp30*tmp340*tmp343*tmp56*tmp71*tmp8+1.*tmp118*tmp147*tmp260*tmp30*tmp340*tmp56*tmp71*tmp88+4.*tmp118*tmp30*tmp340*tmp343*tmp56*tmp71*tmp94+2.*tmp30*tmp365*tmp56*x->data[0]+1.*tmp19*tmp30*tmp56*(-(eta*tmp102*tmp108*tmp19*tmp30*tmp343*tmp44*tmp51)-tmp102*tmp112*tmp113*tmp19*tmp30*tmp343*tmp56-tmp102*tmp108*tmp19*tmp343*tmp56*tmp60+tmp343*tmp358*tmp76-tmp102*tmp108*tmp19*tmp260*tmp30*tmp56*tmp88+tmp260*tmp351*tmp76*tmp88+2.*tmp343*tmp351*x->data[0]-2.*tmp102*tmp108*tmp30*tmp343*tmp56*x->data[0]))-0.25*tmp271*tmp294*tmp297*tmp310*tmp312*tmp336*tmp339*tmp375*tmp66*(1.0/sqrt(tmp263*tmp263*tmp263))-0.25*tmp271*tmp290*tmp303*tmp310*tmp312*tmp336*tmp339*tmp375*tmp66*(1.0/sqrt(tmp266*tmp266*tmp266))-(0.25*tmp154*tmp290*tmp297*tmp310*tmp312*tmp336*tmp375*tmp66*(1.0/(tmp338*tmp338)))/tmp270-tmp26*tmp267*tmp358*tmp445*tmp75*tmp77*tmp78*(1.0/(tmp351*tmp351))+1.*tmp26*tmp267*tmp447*tmp75*tmp77*tmp78*(0.5*tmp416*tmp421*tmp443+tmp312*(-0.5*tmp13*tmp262*tmp267*tmp338*tmp343*tmp427*tmp430-0.5*tmp13*tmp19*tmp267*tmp286*tmp338*tmp343*tmp427*tmp430+0.5*tmp290*tmp294*tmp441+0.5*p->data[2]*tmp262*tmp267*tmp338*tmp343*tmp427*tmp430*tmp8*tmp9-0.25*tmp13*tmp154*tmp262*tmp267*tmp271*tmp343*tmp427*tmp430*x->data[0]-0.25*tmp13*tmp262*tmp297*tmp303*tmp338*tmp343*tmp427*tmp430*x->data[0]-0.5*tmp13*tmp260*tmp262*tmp267*tmp338*tmp427*tmp430*tmp88*x->data[0]-0.5*tmp13*tmp262*tmp267*tmp338*tmp343*tmp427*x->data[0]*(tmp312*tmp402+0.5*tmp332*tmp416*tmp421-2.*eta*tmp19*tmp30*tmp44*tmp51-2.*tmp19*tmp56*tmp60-4.*tmp30*tmp56*x->data[0])+tmp264*(1.*tmp13*tmp262*tmp267*tmp343*tmp437*tmp439+1.*tmp13*tmp19*tmp267*tmp286*tmp343*tmp437*tmp439-p->data[2]*tmp262*tmp267*tmp343*tmp437*tmp439*tmp8*tmp9+1.*tmp13*tmp154*tmp262*tmp267*tmp271*tmp343*tmp437*x->data[0]+0.5*tmp13*tmp262*tmp297*tmp303*tmp343*tmp437*tmp439*x->data[0]+1.*tmp13*tmp260*tmp262*tmp267*tmp437*tmp439*tmp88*x->data[0]+1.*tmp13*tmp262*tmp267*tmp343*tmp439*x->data[0]*(-0.5*eta*tmp134*tmp20*tmp26*tmp435*tmp44*tmp51*tmp71*tmp75+1.*tmp77-0.5*tmp138*tmp20*tmp26*tmp435*tmp60*tmp71*tmp78-0.5*tmp20*tmp26*tmp435*tmp62*tmp66*tmp75*tmp78-tmp20*tmp36*tmp435*tmp71*tmp75*tmp78+1.*tmp435*tmp71*tmp75*tmp78*tmp9-2.*tmp19*tmp95+0.5*tmp20*tmp26*tmp71*tmp75*tmp78*(tmp20*tmp402-12.*tmp19*tmp30*tmp56-4.*eta*tmp30*tmp44*tmp51*tmp94-4.*tmp56*tmp60*tmp94+2.*tmp332*x->data[0])))+0.25*tmp13*tmp262*tmp267*tmp338*tmp343*tmp430*(tmp112*tmp113*tmp147*tmp172*tmp177+2.*eta*tmp108*tmp147*tmp172*tmp44*tmp51*tmp56+2.*tmp108*tmp147*tmp177*tmp30*tmp60+4.*tmp108*tmp172*tmp177*tmp94)*x->data[0]*(1.0/sqrt(tmp426*tmp426*tmp426))))-tmp271*tmp290*tmp297*tmp310*tmp312*tmp336*tmp339*tmp375*tmp62*1./(tmp65*tmp65*tmp65)-1.5*tmp271*tmp290*tmp297*tmp312*tmp336*tmp339*tmp375*tmp66*x->data[0]*pow(tmp76,-2.5)-0.5*tmp132*(-(eta*tmp134*tmp26*tmp44*tmp51*tmp65*tmp75*tmp77)-tmp138*tmp26*tmp60*tmp65*tmp77*tmp78+tmp26*tmp62*tmp75*tmp77*tmp78-2.*tmp36*tmp65*tmp75*tmp77*tmp78-2.*tmp65*tmp75*tmp78*tmp9*tmp95)*(1.0/sqrt(tmp79*tmp79*tmp79))+(0.5*tmp80*(tmp109+tmp114+tmp115+tmp116+tmp117+tmp119+tmp120+tmp121+tmp122-4.*eta*tmp36*tmp82*tmp85+tmp96+tmp97+tmp98))/sqrt(tmp128)))/sqrt(1.+2.*eta*(-1.-0.5*tmp308*tmp36+1.*coeffs->dheffSSv2*eta*tmp34*tmp4+0.5*tmp271*tmp290*tmp297*tmp310*tmp312*tmp336*tmp339*tmp375*tmp66+1.*tmp18*tmp280*tmp71+1.*tmp26*tmp267*tmp445*tmp447*tmp75*tmp77*tmp78+tmp132*tmp80+1.*tmp13*tmp18*tmp71*x->data[0]+1.*tmp13*tmp262*tmp269*tmp271*tmp280*tmp71*x->data[0]));