50#include <lal/LALStdlib.h>
51#include<lal/LALConstants.h>
52#include<lal/LALInspiral.h>
73 hPlus = hPlus + v * v * v * v * v * ( sI*dM*
74 ( (1771./5120. - cI2*1667./5120. + cI4*217./9216. - cI6/9216.)
75 + eta*(681./256. + cI2*13./768. - cI4*35./768. + cI6/2304.)
76 + eta2*(-3451./9216. + cI2*673./3072. - cI4*5./9216. - cI6/3072.) )
77 *cos(phase) +
LAL_PI*( (19./3. + 3.*cI2 - cI4*2./3.)
78 + eta*(-16./3. + cI2*14./3. + 2.*cI4) )*cos(2.*phase)
79 + sI*dM* ( (3537./1024. - cI2*22977./5120. - cI4*15309./5120. + cI6*729./5120.)
80 + eta*(-23829./1280. + cI2*5529./1280. + cI4*7749./1280. - cI6*729./1280.)
81 + eta2*(29127./5120. - cI2*27267./5120. - cI4*1647./5120. + cI6*2187./5120.) )
82 *cos(3.*phase) - (16./3.)*
LAL_PI*(1.+cI2)*sI2*(1.-3.*eta)*cos(4.*phase)
83 + sI*dM*( (-108125./9216. + cI2*40625./9216. + cI4*83125./9216. - cI6*15625./9216.)
84 + eta*(8125./256. - cI2*40625./2304. - cI4*48125./2304. + cI6*15625./2304.)
85 + eta2*(-119375./9216. + cI2*40625./3072. + cI4*44375./9216. - cI6*15625./3072.) )
86 *cos(5.*phase) + dM*sI2*sI2*sI*(1.+cI2)*(117649./46080.)*(1. - 4.*eta + 3.*eta2)
87 *cos(7.*phase) + ( (-9./5. + cI2*14./5. + cI4*7./5.)
88 + eta*(32. + cI2*56./5. - cI4*28./5.) )*sin(2.*phase) + sI2*(1.+cI2)
89 *( (56./5. - log(2.)*32./3.) - eta*(1193./30. - 32.*log(2.)) ) *sin(4.*phase) );
90#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
94 hPlus = hPlus + v * v * v * v * ( (1./120.)*( (22.+396.*cI2 + 145.*cI4 - 5.*cI6)
95 + (5./3.)*eta*(706. - 216.*cI2 - 251.*cI4 + 15.*cI6)
96 - 5.*eta2* (98. - 108.*cI2 + 7.*cI4 + 5.*cI6) )* cos(2.*phase)
97 + (2./15.)*sI2* ( (59. + 35.*cI2 - 8.*cI4) - (5./3.)*eta*(131.+59.*cI2 - 24.*cI4)
98 + 5.*eta2* (21. - 3.*cI2 - 8.*cI4) )*cos(4.*phase)
99 - (81./40.)*(1. - 5.*eta + 5.*eta2)*sI2*sI2*(1.+cI2)* cos(6.*phase)
100 + (sI/40.)*dM* ( (11.+7.*cI2 + 10.*(5. + cI2)*log(2))*sin(phase)
101 - 5.*
LAL_PI*(5. + cI2)*cos(phase) - 27.*(7. - 10.*log(3./2.))*(1 + cI2)*sin(3.*phase)
102 + 135.*
LAL_PI*(1.+cI2)*cos(3.*phase) ) );
103#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
107 hPlus = hPlus + v * v * v * ( (sI/192.)* dM*
108 ( ( (57.+60.*cI2 - cI4) - 2.*eta*(49. - 12.*cI2 - cI4) )* cos(phase)
109 - (27./2.)*( (73. + 40.*cI2 - 9.*cI4) - 2.*eta*(25. - 8.*cI2 - 9.*cI4) )
110 *cos(3.*phase) + (625./2.)*(1. - 2.*eta)*sI2*(1+cI2)* cos(5.*phase) )
111 - 2.*
LAL_PI* (1.+cI2)*cos(2.*phase) );
112#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
116 hPlus = hPlus + v * v *( (1./6.)*( (19.+9.*cI2 - 2.*cI4) - eta*(19. - 11.*cI2 - 6.*cI4) )
117 * cos( 2.*phase) - (4./3.)*sI2*(1.+cI2)*(1. - 3*eta)* cos(4.*phase) );
118#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
122 hPlus = hPlus - v * (sI * dM / 8.) * ( (5. + cI2) * cos(phase) - 9.*(1. + cI2)*cos(3.*phase));
123#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
127 hPlus = hPlus - (1.+cI2) * cos(2.*phase);
129 default:
fprintf(stderr,
"There are no EOB waveforms at order %d in amplitude\n",
params->ampOrder);
139 REAL8 eta2 = eta*eta;
147 switch(
params->ampOrder )
150 hCross = hCross + v * v * v * v * v * ( cI*( (2. - cI2*22./5.) + eta*(-282./5. + cI2*94./5.) )
151 *cos(2.*phase) + cI*sI2*( (-112./5. + log(2.)*64./3.) + eta*(1193./15. - 64.*log(2.)) )
152 *cos(4.*phase) + sI*cI*dM*( (-913./7680. + cI2*1891./11520. - cI4*7./4608.)
153 + eta*(1165./384. - cI2*235./576. + cI4*7./1152.)
154 + eta2*(-1301./4608. + cI2*301./2304. - cI4*7./1536.) )*sin(phase)
155 + cI*
LAL_PI*( (34./3. - cI2*8./3.) - eta*(20./3. - 8.*cI2) )*sin(2.*phase)
156 + sI*cI*dM*( (12501./2560. - cI2*12069./1280. + cI4*1701./2560.)
157 + eta*(-19581./640. + cI2*7821./320. - cI4*1701./640.)
158 + eta2*(18903./2560. - cI2*11403./1280. + cI4*5103./2560.) )*sin(3.*phase)
159 - sI2*cI*
LAL_PI*(32./3.)*(1. - 3.*eta)*sin(4.*phase)
160 + sI*cI*dM*( (-101875./4608. + cI2*6875./256. - cI4*21875./4608.)
161 + eta*(66875./1152. - cI2*44375./576. + cI4*21875./1152.)
162 + eta2*(-100625./4608. + cI2*83125./2304. - cI4*21875./1536.) )*sin(5.*phase)
163 + sI2*sI2*sI*cI*dM*(117649./23040.)*(1. - 4.*eta + 3.*eta2)*sin(7.*phase) );
164#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
168 hCross = hCross + v * v * v * v * ( (cI/60.)*( (68. + 226.*cI2 - 15.*cI4)
169 + (5./3.)*eta*(572. - 490.*cI2 + 45.*cI4)
170 - 5.*eta2*(56. - 70.*cI2 + 15.*cI4) )*sin(2.*phase)
171 + (4./15.)*cI*sI2*( (55. - 12.*cI2) - (5./3.)*eta*(119. - 36.*cI2)
172 + 5.*eta2*(17. - 12.*cI2) )*sin(4.*phase)
173 - (81./20.)*(1. - 5.*eta + 5.*eta2)*cI*sI2*sI2* sin(6.*phase)
174 - (3./20.)*sI*cI*dM* ( (3.+10.*log(2))*cos(phase) + 5.*
LAL_PI*sin(phase)
175 - 9.*(7. - 10.*log(3./2.))*cos(3.*phase) - 45.*
LAL_PI*sin(3.*phase) ) );
176#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
180 hCross = hCross + v * v * v * ( (sI*cI/96.)*dM* ( ( (63. - 5.*cI2) - 2.*eta*(23. - 5.*cI2) )
181 *sin(phase) - (27./2.)*( (67. - 15.*cI2) - 2.*eta*(19. - 15.*cI2) )*sin(3.*phase)
182 + (625./2.)*(1. - 2.*eta)*sI2*sin(5.*phase) ) - 4.*
LAL_PI*cI* sin(2.*phase) ) ;
183#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
187 hCross = hCross + v * v * ( (cI/3.)*( (17. - 4.*cI2) - eta* (13. - 12.*cI2) )
188 *sin(2.*phase) - (8./3.)*(1. - 3.*eta)*cI*sI2*sin(4.*phase) );
189#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
193 hCross = hCross - v * (3./4.)*sI*cI*dM* ( sin(phase) - 3.* sin(3.*phase) );
194#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
198 hCross = hCross - 2.* cI * sin(2.*phase);
200 default:
fprintf(stderr,
"There are no EOB waveforms at order %d in amplitude\n",
params->ampOrder);
REAL4 LALInspiralHCrossPolarization(REAL8 phase, REAL8 v, InspiralTemplate *params)
REAL4 LALInspiralHPlusPolarization(REAL8 phase, REAL8 v, InspiralTemplate *params)
LAL_PNORDER_TWO_POINT_FIVE
LAL_PNORDER_ONE_POINT_FIVE
The inspiral waveform parameter structure containing information about the waveform to be generated.