SEOBNRv2_opt_dx0dx0_sph.h

Go to the documentation of this file.

REAL8 tmp60=1/x->data[0];
REAL8 tmp149=sigmaKerr->data[2]*sigmaKerr->data[2];
REAL8 tmp152=sqrt(tmp149);
REAL8 tmp161=coeffs->KK*eta;
REAL8 tmp164=-1.+tmp161;
REAL8 tmp158=x->data[0]*x->data[0];
REAL8 tmp56=pow(x->data[0],-5.);
REAL8 tmp167=(1.0/(x->data[0]*x->data[0]));
REAL8 tmp191=pow(x->data[0],-6.);
REAL8 tmp196=1./(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp201=1./(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp205=1.*tmp60;
REAL8 tmp206=log(tmp205);
REAL8 tmp207=coeffs->k5l*tmp206;
REAL8 tmp208=coeffs->k5+tmp207;
REAL8 tmp165=(1.0/(tmp164*tmp164));
REAL8 tmp166=1.*tmp165;
REAL8 tmp168=1.*tmp149*tmp167;
REAL8 tmp169=1/tmp164;
REAL8 tmp170=2.*tmp169*tmp60;
REAL8 tmp190=tmp166+tmp168+tmp170;
REAL8 tmp211=1.*coeffs->k4*tmp196;
REAL8 tmp212=1.*coeffs->k3*tmp201;
REAL8 tmp213=1.*coeffs->k2*tmp167;
REAL8 tmp214=1.*coeffs->k1*tmp60;
REAL8 tmp215=1.*tmp208*tmp56;
REAL8 tmp216=1.+tmp211+tmp212+tmp213+tmp214+tmp215;
REAL8 tmp219=coeffs->k0*eta;
REAL8 tmp220=log(tmp216);
REAL8 tmp221=eta*tmp220;
REAL8 tmp222=1.+tmp219+tmp221;
REAL8 tmp159=tmp149+tmp158;
REAL8 tmp58=1.*e3z;
REAL8 tmp59=0.+tmp58;
REAL8 tmp81=p->data[2]*tmp60;
REAL8 tmp82=0.+tmp81;
REAL8 tmp83=tmp59*tmp82;
REAL8 tmp145=0.+tmp83;
REAL8 tmp229=tmp159*tmp159;
REAL8 tmp230=-(tmp149*tmp158*tmp190*tmp222);
REAL8 tmp231=tmp229+tmp230;
REAL8 tmp153=0.;
REAL8 tmp156=2.*tmp152*x->data[0];
REAL8 tmp157=tmp153+tmp156;
REAL8 tmp237=1/tmp231;
REAL8 tmp242=0.+tmp158;
REAL8 tmp253=1.*tmp59;
REAL8 tmp254=0.+tmp253;
REAL8 tmp255=-(p->data[1]*tmp60);
REAL8 tmp256=0.+tmp255;
REAL8 tmp257=tmp254*tmp256;
REAL8 tmp258=0.+tmp257;
REAL8 tmp243=1/tmp242;
REAL8 tmp259=tmp258*tmp258;
REAL8 tmp247=-3.*eta;
REAL8 tmp192=-(coeffs->k5l*tmp191);
REAL8 tmp193=-4.*coeffs->k4*tmp56;
REAL8 tmp198=-3.*coeffs->k3*tmp196;
REAL8 tmp202=-2.*coeffs->k2*tmp201;
REAL8 tmp204=-(coeffs->k1*tmp167);
REAL8 tmp209=-5.*tmp191*tmp208;
REAL8 tmp210=tmp192+tmp193+tmp198+tmp202+tmp204+tmp209;
REAL8 tmp217=1/tmp216;
REAL8 tmp265=0.+p->data[0];
REAL8 tmp266=1.*tmp265;
REAL8 tmp267=0.+tmp266;
REAL8 tmp268=tmp267*tmp267;
REAL8 tmp269=26.+tmp247;
REAL8 tmp270=2.*eta*tmp201*tmp269;
REAL8 tmp271=6.*eta*tmp167;
REAL8 tmp272=1.+tmp270+tmp271;
REAL8 tmp260=(x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp261=(1.0/(tmp242*tmp242));
REAL8 tmp273=log(tmp272);
REAL8 tmp274=1.+tmp273;
REAL8 tmp224=-2.*tmp149*tmp201;
REAL8 tmp225=-2.*tmp167*tmp169;
REAL8 tmp226=tmp224+tmp225;
REAL8 tmp160=4.*tmp159*x->data[0];
REAL8 tmp218=-(eta*tmp149*tmp158*tmp190*tmp210*tmp217);
REAL8 tmp223=-2.*tmp149*tmp190*tmp222*x->data[0];
REAL8 tmp227=-(tmp149*tmp158*tmp222*tmp226);
REAL8 tmp228=tmp160+tmp218+tmp223+tmp227;
REAL8 tmp232=(1.0/(tmp231*tmp231));
REAL8 tmp284=tmp145*tmp145;
REAL8 tmp248=4.+tmp247;
REAL8 tmp249=1.*p->data[0];
REAL8 tmp250=0.+tmp249;
REAL8 tmp251=(tmp250*tmp250*tmp250*tmp250);
REAL8 tmp241=1/tmp190;
REAL8 tmp244=1/tmp222;
REAL8 tmp245=tmp167*tmp231*tmp241*tmp243*tmp244;
REAL8 tmp290=2.*eta*tmp167*tmp248*tmp251;
REAL8 tmp291=1.*tmp158*tmp243*tmp259;
REAL8 tmp292=1.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274;
REAL8 tmp293=1.*tmp158*tmp237*tmp242*tmp284;
REAL8 tmp294=1.+tmp290+tmp291+tmp292+tmp293;
REAL8 tmp276=-6.*eta*tmp196*tmp269;
REAL8 tmp277=-12.*eta*tmp201;
REAL8 tmp278=tmp276+tmp277;
REAL8 tmp279=1/tmp272;
REAL8 tmp313=(x->data[0]*x->data[0]*x->data[0]*x->data[0]);
REAL8 tmp262=-2.*tmp259*tmp260*tmp261;
REAL8 tmp263=2.*p->data[1]*tmp243*tmp254*tmp258;
REAL8 tmp264=2.*tmp243*tmp259*x->data[0];
REAL8 tmp275=1.*eta*tmp158*tmp190*tmp210*tmp217*tmp243*tmp268*tmp274;
REAL8 tmp280=1.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp278*tmp279;
REAL8 tmp281=-2.*tmp190*tmp222*tmp260*tmp261*tmp268*tmp274;
REAL8 tmp282=2.*tmp190*tmp222*tmp243*tmp268*tmp274*x->data[0];
REAL8 tmp283=1.*tmp158*tmp222*tmp226*tmp243*tmp268*tmp274;
REAL8 tmp285=-(tmp158*tmp228*tmp232*tmp242*tmp284);
REAL8 tmp286=2.*tmp237*tmp260*tmp284;
REAL8 tmp287=-2.*p->data[2]*tmp145*tmp237*tmp242*tmp59;
REAL8 tmp288=2.*tmp237*tmp242*tmp284*x->data[0];
REAL8 tmp320=tmp262+tmp263+tmp264+tmp275+tmp280+tmp281+tmp282+tmp283+tmp285+tmp286+tmp287+tmp288;
REAL8 tmp323=0.+tmp291+tmp292+tmp293;
REAL8 tmp5=eta*eta;
REAL8 tmp337=(tmp267*tmp267*tmp267*tmp267);
REAL8 tmp338=tmp190*tmp190;
REAL8 tmp339=pow(x->data[0],6.);
REAL8 tmp340=tmp274*tmp274;
REAL8 tmp342=tmp222*tmp222;
REAL8 tmp350=103.*eta;
REAL8 tmp351=-60.*tmp5;
REAL8 tmp352=tmp350+tmp351;
REAL8 tmp355=3.*eta;
REAL8 tmp356=23.+tmp355;
REAL8 tmp360=-16.*eta;
REAL8 tmp361=21.*tmp5;
REAL8 tmp362=tmp360+tmp361;
REAL8 tmp367=-47.*eta;
REAL8 tmp368=54.*tmp5;
REAL8 tmp369=tmp323*tmp362*x->data[0];
REAL8 tmp370=tmp367+tmp368+tmp369;
REAL8 tmp358=tmp323*tmp323;
REAL8 tmp336=(eta*eta*eta);
REAL8 tmp344=pow(x->data[0],7.);
REAL8 tmp345=1./(tmp242*tmp242*tmp242);
REAL8 tmp347=pow(x->data[0],5.);
REAL8 tmp392=-109.*eta;
REAL8 tmp393=51.*tmp5;
REAL8 tmp394=tmp392+tmp393;
REAL8 tmp399=-6.*eta;
REAL8 tmp400=39.*tmp5;
REAL8 tmp401=tmp399+tmp400;
REAL8 tmp406=16.*eta;
REAL8 tmp407=147.*tmp5;
REAL8 tmp408=tmp323*tmp401*x->data[0];
REAL8 tmp409=tmp406+tmp407+tmp408;
REAL8 tmp427=sqrt(tmp242);
REAL8 tmp308=-3.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp196;
REAL8 tmp309=-36.*eta*sigmaKerr->data[2]*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp310=-30.*eta*sigmaStar->data[2]*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp311=-36.*sigmaKerr->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp278*tmp279;
REAL8 tmp312=-30.*sigmaStar->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp278*tmp279;
REAL8 tmp314=72.*sigmaKerr->data[2]*tmp190*tmp222*tmp261*tmp268*tmp274*tmp313;
REAL8 tmp315=60.*sigmaStar->data[2]*tmp190*tmp222*tmp261*tmp268*tmp274*tmp313;
REAL8 tmp316=-108.*sigmaKerr->data[2]*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274;
REAL8 tmp317=-90.*sigmaStar->data[2]*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274;
REAL8 tmp318=-36.*sigmaKerr->data[2]*tmp222*tmp226*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp319=-30.*sigmaStar->data[2]*tmp222*tmp226*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp321=3.*sigmaKerr->data[2]*tmp320*x->data[0];
REAL8 tmp322=4.*sigmaStar->data[2]*tmp320*x->data[0];
REAL8 tmp324=3.*sigmaKerr->data[2]*tmp323;
REAL8 tmp325=4.*sigmaStar->data[2]*tmp323;
REAL8 tmp326=tmp309+tmp310+tmp311+tmp312+tmp314+tmp315+tmp316+tmp317+tmp318+tmp319+tmp321+tmp322+tmp324+tmp325;
REAL8 tmp327=0.08333333333333333*eta*tmp326*tmp60;
REAL8 tmp328=-8.*sigmaKerr->data[2];
REAL8 tmp329=14.*sigmaStar->data[2];
REAL8 tmp330=-36.*sigmaKerr->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp331=-30.*sigmaStar->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp332=3.*sigmaKerr->data[2]*tmp323*x->data[0];
REAL8 tmp333=4.*sigmaStar->data[2]*tmp323*x->data[0];
REAL8 tmp334=tmp328+tmp329+tmp330+tmp331+tmp332+tmp333;
REAL8 tmp335=-0.08333333333333333*eta*tmp167*tmp334;
REAL8 tmp341=720.*tmp210*tmp217*tmp222*tmp261*tmp336*tmp337*tmp338*tmp339*tmp340;
REAL8 tmp343=720.*tmp261*tmp274*tmp278*tmp279*tmp337*tmp338*tmp339*tmp342*tmp5;
REAL8 tmp346=-1440.*tmp337*tmp338*tmp340*tmp342*tmp344*tmp345*tmp5;
REAL8 tmp348=2160.*tmp261*tmp337*tmp338*tmp340*tmp342*tmp347*tmp5;
REAL8 tmp349=720.*tmp190*tmp226*tmp261*tmp337*tmp339*tmp340*tmp342*tmp5;
REAL8 tmp353=-2.*tmp320*tmp352*x->data[0];
REAL8 tmp354=-2.*tmp323*tmp352;
REAL8 tmp357=-2.*eta*tmp158*tmp320*tmp323*tmp356;
REAL8 tmp359=-2.*eta*tmp356*tmp358*x->data[0];
REAL8 tmp363=tmp320*tmp362*x->data[0];
REAL8 tmp364=tmp323*tmp362;
REAL8 tmp365=tmp363+tmp364;
REAL8 tmp366=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274*tmp365;
REAL8 tmp371=-6.*eta*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp274*tmp370;
REAL8 tmp372=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp278*tmp279*tmp370;
REAL8 tmp373=12.*tmp190*tmp222*tmp261*tmp268*tmp274*tmp313*tmp370;
REAL8 tmp374=-18.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274*tmp370;
REAL8 tmp375=-6.*tmp222*tmp226*tmp243*tmp260*tmp268*tmp274*tmp370;
REAL8 tmp376=tmp341+tmp343+tmp346+tmp348+tmp349+tmp353+tmp354+tmp357+tmp359+tmp366+tmp371+tmp372+tmp373+tmp374+tmp375;
REAL8 tmp377=0.013888888888888888*sigmaStar->data[2]*tmp167*tmp376;
REAL8 tmp378=27.*eta;
REAL8 tmp379=-353.+tmp378;
REAL8 tmp380=-2.*eta*tmp379;
REAL8 tmp381=360.*tmp261*tmp337*tmp338*tmp339*tmp340*tmp342*tmp5;
REAL8 tmp382=-2.*tmp323*tmp352*x->data[0];
REAL8 tmp383=-(eta*tmp158*tmp356*tmp358);
REAL8 tmp384=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274*tmp370;
REAL8 tmp385=tmp380+tmp381+tmp382+tmp383+tmp384;
REAL8 tmp386=-0.027777777777777776*sigmaStar->data[2]*tmp201*tmp385;
REAL8 tmp387=1620.*tmp210*tmp217*tmp222*tmp261*tmp336*tmp337*tmp338*tmp339*tmp340;
REAL8 tmp388=1620.*tmp261*tmp274*tmp278*tmp279*tmp337*tmp338*tmp339*tmp342*tmp5;
REAL8 tmp389=-3240.*tmp337*tmp338*tmp340*tmp342*tmp344*tmp345*tmp5;
REAL8 tmp390=4860.*tmp261*tmp337*tmp338*tmp340*tmp342*tmp347*tmp5;
REAL8 tmp391=1620.*tmp190*tmp226*tmp261*tmp337*tmp339*tmp340*tmp342*tmp5;
REAL8 tmp395=4.*tmp320*tmp394*x->data[0];
REAL8 tmp396=4.*tmp323*tmp394;
REAL8 tmp397=-90.*eta*tmp158*tmp320*tmp323;
REAL8 tmp398=-90.*eta*tmp358*x->data[0];
REAL8 tmp402=tmp320*tmp401*x->data[0];
REAL8 tmp403=tmp323*tmp401;
REAL8 tmp404=tmp402+tmp403;
REAL8 tmp405=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274*tmp404;
REAL8 tmp410=-6.*eta*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp274*tmp409;
REAL8 tmp411=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp278*tmp279*tmp409;
REAL8 tmp412=12.*tmp190*tmp222*tmp261*tmp268*tmp274*tmp313*tmp409;
REAL8 tmp413=-18.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274*tmp409;
REAL8 tmp414=-6.*tmp222*tmp226*tmp243*tmp260*tmp268*tmp274*tmp409;
REAL8 tmp415=tmp387+tmp388+tmp389+tmp390+tmp391+tmp395+tmp396+tmp397+tmp398+tmp405+tmp410+tmp411+tmp412+tmp413+tmp414;
REAL8 tmp416=0.006944444444444444*sigmaKerr->data[2]*tmp167*tmp415;
REAL8 tmp417=8.+tmp355;
REAL8 tmp418=-112.*eta*tmp417;
REAL8 tmp419=810.*tmp261*tmp337*tmp338*tmp339*tmp340*tmp342*tmp5;
REAL8 tmp420=4.*tmp323*tmp394*x->data[0];
REAL8 tmp421=-45.*eta*tmp158*tmp358;
REAL8 tmp422=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274*tmp409;
REAL8 tmp423=tmp418+tmp419+tmp420+tmp421+tmp422;
REAL8 tmp424=-0.013888888888888888*sigmaKerr->data[2]*tmp201*tmp423;
REAL8 tmp425=tmp308+tmp327+tmp335+tmp377+tmp386+tmp416+tmp424;
REAL8 tmp438=1.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp201;
REAL8 tmp439=0.08333333333333333*eta*tmp334*tmp60;
REAL8 tmp440=0.013888888888888888*sigmaStar->data[2]*tmp167*tmp385;
REAL8 tmp441=0.006944444444444444*sigmaKerr->data[2]*tmp167*tmp423;
REAL8 tmp442=0.+sigmaStar->data[2]+tmp438+tmp439+tmp440+tmp441;
REAL8 tmp234=2.*tmp152;
REAL8 tmp235=0.;
REAL8 tmp236=tmp234+tmp235;
REAL8 tmp444=e3z*tmp442;
REAL8 tmp445=0.+tmp444;
REAL8 tmp428=tmp158*tmp190*tmp222;
REAL8 tmp429=sqrt(tmp428);
REAL8 tmp430=-tmp429;
REAL8 tmp431=tmp158*tmp190*tmp222*tmp237*tmp242;
REAL8 tmp432=sqrt(tmp431);
REAL8 tmp433=1.*tmp427*tmp432;
REAL8 tmp434=tmp430+tmp433;
REAL8 tmp435=1.+tmp291+tmp292+tmp293;
REAL8 tmp436=1./sqrt(tmp435);
REAL8 tmp451=1./sqrt(tmp242);
REAL8 tmp455=1./sqrt(tmp428);
REAL8 tmp462=1./sqrt(tmp431);
REAL8 tmp457=2.*tmp190*tmp222*x->data[0];
REAL8 tmp505=tmp242*tmp242;
REAL8 tmp507=tmp254*tmp442;
REAL8 tmp508=0.+tmp507;
REAL8 tmp502=sqrt(tmp435);
REAL8 tmp516=1.+tmp291+tmp292+tmp293+tmp502;
REAL8 tmp528=-(tmp158*tmp190*tmp222*tmp268*tmp274*tmp508);
REAL8 tmp529=tmp242*tmp508*tmp516;
REAL8 tmp530=0.+tmp528+tmp529;
REAL8 tmp475=(1.0/sqrt(tmp242*tmp242*tmp242));
REAL8 tmp476=tmp158*tmp190*tmp222*tmp274;
REAL8 tmp477=sqrt(tmp476);
REAL8 tmp478=0.+tmp234+tmp235;
REAL8 tmp479=tmp231*tmp478;
REAL8 tmp480=-4.*tmp159*x->data[0];
REAL8 tmp481=2.*coeffs->k2;
REAL8 tmp482=3.*coeffs->k3;
REAL8 tmp483=4.*coeffs->k4;
REAL8 tmp484=5.*tmp208*tmp60;
REAL8 tmp485=tmp483+tmp484;
REAL8 tmp486=1.*tmp485*tmp60;
REAL8 tmp487=tmp482+tmp486;
REAL8 tmp488=1.*tmp487*tmp60;
REAL8 tmp489=tmp481+tmp488;
REAL8 tmp490=1.*tmp489*tmp60;
REAL8 tmp491=coeffs->k1+tmp490;
REAL8 tmp492=-(eta*tmp190*tmp217*tmp491);
REAL8 tmp493=1.*tmp169;
REAL8 tmp494=1.*tmp149*tmp60;
REAL8 tmp495=tmp493+tmp494;
REAL8 tmp496=-2.*tmp222*tmp495;
REAL8 tmp497=tmp457+tmp492+tmp496;
REAL8 tmp498=1.*tmp149*tmp497;
REAL8 tmp499=tmp480+tmp498;
REAL8 tmp500=tmp157*tmp499;
REAL8 tmp501=tmp479+tmp500;
REAL8 tmp503=1.+tmp502;
REAL8 tmp504=1/tmp503;
REAL8 tmp538=1.*tmp190*tmp222*tmp237*tmp284*tmp313*tmp505*tmp508;
REAL8 tmp539=1.*tmp158*tmp190*tmp222*tmp530;
REAL8 tmp540=0.+tmp538+tmp539;
REAL8 tmp448=(1.0/sqrt(tmp435*tmp435*tmp435));
REAL8 tmp456=eta*tmp158*tmp190*tmp210*tmp217;
REAL8 tmp458=tmp158*tmp222*tmp226;
REAL8 tmp459=tmp456+tmp457+tmp458;
REAL8 tmp463=-(tmp158*tmp190*tmp222*tmp228*tmp232*tmp242);
REAL8 tmp464=eta*tmp158*tmp190*tmp210*tmp217*tmp237*tmp242;
REAL8 tmp465=2.*tmp190*tmp222*tmp237*tmp260;
REAL8 tmp466=2.*tmp190*tmp222*tmp237*tmp242*x->data[0];
REAL8 tmp467=tmp158*tmp222*tmp226*tmp237*tmp242;
REAL8 tmp468=tmp463+tmp464+tmp465+tmp466+tmp467;
REAL8 tmp522=0.5*tmp320*tmp436;
REAL8 tmp523=tmp262+tmp263+tmp264+tmp275+tmp280+tmp281+tmp282+tmp283+tmp285+tmp286+tmp287+tmp288+tmp522;
REAL8 tmp300=(1.0/(tmp222*tmp222));
REAL8 tmp591=tmp274*tmp313*tmp338*tmp342;
REAL8 tmp592=1./sqrt(tmp591);
REAL8 tmp593=-2.*tmp158*tmp190*tmp222;
REAL8 tmp594=tmp477*tmp497;
REAL8 tmp595=tmp593+tmp594;
REAL8 tmp596=-0.5*tmp145*tmp427*tmp432*tmp503*tmp508*tmp592*tmp595*x->data[0];
REAL8 tmp597=1.*tmp243*x->data[0];
REAL8 tmp598=-4.*tmp190*tmp222*tmp260;
REAL8 tmp599=tmp159*tmp497;
REAL8 tmp600=tmp598+tmp599;
REAL8 tmp601=0.5*tmp159*tmp167*tmp237*tmp241*tmp244*tmp600;
REAL8 tmp602=tmp597+tmp601;
REAL8 tmp603=2.*tmp502;
REAL8 tmp604=1.+tmp603;
REAL8 tmp605=1.*tmp145*tmp427*tmp432*tmp508*tmp602*tmp604*x->data[0];
REAL8 tmp606=0.+tmp605;
REAL8 tmp607=tmp429*tmp606;
REAL8 tmp608=tmp596+tmp607;
REAL8 tmp609=tmp477*tmp608;
REAL8 tmp610=0.+tmp609;
REAL8 tmp612=1/tmp516;
REAL8 tmp304=(1.0/(tmp190*tmp190));
REAL8 tmp581=1./sqrt(tmp476);
REAL8 tmp582=eta*tmp158*tmp190*tmp210*tmp217*tmp274;
REAL8 tmp583=tmp158*tmp190*tmp222*tmp278*tmp279;
REAL8 tmp584=2.*tmp190*tmp222*tmp274*x->data[0];
REAL8 tmp585=tmp158*tmp222*tmp226*tmp274;
REAL8 tmp586=tmp582+tmp583+tmp584+tmp585;
REAL8 tmp546=-2.*eta*tmp210*tmp217*tmp495;
REAL8 tmp547=2.*eta*tmp190*tmp210*tmp217*x->data[0];
REAL8 tmp548=-5.*coeffs->k5l*tmp167;
REAL8 tmp549=-5.*tmp167*tmp208;
REAL8 tmp550=tmp548+tmp549;
REAL8 tmp551=1.*tmp550*tmp60;
REAL8 tmp552=-(tmp167*tmp485);
REAL8 tmp553=tmp551+tmp552;
REAL8 tmp554=1.*tmp553*tmp60;
REAL8 tmp555=-(tmp167*tmp487);
REAL8 tmp556=tmp554+tmp555;
REAL8 tmp557=1.*tmp556*tmp60;
REAL8 tmp558=-(tmp167*tmp489);
REAL8 tmp559=tmp557+tmp558;
REAL8 tmp560=-(eta*tmp190*tmp217*tmp559);
REAL8 tmp561=(1.0/(tmp216*tmp216));
REAL8 tmp562=eta*tmp190*tmp210*tmp491*tmp561;
REAL8 tmp563=-(eta*tmp217*tmp226*tmp491);
REAL8 tmp564=2.*tmp190*tmp222;
REAL8 tmp565=2.*tmp149*tmp167*tmp222;
REAL8 tmp566=2.*tmp222*tmp226*x->data[0];
REAL8 tmp567=tmp546+tmp547+tmp560+tmp562+tmp563+tmp564+tmp565+tmp566;
REAL8 tmp48=s1Vec->data[2]*s1Vec->data[2];
REAL8 tmp53=s2Vec->data[2]*s2Vec->data[2];
REAL8 tmp54=tmp48+tmp53;
REAL8 tmp246=1./sqrt(tmp245);
REAL8 tmp298=sqrt(tmp294);
REAL8 tmp472=tmp442*tmp442;
REAL8 tmp473=0.+tmp472;
REAL8 tmp575=1./(tmp231*tmp231*tmp231);
REAL8 tmp690=pow(x->data[0],-7.);
REAL8 tmp252=-4.*eta*tmp201*tmp248*tmp251;
REAL8 tmp289=tmp252+tmp262+tmp263+tmp264+tmp275+tmp280+tmp281+tmp282+tmp283+tmp285+tmp286+tmp287+tmp288;
REAL8 tmp688=tmp210*tmp210;
REAL8 tmp691=11.*coeffs->k5l*tmp690;
REAL8 tmp692=20.*coeffs->k4*tmp191;
REAL8 tmp693=12.*coeffs->k3*tmp56;
REAL8 tmp694=6.*coeffs->k2*tmp196;
REAL8 tmp695=2.*coeffs->k1*tmp201;
REAL8 tmp696=30.*tmp208*tmp690;
REAL8 tmp697=tmp691+tmp692+tmp693+tmp694+tmp695+tmp696;
REAL8 tmp703=6.*tmp149*tmp196;
REAL8 tmp704=4.*tmp169*tmp201;
REAL8 tmp705=tmp703+tmp704;
REAL8 tmp684=tmp228*tmp228;
REAL8 tmp686=8.*tmp158;
REAL8 tmp687=4.*tmp159;
REAL8 tmp689=1.*eta*tmp149*tmp158*tmp190*tmp561*tmp688;
REAL8 tmp698=-(eta*tmp149*tmp158*tmp190*tmp217*tmp697);
REAL8 tmp699=-4.*eta*tmp149*tmp190*tmp210*tmp217*x->data[0];
REAL8 tmp700=-2.*eta*tmp149*tmp158*tmp210*tmp217*tmp226;
REAL8 tmp701=-2.*tmp149*tmp190*tmp222;
REAL8 tmp702=-4.*tmp149*tmp222*tmp226*x->data[0];
REAL8 tmp706=-(tmp149*tmp158*tmp222*tmp705);
REAL8 tmp707=tmp686+tmp687+tmp689+tmp698+tmp699+tmp700+tmp701+tmp702+tmp706;
REAL8 tmp295=1./sqrt(tmp294);
REAL8 tmp297=(1.0/sqrt(tmp245*tmp245*tmp245));
REAL8 tmp299=tmp167*tmp228*tmp241*tmp243*tmp244;
REAL8 tmp301=-(eta*tmp167*tmp210*tmp217*tmp231*tmp241*tmp243*tmp300);
REAL8 tmp302=-2.*tmp231*tmp241*tmp244*tmp261*tmp60;
REAL8 tmp303=-2.*tmp201*tmp231*tmp241*tmp243*tmp244;
REAL8 tmp305=-(tmp167*tmp226*tmp231*tmp243*tmp244*tmp304);
REAL8 tmp306=tmp299+tmp301+tmp302+tmp303+tmp305;
REAL8 tmp736=tmp278*tmp278;
REAL8 tmp737=(1.0/(tmp272*tmp272));
REAL8 tmp739=24.*eta*tmp269*tmp56;
REAL8 tmp740=36.*eta*tmp196;
REAL8 tmp741=tmp739+tmp740;
REAL8 tmp719=8.*tmp259*tmp313*tmp345;
REAL8 tmp720=-8.*p->data[1]*tmp254*tmp258*tmp261*x->data[0];
REAL8 tmp721=-10.*tmp158*tmp259*tmp261;
REAL8 tmp722=2.*tmp243*tmp259;
REAL8 tmp723=tmp254*tmp254;
REAL8 tmp724=p->data[1]*p->data[1];
REAL8 tmp725=2.*tmp167*tmp243*tmp723*tmp724;
REAL8 tmp726=4.*p->data[1]*tmp243*tmp254*tmp258*tmp60;
REAL8 tmp727=-(eta*tmp158*tmp190*tmp243*tmp268*tmp274*tmp561*tmp688);
REAL8 tmp728=1.*eta*tmp158*tmp190*tmp217*tmp243*tmp268*tmp274*tmp697;
REAL8 tmp729=2.*eta*tmp158*tmp190*tmp210*tmp217*tmp243*tmp268*tmp278*tmp279;
REAL8 tmp730=-4.*eta*tmp190*tmp210*tmp217*tmp260*tmp261*tmp268*tmp274;
REAL8 tmp731=4.*eta*tmp190*tmp210*tmp217*tmp243*tmp268*tmp274*x->data[0];
REAL8 tmp732=2.*eta*tmp158*tmp210*tmp217*tmp226*tmp243*tmp268*tmp274;
REAL8 tmp733=-4.*tmp190*tmp222*tmp260*tmp261*tmp268*tmp278*tmp279;
REAL8 tmp734=4.*tmp190*tmp222*tmp243*tmp268*tmp278*tmp279*x->data[0];
REAL8 tmp735=2.*tmp158*tmp222*tmp226*tmp243*tmp268*tmp278*tmp279;
REAL8 tmp738=-(tmp158*tmp190*tmp222*tmp243*tmp268*tmp736*tmp737);
REAL8 tmp742=1.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp279*tmp741;
REAL8 tmp743=8.*tmp190*tmp222*tmp268*tmp274*tmp313*tmp345;
REAL8 tmp744=-10.*tmp158*tmp190*tmp222*tmp261*tmp268*tmp274;
REAL8 tmp745=-4.*tmp222*tmp226*tmp260*tmp261*tmp268*tmp274;
REAL8 tmp746=2.*tmp190*tmp222*tmp243*tmp268*tmp274;
REAL8 tmp747=4.*tmp222*tmp226*tmp243*tmp268*tmp274*x->data[0];
REAL8 tmp748=1.*tmp158*tmp222*tmp243*tmp268*tmp274*tmp705;
REAL8 tmp749=2.*tmp158*tmp242*tmp284*tmp575*tmp684;
REAL8 tmp750=-(tmp158*tmp232*tmp242*tmp284*tmp707);
REAL8 tmp751=-4.*tmp228*tmp232*tmp260*tmp284;
REAL8 tmp752=4.*p->data[2]*tmp145*tmp228*tmp232*tmp242*tmp59;
REAL8 tmp753=-4.*tmp228*tmp232*tmp242*tmp284*x->data[0];
REAL8 tmp754=-8.*p->data[2]*tmp145*tmp237*tmp59*x->data[0];
REAL8 tmp755=10.*tmp158*tmp237*tmp284;
REAL8 tmp756=2.*tmp237*tmp242*tmp284;
REAL8 tmp757=tmp59*tmp59;
REAL8 tmp758=p->data[2]*p->data[2];
REAL8 tmp759=2.*tmp167*tmp237*tmp242*tmp757*tmp758;
REAL8 tmp760=-4.*p->data[2]*tmp145*tmp237*tmp242*tmp59*tmp60;
REAL8 tmp825=tmp719+tmp720+tmp721+tmp722+tmp725+tmp726+tmp727+tmp728+tmp729+tmp730+tmp731+tmp732+tmp733+tmp734+tmp735+tmp738+tmp742+tmp743+tmp744+tmp745+tmp746+tmp747+tmp748+tmp749+tmp750+tmp751+tmp752+tmp753+tmp754+tmp755+tmp756+tmp759+tmp760;
REAL8 tmp780=tmp226*tmp226;
REAL8 tmp834=(eta*eta*eta*eta);
REAL8 tmp848=pow(x->data[0],8.);
REAL8 tmp849=1./(tmp242*tmp242*tmp242*tmp242);
REAL8 tmp859=tmp320*tmp320;
REAL8 tmp790=12.*coeffs->d1v2*eta*sigmaKerr->data[2]*tmp56;
REAL8 tmp791=36.*eta*sigmaKerr->data[2]*tmp190*tmp243*tmp260*tmp268*tmp274*tmp561*tmp688;
REAL8 tmp792=30.*eta*sigmaStar->data[2]*tmp190*tmp243*tmp260*tmp268*tmp274*tmp561*tmp688;
REAL8 tmp793=-36.*eta*sigmaKerr->data[2]*tmp190*tmp217*tmp243*tmp260*tmp268*tmp274*tmp697;
REAL8 tmp794=-30.*eta*sigmaStar->data[2]*tmp190*tmp217*tmp243*tmp260*tmp268*tmp274*tmp697;
REAL8 tmp795=-72.*eta*sigmaKerr->data[2]*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp278*tmp279;
REAL8 tmp796=-60.*eta*sigmaStar->data[2]*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp278*tmp279;
REAL8 tmp797=144.*eta*sigmaKerr->data[2]*tmp190*tmp210*tmp217*tmp261*tmp268*tmp274*tmp313;
REAL8 tmp798=120.*eta*sigmaStar->data[2]*tmp190*tmp210*tmp217*tmp261*tmp268*tmp274*tmp313;
REAL8 tmp799=-216.*eta*sigmaKerr->data[2]*tmp158*tmp190*tmp210*tmp217*tmp243*tmp268*tmp274;
REAL8 tmp800=-180.*eta*sigmaStar->data[2]*tmp158*tmp190*tmp210*tmp217*tmp243*tmp268*tmp274;
REAL8 tmp801=-72.*eta*sigmaKerr->data[2]*tmp210*tmp217*tmp226*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp802=-60.*eta*sigmaStar->data[2]*tmp210*tmp217*tmp226*tmp243*tmp260*tmp268*tmp274;
REAL8 tmp803=144.*sigmaKerr->data[2]*tmp190*tmp222*tmp261*tmp268*tmp278*tmp279*tmp313;
REAL8 tmp804=120.*sigmaStar->data[2]*tmp190*tmp222*tmp261*tmp268*tmp278*tmp279*tmp313;
REAL8 tmp805=-216.*sigmaKerr->data[2]*tmp158*tmp190*tmp222*tmp243*tmp268*tmp278*tmp279;
REAL8 tmp806=-180.*sigmaStar->data[2]*tmp158*tmp190*tmp222*tmp243*tmp268*tmp278*tmp279;
REAL8 tmp807=-72.*sigmaKerr->data[2]*tmp222*tmp226*tmp243*tmp260*tmp268*tmp278*tmp279;
REAL8 tmp808=-60.*sigmaStar->data[2]*tmp222*tmp226*tmp243*tmp260*tmp268*tmp278*tmp279;
REAL8 tmp809=36.*sigmaKerr->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp736*tmp737;
REAL8 tmp810=30.*sigmaStar->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp736*tmp737;
REAL8 tmp811=-36.*sigmaKerr->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp279*tmp741;
REAL8 tmp812=-30.*sigmaStar->data[2]*tmp190*tmp222*tmp243*tmp260*tmp268*tmp279*tmp741;
REAL8 tmp813=-288.*sigmaKerr->data[2]*tmp190*tmp222*tmp268*tmp274*tmp345*tmp347;
REAL8 tmp814=-240.*sigmaStar->data[2]*tmp190*tmp222*tmp268*tmp274*tmp345*tmp347;
REAL8 tmp815=504.*sigmaKerr->data[2]*tmp190*tmp222*tmp260*tmp261*tmp268*tmp274;
REAL8 tmp816=420.*sigmaStar->data[2]*tmp190*tmp222*tmp260*tmp261*tmp268*tmp274;
REAL8 tmp817=144.*sigmaKerr->data[2]*tmp222*tmp226*tmp261*tmp268*tmp274*tmp313;
REAL8 tmp818=120.*sigmaStar->data[2]*tmp222*tmp226*tmp261*tmp268*tmp274*tmp313;
REAL8 tmp819=-216.*sigmaKerr->data[2]*tmp190*tmp222*tmp243*tmp268*tmp274*x->data[0];
REAL8 tmp820=-180.*sigmaStar->data[2]*tmp190*tmp222*tmp243*tmp268*tmp274*x->data[0];
REAL8 tmp821=-216.*sigmaKerr->data[2]*tmp158*tmp222*tmp226*tmp243*tmp268*tmp274;
REAL8 tmp822=-180.*sigmaStar->data[2]*tmp158*tmp222*tmp226*tmp243*tmp268*tmp274;
REAL8 tmp823=-36.*sigmaKerr->data[2]*tmp222*tmp243*tmp260*tmp268*tmp274*tmp705;
REAL8 tmp824=-30.*sigmaStar->data[2]*tmp222*tmp243*tmp260*tmp268*tmp274*tmp705;
REAL8 tmp826=3.*sigmaKerr->data[2]*tmp825*x->data[0];
REAL8 tmp827=4.*sigmaStar->data[2]*tmp825*x->data[0];
REAL8 tmp828=6.*sigmaKerr->data[2]*tmp320;
REAL8 tmp829=8.*sigmaStar->data[2]*tmp320;
REAL8 tmp830=tmp791+tmp792+tmp793+tmp794+tmp795+tmp796+tmp797+tmp798+tmp799+tmp800+tmp801+tmp802+tmp803+tmp804+tmp805+tmp806+tmp807+tmp808+tmp809+tmp810+tmp811+tmp812+tmp813+tmp814+tmp815+tmp816+tmp817+tmp818+tmp819+tmp820+tmp821+tmp822+tmp823+tmp824+tmp826+tmp827+tmp828+tmp829;
REAL8 tmp831=0.08333333333333333*eta*tmp60*tmp830;
REAL8 tmp832=-0.16666666666666666*eta*tmp167*tmp326;
REAL8 tmp833=0.16666666666666666*eta*tmp201*tmp334;
REAL8 tmp835=720.*tmp261*tmp337*tmp338*tmp339*tmp340*tmp561*tmp688*tmp834;
REAL8 tmp836=-720.*tmp222*tmp261*tmp336*tmp337*tmp338*tmp339*tmp340*tmp561*tmp688;
REAL8 tmp837=720.*tmp217*tmp222*tmp261*tmp336*tmp337*tmp338*tmp339*tmp340*tmp697;
REAL8 tmp838=2880.*tmp210*tmp217*tmp222*tmp261*tmp274*tmp278*tmp279*tmp336*tmp337*tmp338*tmp339;
REAL8 tmp839=-5760.*tmp210*tmp217*tmp222*tmp336*tmp337*tmp338*tmp340*tmp344*tmp345;
REAL8 tmp840=8640.*tmp210*tmp217*tmp222*tmp261*tmp336*tmp337*tmp338*tmp340*tmp347;
REAL8 tmp841=2880.*tmp190*tmp210*tmp217*tmp222*tmp226*tmp261*tmp336*tmp337*tmp339*tmp340;
REAL8 tmp842=720.*tmp261*tmp337*tmp338*tmp339*tmp342*tmp5*tmp736*tmp737;
REAL8 tmp843=-5760.*tmp274*tmp278*tmp279*tmp337*tmp338*tmp342*tmp344*tmp345*tmp5;
REAL8 tmp844=8640.*tmp261*tmp274*tmp278*tmp279*tmp337*tmp338*tmp342*tmp347*tmp5;
REAL8 tmp845=2880.*tmp190*tmp226*tmp261*tmp274*tmp278*tmp279*tmp337*tmp339*tmp342*tmp5;
REAL8 tmp846=-720.*tmp261*tmp274*tmp337*tmp338*tmp339*tmp342*tmp5*tmp736*tmp737;
REAL8 tmp847=720.*tmp261*tmp274*tmp279*tmp337*tmp338*tmp339*tmp342*tmp5*tmp741;
REAL8 tmp850=8640.*tmp337*tmp338*tmp340*tmp342*tmp5*tmp848*tmp849;
REAL8 tmp851=-18720.*tmp337*tmp338*tmp339*tmp340*tmp342*tmp345*tmp5;
REAL8 tmp852=-5760.*tmp190*tmp226*tmp337*tmp340*tmp342*tmp344*tmp345*tmp5;
REAL8 tmp853=10800.*tmp261*tmp313*tmp337*tmp338*tmp340*tmp342*tmp5;
REAL8 tmp854=8640.*tmp190*tmp226*tmp261*tmp337*tmp340*tmp342*tmp347*tmp5;
REAL8 tmp855=720.*tmp261*tmp337*tmp339*tmp340*tmp342*tmp5*tmp780;
REAL8 tmp856=720.*tmp190*tmp261*tmp337*tmp339*tmp340*tmp342*tmp5*tmp705;
REAL8 tmp857=-2.*tmp352*tmp825*x->data[0];
REAL8 tmp858=-4.*tmp320*tmp352;
REAL8 tmp860=-2.*eta*tmp158*tmp356*tmp859;
REAL8 tmp861=-2.*eta*tmp158*tmp323*tmp356*tmp825;
REAL8 tmp862=-8.*eta*tmp320*tmp323*tmp356*x->data[0];
REAL8 tmp863=-2.*eta*tmp356*tmp358;
REAL8 tmp864=tmp362*tmp825*x->data[0];
REAL8 tmp865=2.*tmp320*tmp362;
REAL8 tmp866=tmp864+tmp865;
REAL8 tmp867=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274*tmp866;
REAL8 tmp868=-12.*eta*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp274*tmp365;
REAL8 tmp869=-12.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp278*tmp279*tmp365;
REAL8 tmp870=24.*tmp190*tmp222*tmp261*tmp268*tmp274*tmp313*tmp365;
REAL8 tmp871=-36.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274*tmp365;
REAL8 tmp872=-12.*tmp222*tmp226*tmp243*tmp260*tmp268*tmp274*tmp365;
REAL8 tmp873=6.*eta*tmp190*tmp243*tmp260*tmp268*tmp274*tmp370*tmp561*tmp688;
REAL8 tmp874=-6.*eta*tmp190*tmp217*tmp243*tmp260*tmp268*tmp274*tmp370*tmp697;
REAL8 tmp875=-12.*eta*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp278*tmp279*tmp370;
REAL8 tmp876=24.*eta*tmp190*tmp210*tmp217*tmp261*tmp268*tmp274*tmp313*tmp370;
REAL8 tmp877=-36.*eta*tmp158*tmp190*tmp210*tmp217*tmp243*tmp268*tmp274*tmp370;
REAL8 tmp878=-12.*eta*tmp210*tmp217*tmp226*tmp243*tmp260*tmp268*tmp274*tmp370;
REAL8 tmp879=24.*tmp190*tmp222*tmp261*tmp268*tmp278*tmp279*tmp313*tmp370;
REAL8 tmp880=-36.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp278*tmp279*tmp370;
REAL8 tmp881=-12.*tmp222*tmp226*tmp243*tmp260*tmp268*tmp278*tmp279*tmp370;
REAL8 tmp882=6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp370*tmp736*tmp737;
REAL8 tmp883=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp279*tmp370*tmp741;
REAL8 tmp884=-48.*tmp190*tmp222*tmp268*tmp274*tmp345*tmp347*tmp370;
REAL8 tmp885=84.*tmp190*tmp222*tmp260*tmp261*tmp268*tmp274*tmp370;
REAL8 tmp886=24.*tmp222*tmp226*tmp261*tmp268*tmp274*tmp313*tmp370;
REAL8 tmp887=-36.*tmp190*tmp222*tmp243*tmp268*tmp274*tmp370*x->data[0];
REAL8 tmp888=-36.*tmp158*tmp222*tmp226*tmp243*tmp268*tmp274*tmp370;
REAL8 tmp889=-6.*tmp222*tmp243*tmp260*tmp268*tmp274*tmp370*tmp705;
REAL8 tmp890=tmp835+tmp836+tmp837+tmp838+tmp839+tmp840+tmp841+tmp842+tmp843+tmp844+tmp845+tmp846+tmp847+tmp850+tmp851+tmp852+tmp853+tmp854+tmp855+tmp856+tmp857+tmp858+tmp860+tmp861+tmp862+tmp863+tmp867+tmp868+tmp869+tmp870+tmp871+tmp872+tmp873+tmp874+tmp875+tmp876+tmp877+tmp878+tmp879+tmp880+tmp881+tmp882+tmp883+tmp884+tmp885+tmp886+tmp887+tmp888+tmp889;
REAL8 tmp891=0.013888888888888888*sigmaStar->data[2]*tmp167*tmp890;
REAL8 tmp892=-0.05555555555555555*sigmaStar->data[2]*tmp201*tmp376;
REAL8 tmp893=0.08333333333333333*sigmaStar->data[2]*tmp196*tmp385;
REAL8 tmp894=1620.*tmp261*tmp337*tmp338*tmp339*tmp340*tmp561*tmp688*tmp834;
REAL8 tmp895=-1620.*tmp222*tmp261*tmp336*tmp337*tmp338*tmp339*tmp340*tmp561*tmp688;
REAL8 tmp896=1620.*tmp217*tmp222*tmp261*tmp336*tmp337*tmp338*tmp339*tmp340*tmp697;
REAL8 tmp897=6480.*tmp210*tmp217*tmp222*tmp261*tmp274*tmp278*tmp279*tmp336*tmp337*tmp338*tmp339;
REAL8 tmp898=-12960.*tmp210*tmp217*tmp222*tmp336*tmp337*tmp338*tmp340*tmp344*tmp345;
REAL8 tmp899=19440.*tmp210*tmp217*tmp222*tmp261*tmp336*tmp337*tmp338*tmp340*tmp347;
REAL8 tmp900=6480.*tmp190*tmp210*tmp217*tmp222*tmp226*tmp261*tmp336*tmp337*tmp339*tmp340;
REAL8 tmp901=1620.*tmp261*tmp337*tmp338*tmp339*tmp342*tmp5*tmp736*tmp737;
REAL8 tmp902=-12960.*tmp274*tmp278*tmp279*tmp337*tmp338*tmp342*tmp344*tmp345*tmp5;
REAL8 tmp903=19440.*tmp261*tmp274*tmp278*tmp279*tmp337*tmp338*tmp342*tmp347*tmp5;
REAL8 tmp904=6480.*tmp190*tmp226*tmp261*tmp274*tmp278*tmp279*tmp337*tmp339*tmp342*tmp5;
REAL8 tmp905=-1620.*tmp261*tmp274*tmp337*tmp338*tmp339*tmp342*tmp5*tmp736*tmp737;
REAL8 tmp906=1620.*tmp261*tmp274*tmp279*tmp337*tmp338*tmp339*tmp342*tmp5*tmp741;
REAL8 tmp907=19440.*tmp337*tmp338*tmp340*tmp342*tmp5*tmp848*tmp849;
REAL8 tmp908=-42120.*tmp337*tmp338*tmp339*tmp340*tmp342*tmp345*tmp5;
REAL8 tmp909=-12960.*tmp190*tmp226*tmp337*tmp340*tmp342*tmp344*tmp345*tmp5;
REAL8 tmp910=24300.*tmp261*tmp313*tmp337*tmp338*tmp340*tmp342*tmp5;
REAL8 tmp911=19440.*tmp190*tmp226*tmp261*tmp337*tmp340*tmp342*tmp347*tmp5;
REAL8 tmp912=1620.*tmp261*tmp337*tmp339*tmp340*tmp342*tmp5*tmp780;
REAL8 tmp913=1620.*tmp190*tmp261*tmp337*tmp339*tmp340*tmp342*tmp5*tmp705;
REAL8 tmp914=4.*tmp394*tmp825*x->data[0];
REAL8 tmp915=8.*tmp320*tmp394;
REAL8 tmp916=-90.*eta*tmp158*tmp859;
REAL8 tmp917=-90.*eta*tmp158*tmp323*tmp825;
REAL8 tmp918=-360.*eta*tmp320*tmp323*x->data[0];
REAL8 tmp919=-90.*eta*tmp358;
REAL8 tmp920=tmp401*tmp825*x->data[0];
REAL8 tmp921=2.*tmp320*tmp401;
REAL8 tmp922=tmp920+tmp921;
REAL8 tmp923=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp274*tmp922;
REAL8 tmp924=-12.*eta*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp274*tmp404;
REAL8 tmp925=-12.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp278*tmp279*tmp404;
REAL8 tmp926=24.*tmp190*tmp222*tmp261*tmp268*tmp274*tmp313*tmp404;
REAL8 tmp927=-36.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp274*tmp404;
REAL8 tmp928=-12.*tmp222*tmp226*tmp243*tmp260*tmp268*tmp274*tmp404;
REAL8 tmp929=6.*eta*tmp190*tmp243*tmp260*tmp268*tmp274*tmp409*tmp561*tmp688;
REAL8 tmp930=-6.*eta*tmp190*tmp217*tmp243*tmp260*tmp268*tmp274*tmp409*tmp697;
REAL8 tmp931=-12.*eta*tmp190*tmp210*tmp217*tmp243*tmp260*tmp268*tmp278*tmp279*tmp409;
REAL8 tmp932=24.*eta*tmp190*tmp210*tmp217*tmp261*tmp268*tmp274*tmp313*tmp409;
REAL8 tmp933=-36.*eta*tmp158*tmp190*tmp210*tmp217*tmp243*tmp268*tmp274*tmp409;
REAL8 tmp934=-12.*eta*tmp210*tmp217*tmp226*tmp243*tmp260*tmp268*tmp274*tmp409;
REAL8 tmp935=24.*tmp190*tmp222*tmp261*tmp268*tmp278*tmp279*tmp313*tmp409;
REAL8 tmp936=-36.*tmp158*tmp190*tmp222*tmp243*tmp268*tmp278*tmp279*tmp409;
REAL8 tmp937=-12.*tmp222*tmp226*tmp243*tmp260*tmp268*tmp278*tmp279*tmp409;
REAL8 tmp938=6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp409*tmp736*tmp737;
REAL8 tmp939=-6.*tmp190*tmp222*tmp243*tmp260*tmp268*tmp279*tmp409*tmp741;
REAL8 tmp940=-48.*tmp190*tmp222*tmp268*tmp274*tmp345*tmp347*tmp409;
REAL8 tmp941=84.*tmp190*tmp222*tmp260*tmp261*tmp268*tmp274*tmp409;
REAL8 tmp942=24.*tmp222*tmp226*tmp261*tmp268*tmp274*tmp313*tmp409;
REAL8 tmp943=-36.*tmp190*tmp222*tmp243*tmp268*tmp274*tmp409*x->data[0];
REAL8 tmp944=-36.*tmp158*tmp222*tmp226*tmp243*tmp268*tmp274*tmp409;
REAL8 tmp945=-6.*tmp222*tmp243*tmp260*tmp268*tmp274*tmp409*tmp705;
REAL8 tmp946=tmp894+tmp895+tmp896+tmp897+tmp898+tmp899+tmp900+tmp901+tmp902+tmp903+tmp904+tmp905+tmp906+tmp907+tmp908+tmp909+tmp910+tmp911+tmp912+tmp913+tmp914+tmp915+tmp916+tmp917+tmp918+tmp919+tmp923+tmp924+tmp925+tmp926+tmp927+tmp928+tmp929+tmp930+tmp931+tmp932+tmp933+tmp934+tmp935+tmp936+tmp937+tmp938+tmp939+tmp940+tmp941+tmp942+tmp943+tmp944+tmp945;
REAL8 tmp947=0.006944444444444444*sigmaKerr->data[2]*tmp167*tmp946;
REAL8 tmp948=-0.027777777777777776*sigmaKerr->data[2]*tmp201*tmp415;
REAL8 tmp949=0.041666666666666664*sigmaKerr->data[2]*tmp196*tmp423;
REAL8 tmp950=tmp790+tmp831+tmp832+tmp833+tmp891+tmp892+tmp893+tmp947+tmp948+tmp949;
REAL8 tmp460=-0.5*tmp455*tmp459;
REAL8 tmp461=1.*tmp432*tmp451*x->data[0];
REAL8 tmp469=0.5*tmp427*tmp462*tmp468;
REAL8 tmp470=tmp460+tmp461+tmp469;
REAL8 tmp579=(1.0/sqrt(tmp428*tmp428*tmp428));
REAL8 tmp588=(1.0/sqrt(tmp431*tmp431*tmp431));
REAL8 tmp515=-(tmp158*tmp190*tmp222*tmp254*tmp268*tmp274*tmp425);
REAL8 tmp517=tmp242*tmp254*tmp425*tmp516;
REAL8 tmp518=-(eta*tmp158*tmp190*tmp210*tmp217*tmp268*tmp274*tmp508);
REAL8 tmp519=-(tmp158*tmp190*tmp222*tmp268*tmp278*tmp279*tmp508);
REAL8 tmp520=-2.*tmp190*tmp222*tmp268*tmp274*tmp508*x->data[0];
REAL8 tmp521=-(tmp158*tmp222*tmp226*tmp268*tmp274*tmp508);
REAL8 tmp524=tmp242*tmp508*tmp523;
REAL8 tmp525=2.*tmp508*tmp516*x->data[0];
REAL8 tmp526=tmp515+tmp517+tmp518+tmp519+tmp520+tmp521+tmp524+tmp525;
REAL8 tmp536=1/tmp435;
REAL8 tmp537=(1.0/(tmp503*tmp503));
REAL8 tmp506=1.*tmp190*tmp222*tmp237*tmp254*tmp284*tmp313*tmp425*tmp505;
REAL8 tmp509=-(tmp190*tmp222*tmp228*tmp232*tmp284*tmp313*tmp505*tmp508);
REAL8 tmp510=1.*eta*tmp190*tmp210*tmp217*tmp237*tmp284*tmp313*tmp505*tmp508;
REAL8 tmp511=4.*tmp190*tmp222*tmp237*tmp242*tmp284*tmp347*tmp508;
REAL8 tmp512=-2.*p->data[2]*tmp145*tmp158*tmp190*tmp222*tmp237*tmp505*tmp508*tmp59;
REAL8 tmp513=4.*tmp190*tmp222*tmp237*tmp260*tmp284*tmp505*tmp508;
REAL8 tmp514=1.*tmp222*tmp226*tmp237*tmp284*tmp313*tmp505*tmp508;
REAL8 tmp527=1.*tmp158*tmp190*tmp222*tmp526;
REAL8 tmp531=1.*eta*tmp158*tmp190*tmp210*tmp217*tmp530;
REAL8 tmp532=2.*tmp190*tmp222*tmp530*x->data[0];
REAL8 tmp533=1.*tmp158*tmp222*tmp226*tmp530;
REAL8 tmp534=tmp506+tmp509+tmp510+tmp511+tmp512+tmp513+tmp514+tmp527+tmp531+tmp532+tmp533;
REAL8 tmp542=tmp228*tmp478;
REAL8 tmp543=0.;
REAL8 tmp544=-8.*tmp158;
REAL8 tmp545=-4.*tmp159;
REAL8 tmp568=1.*tmp149*tmp567;
REAL8 tmp569=tmp544+tmp545+tmp568;
REAL8 tmp570=tmp157*tmp569;
REAL8 tmp571=tmp236*tmp499;
REAL8 tmp572=tmp542+tmp543+tmp570+tmp571;
REAL8 tmp577=pow(tmp242,-2.5);
REAL8 tmp969=pow(tmp435,-2.5);
REAL8 tmp989=-(eta*tmp158*tmp190*tmp561*tmp688);
REAL8 tmp990=eta*tmp158*tmp190*tmp217*tmp697;
REAL8 tmp991=4.*eta*tmp190*tmp210*tmp217*x->data[0];
REAL8 tmp992=2.*eta*tmp158*tmp210*tmp217*tmp226;
REAL8 tmp993=4.*tmp222*tmp226*x->data[0];
REAL8 tmp994=tmp158*tmp222*tmp705;
REAL8 tmp995=tmp564+tmp989+tmp990+tmp991+tmp992+tmp993+tmp994;
REAL8 tmp997=tmp459*tmp459;
REAL8 tmp1001=2.*tmp158*tmp190*tmp222*tmp242*tmp575*tmp684;
REAL8 tmp1002=-(tmp158*tmp190*tmp222*tmp232*tmp242*tmp707);
REAL8 tmp1003=-2.*eta*tmp158*tmp190*tmp210*tmp217*tmp228*tmp232*tmp242;
REAL8 tmp1004=-4.*tmp190*tmp222*tmp228*tmp232*tmp260;
REAL8 tmp1005=-4.*tmp190*tmp222*tmp228*tmp232*tmp242*x->data[0];
REAL8 tmp1006=-2.*tmp158*tmp222*tmp226*tmp228*tmp232*tmp242;
REAL8 tmp1007=-(eta*tmp158*tmp190*tmp237*tmp242*tmp561*tmp688);
REAL8 tmp1008=eta*tmp158*tmp190*tmp217*tmp237*tmp242*tmp697;
REAL8 tmp1009=4.*eta*tmp190*tmp210*tmp217*tmp237*tmp260;
REAL8 tmp1010=4.*eta*tmp190*tmp210*tmp217*tmp237*tmp242*x->data[0];
REAL8 tmp1011=2.*eta*tmp158*tmp210*tmp217*tmp226*tmp237*tmp242;
REAL8 tmp1012=10.*tmp158*tmp190*tmp222*tmp237;
REAL8 tmp1013=4.*tmp222*tmp226*tmp237*tmp260;
REAL8 tmp1014=2.*tmp190*tmp222*tmp237*tmp242;
REAL8 tmp1015=4.*tmp222*tmp226*tmp237*tmp242*x->data[0];
REAL8 tmp1016=tmp158*tmp222*tmp237*tmp242*tmp705;
REAL8 tmp1017=tmp1001+tmp1002+tmp1003+tmp1004+tmp1005+tmp1006+tmp1007+tmp1008+tmp1009+tmp1010+tmp1011+tmp1012+tmp1013+tmp1014+tmp1015+tmp1016;
REAL8 tmp1020=tmp468*tmp468;
REAL8 tmp1075=-0.25*tmp448*tmp859;
REAL8 tmp1076=0.5*tmp436*tmp825;
REAL8 tmp1077=tmp1075+tmp1076+tmp719+tmp720+tmp721+tmp722+tmp725+tmp726+tmp727+tmp728+tmp729+tmp730+tmp731+tmp732+tmp733+tmp734+tmp735+tmp738+tmp742+tmp743+tmp744+tmp745+tmp746+tmp747+tmp748+tmp749+tmp750+tmp751+tmp752+tmp753+tmp754+tmp755+tmp756+tmp759+tmp760;
REAL8 tmp590=(1.0/(tmp516*tmp516));
REAL8 tmp768=1./(tmp222*tmp222*tmp222);
REAL8 tmp781=1./(tmp190*tmp190*tmp190);
REAL8 tmp618=0.5*tmp581*tmp586*tmp608;
REAL8 tmp619=-0.5*tmp145*tmp254*tmp425*tmp427*tmp432*tmp503*tmp592*tmp595*x->data[0];
REAL8 tmp620=-0.25*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp592*tmp595*x->data[0];
REAL8 tmp621=-0.5*tmp145*tmp158*tmp432*tmp451*tmp503*tmp508*tmp592*tmp595;
REAL8 tmp622=-0.5*tmp145*tmp427*tmp432*tmp503*tmp508*tmp592*tmp595;
REAL8 tmp623=0.5*p->data[2]*tmp427*tmp432*tmp503*tmp508*tmp59*tmp592*tmp595*tmp60;
REAL8 tmp624=(1.0/sqrt(tmp591*tmp591*tmp591));
REAL8 tmp625=2.*eta*tmp210*tmp217*tmp222*tmp274*tmp313*tmp338;
REAL8 tmp626=tmp278*tmp279*tmp313*tmp338*tmp342;
REAL8 tmp627=4.*tmp260*tmp274*tmp338*tmp342;
REAL8 tmp628=2.*tmp190*tmp226*tmp274*tmp313*tmp342;
REAL8 tmp629=tmp625+tmp626+tmp627+tmp628;
REAL8 tmp630=0.25*tmp145*tmp427*tmp432*tmp503*tmp508*tmp595*tmp624*tmp629*x->data[0];
REAL8 tmp631=-0.25*tmp145*tmp427*tmp462*tmp468*tmp503*tmp508*tmp592*tmp595*x->data[0];
REAL8 tmp632=-2.*eta*tmp158*tmp190*tmp210*tmp217;
REAL8 tmp633=-4.*tmp190*tmp222*x->data[0];
REAL8 tmp634=-2.*tmp158*tmp222*tmp226;
REAL8 tmp635=tmp477*tmp567;
REAL8 tmp636=0.5*tmp497*tmp581*tmp586;
REAL8 tmp637=tmp632+tmp633+tmp634+tmp635+tmp636;
REAL8 tmp638=-0.5*tmp145*tmp427*tmp432*tmp503*tmp508*tmp592*tmp637*x->data[0];
REAL8 tmp639=0.5*tmp455*tmp459*tmp606;
REAL8 tmp640=1.*tmp145*tmp254*tmp425*tmp427*tmp432*tmp602*tmp604*x->data[0];
REAL8 tmp641=1.*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp602*x->data[0];
REAL8 tmp642=-2.*tmp158*tmp261;
REAL8 tmp643=1.*tmp243;
REAL8 tmp644=-4.*eta*tmp190*tmp210*tmp217*tmp260;
REAL8 tmp645=-12.*tmp158*tmp190*tmp222;
REAL8 tmp646=-4.*tmp222*tmp226*tmp260;
REAL8 tmp647=tmp159*tmp567;
REAL8 tmp648=2.*tmp497*x->data[0];
REAL8 tmp649=tmp644+tmp645+tmp646+tmp647+tmp648;
REAL8 tmp650=0.5*tmp159*tmp167*tmp237*tmp241*tmp244*tmp649;
REAL8 tmp651=-0.5*tmp159*tmp167*tmp228*tmp232*tmp241*tmp244*tmp600;
REAL8 tmp652=-0.5*eta*tmp159*tmp167*tmp210*tmp217*tmp237*tmp241*tmp300*tmp600;
REAL8 tmp653=1.*tmp237*tmp241*tmp244*tmp60*tmp600;
REAL8 tmp654=-(tmp159*tmp201*tmp237*tmp241*tmp244*tmp600);
REAL8 tmp655=-0.5*tmp159*tmp167*tmp226*tmp237*tmp244*tmp304*tmp600;
REAL8 tmp656=tmp642+tmp643+tmp650+tmp651+tmp652+tmp653+tmp654+tmp655;
REAL8 tmp657=1.*tmp145*tmp427*tmp432*tmp508*tmp604*tmp656*x->data[0];
REAL8 tmp658=1.*tmp145*tmp158*tmp432*tmp451*tmp508*tmp602*tmp604;
REAL8 tmp659=1.*tmp145*tmp427*tmp432*tmp508*tmp602*tmp604;
REAL8 tmp660=-(p->data[2]*tmp427*tmp432*tmp508*tmp59*tmp60*tmp602*tmp604);
REAL8 tmp661=0.5*tmp145*tmp427*tmp462*tmp468*tmp508*tmp602*tmp604*x->data[0];
REAL8 tmp662=tmp640+tmp641+tmp657+tmp658+tmp659+tmp660+tmp661;
REAL8 tmp663=tmp429*tmp662;
REAL8 tmp664=tmp619+tmp620+tmp621+tmp622+tmp623+tmp630+tmp631+tmp638+tmp639+tmp663;
REAL8 tmp665=tmp477*tmp664;
REAL8 tmp666=tmp618+tmp665;
REAL8 tmp1183=-(eta*tmp158*tmp190*tmp274*tmp561*tmp688);
REAL8 tmp1184=eta*tmp158*tmp190*tmp217*tmp274*tmp697;
REAL8 tmp1185=2.*eta*tmp158*tmp190*tmp210*tmp217*tmp278*tmp279;
REAL8 tmp1186=4.*eta*tmp190*tmp210*tmp217*tmp274*x->data[0];
REAL8 tmp1187=2.*eta*tmp158*tmp210*tmp217*tmp226*tmp274;
REAL8 tmp1188=4.*tmp190*tmp222*tmp278*tmp279*x->data[0];
REAL8 tmp1189=2.*tmp158*tmp222*tmp226*tmp278*tmp279;
REAL8 tmp1190=-(tmp158*tmp190*tmp222*tmp736*tmp737);
REAL8 tmp1191=tmp158*tmp190*tmp222*tmp279*tmp741;
REAL8 tmp1192=2.*tmp190*tmp222*tmp274;
REAL8 tmp1193=4.*tmp222*tmp226*tmp274*x->data[0];
REAL8 tmp1194=tmp158*tmp222*tmp274*tmp705;
REAL8 tmp1195=tmp1183+tmp1184+tmp1185+tmp1186+tmp1187+tmp1188+tmp1189+tmp1190+tmp1191+tmp1192+tmp1193+tmp1194;
REAL8 tmp1200=(1.0/sqrt(tmp476*tmp476*tmp476));
REAL8 tmp1201=tmp586*tmp586;
REAL8 tmp1124=2.*eta*tmp495*tmp561*tmp688;
REAL8 tmp1125=-2.*eta*tmp190*tmp561*tmp688*x->data[0];
REAL8 tmp1126=-2.*eta*tmp217*tmp495*tmp697;
REAL8 tmp1127=2.*eta*tmp190*tmp217*tmp697*x->data[0];
REAL8 tmp1128=4.*eta*tmp190*tmp210*tmp217;
REAL8 tmp1129=4.*eta*tmp149*tmp167*tmp210*tmp217;
REAL8 tmp1130=4.*eta*tmp210*tmp217*tmp226*x->data[0];
REAL8 tmp1131=15.*coeffs->k5l*tmp201;
REAL8 tmp1132=10.*tmp201*tmp208;
REAL8 tmp1133=tmp1131+tmp1132;
REAL8 tmp1134=1.*tmp1133*tmp60;
REAL8 tmp1135=-2.*tmp167*tmp550;
REAL8 tmp1136=2.*tmp201*tmp485;
REAL8 tmp1137=tmp1134+tmp1135+tmp1136;
REAL8 tmp1138=1.*tmp1137*tmp60;
REAL8 tmp1139=-2.*tmp167*tmp553;
REAL8 tmp1140=2.*tmp201*tmp487;
REAL8 tmp1141=tmp1138+tmp1139+tmp1140;
REAL8 tmp1142=1.*tmp1141*tmp60;
REAL8 tmp1143=-2.*tmp167*tmp556;
REAL8 tmp1144=2.*tmp201*tmp489;
REAL8 tmp1145=tmp1142+tmp1143+tmp1144;
REAL8 tmp1146=-(eta*tmp1145*tmp190*tmp217);
REAL8 tmp1147=2.*eta*tmp190*tmp210*tmp559*tmp561;
REAL8 tmp1148=-2.*eta*tmp217*tmp226*tmp559;
REAL8 tmp1149=1./(tmp216*tmp216*tmp216);
REAL8 tmp1150=-2.*eta*tmp1149*tmp190*tmp491*tmp688;
REAL8 tmp1151=eta*tmp190*tmp491*tmp561*tmp697;
REAL8 tmp1152=2.*eta*tmp210*tmp226*tmp491*tmp561;
REAL8 tmp1153=-(eta*tmp217*tmp491*tmp705);
REAL8 tmp1154=4.*tmp222*tmp226;
REAL8 tmp1155=-4.*tmp149*tmp201*tmp222;
REAL8 tmp1156=2.*tmp222*tmp705*x->data[0];
REAL8 tmp1157=tmp1124+tmp1125+tmp1126+tmp1127+tmp1128+tmp1129+tmp1130+tmp1146+tmp1147+tmp1148+tmp1150+tmp1151+tmp1152+tmp1153+tmp1154+tmp1155+tmp1156;
REAL8 tmp670=1.*coeffs->dheffSSv2*eta*tmp196*tmp54;
REAL8 tmp671=1.*tmp145*tmp157*tmp237*x->data[0];
REAL8 tmp672=tmp246*tmp298;
REAL8 tmp673=1.*tmp157*tmp237*tmp445;
REAL8 tmp674=1.*tmp145*tmp237*tmp427*tmp434*tmp436*tmp445*x->data[0];
REAL8 tmp675=-0.5*tmp201*tmp473;
REAL8 tmp676=0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540;
REAL8 tmp677=1.*tmp167*tmp241*tmp243*tmp244*tmp432*tmp610*tmp612;
REAL8 tmp678=-1.+tmp670+tmp671+tmp672+tmp673+tmp674+tmp675+tmp676+tmp677;
REAL8 tmp679=2.*eta*tmp678;
REAL8 tmp680=1.+tmp679;
REAL8 d2Hdx0dx0=-(tmp5*(1.0/sqrt(tmp680*tmp680*tmp680))*pow(1.*tmp145*tmp157*tmp237+0.5*tmp246*tmp289*tmp295-0.5*tmp297*tmp298*tmp306+1.*e3z*tmp157*tmp237*tmp425-tmp201*tmp425*tmp442-tmp157*tmp228*tmp232*tmp445+1.*tmp236*tmp237*tmp445+1.*tmp145*tmp237*tmp427*tmp434*tmp436*tmp445+1.*tmp145*tmp158*tmp237*tmp434*tmp436*tmp445*tmp451+1.5*tmp196*tmp473+0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp534-0.25*tmp232*tmp320*tmp448*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540-0.25*tmp232*tmp320*tmp455*tmp462*tmp475*tmp477*tmp501*tmp536*tmp537*tmp540-4.*coeffs->dheffSSv2*eta*tmp54*tmp56+0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp477*tmp504*tmp540*tmp572-tmp228*tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp575-0.25*tmp232*tmp436*tmp459*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp579+0.25*tmp232*tmp436*tmp455*tmp462*tmp475*tmp501*tmp504*tmp540*tmp581*tmp586-0.25*tmp232*tmp436*tmp455*tmp468*tmp475*tmp477*tmp501*tmp504*tmp540*tmp588-p->data[2]*tmp157*tmp237*tmp59*tmp60-p->data[2]*tmp237*tmp427*tmp434*tmp436*tmp445*tmp59*tmp60-tmp167*tmp241*tmp243*tmp244*tmp432*tmp523*tmp590*tmp610-2.*tmp201*tmp241*tmp243*tmp244*tmp432*tmp610*tmp612-eta*tmp167*tmp210*tmp217*tmp241*tmp243*tmp300*tmp432*tmp610*tmp612-tmp167*tmp226*tmp243*tmp244*tmp304*tmp432*tmp610*tmp612+0.5*tmp167*tmp241*tmp243*tmp244*tmp462*tmp468*tmp610*tmp612-2.*tmp241*tmp244*tmp261*tmp432*tmp60*tmp610*tmp612+1.*tmp167*tmp241*tmp243*tmp244*tmp432*tmp612*tmp666-tmp145*tmp157*tmp228*tmp232*x->data[0]+1.*tmp145*tmp236*tmp237*x->data[0]+1.*e3z*tmp145*tmp237*tmp425*tmp427*tmp434*tmp436*x->data[0]-tmp145*tmp228*tmp232*tmp427*tmp434*tmp436*tmp445*x->data[0]-0.5*tmp145*tmp237*tmp320*tmp427*tmp434*tmp445*tmp448*x->data[0]+1.*tmp145*tmp237*tmp427*tmp436*tmp445*tmp470*x->data[0]-1.5*tmp232*tmp436*tmp455*tmp462*tmp477*tmp501*tmp504*tmp540*tmp577*x->data[0],2.))+(1.*eta*(0.-2.*tmp145*tmp157*tmp228*tmp232+2.*tmp145*tmp236*tmp237-0.5*tmp289*tmp295*tmp297*tmp306-2.*e3z*tmp157*tmp228*tmp232*tmp425+2.*e3z*tmp236*tmp237*tmp425+2.*e3z*tmp145*tmp237*tmp425*tmp427*tmp434*tmp436+6.*tmp196*tmp425*tmp442-2.*tmp228*tmp232*tmp236*tmp445-2.*tmp145*tmp228*tmp232*tmp427*tmp434*tmp436*tmp445-tmp145*tmp237*tmp320*tmp427*tmp434*tmp445*tmp448+2.*e3z*tmp145*tmp158*tmp237*tmp425*tmp434*tmp436*tmp451-2.*tmp145*tmp158*tmp228*tmp232*tmp434*tmp436*tmp445*tmp451-tmp145*tmp158*tmp237*tmp320*tmp434*tmp445*tmp448*tmp451+2.*tmp145*tmp237*tmp427*tmp436*tmp445*tmp470+2.*tmp145*tmp158*tmp237*tmp436*tmp445*tmp451*tmp470-tmp145*tmp237*tmp260*tmp434*tmp436*tmp445*tmp475-0.5*tmp232*tmp320*tmp448*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp534-0.5*tmp232*tmp320*tmp455*tmp462*tmp475*tmp477*tmp501*tmp534*tmp536*tmp537+20.*coeffs->dheffSSv2*eta*tmp191*tmp54-0.125*tmp1200*tmp1201*tmp232*tmp436*tmp455*tmp462*tmp475*tmp501*tmp504*tmp540-6.*tmp473*tmp56+1.*tmp232*tmp436*tmp455*tmp462*tmp475*tmp477*tmp504*tmp534*tmp572-0.5*tmp232*tmp320*tmp448*tmp455*tmp462*tmp475*tmp477*tmp504*tmp540*tmp572-0.5*tmp232*tmp320*tmp455*tmp462*tmp475*tmp477*tmp536*tmp537*tmp540*tmp572-2.*tmp228*tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp534*tmp575+1.*tmp228*tmp320*tmp448*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp575+1.*tmp228*tmp320*tmp455*tmp462*tmp475*tmp477*tmp501*tmp536*tmp537*tmp540*tmp575-2.*tmp228*tmp436*tmp455*tmp462*tmp475*tmp477*tmp504*tmp540*tmp572*tmp575-1.5*tmp232*tmp436*tmp455*tmp462*tmp477*tmp501*tmp504*tmp540*tmp577-0.5*tmp232*tmp436*tmp459*tmp462*tmp475*tmp477*tmp501*tmp504*tmp534*tmp579+0.25*tmp232*tmp320*tmp448*tmp459*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp579+0.25*tmp232*tmp320*tmp459*tmp462*tmp475*tmp477*tmp501*tmp536*tmp537*tmp540*tmp579-0.5*tmp232*tmp436*tmp459*tmp462*tmp475*tmp477*tmp504*tmp540*tmp572*tmp579+1.*tmp228*tmp436*tmp459*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp575*tmp579+0.25*tmp1195*tmp232*tmp436*tmp455*tmp462*tmp475*tmp501*tmp504*tmp540*tmp581+0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp501*tmp504*tmp534*tmp581*tmp586-0.25*tmp232*tmp320*tmp448*tmp455*tmp462*tmp475*tmp501*tmp504*tmp540*tmp581*tmp586-0.25*tmp232*tmp320*tmp455*tmp462*tmp475*tmp501*tmp536*tmp537*tmp540*tmp581*tmp586+0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp504*tmp540*tmp572*tmp581*tmp586-tmp228*tmp436*tmp455*tmp462*tmp475*tmp501*tmp504*tmp540*tmp575*tmp581*tmp586-0.25*tmp232*tmp436*tmp459*tmp462*tmp475*tmp501*tmp504*tmp540*tmp579*tmp581*tmp586-0.5*tmp232*tmp436*tmp455*tmp468*tmp475*tmp477*tmp501*tmp504*tmp534*tmp588-0.25*tmp1017*tmp232*tmp436*tmp455*tmp475*tmp477*tmp501*tmp504*tmp540*tmp588+0.25*tmp232*tmp320*tmp448*tmp455*tmp468*tmp475*tmp477*tmp501*tmp504*tmp540*tmp588+0.25*tmp232*tmp320*tmp455*tmp468*tmp475*tmp477*tmp501*tmp536*tmp537*tmp540*tmp588-0.5*tmp232*tmp436*tmp455*tmp468*tmp475*tmp477*tmp504*tmp540*tmp572*tmp588+1.*tmp228*tmp436*tmp455*tmp468*tmp475*tmp477*tmp501*tmp504*tmp540*tmp575*tmp588+0.25*tmp232*tmp436*tmp459*tmp468*tmp475*tmp477*tmp501*tmp504*tmp540*tmp579*tmp588-0.25*tmp232*tmp436*tmp455*tmp468*tmp475*tmp501*tmp504*tmp540*tmp581*tmp586*tmp588-2.*p->data[2]*tmp237*tmp434*tmp436*tmp445*tmp451*tmp59+2.*p->data[2]*tmp157*tmp228*tmp232*tmp59*tmp60-2.*p->data[2]*tmp236*tmp237*tmp59*tmp60-2.*e3z*p->data[2]*tmp237*tmp425*tmp427*tmp434*tmp436*tmp59*tmp60+2.*p->data[2]*tmp228*tmp232*tmp427*tmp434*tmp436*tmp445*tmp59*tmp60+1.*p->data[2]*tmp237*tmp320*tmp427*tmp434*tmp445*tmp448*tmp59*tmp60-2.*p->data[2]*tmp237*tmp427*tmp436*tmp445*tmp470*tmp59*tmp60-tmp1077*tmp167*tmp241*tmp243*tmp244*tmp432*tmp590*tmp610+4.*tmp201*tmp241*tmp243*tmp244*tmp432*tmp523*tmp590*tmp610+2.*eta*tmp167*tmp210*tmp217*tmp241*tmp243*tmp300*tmp432*tmp523*tmp590*tmp610+2.*tmp167*tmp226*tmp243*tmp244*tmp304*tmp432*tmp523*tmp590*tmp610-tmp167*tmp241*tmp243*tmp244*tmp462*tmp468*tmp523*tmp590*tmp610+4.*tmp241*tmp244*tmp261*tmp432*tmp523*tmp590*tmp60*tmp610+6.*tmp196*tmp241*tmp243*tmp244*tmp432*tmp610*tmp612+6.*tmp167*tmp241*tmp244*tmp261*tmp432*tmp610*tmp612+4.*eta*tmp201*tmp210*tmp217*tmp241*tmp243*tmp300*tmp432*tmp610*tmp612+4.*tmp201*tmp226*tmp243*tmp244*tmp304*tmp432*tmp610*tmp612+2.*eta*tmp167*tmp210*tmp217*tmp226*tmp243*tmp300*tmp304*tmp432*tmp610*tmp612+8.*tmp241*tmp244*tmp345*tmp432*tmp610*tmp612+0.5*tmp1017*tmp167*tmp241*tmp243*tmp244*tmp462*tmp610*tmp612-2.*tmp201*tmp241*tmp243*tmp244*tmp462*tmp468*tmp610*tmp612-eta*tmp167*tmp210*tmp217*tmp241*tmp243*tmp300*tmp462*tmp468*tmp610*tmp612-tmp167*tmp226*tmp243*tmp244*tmp304*tmp462*tmp468*tmp610*tmp612-0.25*tmp1020*tmp167*tmp241*tmp243*tmp244*tmp588*tmp610*tmp612+4.*eta*tmp210*tmp217*tmp241*tmp261*tmp300*tmp432*tmp60*tmp610*tmp612+4.*tmp226*tmp244*tmp261*tmp304*tmp432*tmp60*tmp610*tmp612-2.*tmp241*tmp244*tmp261*tmp462*tmp468*tmp60*tmp610*tmp612-2.*tmp167*tmp241*tmp243*tmp244*tmp432*tmp523*tmp590*tmp666-4.*tmp201*tmp241*tmp243*tmp244*tmp432*tmp612*tmp666-2.*eta*tmp167*tmp210*tmp217*tmp241*tmp243*tmp300*tmp432*tmp612*tmp666-2.*tmp167*tmp226*tmp243*tmp244*tmp304*tmp432*tmp612*tmp666+1.*tmp167*tmp241*tmp243*tmp244*tmp462*tmp468*tmp612*tmp666-4.*tmp241*tmp244*tmp261*tmp432*tmp60*tmp612*tmp666+2.*tmp157*tmp445*tmp575*tmp684+1.*eta*tmp167*tmp241*tmp243*tmp300*tmp432*tmp561*tmp610*tmp612*tmp688-eta*tmp167*tmp217*tmp241*tmp243*tmp300*tmp432*tmp610*tmp612*tmp697-tmp167*tmp243*tmp244*tmp304*tmp432*tmp610*tmp612*tmp705-tmp157*tmp232*tmp445*tmp707-tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp575*tmp707+0.5*tmp246*tmp295*(12.*eta*tmp196*tmp248*tmp251+tmp719+tmp720+tmp721+tmp722+tmp725+tmp726+tmp727+tmp728+tmp729+tmp730+tmp731+tmp732+tmp733+tmp734+tmp735+tmp738+tmp742+tmp743+tmp744+tmp745+tmp746+tmp747+tmp748+tmp749+tmp750+tmp751+tmp752+tmp753+tmp754+tmp755+tmp756+tmp759+tmp760)+2.*tmp167*tmp241*tmp243*tmp432*tmp5*tmp561*tmp610*tmp612*tmp688*tmp768+2.*tmp167*tmp243*tmp244*tmp432*tmp610*tmp612*tmp780*tmp781-0.5*tmp297*tmp298*(-4.*tmp201*tmp228*tmp241*tmp243*tmp244+6.*tmp196*tmp231*tmp241*tmp243*tmp244+6.*tmp167*tmp231*tmp241*tmp244*tmp261-2.*eta*tmp167*tmp210*tmp217*tmp228*tmp241*tmp243*tmp300+4.*eta*tmp201*tmp210*tmp217*tmp231*tmp241*tmp243*tmp300-2.*tmp167*tmp226*tmp228*tmp243*tmp244*tmp304+4.*tmp201*tmp226*tmp231*tmp243*tmp244*tmp304+2.*eta*tmp167*tmp210*tmp217*tmp226*tmp231*tmp243*tmp300*tmp304+8.*tmp231*tmp241*tmp244*tmp345-4.*tmp228*tmp241*tmp244*tmp261*tmp60+4.*eta*tmp210*tmp217*tmp231*tmp241*tmp261*tmp300*tmp60+4.*tmp226*tmp231*tmp244*tmp261*tmp304*tmp60+eta*tmp167*tmp231*tmp241*tmp243*tmp300*tmp561*tmp688-eta*tmp167*tmp217*tmp231*tmp241*tmp243*tmp300*tmp697-tmp167*tmp231*tmp243*tmp244*tmp304*tmp705+tmp167*tmp241*tmp243*tmp244*tmp707+2.*tmp167*tmp231*tmp241*tmp243*tmp5*tmp561*tmp688*tmp768+2.*tmp167*tmp231*tmp243*tmp244*tmp780*tmp781)-0.25*tmp232*tmp448*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp825-0.25*tmp232*tmp455*tmp462*tmp475*tmp477*tmp501*tmp536*tmp537*tmp540*tmp825+1.*e3z*tmp157*tmp237*tmp950-tmp201*tmp442*tmp950+0.375*tmp232*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp859*tmp969-0.25*tmp232*tmp436*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp579*tmp995+0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp477*tmp504*tmp540*(0.+2.*tmp236*tmp569+tmp478*tmp707+tmp157*(1.*tmp1157*tmp149-24.*x->data[0]))-2.*tmp145*tmp228*tmp232*tmp236*x->data[0]-2.*e3z*tmp145*tmp228*tmp232*tmp425*tmp427*tmp434*tmp436*x->data[0]-e3z*tmp145*tmp237*tmp320*tmp425*tmp427*tmp434*tmp448*x->data[0]+1.*tmp145*tmp228*tmp232*tmp320*tmp427*tmp434*tmp445*tmp448*x->data[0]+3.*tmp145*tmp237*tmp434*tmp436*tmp445*tmp451*x->data[0]+2.*e3z*tmp145*tmp237*tmp425*tmp427*tmp436*tmp470*x->data[0]-2.*tmp145*tmp228*tmp232*tmp427*tmp436*tmp445*tmp470*x->data[0]-tmp145*tmp237*tmp320*tmp427*tmp445*tmp448*tmp470*x->data[0]-3.*tmp232*tmp436*tmp455*tmp462*tmp477*tmp501*tmp504*tmp534*tmp577*x->data[0]+1.5*tmp232*tmp320*tmp448*tmp455*tmp462*tmp477*tmp501*tmp504*tmp540*tmp577*x->data[0]+1.5*tmp232*tmp320*tmp455*tmp462*tmp477*tmp501*tmp536*tmp537*tmp540*tmp577*x->data[0]-3.*tmp232*tmp436*tmp455*tmp462*tmp477*tmp504*tmp540*tmp572*tmp577*x->data[0]+6.*tmp228*tmp436*tmp455*tmp462*tmp477*tmp501*tmp504*tmp540*tmp575*tmp577*x->data[0]+1.5*tmp232*tmp436*tmp459*tmp462*tmp477*tmp501*tmp504*tmp540*tmp577*tmp579*x->data[0]-1.5*tmp232*tmp436*tmp455*tmp462*tmp501*tmp504*tmp540*tmp577*tmp581*tmp586*x->data[0]+1.5*tmp232*tmp436*tmp455*tmp468*tmp477*tmp501*tmp504*tmp540*tmp577*tmp588*x->data[0]+2.*tmp145*tmp157*tmp575*tmp684*x->data[0]+2.*tmp145*tmp427*tmp434*tmp436*tmp445*tmp575*tmp684*x->data[0]-tmp145*tmp157*tmp232*tmp707*x->data[0]-tmp145*tmp232*tmp427*tmp434*tmp436*tmp445*tmp707*x->data[0]-0.5*tmp145*tmp237*tmp427*tmp434*tmp445*tmp448*tmp825*x->data[0]+1.*e3z*tmp145*tmp237*tmp427*tmp434*tmp436*tmp950*x->data[0]+0.75*tmp145*tmp237*tmp427*tmp434*tmp445*tmp859*tmp969*x->data[0]+1.*tmp145*tmp237*tmp427*tmp436*tmp445*x->data[0]*(1.*tmp432*tmp451+0.5*tmp1017*tmp427*tmp462-tmp158*tmp432*tmp475-0.25*tmp1020*tmp427*tmp588-0.5*tmp455*tmp995+0.25*tmp579*tmp997+1.*tmp451*tmp462*tmp468*x->data[0])+0.5*tmp232*tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*(8.*tmp190*tmp222*tmp237*tmp242*tmp254*tmp284*tmp347*tmp425+8.*tmp190*tmp222*tmp237*tmp254*tmp260*tmp284*tmp425*tmp505-2.*tmp190*tmp222*tmp228*tmp232*tmp254*tmp284*tmp313*tmp425*tmp505+2.*eta*tmp190*tmp210*tmp217*tmp237*tmp254*tmp284*tmp313*tmp425*tmp505+2.*tmp222*tmp226*tmp237*tmp254*tmp284*tmp313*tmp425*tmp505+36.*tmp190*tmp222*tmp237*tmp242*tmp284*tmp313*tmp508+8.*tmp190*tmp222*tmp237*tmp284*tmp339*tmp508-8.*tmp190*tmp222*tmp228*tmp232*tmp242*tmp284*tmp347*tmp508+8.*eta*tmp190*tmp210*tmp217*tmp237*tmp242*tmp284*tmp347*tmp508+8.*tmp222*tmp226*tmp237*tmp242*tmp284*tmp347*tmp508+12.*tmp158*tmp190*tmp222*tmp237*tmp284*tmp505*tmp508-8.*tmp190*tmp222*tmp228*tmp232*tmp260*tmp284*tmp505*tmp508+8.*eta*tmp190*tmp210*tmp217*tmp237*tmp260*tmp284*tmp505*tmp508+8.*tmp222*tmp226*tmp237*tmp260*tmp284*tmp505*tmp508-2.*eta*tmp190*tmp210*tmp217*tmp228*tmp232*tmp284*tmp313*tmp505*tmp508-2.*tmp222*tmp226*tmp228*tmp232*tmp284*tmp313*tmp505*tmp508+2.*eta*tmp210*tmp217*tmp226*tmp237*tmp284*tmp313*tmp505*tmp508+2.*eta*tmp158*tmp190*tmp210*tmp217*tmp526+2.*tmp158*tmp222*tmp226*tmp526+2.*tmp190*tmp222*tmp530+2.*eta*tmp158*tmp210*tmp217*tmp226*tmp530-4.*p->data[2]*tmp145*tmp158*tmp190*tmp222*tmp237*tmp254*tmp425*tmp505*tmp59-16.*p->data[2]*tmp145*tmp190*tmp222*tmp237*tmp242*tmp260*tmp508*tmp59+4.*p->data[2]*tmp145*tmp158*tmp190*tmp222*tmp228*tmp232*tmp505*tmp508*tmp59-4.*eta*p->data[2]*tmp145*tmp158*tmp190*tmp210*tmp217*tmp237*tmp505*tmp508*tmp59-4.*p->data[2]*tmp145*tmp158*tmp222*tmp226*tmp237*tmp505*tmp508*tmp59+2.*tmp190*tmp222*tmp284*tmp313*tmp505*tmp508*tmp575*tmp684-eta*tmp190*tmp237*tmp284*tmp313*tmp505*tmp508*tmp561*tmp688-eta*tmp158*tmp190*tmp530*tmp561*tmp688+1.*eta*tmp190*tmp217*tmp237*tmp284*tmp313*tmp505*tmp508*tmp697+1.*eta*tmp158*tmp190*tmp217*tmp530*tmp697+1.*tmp222*tmp237*tmp284*tmp313*tmp505*tmp508*tmp705+1.*tmp158*tmp222*tmp530*tmp705-tmp190*tmp222*tmp232*tmp284*tmp313*tmp505*tmp508*tmp707+2.*tmp190*tmp222*tmp237*tmp505*tmp508*tmp757*tmp758+1.*tmp190*tmp222*tmp237*tmp254*tmp284*tmp313*tmp505*tmp950+4.*tmp190*tmp222*tmp526*x->data[0]+4.*eta*tmp190*tmp210*tmp217*tmp530*x->data[0]+4.*tmp222*tmp226*tmp530*x->data[0]-12.*p->data[2]*tmp145*tmp190*tmp222*tmp237*tmp505*tmp508*tmp59*x->data[0]+1.*tmp158*tmp190*tmp222*(-2.*eta*tmp158*tmp190*tmp210*tmp217*tmp254*tmp268*tmp274*tmp425-2.*tmp158*tmp222*tmp226*tmp254*tmp268*tmp274*tmp425-2.*tmp158*tmp190*tmp222*tmp254*tmp268*tmp278*tmp279*tmp425+tmp1077*tmp242*tmp508-2.*tmp190*tmp222*tmp268*tmp274*tmp508-2.*eta*tmp158*tmp210*tmp217*tmp226*tmp268*tmp274*tmp508-2.*eta*tmp158*tmp190*tmp210*tmp217*tmp268*tmp278*tmp279*tmp508-2.*tmp158*tmp222*tmp226*tmp268*tmp278*tmp279*tmp508+2.*tmp508*tmp516+2.*tmp242*tmp254*tmp425*tmp523+eta*tmp158*tmp190*tmp268*tmp274*tmp508*tmp561*tmp688-eta*tmp158*tmp190*tmp217*tmp268*tmp274*tmp508*tmp697-tmp158*tmp222*tmp268*tmp274*tmp508*tmp705+tmp158*tmp190*tmp222*tmp268*tmp508*tmp736*tmp737-tmp158*tmp190*tmp222*tmp268*tmp279*tmp508*tmp741-tmp158*tmp190*tmp222*tmp254*tmp268*tmp274*tmp950+tmp242*tmp254*tmp516*tmp950-4.*tmp190*tmp222*tmp254*tmp268*tmp274*tmp425*x->data[0]-4.*eta*tmp190*tmp210*tmp217*tmp268*tmp274*tmp508*x->data[0]-4.*tmp222*tmp226*tmp268*tmp274*tmp508*x->data[0]-4.*tmp190*tmp222*tmp268*tmp278*tmp279*tmp508*x->data[0]+4.*tmp254*tmp425*tmp516*x->data[0]+4.*tmp508*tmp523*x->data[0]))-tmp201*(tmp425*tmp425)+3.*tmp436*tmp455*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp684*1./(tmp231*tmp231*tmp231*tmp231)+7.5*tmp158*tmp232*tmp436*tmp455*tmp462*tmp477*tmp501*tmp504*tmp540*pow(tmp242,-3.5)+0.75*tmp298*(tmp306*tmp306)*pow(tmp245,-2.5)-0.25*tmp246*(tmp289*tmp289)*(1.0/sqrt(tmp294*tmp294*tmp294))+0.375*tmp232*tmp436*tmp462*tmp475*tmp477*tmp501*tmp504*tmp540*tmp997*pow(tmp428,-2.5)+0.375*tmp1020*tmp232*tmp436*tmp455*tmp475*tmp477*tmp501*tmp504*tmp540*pow(tmp431,-2.5)+0.375*tmp232*tmp455*tmp462*tmp475*tmp477*tmp501*tmp537*tmp540*tmp859*(1.0/(tmp435*tmp435))+0.25*tmp232*tmp448*tmp455*tmp462*tmp475*tmp477*tmp501*tmp540*tmp859*1./(tmp503*tmp503*tmp503)+2.*tmp167*tmp241*tmp243*tmp244*tmp432*tmp610*(tmp523*tmp523)*1./(tmp516*tmp516*tmp516)+1.*tmp167*tmp241*tmp243*tmp244*tmp432*tmp612*(-0.25*tmp1200*tmp1201*tmp608+0.5*tmp1195*tmp581*tmp608+tmp581*tmp586*tmp664+tmp477*(0.-tmp145*tmp254*tmp425*tmp427*tmp432*tmp503*tmp592*tmp595-tmp145*tmp158*tmp254*tmp425*tmp432*tmp451*tmp503*tmp592*tmp595-0.5*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp592*tmp595-0.5*tmp145*tmp158*tmp320*tmp432*tmp436*tmp451*tmp508*tmp592*tmp595-0.5*tmp145*tmp427*tmp462*tmp468*tmp503*tmp508*tmp592*tmp595-0.5*tmp145*tmp158*tmp451*tmp462*tmp468*tmp503*tmp508*tmp592*tmp595+0.5*tmp145*tmp260*tmp432*tmp475*tmp503*tmp508*tmp592*tmp595+1.*p->data[2]*tmp432*tmp451*tmp503*tmp508*tmp59*tmp592*tmp595+1.*p->data[2]*tmp254*tmp425*tmp427*tmp432*tmp503*tmp59*tmp592*tmp595*tmp60+0.5*p->data[2]*tmp320*tmp427*tmp432*tmp436*tmp508*tmp59*tmp592*tmp595*tmp60+0.5*p->data[2]*tmp427*tmp462*tmp468*tmp503*tmp508*tmp59*tmp592*tmp595*tmp60+0.5*tmp145*tmp427*tmp432*tmp503*tmp508*tmp595*tmp624*tmp629+0.5*tmp145*tmp158*tmp432*tmp451*tmp503*tmp508*tmp595*tmp624*tmp629-0.5*p->data[2]*tmp427*tmp432*tmp503*tmp508*tmp59*tmp595*tmp60*tmp624*tmp629-tmp145*tmp427*tmp432*tmp503*tmp508*tmp592*tmp637-tmp145*tmp158*tmp432*tmp451*tmp503*tmp508*tmp592*tmp637+1.*p->data[2]*tmp427*tmp432*tmp503*tmp508*tmp59*tmp592*tmp60*tmp637+tmp455*tmp459*tmp662+0.5*tmp455*tmp606*tmp995-0.25*tmp579*tmp606*tmp997-0.5*tmp145*tmp254*tmp320*tmp425*tmp427*tmp432*tmp436*tmp592*tmp595*x->data[0]-0.5*tmp145*tmp254*tmp425*tmp427*tmp462*tmp468*tmp503*tmp592*tmp595*x->data[0]-0.25*tmp145*tmp320*tmp427*tmp436*tmp462*tmp468*tmp508*tmp592*tmp595*x->data[0]-1.5*tmp145*tmp432*tmp451*tmp503*tmp508*tmp592*tmp595*x->data[0]-0.25*tmp1017*tmp145*tmp427*tmp462*tmp503*tmp508*tmp592*tmp595*x->data[0]+0.125*tmp1020*tmp145*tmp427*tmp503*tmp508*tmp588*tmp592*tmp595*x->data[0]+0.5*tmp145*tmp254*tmp425*tmp427*tmp432*tmp503*tmp595*tmp624*tmp629*x->data[0]+0.25*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp595*tmp624*tmp629*x->data[0]+0.25*tmp145*tmp427*tmp462*tmp468*tmp503*tmp508*tmp595*tmp624*tmp629*x->data[0]-tmp145*tmp254*tmp425*tmp427*tmp432*tmp503*tmp592*tmp637*x->data[0]-0.5*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp592*tmp637*x->data[0]-0.5*tmp145*tmp427*tmp462*tmp468*tmp503*tmp508*tmp592*tmp637*x->data[0]+0.5*tmp145*tmp427*tmp432*tmp503*tmp508*tmp624*tmp629*tmp637*x->data[0]+0.25*tmp145*tmp427*tmp432*tmp503*tmp508*tmp595*tmp624*(8.*eta*tmp190*tmp210*tmp217*tmp222*tmp226*tmp274*tmp313+16.*eta*tmp210*tmp217*tmp222*tmp260*tmp274*tmp338+4.*eta*tmp210*tmp217*tmp222*tmp278*tmp279*tmp313*tmp338+16.*tmp190*tmp226*tmp260*tmp274*tmp342+4.*tmp190*tmp226*tmp278*tmp279*tmp313*tmp342+12.*tmp158*tmp274*tmp338*tmp342+8.*tmp260*tmp278*tmp279*tmp338*tmp342-2.*eta*tmp222*tmp274*tmp313*tmp338*tmp561*tmp688+2.*tmp274*tmp313*tmp338*tmp5*tmp561*tmp688+2.*eta*tmp217*tmp222*tmp274*tmp313*tmp338*tmp697+2.*tmp190*tmp274*tmp313*tmp342*tmp705-tmp313*tmp338*tmp342*tmp736*tmp737+tmp279*tmp313*tmp338*tmp342*tmp741+2.*tmp274*tmp313*tmp342*tmp780)*x->data[0]-0.25*tmp145*tmp427*tmp432*tmp436*tmp508*tmp592*tmp595*tmp825*x->data[0]+0.125*tmp145*tmp427*tmp432*tmp448*tmp508*tmp592*tmp595*tmp859*x->data[0]-0.5*tmp145*tmp254*tmp427*tmp432*tmp503*tmp592*tmp595*tmp950*x->data[0]-0.5*tmp145*tmp427*tmp432*tmp503*tmp508*tmp592*x->data[0]*(-4.*tmp190*tmp222-4.*eta*tmp158*tmp210*tmp217*tmp226+tmp1157*tmp477-0.25*tmp1200*tmp1201*tmp497+0.5*tmp1195*tmp497*tmp581+tmp567*tmp581*tmp586+2.*eta*tmp158*tmp190*tmp561*tmp688-2.*eta*tmp158*tmp190*tmp217*tmp697-2.*tmp158*tmp222*tmp705-8.*eta*tmp190*tmp210*tmp217*x->data[0]-8.*tmp222*tmp226*x->data[0])+tmp429*(0.+2.*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp602+2.*tmp145*tmp158*tmp320*tmp432*tmp436*tmp451*tmp508*tmp602-2.*p->data[2]*tmp320*tmp427*tmp432*tmp436*tmp508*tmp59*tmp60*tmp602+2.*tmp145*tmp254*tmp425*tmp427*tmp432*tmp602*tmp604+2.*tmp145*tmp158*tmp254*tmp425*tmp432*tmp451*tmp602*tmp604+1.*tmp145*tmp427*tmp462*tmp468*tmp508*tmp602*tmp604+1.*tmp145*tmp158*tmp451*tmp462*tmp468*tmp508*tmp602*tmp604-tmp145*tmp260*tmp432*tmp475*tmp508*tmp602*tmp604-2.*p->data[2]*tmp432*tmp451*tmp508*tmp59*tmp602*tmp604-2.*p->data[2]*tmp254*tmp425*tmp427*tmp432*tmp59*tmp60*tmp602*tmp604-p->data[2]*tmp427*tmp462*tmp468*tmp508*tmp59*tmp60*tmp602*tmp604+2.*tmp145*tmp427*tmp432*tmp508*tmp604*tmp656+2.*tmp145*tmp158*tmp432*tmp451*tmp508*tmp604*tmp656-2.*p->data[2]*tmp427*tmp432*tmp508*tmp59*tmp60*tmp604*tmp656+2.*tmp145*tmp254*tmp320*tmp425*tmp427*tmp432*tmp436*tmp602*x->data[0]+1.*tmp145*tmp320*tmp427*tmp436*tmp462*tmp468*tmp508*tmp602*x->data[0]+1.*tmp145*tmp254*tmp425*tmp427*tmp462*tmp468*tmp602*tmp604*x->data[0]+3.*tmp145*tmp432*tmp451*tmp508*tmp602*tmp604*x->data[0]+0.5*tmp1017*tmp145*tmp427*tmp462*tmp508*tmp602*tmp604*x->data[0]-0.25*tmp1020*tmp145*tmp427*tmp508*tmp588*tmp602*tmp604*x->data[0]+2.*tmp145*tmp320*tmp427*tmp432*tmp436*tmp508*tmp656*x->data[0]+2.*tmp145*tmp254*tmp425*tmp427*tmp432*tmp604*tmp656*x->data[0]+1.*tmp145*tmp427*tmp462*tmp468*tmp508*tmp604*tmp656*x->data[0]+1.*tmp145*tmp427*tmp432*tmp436*tmp508*tmp602*tmp825*x->data[0]-0.5*tmp145*tmp427*tmp432*tmp448*tmp508*tmp602*tmp859*x->data[0]+1.*tmp145*tmp254*tmp427*tmp432*tmp602*tmp604*tmp950*x->data[0]+1.*tmp145*tmp427*tmp432*tmp508*tmp604*x->data[0]*(8.*tmp260*tmp345+2.*tmp159*tmp201*tmp228*tmp232*tmp241*tmp244*tmp600-3.*tmp167*tmp237*tmp241*tmp244*tmp600+3.*tmp159*tmp196*tmp237*tmp241*tmp244*tmp600+1.*eta*tmp159*tmp167*tmp210*tmp217*tmp228*tmp232*tmp241*tmp300*tmp600+2.*eta*tmp159*tmp201*tmp210*tmp217*tmp237*tmp241*tmp300*tmp600+1.*tmp159*tmp167*tmp226*tmp228*tmp232*tmp244*tmp304*tmp600+2.*tmp159*tmp201*tmp226*tmp237*tmp244*tmp304*tmp600+1.*eta*tmp159*tmp167*tmp210*tmp217*tmp226*tmp237*tmp300*tmp304*tmp600-2.*tmp228*tmp232*tmp241*tmp244*tmp60*tmp600-2.*eta*tmp210*tmp217*tmp237*tmp241*tmp300*tmp60*tmp600-2.*tmp226*tmp237*tmp244*tmp304*tmp60*tmp600-tmp159*tmp167*tmp228*tmp232*tmp241*tmp244*tmp649-2.*tmp159*tmp201*tmp237*tmp241*tmp244*tmp649-eta*tmp159*tmp167*tmp210*tmp217*tmp237*tmp241*tmp300*tmp649-tmp159*tmp167*tmp226*tmp237*tmp244*tmp304*tmp649+2.*tmp237*tmp241*tmp244*tmp60*tmp649+1.*tmp159*tmp167*tmp241*tmp244*tmp575*tmp600*tmp684+0.5*eta*tmp159*tmp167*tmp237*tmp241*tmp300*tmp561*tmp600*tmp688-0.5*eta*tmp159*tmp167*tmp217*tmp237*tmp241*tmp300*tmp600*tmp697-0.5*tmp159*tmp167*tmp237*tmp244*tmp304*tmp600*tmp705-0.5*tmp159*tmp167*tmp232*tmp241*tmp244*tmp600*tmp707+1.*tmp159*tmp167*tmp237*tmp241*tmp5*tmp561*tmp600*tmp688*tmp768+1.*tmp159*tmp167*tmp237*tmp244*tmp600*tmp780*tmp781-6.*tmp261*x->data[0]+0.5*tmp159*tmp167*tmp237*tmp241*tmp244*(tmp1157*tmp159-24.*eta*tmp158*tmp190*tmp210*tmp217-24.*tmp158*tmp222*tmp226-8.*eta*tmp210*tmp217*tmp226*tmp260+2.*tmp497+4.*eta*tmp190*tmp260*tmp561*tmp688-4.*eta*tmp190*tmp217*tmp260*tmp697-4.*tmp222*tmp260*tmp705-24.*tmp190*tmp222*x->data[0]+4.*tmp567*x->data[0])))-0.375*tmp145*tmp427*tmp432*tmp503*tmp508*tmp595*x->data[0]*(tmp629*tmp629)*pow(tmp591,-2.5)))))/sqrt(tmp680);