Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALSimulation 6.2.0.1-5e288d3
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
exact_derivatives-z.c
Go to the documentation of this file.
1#define Power(A,B) pow(A,B)
2#define Sqrt(A) sqrt(A)
3#define Log(A) log(A)
4const double r2prm = 2*x->data[2];
5const double rprm = r2prm/(2.*Sqrt(r2));
6const double uprm = (-rprm)/((r)*(r));
7const double u2prm = 2*u*uprm;
8const double u3prm = u2*uprm + u*u2prm;
9const double u4prm = 2*u2*u2prm;
10const double u5prm = u4*uprm + u*u4prm;
11const double etau3prm = eta*u3prm;
12const double etau4prm = eta*u4prm;
13const double nxprm = x->data[0]*uprm;
14const double nyprm = x->data[1]*uprm;
15const double nzprm = x->data[2]*uprm + u;
17const double xi2prm = -2*costheta*costhetaprm;
24const double w2prm = r2prm;
25const double rho2prm = 2*a2*costheta*costhetaprm + r2prm;
26const double bulkprm = 2.*invm1PlusetaKK*uprm + a2*u2prm;
27const double loguprm = uprm/u;
28const double logargprm = coeffs->k5l*u5*loguprm + coeffs->k1*uprm + coeffs->k2*u2prm + coeffs->k3*u3prm + coeffs->k4*u4prm + coeffs->k5*u5prm + coeffs->k5l*logu*u5prm;
34const double deltaUUSCOREupt7prm = coeffs->k5l*loguprm;
43const double Lambdaprm = -(a2*xi2*deltaTprm) + 2*w2*w2prm - a2*deltaT*xi2prm;
52const double expMUprm = rho2prm/(2.*Sqrt(rho2));
54const double expMUsqprm = 2*expMU*expMUprm;
55const double expnusqprm = 2*expnu*expnuprm;
62const double DDprm = ((-8.666666666666666 + eta)*etau3prm - eta*u2prm)/(-0.16666666666666666 + (-8.666666666666666 + eta)*etau3 - eta*u2);
64const double wwprm = a*(2.*rprm + (a2*coeffs->b3 + coeffs->bb3)*eta*uprm);
65const double Bprm = deltaTprm/(2.*Sqrt(deltaT));
66const double sqrtdeltaTprm = Bprm;
67const double sqrtdeltaRprm = deltaRprm/(2.*Sqrt(deltaR));
76const double wwUSCORErprm = -(a*(a2*coeffs->b3 + coeffs->bb3)*eta*u2prm);
92const double etaover12rprm = 0.08333333333333333*eta*uprm;
93const double csiprm = (deltaR*w2*deltaTprm + deltaT*(w2*deltaRprm - 2*deltaR*w2prm))/(2.*Sqrt(deltaR*deltaT)*((w2)*(w2)));
94const double csi1prm = (1.-fabs(1.-tortoise)) * (csiprm);
95const double csi2prm = (0.5-copysign(0.5, 1.5-tortoise)) * (csiprm);
96const double prTprm = (p->data[0]*nx + p->data[1]*ny + p->data[2]*nz)*csi2prm + csi2*(p->data[0]*nxprm + p->data[1]*nyprm + p->data[2]*nzprm);
97const double prTtimesoneminuscsi1invprm = (prT*csi1prm)/((csi1)*(csi1)) + (1. - 1./csi1)*prTprm;
103const double pvrsqprm = 2*pvr*pvrprm;
105const double pnsqprm = 2*pn*pnprm;
106const double prprm = pnprm;
107const double prsqprm = 2*pr*prprm;
108const double pfprm = pxirprm;
109const double pxirsqprm = 2*pxir*pxirprm;
111const double prT4prm = 4*((prT)*(prT)*(prT))*prTprm;
114const double Hnspt5prm = qq*u2prm;
118const double Hnspt1prm = ww*pfprm + pf*wwprm;
125const double ppprm = Qprm;
126const double sKerrmultfactprm = r*(-36.*pn2prm + 3.*ppprm) + (-36.*pn2 + 3.*pp)*rprm;
127const double sStarmultfactprm = r*(-30.*pn2prm + 4.*ppprm) + (-30.*pn2 + 4.*pp)*rprm;
131const double pn2ppprm = pp*pn2prm + pn2*ppprm;
132const double pp2prm = 2*pp*ppprm;
133const double pn2u2prm = u2*pn2prm + pn2*u2prm;
134const double ppu2prm = u2*ppprm + pp*u2prm;
136const double sMultiplier1pt6prm = -720.*pn2*pn2prm + 126.*pn2ppprm + 3.*pp2prm;
137const double sMultiplier1pt5prm = -96.*pn2ppprm + 23.*pp2prm;
142const double sMultiplier1prm = -0.013888888888888888*eta*u2*sMultiplier1pt1prm - 0.013888888888888888*eta*sMultiplier1pt1*u2prm;
143const double sMultiplier2pt6prm = 5.625*pn2u2*pn2prm - 1.625*pn2ppu2prm + 5.625*pn2*pn2u2prm;
144const double sMultiplier2pt5prm = 0.25*pn2ppu2prm - 0.3125*u2*pp2prm - 0.3125*pp2*u2prm;
145const double sMultiplier2pt4prm = -6.125*pn2u2prm + 1.4166666666666665*ppu2prm + sMultiplier2pt6*rprm + r*sMultiplier2pt6prm;
146const double sMultiplier2pt3prm = -0.6666666666666666*pn2u2prm - 3.0277777777777777*ppu2prm + sMultiplier2pt5*rprm + r*sMultiplier2pt5prm;
150const double deltaSigmaStarUSCOREx2prm = deltaSigmaStarUSCOREx1prm + sigmaStar->data[0]*sMultiplier1prm + sigmaKerr->data[0]*sMultiplier2prm;
151const double deltaSigmaStarUSCOREy2prm = deltaSigmaStarUSCOREy1prm + sigmaStar->data[1]*sMultiplier1prm + sigmaKerr->data[1]*sMultiplier2prm;
152const double deltaSigmaStarUSCOREz2prm = deltaSigmaStarUSCOREz1prm + sigmaStar->data[2]*sMultiplier1prm + sigmaKerr->data[2]*sMultiplier2prm;
153const double deltaSigmaStarUSCOREx3prm = deltaSigmaStarUSCOREx2prm + coeffs->d1*sigmaStar->data[0]*etau3prm;
154const double deltaSigmaStarUSCOREy3prm = deltaSigmaStarUSCOREy2prm + coeffs->d1*sigmaStar->data[1]*etau3prm;
155const double deltaSigmaStarUSCOREz3prm = deltaSigmaStarUSCOREz2prm + coeffs->d1*sigmaStar->data[2]*etau3prm;
156const double deltaSigmaStarUSCORExprm = deltaSigmaStarUSCOREx3prm + coeffs->d1v2*sigmaKerr->data[0]*etau3prm;
157const double deltaSigmaStarUSCOREyprm = deltaSigmaStarUSCOREy3prm + coeffs->d1v2*sigmaKerr->data[1]*etau3prm;
158const double deltaSigmaStarUSCOREzprm = deltaSigmaStarUSCOREz3prm + coeffs->d1v2*sigmaKerr->data[2]*etau3prm;
163const double svprm = vx*sxprm + vy*syprm + vz*szprm + sx*vxprm + sy*vyprm + sz*vzprm;
164const double snprm = sx*nxprm + sy*nyprm + sz*nzprm + nx*sxprm + ny*syprm + nz*szprm;
166const double sqrtQprm = Qprm/(2.*Sqrt(Q));
167const double oneplus2sqrtQprm = 2.*sqrtQprm;
191const double Hwcospt8prm = sv*pvrprm + pvr*svprm;
212const double HSONLpt1bprm = xi2*pnprm + pn*xi2prm;
217const double HSONLpt2bprm = sxi*Bprm + B*sxiprm;
224const double HSONLpt4eprm = xi2*snprm + sn*xi2prm;
232const double Hsspt1prm = -0.5*(-6.*sn*snprm + 2*(sx*sxprm + sy*syprm + sz*szprm));
233const double Hssprm = u3*Hsspt1prm + Hsspt1*u3prm;
234const double Hpt1prm = (s1dots1 + s2dots2)*etau4prm;
235const double Hprm = Hnsprm + (coeffs->dheffSSv2 + coeffs->dheffSS*sKerrdotsStar)*Hpt1prm + Hsprm + Hssprm;
236const double Hrealprm = (eta*Hprm)*invHreal;
const double HSONLpt4b
const double sKerrmultfact
const double Hwrpt3b
const double pn2u2
const double Hwcospt3
const double HSONLpt2a
const double Hnspt1
const double pn2
const double sMultiplier2pt6
const double sx
const double pp
const double Hsspt1
const double Hwrpt2d
const double Hwrpt1
const double pvr
const double Hwrpt2e
const double Hnspt2
const double HSONLpt1
const double sxi
const double invdeltaTsqrtQplusQ
const double expMUsqsqrtQplusQ
const double csi1
const double HSONLmult
const double oneplus2sqrtQ
const double Hwcospt9
const double pvrsq
const double sStarmultfact
const double Hnspt7
const double prT
const double Hwrpt2
const double HSONLpt1a
const double sMultiplier1pt1
const double HSOLpt3
const double sz
const double tmpP0
const double prsq
const double Q
const double sMultiplier1pt5
const double invHreal
const double Hwrpt3a
const double Hnspt6
const double deltaTsqrtQ
const double Hwcospt5
const double HSONLmult2
const double Hwcospt7
const double HSONLpt2
const double prT4
const double pnsq
const double HSONLpt3b
const double twoB1psqrtQsqrtQ
const double Hwrpt4
const double pn
const double sy
const double oneplus1sqrtQ
const double csi2
const double HSONLpt3c
const double Hwcospt8
const double Hwcospt6
const double pxirsq
const double HSOLpt4
const double Hwrpt1a
const double Hwcospt1
const double HSONLpt2d
const double invdeltatTsqrtQ
const double Hwr
const double Hnspt4
const double Hwrpt2b
const double tmpP1
const double pn2pp
const double sMultiplier1pt6
const double pp2
const double Hwcospt2
const double sv
const double Hwrpt2f
const double HSOLpt1
const double Hnspt3
const double HSONLpt2b
const double HSOLpt2
const double Hwrpt4a
const double sMultiplier1pt3
const double Hwrpt2c
const double HSONLpt3a
const double sMultiplier1pt4
const double sqrtQ
const double HSOLpt5
const double Hwrpt3c
const double HSONLpt3
const double Qpt1
const double sKerrdotsStar
const double pxir
const double Hnspt5
const double Qpt3
const double HSONLpt4d
const double sMultiplier2pt4
const double Qpt2
const double HSONLpt4e
const double sMultiplier2pt3
const double tmpP2
const double sMultiplier2pt5
const double HSONLpt4c
const double pr
const double pf
const double s3
const double deltaTsqrtQplusQ
const double HSONLpt4
const double HSONLpt1b
const double invtwoB1psqrtQsqrtQ
const double HSONLpt4a
const double sn
const double Hwrpt1b
const double HSONLpt2c
const double Hwcos
const double ptheta2
const double Hwrpt2g
const double Hwcospt4
const double prTtimesoneminuscsi1inv
const double Hwrpt3
const double Hwrpt2a
const double sStarUSCOREy
const double ww
const double deltaTUSCOREr
const double deltaUUSCOREupt5
const double deltaU
const double nur
const double deltaUUSCOREupt4
const double etau3
const double xiUSCOREy
const double mucos
const double expnu
const double invLambda
const double invdeltaT
const double e3USCOREz
const double deltaR
const double e3USCOREx
const double invxi2
const double wwUSCOREr
const double invexpMU
const double deltaUUSCOREupt3
const double LambdaUSCOREr
const double u
const double s2dots2
const double w2
const double sqrtdeltaT
const double nucospt2
const double wcospt2
const double logu
const double deltaUUSCOREupt7
const double rho2invLambda
const double expMUsq
const double nucospt3
const double u3
const double invexpMUsq
const double nucospt1
const double bulk
const double deltaUUSCOREupt2
const double expnusq
const double wr
const double r2
const double invsqrtdeltaT
const double ny
const double invexpMUcubinvexpnu
const double wcospt1
const double invrho2xi2Lambda
const double sStarUSCOREz
const double nucos
const double vy
const double wcos
const double invm1PlusetaKK
const double deltaUUSCOREupt6
const double BR
const double u5
const double nurpt2
const double qq
const double nurpt1
const double invexpnuexpMU
const double a2
const double s1dots1
const double expMU
const double deltaUUSCOREupt1
const double xi2
const double sqrtdeltaTdeltaTsqrtdeltaR
const double w
const double sKerrUSCOREy
const double sStarUSCOREx
const double invrho2
const double sKerrUSCOREx
const double expMUexpnu
const double rho2xi2Lambda
const double rho2
const double sKerrUSCOREz
const double expMUsqexpnusq
const double u2
const double a2costheta
const double e3USCOREy
const double xiUSCOREx
const double vz
const double invLambdasq
const double etaover12r
const double xiUSCOREz
const double u4
const double mur
const double expnuinvexpMU2
const double invonepluslogarg
const double DD
const double Lambda
const double costheta
const double nz
const double invdeltaTsqrtdeltaTsqrtdeltaR
const double B
const double vx
const double deltaTsqrtdeltaR
const double invsqrtdeltaR
const double nx
const double onepluslogarg
const double logTerms
const double sqrtdeltaR
const double Hwrpt2eprm
const double HSOLpt4prm
const double Hwrpt4aprm
const double Hwrpt2bprm
const double DDprm
const double deltaSigmaStarUSCOREx3prm
const double sMultiplier1prm
const double mucosprm
const double expMUexpnuprm
const double HSOLpt2prm
const double wprm
const double nzprm
const double Hwrpt3cprm
const double sxiprm
const double xiUSCOREyprm
const double HSONLpt3bprm
const double deltaTsqrtdeltaRprm
const double sMultiplier1pt4prm
const double HSONLpt2prm
const double Lambdaprm
const double Hnspt5prm
const double deltaSigmaStarUSCOREzprm
const double Hwrpt2aprm
const double twoB1psqrtQsqrtQprm
const double HSONLpt2cprm
const double s3prm
const double Hssprm
const double sMultiplier2pt1prm
const double deltaSigmaStarUSCOREx1prm
const double Hwrpt3bprm
const double rho2xi2Lambdaprm
const double invonepluslogargprm
const double invxi2prm
const double pvrsqprm
const double HSONLpt3prm
const double logTermsprm
const double sMultiplier2pt4prm
const double expMUsqprm
const double etau4prm
const double wrprm
const double invexpMUcubinvexpnuprm
const double Hwcospt9prm
const double sqrtdeltaTdeltaTsqrtdeltaRprm
const double sMultiplier1pt3prm
const double Hwrpt1aprm
const double sMultiplier2pt5prm
const double invdeltaTsqrtQplusQprm
const double Hwrpt2gprm
const double invtwoB1psqrtQsqrtQprm
const double pnsqprm
const double HSONLpt4prm
const double HSONLpt2aprm
const double sMultiplier2pt3prm
const double deltaTsqrtQplusQprm
const double sqrtdeltaRprm
const double w2prm
#define Sqrt(A)
const double Hsprm
const double murprm
const double ppu2prm
const double deltaSigmaStarUSCOREyprm
const double etau3prm
const double deltaTUSCORErprm
const double prTprm
const double HSOLpt5prm
const double Hwcospt2prm
const double Hwrpt3aprm
const double rho2prm
const double Hwrpt2cprm
const double sxprm
const double Hwcospt5prm
const double HSONLmultprm
const double deltaTprm
const double HSOLprm
const double deltaUUSCOREuprm
const double csi1prm
const double invrho2prm
const double u5prm
const double deltaSigmaStarUSCORExprm
const double syprm
const double HSONLpt4dprm
const double nxprm
const double prsqprm
const double nucospt2prm
const double r2prm
const double HSONLpt3aprm
const double invexpMUprm
const double Hwrpt1prm
const double deltaUUSCOREupt5prm
const double wwUSCORErprm
const double uprm
const double bulkprm
const double expMUsqexpnusqprm
const double deltaSigmaStarUSCOREx2prm
const double szprm
const double wcosprm
const double pvrprm
const double Hwrpt2fprm
const double BRprm
const double xiUSCORExprm
const double nucosprm
const double prT4prm
const double Qpt2prm
const double csiprm
const double pxirsqprm
const double HSONLpt2bprm
const double sMultiplier1pt2prm
const double HSONLmult2prm
const double Bprm
const double deltaTsqrtQprm
const double invexpnuexpMUprm
const double sKerrmultfactprm
const double etaover12rprm
const double HSONLpt1bprm
const double sMultiplier1pt6prm
const double ppprm
const double expMUsqsqrtQplusQprm
const double tmpP0prm
const double a2costhetaprm
const double invsqrtdeltaTprm
const double tmpP1prm
const double snprm
const double Hwrpt3prm
const double HSONLpt2dprm
const double wcospt1prm
const double expnuinvexpMU2prm
const double Hsspt1prm
const double Hwrpt1bprm
const double Qpt3prm
const double tmpP2prm
const double oneplus2sqrtQprm
const double Hnspt4prm
const double Hrealprm
const double expMUprm
const double HSONLpt4cprm
const double nucospt1prm
const double deltaUUSCOREupt7prm
const double Hwcospt8prm
const double HSONLpt3cprm
const double expnusqprm
const double deltaUUSCOREupt3prm
const double invsqrtdeltaRprm
const double sqrtQprm
const double Hpt1prm
const double HSONLpt1prm
const double Hnspt7prm
const double Hwcospt6prm
const double Hwrpt2prm
const double pn2u2prm
const double deltaSigmaStarUSCOREy1prm
const double svprm
const double invdeltaTsqrtdeltaTsqrtdeltaRprm
const double deltaSigmaStarUSCOREy2prm
const double HSOLpt1prm
const double expnuprm
const double Qprm
const double sStarmultfactprm
const double nyprm
const double Hwcospt1prm
const double invrho2xi2Lambdaprm
const double invLambdasqprm
const double u2prm
const double vzprm
const double Hwrpt4prm
const double u4prm
const double wcospt2prm
const double u3prm
const double Hprm
const double HSONLpt4eprm
const double Hnspt6prm
const double sMultiplier1pt5prm
const double deltaSigmaStarUSCOREz3prm
const double sMultiplier1pt1prm
const double LambdaUSCORErprm
const double sMultiplier2prm
const double pn2prm
const double HSOLpt3prm
const double prTtimesoneminuscsi1invprm
const double sMultiplier2pt6prm
const double invexpMUsqprm
const double vyprm
const double Qpt1prm
const double nurpt1prm
const double deltaUUSCOREupt4prm
const double deltaRprm
const double HSONLpt4aprm
const double nurpt2prm
const double deltaUUSCOREupt2prm
const double ptheta2prm
const double sqrtdeltaTprm
const double Hnspt2prm
const double onepluslogargprm
const double invLambdaprm
const double HSONLpt1aprm
const double csi2prm
const double nucospt3prm
const double costhetaprm
const double vxprm
const double prprm
const double oneplus1sqrtQprm
const double deltaSigmaStarUSCOREz2prm
const double xi2prm
const double loguprm
const double deltaUUSCOREupt1prm
const double rprm
const double xiUSCOREzprm
const double Hwcospt3prm
const double deltaSigmaStarUSCOREz1prm
const double HSONLpt4bprm
const double sMultiplier2pt2prm
const double pxirprm
const double Hnspt3prm
const double Hwcosprm
const double Hwrprm
const double wwprm
const double HSONLprm
const double deltaUUSCOREupt6prm
const double Hwrpt2dprm
const double pn2ppprm
const double pn2ppu2prm
const double pnprm
const double Hnsprm
const double logargprm
const double Hnspt1prm
const double pp2prm
const double Hwcospt7prm
const double rho2invLambdaprm
const double deltaUprm
const double invdeltaTprm
const double nurprm
const double Hwcospt4prm
const double deltaSigmaStarUSCOREy3prm
const double invdeltatTsqrtQprm
const double pfprm
static const INT4 r
static const INT4 a
p
x
double deltaT
Definition: unicorn.c:24