exact_derivatives-px.c

Go to the documentation of this file.

#define Power(A,B) pow(A,B)
#define Sqrt(A) sqrt(A)
#define Log(A) log(A)
const double prTprm  =  csi2*nx;
const double prTtimesoneminuscsi1invprm  =  (1. - 1./csi1)*prTprm;
const double tmpP0prm  =  1.0 - nx*prTtimesoneminuscsi1invprm;
const double tmpP1prm  =  -(ny*prTtimesoneminuscsi1invprm);
const double tmpP2prm  =  -(nz*prTtimesoneminuscsi1invprm);
const double pxirprm  =  r*(xiUSCOREx*tmpP0prm + xiUSCOREy*tmpP1prm + xiUSCOREz*tmpP2prm);
const double pvrprm  =  r*(vx*tmpP0prm + vy*tmpP1prm + vz*tmpP2prm);
const double pvrsqprm  =  2*pvr*pvrprm;
const double pnprm  =  nx*tmpP0prm + ny*tmpP1prm + nz*tmpP2prm;
const double pnsqprm  =  2*pn*pnprm;
const double prprm  =  pnprm;
const double prsqprm  =  2*pr*prprm;
const double pfprm  =  pxirprm;
const double pxirsqprm  =  2*pxir*pxirprm;
const double ptheta2prm  =  invxi2*pvrsqprm;
const double prT4prm  =  4*((prT)*(prT)*(prT))*prTprm;
const double Hnspt4prm  =  2*Hnspt6*pf*pfprm + Hnspt7*prsqprm + Hnspt5*prT4prm + invrho2*ptheta2prm;
const double Hnspt3prm  =  deltaT*Hnspt4prm;
const double Hnspt2prm  =  rho2*Hnspt3prm;
const double Hnspt1prm  =  ww*pfprm;
const double Hnsprm  =  invLambda*Hnspt1prm + (invLambda*Hnspt2prm)/(2.*Sqrt(invLambda*Hnspt2));
const double Qprm  =  Qpt3*pnsqprm + Qpt1*pvrsqprm + Qpt2*pxirsqprm;
const double pn2prm  =  deltaR*invrho2*prsqprm;
const double ppprm  =  Qprm;
const double sKerrmultfactprm  =  -36.*r*pn2prm + 3.*r*ppprm;
const double sStarmultfactprm  =  r*(-30.*pn2prm + 4.*ppprm);
const double deltaSigmaStarUSCOREx1prm  =  etaover12r*(sKerrUSCOREx*sKerrmultfactprm + sStarUSCOREx*sStarmultfactprm);
const double deltaSigmaStarUSCOREy1prm  =  etaover12r*(sKerrUSCOREy*sKerrmultfactprm + sStarUSCOREy*sStarmultfactprm);
const double deltaSigmaStarUSCOREz1prm  =  etaover12r*(sKerrUSCOREz*sKerrmultfactprm + sStarUSCOREz*sStarmultfactprm);
const double pn2ppprm  =  pp*pn2prm + pn2*ppprm;
const double pp2prm  =  2*pp*ppprm;
const double pn2u2prm  =  u2*pn2prm;
const double ppu2prm  =  u2*ppprm;
const double pn2ppu2prm  =  u2*pn2ppprm;
const double sMultiplier1pt6prm  =  -720.*pn2*pn2prm + 126.*pn2ppprm + 3.*pp2prm;
const double sMultiplier1pt5prm  =  -96.*pn2ppprm + 23.*pp2prm;
const double sMultiplier1pt4prm  =  324.*pn2prm - 120.*ppprm + r*sMultiplier1pt6prm;
const double sMultiplier1pt3prm  =  -282.*pn2prm + 206.*ppprm + r*sMultiplier1pt5prm;
const double sMultiplier1pt2prm  =  r*sMultiplier1pt4prm;
const double sMultiplier1pt1prm  =  eta*sMultiplier1pt2prm + r*sMultiplier1pt3prm;
const double sMultiplier1prm  =  -0.013888888888888888*eta*u2*sMultiplier1pt1prm;
const double sMultiplier2pt6prm  =  5.625*pn2u2*pn2prm - 1.625*pn2ppu2prm + 5.625*pn2*pn2u2prm;
const double sMultiplier2pt5prm  =  0.25*pn2ppu2prm - 0.3125*u2*pp2prm;
const double sMultiplier2pt4prm  =  -6.125*pn2u2prm + 1.4166666666666665*ppu2prm + r*sMultiplier2pt6prm;
const double sMultiplier2pt3prm  =  -0.6666666666666666*pn2u2prm - 3.0277777777777777*ppu2prm + r*sMultiplier2pt5prm;
const double sMultiplier2pt2prm  =  r*sMultiplier2pt4prm;
const double sMultiplier2pt1prm  =  eta*sMultiplier2pt2prm + r*sMultiplier2pt3prm;
const double sMultiplier2prm  =  eta*sMultiplier2pt1prm;
const double deltaSigmaStarUSCOREx2prm  =  deltaSigmaStarUSCOREx1prm + sigmaStar->data[0]*sMultiplier1prm + sigmaKerr->data[0]*sMultiplier2prm;
const double deltaSigmaStarUSCOREy2prm  =  deltaSigmaStarUSCOREy1prm + sigmaStar->data[1]*sMultiplier1prm + sigmaKerr->data[1]*sMultiplier2prm;
const double deltaSigmaStarUSCOREz2prm  =  deltaSigmaStarUSCOREz1prm + sigmaStar->data[2]*sMultiplier1prm + sigmaKerr->data[2]*sMultiplier2prm;
const double deltaSigmaStarUSCOREx3prm  =  deltaSigmaStarUSCOREx2prm;
const double deltaSigmaStarUSCOREy3prm  =  deltaSigmaStarUSCOREy2prm;
const double deltaSigmaStarUSCOREz3prm  =  deltaSigmaStarUSCOREz2prm;
const double deltaSigmaStarUSCORExprm  =  deltaSigmaStarUSCOREx3prm;
const double deltaSigmaStarUSCOREyprm  =  deltaSigmaStarUSCOREy3prm;
const double deltaSigmaStarUSCOREzprm  =  deltaSigmaStarUSCOREz3prm;
const double sxprm  =  deltaSigmaStarUSCORExprm;
const double syprm  =  deltaSigmaStarUSCOREyprm;
const double szprm  =  deltaSigmaStarUSCOREzprm;
const double sxiprm  =  xiUSCOREx*sxprm + xiUSCOREy*syprm + xiUSCOREz*szprm;
const double svprm  =  vx*sxprm + vy*syprm + vz*szprm;
const double snprm  =  nx*sxprm + ny*syprm + nz*szprm;
const double s3prm  =  e3USCOREx*sxprm + e3USCOREy*syprm + e3USCOREz*szprm;
const double sqrtQprm  =  Qprm/(2.*Sqrt(Q));
const double oneplus2sqrtQprm  =  2.*sqrtQprm;
const double oneplus1sqrtQprm  =  oneplus2sqrtQprm - sqrtQprm;
const double twoB1psqrtQsqrtQprm  =  2.*B*sqrtQ*oneplus1sqrtQprm + 2.*B*oneplus1sqrtQ*sqrtQprm;
const double invtwoB1psqrtQsqrtQprm  =  (-twoB1psqrtQsqrtQprm)/((twoB1psqrtQsqrtQ)*(twoB1psqrtQsqrtQ));
const double expMUsqsqrtQplusQprm  =  expMUsq*(Qprm + sqrtQprm);
const double Hwrpt4aprm  =  sv*pxirsqprm + pxirsq*svprm;
const double Hwrpt4prm  =  expMUsqexpnusq*Hwrpt4aprm;
const double Hwrpt3cprm  =  sxi*pxirprm + pxir*sxiprm;
const double Hwrpt3bprm  =  pvr*Hwrpt3cprm + Hwrpt3c*pvrprm;
const double Hwrpt3aprm  =  expMUexpnu*Hwrpt3bprm;
const double Hwrpt3prm  =  B*Hwrpt3aprm;
const double Hwrpt2gprm  =  deltaR*svprm;
const double Hwrpt2fprm  =  sqrtdeltaR*snprm;
const double Hwrpt2eprm  =  pvr*Hwrpt2fprm + Hwrpt2f*pvrprm;
const double Hwrpt2dprm  =  pnsq*Hwrpt2gprm + Hwrpt2g*pnsqprm;
const double Hwrpt2cprm  =  pn*Hwrpt2eprm + Hwrpt2e*pnprm;
const double Hwrpt2bprm  =  sv*expMUsqsqrtQplusQprm + expMUsqsqrtQplusQ*svprm;
const double Hwrpt2aprm  =  xi2*(Hwrpt2bprm + Hwrpt2cprm - Hwrpt2dprm);
const double Hwrpt2prm  =  deltaT*Hwrpt2aprm;
const double Hwrpt1bprm  =  invxi2*invtwoB1psqrtQsqrtQprm;
const double Hwrpt1aprm  =  sqrtdeltaR*Hwrpt1bprm;
const double Hwrpt1prm  =  invexpMUcubinvexpnu*Hwrpt1aprm;
const double Hwrprm  =  (Hwrpt2 - Hwrpt3 + Hwrpt4)*Hwrpt1prm + Hwrpt1*(Hwrpt2prm - Hwrpt3prm + Hwrpt4prm);
const double Hwcospt9prm  =  sxi*pxirprm + pxir*sxiprm;
const double Hwcospt8prm  =  sv*pvrprm + pvr*svprm;
const double Hwcospt7prm  =  B*Hwcospt8prm - expMUexpnu*Hwcospt9prm;
const double Hwcospt6prm  =  sqrtdeltaR*Hwcospt7prm;
const double Hwcospt5prm  =  -(xi2*expMUsqsqrtQplusQprm) + pvrsqprm;
const double Hwcospt4prm  =  pn*Hwcospt6prm + Hwcospt6*pnprm;
const double Hwcospt3prm  =  deltaT*Hwcospt5prm - expMUsqexpnusq*pxirsqprm;
const double Hwcospt2prm  =  sn*Hwcospt3prm - B*Hwcospt4prm + Hwcospt3*snprm;
const double Hwcospt1prm  =  invexpMUcubinvexpnu*Hwcospt2prm;
const double Hwcosprm  =  invtwoB1psqrtQsqrtQ*Hwcospt1prm + Hwcospt1*invtwoB1psqrtQsqrtQprm;
const double deltaTsqrtQprm  =  deltaT*sqrtQprm;
const double invdeltatTsqrtQprm  =  (-deltaTsqrtQprm)/((deltaTsqrtQ)*(deltaTsqrtQ));
const double HSOLpt5prm  =  (-B + expMUexpnu)*pxirprm;
const double HSOLpt4prm  =  invexpMU*HSOLpt5prm;
const double HSOLpt3prm  =  expnusq*HSOLpt4prm;
const double HSOLpt2prm  =  s3*HSOLpt3prm + HSOLpt3*s3prm;
const double HSOLpt1prm  =  invxi2*HSOLpt2prm;
const double HSOLprm  =  invdeltatTsqrtQ*HSOLpt1prm + HSOLpt1*invdeltatTsqrtQprm;
const double deltaTsqrtQplusQprm  =  deltaT*(Qprm + sqrtQprm);
const double invdeltaTsqrtQplusQprm  =  (-deltaTsqrtQplusQprm)/((deltaTsqrtQplusQ)*(deltaTsqrtQplusQ));
const double HSONLmult2prm  =  invxi2*invdeltaTsqrtQplusQprm;
const double HSONLmultprm  =  expnuinvexpMU2*HSONLmult2prm;
const double HSONLpt1bprm  =  xi2*pnprm;
const double HSONLpt1aprm  =  (-mucos + nucos)*HSONLpt1bprm + (mur - nur)*pvrprm;
const double HSONLpt1prm  =  sqrtQ*HSONLpt1aprm - mucos*HSONLpt1bprm + mur*pvrprm + HSONLpt1a*sqrtQprm;
const double HSONLpt2dprm  =  nur*pxirprm;
const double HSONLpt2cprm  =  oneplus2sqrtQ*HSONLpt2dprm + HSONLpt2d*oneplus2sqrtQprm;
const double HSONLpt2bprm  =  B*sxiprm;
const double HSONLpt2aprm  =  expMUexpnu*HSONLpt2cprm;
const double HSONLpt2prm  =  HSONLpt2b*HSONLpt1prm + sv*HSONLpt2aprm + HSONLpt1*HSONLpt2bprm + HSONLpt2a*svprm;
const double HSONLpt3cprm  =  sv*pxirprm + pxir*svprm;
const double HSONLpt3bprm  =  oneplus1sqrtQ*HSONLpt3cprm + HSONLpt3c*oneplus1sqrtQprm;
const double HSONLpt3aprm  =  expMUexpnu*HSONLpt3bprm;
const double HSONLpt3prm  =  B*HSONLpt2prm - BR*HSONLpt3aprm;
const double HSONLpt4eprm  =  xi2*snprm;
const double HSONLpt4dprm  =  oneplus2sqrtQ*HSONLpt4eprm + HSONLpt4e*oneplus2sqrtQprm;
const double HSONLpt4cprm  =  pxir*HSONLpt4dprm + HSONLpt4d*pxirprm;
const double HSONLpt4bprm  =  nucos*HSONLpt4cprm;
const double HSONLpt4aprm  =  expMUexpnu*HSONLpt4bprm;
const double HSONLpt4prm  =  sqrtdeltaR*HSONLpt3prm - B*HSONLpt4aprm;
const double HSONLprm  =  HSONLpt4*HSONLmultprm + HSONLmult*HSONLpt4prm;
const double Hsprm  =  HSOLprm + HSONLprm + wcos*Hwcosprm + wr*Hwrprm + w*s3prm;
const double Hsspt1prm  =  -0.5*(-6.*sn*snprm + 2*(sx*sxprm + sy*syprm + sz*szprm));
const double Hssprm  =  u3*Hsspt1prm;
const double Hprm  =  Hnsprm + Hsprm + Hssprm;
const double Hrealprm  =  (eta*Hprm)*invHreal;