LALSimulation  5.4.0.1-fe68b98
LALSimInspiralOptimizedCoefficientsEccentricityFD.c
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2014 Eliu Huerta
3  *
4  * This program is free software; you can redistribute it and/or modify
5  * it under the terms of the GNU General Public License as published by
6  * the Free Software Foundation; either version 2 of the License, or
7  * (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with with program; see the file COPYING. If not, write to the
16  * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
17  * MA 02110-1301 USA
18  */
19 
20 #include <lal/LALConstants.h>
21 #include <lal/LALAtomicDatatypes.h>
22 
23 #include <math.h>
24 
25 #ifdef __GNUC__
26 #define UNUSED __attribute__ ((unused))
27 #else
28 #define UNUSED
29 #endif
30 
31 
32 #define gamma (0.577215664901532860606512090)
33 
34 
35 static REAL8 UNUSED
36 C0(REAL8 eta)
37 {
38  return ( 3./(128.*eta));
39 }
40 
41 static REAL8 UNUSED
42 C1(REAL8 Mtotal)
43 {
44  return(LAL_PI*Mtotal);
45 }
46 
47 /*case 1*/
48 
49 static REAL8 UNUSED
50 C2(REAL8 e0, REAL8 f0)
51 {
52  REAL8 e2,e4,e6,e8,p1;
53  e2=e0*e0;
54  e4=e2*e2;
55  e6=e4*e2;
56  e8=e6*e2;
57  p1= pow(f0, 19./9.);
58  return(p1*((-2355.*e2)/1462. - (2608555.*e4)/444448. - (1326481225.*e6)/1.01334144e8 - (6505217202575.*e8)/2.77250217984e11) );
59 }
60 
61 static REAL8 UNUSED
62 C3(REAL8 e0, REAL8 f0)
63 {
64  REAL8 e2,e4,e6,e8,p2;
65  e2=e0*e0;
66  e4=e2*e2;
67  e6=e4*e2;
68  e8=e6*e2;
69  p2= pow(f0, 19./3.);
70  return( p2*((-75356125.*e6)/3.326976e6 - (250408403375.*e8)/1.011400704e9) );
71 }
72 
73 static REAL8 UNUSED
74 C4(REAL8 e0, REAL8 f0)
75 {
76  REAL8 e2,e4,e6,e8,p3;
77  e2=e0*e0;
78  e4=e2*e2;
79  e6=e4*e2;
80  e8=e6*e2;
81  p3= pow(f0, 38./9.);
82  return( p3*((5222765.*e4)/998944. + (17355248095.*e6)/4.55518464e8 +
83  (128274289063885.*e8)/8.30865678336e11) );
84 }
85 
86 static REAL8 UNUSED
87 C5(REAL8 e0, REAL8 f0)
88 {
89  REAL8 e2,e4,e6,e8,p4;
90  e2=e0*e0;
91  e4=e2*e2;
92  e6=e4*e2;
93  e8=e6*e2;
94  p4= pow(f0, 76./9.);
95  return( (4537813337273.*e8*p4)/(39444627456.) );
96 }
97 
98 /*case 2*/
99 
100 static REAL8 UNUSED
101 C6(REAL8 eta)
102 {
103  return( 3715./756. + (55.*eta)/9.);
104 }
105 
106 static REAL8 UNUSED
107 C7(REAL8 eta, REAL8 e0, REAL8 f0)
108 {
109  REAL8 e2,e4,e6,e8,p1;
110  e2=e0*e0;
111  e4=e2*e2;
112  e6=e4*e2;
113  e8=e6*e2;
114  p1= pow(f0, 19./9.);
115  return(p1*((-583255.*e2)/122808. - (1938156365.*e4)/1.12000896e8 -(985575550175*e6)/2.5536204288e10 - (690482340216175.*e8)/9.981007847424e12 - (8635.*e2*eta)/1462. - (28694105.*e4*eta)/1.333344e6 - (14591293475.*e6*eta)/3.04002432e8 - (71557389228325.*e8*eta)/8.31750653952e11) );
116 
117 }
118 
119 static REAL8 UNUSED
120 C8(REAL8 eta,REAL8 e0, REAL8 f0)
121 {
122  REAL8 e2,e4,e6,e8,p2;
123  e2=e0*e0;
124  e4=e2*e2;
125  e6=e4*e2;
126  e8=e6*e2;
127  p2= pow(f0, 19./3.);
128  return(p2*((-31207858707289299381995.*e6)/5.3807139303898e20 -(103703714484322341846369385.*e8)/1.635737034838499e23 -(462027517873731215615.*e6*eta)/6.405611821892619e18 -(1535317441894408829488645.*e8*eta)/1.9473059938553564e21) );
129 
130 }
131 
132 static REAL8 UNUSED
133 C9(REAL8 eta,REAL8 e0, REAL8 f0)
134 {
135  REAL8 e2,e4,e6,e8,p3;
136  e2=e0*e0;
137  e4=e2*e2;
138  e6=e4*e2;
139  e8=e6*e2;
140  p3= pow(f0, 38./9.);
141  return(p3*((1867778296655755.*e4)/1.34516772125568e14 +(6206627279787073865.*e6)/6.133964808925901e16 +(45873772442848006154795.*e8)/1.1188351811480843e20 +(27652168591135.*e4*eta)/1.601390144352e12 +(91888156228341605.*e6*eta)/7.30233905824512e14 +(679154100768947601215.*e8*eta)/1.33194664422391e18 ) );
142 }
143 
144 
145 static REAL8 UNUSED
146 C10(REAL8 eta,REAL8 e0, REAL8 f0)
147 {
148  REAL8 e2,e4,e6,e8,p4;
149  e2=e0*e0;
150  e4=e2*e2;
151  e6=e4*e2;
152  e8=e6*e2;
153  p4= pow(f0, 76./9.);
154  return(p4*((2222952136716664456903429135345831.*e8)/7.66183921683643e30 +
155  (32910462320165961003953863376587.*e8*eta)/9.121237162900513e28) );
156 }
157 
158 /*case 3*/
159 
160 static REAL8 UNUSED
161 C11(REAL8 UNUSED eta)
162 {
163  return( -16.*LAL_PI);
164 }
165 
166 
167 static REAL8 UNUSED
168 C12(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
169 {
170  REAL8 e2,e4,e6,e8,p1;
171  e2=e0*e0;
172  e4=e2*e2;
173  e6=e4*e2;
174  e8=e6*e2;
175  p1=pow(f0, 19./9.);
176  return(LAL_PI*p1*((7536.*e2)/731. + (521711.*e4)/13889. + (265296245.*e6)/3.166692e6 + (1301043440515.*e8)/8.664069312e9 ) );
177 }
178 
179 
180 static REAL8 UNUSED
181 C13(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
182 {
183  REAL8 e2,e4,e6,e8,p2;
184  e2=e0*e0;
185  e4=e2*e2;
186  e6=e4*e2;
187  e8=e6*e2;
188  p2=pow(f0, 19./3.);
189  return(LAL_PI*p2*((7800202806596945705.*e6)/6.672512314471478e16 +(25920073926321650577715.*e8)/2.0284437435993293e19) );
190 }
191 
192 static REAL8 UNUSED
193 C14(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
194 {
195  REAL8 e2,e4,e6,e8,p3;
196  e2=e0*e0;
197  e4=e2*e2;
198  e6=e4*e2;
199  e8=e6*e2;
200  p3=pow(f0, 38./9.);
201  return(LAL_PI*p3*((-475065859669.*e4)/1.6681147337e10 - (1578643851680087.*e6)/7.606603185672e12 -(11667906828579178421.*e8)/1.3874444210665728e16) );
202 }
203 
204 static REAL8 UNUSED
205 C15(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
206 {
207  REAL8 e2,e4,e6,e8,p4;
208  e2=e0*e0;
209  e4=e2*e2;
210  e6=e4*e2;
211  e8=e6*e2;
212  p4=pow(f0, 76./9.);
213  return(-((2759558801317902042155334469457.*e8*p4*LAL_PI)/4750644355677350183868572160.));
214 }
215 
216 /*case 4*/
217 
218 static REAL8 UNUSED
219 C16(REAL8 eta)
220 {
221  return(15293365./508032. + (27145.*eta)/504. + (3085.*eta*eta)/72.);
222 }
223 
224 static REAL8 UNUSED
225 C17(REAL8 eta,REAL8 e0, REAL8 f0)
226 {
227  REAL8 e2,e4,e6,e8,p1;
228  e2=e0*e0;
229  e4=e2*e2;
230  e6=e4*e2;
231  e8=e6*e2;
232  p1=pow(f0, 19./9.);
233  return(p1*((-2401058305.*e2)/2.47580928e8 - (7978716747515.*e4)/2.25793806336e11 -(4057272307914425.*e6)/5.1480987844608e13 -(2842476030950240425.*e8)/2.0121711820406784e16 - (4261765.*e2*eta)/245616. -(14161845095.*e4*eta)/2.24001792e8 - (7201466570525.*e6*eta)/5.1072408576e10 -(5045260598968525.*e8*eta)/1.9962015694848e13 -(484345.*e2*eta*eta)/35088. - (1609478435.*e4*eta*eta)/3.2000256e7 -(818438915825.*e6*eta*eta)/7.296058368e9 -(4013719013988775.*e8*eta*eta)/1.9962015694848e13) );
234 }
235 
236 static REAL8 UNUSED
237 C18(REAL8 eta,REAL8 e0, REAL8 f0)
238 {
239  REAL8 e2,e4,e6,e8,p2;
240  e2=e0*e0;
241  e4=e2*e2;
242  e6=e4*e2;
243  e8=e6*e2;
244  p2=pow(f0, 19./3.);
245  return(p2*( (-110474729572613560552840685.*e6)/1.0847519283665837e24 -(367107526369794861717089596255.*e8)/3.297645862234414e26 -(196087423156943883913505.*e6*eta)/1.07614278607796e21 -(651598507150524526244577115.*e8*eta)/3.271474069676998e23 -(22285124348468295519365.*e6*eta*eta)/1.5373468372542285e20 -(74053468209960146010849895.*e8*eta*eta)/4.673534385252855e22 ) );
246 }
247 
248 static REAL8 UNUSED
249 C19(REAL8 eta,REAL8 e0, REAL8 f0)
250 {
251  REAL8 e2,e4,e6,e8,p3;
252  e2=e0*e0;
253  e4=e2*e2;
254  e6=e4*e2;
255  e8=e6*e2;
256  p3=pow(f0, 38./9.);
257  return(p3*((6841716503626986565.*e4)/2.711858126051451e17 +(22735023941552476355495.*e6)/1.2366073054794616e20 +(168036723934429498871218085.*e8)/2.255571725194538e23 +(12143723404950745.*e4*eta)/2.69033544251136e14 +(40353592874651325635.*e6*eta)/1.2267929617851802e17 +(298257242353143912203705.*e8*eta)/2.2376703622961686e20 +(1380121079545885.*e4*eta*eta)/3.8433363464448e13 +(4586142347330975855.*e6*eta*eta)/1.7525613739788288e16 +(33896614207384379043965.*e8*eta*eta)/3.1966719461373837e19 ) );
258 }
259 
260 
261 static REAL8 UNUSED
262 C20(REAL8 eta,REAL8 e0, REAL8 f0)
263 {
264  REAL8 e2,e4,e6,e8,p4;
265  e2=e0*e0;
266  e4=e2*e2;
267  e6=e4*e2;
268  e8=e6*e2;
269  p4=pow(f0, 76./9.);
270  return(p4*( (7770499311332403684925585429305516241.*e8)/1.5446267861142245e34 +(13792268987637324946295665896844693.*e8*eta)/1.532367843367286e31 +(1567476508633676458254637292015689.*e8*eta*eta)/2.189096919096123e30 ));
271 }
272 
273 /*case 5*/
274 
275 static REAL8 UNUSED
276 C21(REAL8 eta)
277 {
278  return( LAL_PI*(38645./756. -65.*eta/9.));
279 }
280 
281 static REAL8 UNUSED
282 C22(REAL8 eta,REAL8 e0, REAL8 f0)
283 {
284  REAL8 e2,e4,e6,e8,p1;
285  e2=e0*e0;
286  e4=e2*e2;
287  e6=e4*e2;
288  e8=e6*e2;
289  p1=pow(f0, 19./9.);
290  return(LAL_PI*p1*( (6067265.*e2)/122808. + (20161521595.*e4)/1.12000896e8 +(10252373388025.*e6)/2.5536204288e10 +(7182689108386025.*e8)/9.981007847424e12 - (10205.*e2*eta)/1462. -(33911215.*e4*eta)/1.333344e6 - (17244255925.*e6*eta)/3.04002432e8 -(84567823633475.*e8*eta)/8.31750653952e11) );
291 }
292 
293 static REAL8 UNUSED
294 C23(REAL8 eta,REAL8 e0, REAL8 f0)
295 {
296  REAL8 e2,e4,e6,e8,p2;
297  e2=e0*e0;
298  e4=e2*e2;
299  e6=e4*e2;
300  e8=e6*e2;
301  p2=pow(f0, 19./3.);
302  return(LAL_PI*p2*((36828693183369973116475.*e6)/7.686734186271143e19 +(122381747448338420666046425.*e8)/2.3367671926264273e22 -(433615096349678814025.*e6*eta)/6.405611821892619e18 -(1440902965169982699005075.*e8*eta)/1.9473059938553564e21 ) );
303 }
304 
305 static REAL8 UNUSED
306 C24(REAL8 eta,REAL8 e0, REAL8 f0)
307 {
308  REAL8 e2,e4,e6,e8,p3;
309  e2=e0*e0;
310  e4=e2*e2;
311  e6=e4*e2;
312  e8=e6*e2;
313  p3=pow(f0, 38./9.);
314  return( LAL_PI*p3*( (-2316846009950855.*e4)/1.9216681732224e13 -(7698879291066691165.*e6)/8.762806869894144e15 -(56903148963613058870695.*e8)/1.5983359730686919e19 +(27278171420045.*e4*eta)/1.601390144352e12 +(90645363628809535.*e6*eta)/7.30233905824512e14 +(669968502482700007405.*e8*eta)/1.33194664422391e18));
315 }
316 
317 static REAL8 UNUSED
318 C25(REAL8 eta,REAL8 e0, REAL8 f0)
319 {
320  REAL8 e2,e4,e6,e8,p4;
321  e2=e0*e0;
322  e4=e2*e2;
323  e6=e4*e2;
324  e8=e6*e2;
325  p4=pow(f0, 76./9.);
326  return(LAL_PI*p4*( (-18041156334765213110425249225007393.*e8)/7.66183921683643e30 + (2334216112662079584736091244017.*e8*eta)/7.016336279154241e27));
327 }
328 
329 /*case 6*/
330 
331 static REAL8 UNUSED
332 C26(REAL8 eta)
333 {
334  return( 11583231236531./4694215680.- (15737765635.*eta)/3.048192e6 +
335  (76055.*eta*eta)/1728. - (127825.*eta*eta*eta)/1296. -
336  (6848.*gamma)/21. -(640.*LAL_PI*LAL_PI)/3. +(2255.*eta*LAL_PI*LAL_PI)/12.);
337 }
338 
339 /* The following cases should be multiplied by log(4.*pow(LAL_PI*Mtotal*f, 1./3.)) */
340 static REAL8 UNUSED
341 C27(REAL8 UNUSED eta)
342 {
343  return( -6848./21. );/* this term is multiplied by log(4x)*/
344 }
345 
346 static REAL8 UNUSED
347 C28(REAL8 UNUSED eta,REAL8 e0, REAL8 f0)
348 {
349  REAL8 e2,e4,e6,e8,p1;
350  e2=e0*e0;
351  e4=e2*e2;
352  e6=e4*e2;
353  e8=e6*e2;
354  p1=pow(f0, 19./9.);
355  return(p1*((-537568.*e2)/5117. - (111646154.*e4)/291669. - (28386698215.*e6)/3.3250266e7 -(19887378305015.*e8)/1.2996103968e10) );/* this term is multiplied by log(4x)*/
356 }
357 
358 
359 static REAL8 UNUSED
360 C29(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
361 {
362  REAL8 e2,e4,e6,e8,p2;
363  e2=e0*e0;
364  e4=e2*e2;
365  e6=e4*e2;
366  e8=e6*e2;
367  p2=pow(f0, 19./3.);
368  return(p2*( (-657204451373573967571.*e6)/7.006137930195053e17 -(2183890391914386294238433.*e8)/2.1298659307792957e20) );/* this term is multiplied by log(4x)*/
369 }
370 
371 static REAL8 UNUSED
372 C30(REAL8 UNUSED eta,REAL8 e0, REAL8 f0)
373 {
374  REAL8 e2,e4,e6,e8,p3;
375  e2=e0*e0;
376  e4=e2*e2;
377  e6=e4*e2;
378  e8=e6*e2;
379  p3=pow(f0, 38./9.);
380  return(p3*((83880153412870.*e4)/3.50304094077e11 +(139366874895483505.*e6)/7.9869333449556e13 + (1030073825416757818915.*e8)/1.4568166421199014e17 ) );/* this term is multiplied by log(4x)*/
381 }
382 
383 static REAL8 UNUSED
384 C31(REAL8 UNUSED eta,REAL8 e0, REAL8 f0)
385 {
386  REAL8 e2,e4,e6,e8,p4;
387  e2=e0*e0;
388  e4=e2*e2;
389  e6=e4*e2;
390  e8=e6*e2;
391  p4=pow(f0, 76./9.);
392  return((229018097854706671787019720252139.*e8*p4)/49881765734612176930620007680.);/* this term is multiplied by log(4x)*/
393 }
394 
395 static REAL8 UNUSED
396 C32(REAL8 eta, REAL8 e0, REAL8 f0)
397 {
398  REAL8 e2,e4,e6,e8,p1;
399  e2=e0*e0;
400  e4=e2*e2;
401  e6=e4*e2;
402  e8=e6*e2;
403  p1=pow(f0, 19./9.);
404  return(p1*( (1578237767500487.*e2)/2.28764777472e12 +(5244484101404118301.*e4)/2.08633477054464e15 +(533376501191162085059.*e6)/9.513686553683558e16 +(373677141943502503806739.*e8)/3.718492344411174e19 -(2470829204695.*e2*eta)/1.485485568e9 -(8210565447201485.*e4*eta)/1.354762838016e12 -(4175170127655540575.*e6*eta)/3.08885927067648e14 -(2925073821111304814575.*e8*eta)/1.207302709224407e17 +(11940635.*e2*eta*eta)/842112. +(39678730105.*e4*eta*eta)/7.68006144e8 +(20177105913475.*e6*eta*eta)/1.75105400832e11 +(98950858868368325.*e8*eta*eta)/4.79088376676352e14 -(20068525.*e2*eta*eta*eta)/631584. -(66687708575.*e4*eta*eta*eta)/5.76004608e8 -(33911492517125.*e6*eta*eta*eta)/1.31329050624e11 -(166305877783829875.*e8*eta*eta*eta)/3.59316282507264e14 -(537568.*e2*gamma)/5117. - (111646154.*e4*gamma)/291669. -(28386698215.*e6*gamma)/3.3250266e7 -(19887378305015.*e8*gamma)/1.2996103968e10 -(50240.*e2*LAL_PI*LAL_PI)/731. -(10434220.*e4*LAL_PI*LAL_PI)/41667. -(1326481225.*e6*LAL_PI*LAL_PI)/2.375019e6 -(6505217202575.*e8*LAL_PI*LAL_PI)/6.498051984e9 +(354035.*e2*eta*LAL_PI*LAL_PI)/5848. +(1176458305.*e4*eta*LAL_PI*LAL_PI)/5.333376e6 +(598243032475.*e6*eta*LAL_PI*LAL_PI)/1.216009728e9 +(2933852958361325.*e8*eta*LAL_PI*LAL_PI)/3.327002615808e12 ) );
405 }
406 
407 static REAL8 UNUSED
408 C33(REAL8 eta, REAL8 e0, REAL8 f0)
409 {
410  REAL8 e2,e4,e6,e8,p2;
411  e2=e0*e0;
412  e4=e2*e2;
413  e6=e4*e2;
414  e8=e6*e2;
415  p2=pow(f0, 19./3.);
416  return(p2*( (62534662812190508908451288422363.*e6)/1.0023107818107234e28 +(207802684524909061102783631427512249.*e8)/3.047024776704599e30 -(96662893738280943308126058155.*e6*eta)/6.508511570199502e24 -(321210795892307574612902891249065.*e8*eta)/1.9785875173406487e27 +(8813910165617557415555.*e6*eta*eta)/6.961570583792733e19 +(29288623480347143291889265.*e8*eta*eta)/2.116317457472991e22 -(785113635484365349577225.*e6*eta*eta*eta)/2.7672243070576115e21 -(2608932610714546056645118675.*e8*eta*eta*eta)/8.412361893455138e23 -(657204451373573967571.*e6*gamma)/7.006137930195053e17 -(2183890391914386294238433.*e8*gamma)/2.1298659307792957e20 -(30710488381942708765.*e6*LAL_PI*LAL_PI)/5.004384235853609e16 - (102050952893195621226095.*e8*LAL_PI*LAL_PI)/1.5213328076994972e19 +(13850430260256161653015.*e6*eta*LAL_PI*LAL_PI)/2.5622447287570477e19 + (46024979754831225172968845.*e8*eta*LAL_PI*LAL_PI)/7.789223975421425e21) );
417 }
418 
419 static REAL8 UNUSED
420 C34(REAL8 eta, REAL8 e0, REAL8 f0)
421 {
422  REAL8 e2,e4,e6,e8,p3;
423  e2=e0*e0;
424  e4=e2*e2;
425  e6=e4*e2;
426  e8=e6*e2;
427  p3=pow(f0, 38./9.);
428  return(p3*( (-796520595220381729758415.*e4)/5.0115138169430816e20 -(2646837937917328487987213045.*e6)/2.285250300526045e23 -(19563030899655064495340095274735.*e8)/4.168296548159506e26 +(6168627083362586227675.*e4*eta)/1.6271148756308705e18 +(20498347798013874034564025.*e6*eta)/7.41964383287677e20 +(151505237861278885566710654075.*e8*eta)/1.3533430351167228e24 -(562467383866675.*e4*eta*eta)/1.7403787229184e13 -(1869079116588961025.*e6*eta*eta)/7.936126976507904e15 -(13814541490402312805075.*e8*eta*eta)/1.4475495605150417e19 +(50102713130841625.*e4*eta*eta*eta)/6.91800542360064e14 +(166491315733786719875.*e6*eta*eta*eta)/3.154610473161892e17 +(1230553147045766992549625.*e8*eta*eta*eta)/5.754009503047291e20 +(83880153412870.*e4*gamma)/3.50304094077e11 +(139366874895483505.*e6*gamma)/7.9869333449556e13 +(1030073825416757818915.*e8*gamma)/1.4568166421199014e17 +(7839266674100.*e4*LAL_PI*LAL_PI)/5.0043442011e10 +(6512470789508575.*e6*LAL_PI*LAL_PI)/5.704952389254e12 +(48134290907325131725.*e8*LAL_PI*LAL_PI)/1.0405833157999296e16 - (883877317504775.*e4*eta*LAL_PI*LAL_PI)/6.405560577408e12 -(2937124326068367325.*e6*eta*LAL_PI*LAL_PI)/2.920935623298048e15 - (21708565199203634407975.*e8*eta*LAL_PI*LAL_PI)/5.32778657689564e18) );
429 }
430 
431 static REAL8 UNUSED
432 C35(REAL8 eta, REAL8 e0, REAL8 f0)
433 {
434  REAL8 e2,e4,e6,e8,p4;
435  e2=e0*e0;
436  e4=e2*e2;
437  e6=e4*e2;
438  e8=e6*e2;
439  p4=pow(f0, 76./9.);
440  return(p4*( (-21803333020968369933785155127082363071942707.*e8)/7.136175751847717e38 + (6736884392917513798607444845703601288679.*e8*eta)/9.267760716685346e34 -(614282187703745932737146352451999.*e8*eta*eta)/9.91289170911452e29 + (54718202538837159478833263067719005.*e8*eta*eta*eta)/3.940374454373021e31 +(229018097854706671787019720252139.*e8*gamma)/4.988176573461218e28 +(2140356054716884783056259067777.*e8*LAL_PI*LAL_PI)/7.125966533516026e26 - (965300580677315037158372839567427.*e8*eta*LAL_PI*LAL_PI)/3.648494865160205e29));
441 }
442 
443 /* case 7*/
444 
445 static REAL8 UNUSED
446 C36(REAL8 eta)
447 {
448  return( LAL_PI*(77096675./254016. + (378515.*eta)/1512. - (74045.*eta*eta)/756.));
449 }
450 
451 static REAL8 UNUSED
452 C37(REAL8 eta, REAL8 e0, REAL8 f0)
453 {
454  REAL8 e2,e4,e6,e8,p1;
455  e2=e0*e0;
456  e4=e2*e2;
457  e6=e4*e2;
458  e8=e6*e2;
459  p1=pow(f0, 19./9.);
460  return(LAL_PI*p1*((12104177975.*e2)/6.1895232e7 + (40222183410925.*e4)/5.6448451584e10 +(20453458379485375.*e6)/1.2870246961152e13 +(14329446184895255375.*e8)/5.030427955101696e15 + (59426855*e2*eta)/368424. +(197475439165.*e4*eta)/3.36002688e8 +(100418608176175.*e6*eta)/7.6608612864e10 +(70352065412362175.*e8*eta)/2.9943023542272e13 -(11625065.*e2*eta*eta)/184212. -(38630090995.*e4*eta*eta)/1.68001344e8 -(19643860461025.*e6*eta*eta)/3.8304306432e10 -(13762251650419025.*e8*eta*eta)/1.4971511771136e13 ) );
461 
462 }
463 
464 
465 static REAL8 UNUSED
466 C38(REAL8 eta, REAL8 e0, REAL8 f0)
467 {
468  REAL8 e2,e4,e6,e8,p2;
469  e2=e0*e0;
470  e4=e2*e2;
471  e6=e4*e2;
472  e8=e6*e2;
473  p2=pow(f0, 19./3.);
474  return(LAL_PI*p2*((434593322826290843757317275.*e6)/2.7118798209164593e23 +(1444153611751764473805565304825.*e8)/8.244114655586036e25 +(2133685941573919740699595.*e6*eta)/1.61421417911694e21 +(7090238383850135298344754185.*e8*eta)/4.907211104515498e23 -(417391055952448085809285.*e6*eta*eta)/8.0710708955847e20 -(1386990478929984989144254055.*e8*eta*eta)/2.453605552257749e23 ) );
475 }
476 
477 static REAL8 UNUSED
478 C39(REAL8 eta, REAL8 e0, REAL8 f0)
479 {
480  REAL8 e2,e4,e6,e8,p3;
481  e2=e0*e0;
482  e4=e2*e2;
483  e6=e4*e2;
484  e8=e6*e2;
485  p3=pow(f0, 38./9.);
486  return(LAL_PI*p3*( (-28083426938595740975.*e4)/6.779645315128627e16 -(93321227716953647259925.*e6)/3.091518263698654e19 -(689746068418917003152253775.*e8)/5.638929312986345e22 -(137878817052260255.*e4*eta)/4.03550316376704e14 -(458171309064660827365.*e6*eta)/1.8401894426777702e17 -(3386387715003096689295295.*e8*eta)/3.3565055434442526e20 +(26971816199185265.*e4*eta*eta)/2.01775158188352e14 +(89627345229892635595.*e6*eta*eta)/9.200947213388851e16 +(662444231688055412226385.*e8*eta*eta)/1.6782527717221263e20) );
487 }
488 
489 static REAL8 UNUSED
490 C40(REAL8 eta, REAL8 e0, REAL8 f0)
491 {
492  REAL8 e2,e4,e6,e8,p4;
493  e2=e0*e0;
494  e4=e2*e2;
495  e6=e4*e2;
496  e8=e6*e2;
497  p4=pow(f0, 76./9.);
498  return(LAL_PI*p4*((-1776389932392373211125413523490219335.*e8)/2.2715099795797417e32 -(8721390841557034022662273046715703.*e8*eta)/1.3520892735593702e30 +(1706076073241722479183197515934809.*e8*eta*eta)/6.760446367796851e29 ) );
499 }
500 
501  /* The following coefficients are for the real part of zeta_l's*/
502 
503  /*case 1 of zeta_re*/
504 static REAL8 UNUSED
505 z1(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
506 {
507  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p5, s2i,c2i;
508  e2=e0*e0;
509  e3=e2*e0;
510  e4=e3*e0;
511  e5=e4*e0;
512  e6=e5*e0;
513  e7=e6*e0;
514  c2b=cos(2.*bet);
515  s2b=sin(2.*bet);
516  si=sin(inc);
517  ci=cos(inc);
518  s2i=si*si;
519  c2i=ci*ci;
520  P=FPlus;
521  Q=FCross;
522  p5=pow(f0, 133./18.);
523  return( p5*( (46031066168471.*c2b*e7*P)/(1.2136808448e11) + (46031066168471.*c2b*c2i*e7*P)/(1.2136808448e11) - (46031066168471.*ci*e7*Q*s2b)/(6.068404224e10) - (8391437082143.*e7*P*s2i)/(3.6410425344e10)));
524 }
525 
526 static REAL8 UNUSED
527 z2(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
528 {
529  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p6, s2i,c2i;
530  e2=e0*e0;
531  e3=e2*e0;
532  e4=e3*e0;
533  e5=e4*e0;
534  e6=e5*e0;
535  e7=e6*e0;
536  c2b=cos(2.*bet);
537  s2b=sin(2.*bet);
538  si=sin(inc);
539  ci=cos(inc);
540  s2i=si*si;
541  c2i=ci*ci;
542  P=FPlus;
543  Q=FCross;
544  p6=pow(f0, 95./18.);
545  return(p6*( -(717415013.*c2b*e5*P)/(1.3307904e7) -(717415013.*c2b*c2i*e5*P)/(1.3307904e7) -(11919850440995.*c2b*e7*P)/(2.4273616896e10) -(11919850440995.*c2b*c2i*e7*P)/(2.4273616896e10) +(717415013.*ci*e5*Q*s2b)/(6.653952e6) +(11919850440995.*ci*e7*Q*s2b)/(1.2136808448e10)+(220389695.*e5*P*s2i)/(6.653952e6) +(3661774782425.*e7*P*s2i)/(1.2136808448e10)) );
546 }
547 
548 static REAL8 UNUSED
549 z3(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
550 {
551  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p7,s2i,c2i;
552  e2=e0*e0;
553  e3=e2*e0;
554  e4=e3*e0;
555  e5=e4*e0;
556  e6=e5*e0;
557  e7=e6*e0;
558  c2b=cos(2.*bet);
559  s2b=sin(2.*bet);
560  si=sin(inc);
561  ci=cos(inc);
562  s2i=si*si;
563  c2i=ci*ci;
564  P=FPlus;
565  Q=FCross;
566  p7=pow(f0, 19./6.);
567  return(p7*( (30299.*c2b*e3*P)/(3648.) +(30299.*c2b*c2i*e3*P)/(3648.) +(100683577.*c2b*e5*P)/(2.217984e6) +(100683577.*c2b*c2i*e5*P)/(2.217984e6) +(384584085937.*c2b*e7*P)/(2.697068544e9) +(384584085937.*c2b*c2i*e7*P)/(2.697068544e9)-(30299.*ci*e3*Q*s2b)/(1824.) -(100683577.*ci*e5*Q*s2b)/(1.108992e6) -(384584085937.*ci*e7*Q*s2b)/(1.348534272e9)-(9517.*e3*P*s2i)/(1824.) -(31624991.*e5*P*s2i)/(1.108992e6) -(120798928871.*e7*P*s2i)/(1.348534272e9)) );
568 }
569 
570 static REAL8 UNUSED
571 z4(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
572 {
573  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p8,s2i,c2i;
574  e2=e0*e0;
575  e3=e2*e0;
576  e4=e3*e0;
577  e5=e4*e0;
578  e6=e5*e0;
579  e7=e6*e0;
580  c2b=cos(2.*bet);
581  s2b=sin(2.*bet);
582  si=sin(inc);
583  ci=cos(inc);
584  s2i=si*si;
585  c2i=ci*ci;
586  P=FPlus;
587  Q=FCross;
588  p8=pow(f0, 19./18.);
589  return(p8*(-(3.*c2b*e0*P)/(2.) -(3.*c2b*c2i*e0*P)/(2.) -(3323.*c2b*e3*P)/(1216.) -(3323.*c2b*c2i*e3*P)/(1216.) -(15994231.*c2b*e5*P)/(4.435968e6) -(15994231.*c2b*c2i*e5*P)/(4.435968e6) -(105734339801.*c2b*e7*P)/(2.4273616896e10) -(105734339801.*c2b*c2i*e7*P)/(2.4273616896e10)+3.*ci*e0*Q*s2b +(3323.*ci*e3*Q*s2b)/(608.) +(15994231.*ci*e5*Q*s2b)/(2.217984e6) +(105734339801.*ci*e7*Q*s2b)/(1.2136808448e10)+e0*P*s2i +(3323.*e3*P*s2i)/(1824.) +(15994231.*e5*P*s2i)/(6.653952e6) +(105734339801.*e7*P*s2i)/(3.6410425344e10)) );
590 }
591 
592  /* case 2 of zeta_real*/
593 
594 static REAL8 UNUSED
595 z5(REAL8 UNUSED f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
596 {
597  REAL8 c2b,s2b,ci,P,Q,c2i;
598  c2b=cos(2.*bet);
599  s2b=sin(2.*bet);
600  ci=cos(inc);
601  c2i=ci*ci;
602  P=FPlus;
603  Q=FCross;
604  return(2.*c2b*P + 2.*c2b*c2i*P - 4.*ci*Q*s2b);
605 }
606 
607 static REAL8 UNUSED
608 z6(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
609 {
610  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p1,s2i,c2i;
611  e2=e0*e0;
612  e3=e2*e0;
613  e4=e3*e0;
614  e5=e4*e0;
615  e6=e5*e0;
616  e7=e6*e0;
617  e8=e7*e0;
618  c2b=cos(2.*bet);
619  s2b=sin(2.*bet);
620  si=sin(inc);
621  ci=cos(inc);
622  s2i=si*si;
623  c2i=ci*ci;
624  P=FPlus;
625  Q=FCross;
626  p1=pow(f0, 19./9.);
627  return(p1*(-(277.*c2b*e2*P)/(24.) -(277.*c2b*c2i*e2*P)/(24.) -(920471.*c2b*e4*P)/(21888.) -(920471.*c2b*c2i*e4*P)/(21888.) -(468070445.*c2b*e6*P)/(4.990464e6) -(468070445.*c2b*c2i*e6*P)/(4.990464e6) -(2295471547915.*c2b*e8*P)/(1.3653909504e10) -(2295471547915.*c2b*c2i*e8*P)/(1.3653909504e10)+e2*P*s2i +(3323.*e4*P*s2i)/(912.) +(1689785.*e6*P*s2i)/(207936.) +(8286900895.*e8*P*s2i)/(5.68912896e8)+(277.*ci*e2*Q*s2b)/(12.) +(920471.*ci*e4*Q*s2b)/(10944.) +(468070445.*ci*e6*Q*s2b)/(2.495232e6) +(2295471547915.*ci*e8*Q*s2b)/(6.826954752e9)));
628 }
629 
630 static REAL8 UNUSED
631 z7(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
632 {
633  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p2,s2i,c2i;
634  e2=e0*e0;
635  e3=e2*e0;
636  e4=e3*e0;
637  e5=e4*e0;
638  e6=e5*e0;
639  e7=e6*e0;
640  e8=e7*e0;
641  c2b=cos(2.*bet);
642  s2b=sin(2.*bet);
643  si=sin(inc);
644  ci=cos(inc);
645  s2i=si*si;
646  c2i=ci*ci;
647  P=FPlus;
648  Q=FCross;
649  p2=pow(f0, 19./3.);
650  return(p2*(-(103729904239.*c2b*e6*P)/(1.9961856e8) -(103729904239.*c2b*c2i*e6*P)/(1.9961856e8) -(344694471786197.*c2b*e8*P)/(6.068404224e10) -(344694471786197.*c2b*c2i*e8*P)/(6.068404224e10)+(103729904239.*ci*e6*Q*s2b)/(9.980928e7) +(344694471786197.*ci*e8*Q*s2b)/(3.034202112e10)+(29064841.*e6*P*s2i)/(554496.) +(96582466643.*e8*P*s2i)/(1.68566784e8)) );
651 }
652 
653 static REAL8 UNUSED
654 z8(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
655 {
656  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p3,s2i,c2i;
657  e2=e0*e0;
658  e3=e2*e0;
659  e4=e3*e0;
660  e5=e4*e0;
661  e6=e5*e0;
662  e7=e6*e0;
663  e8=e7*e0;
664  c2b=cos(2.*bet);
665  s2b=sin(2.*bet);
666  si=sin(inc);
667  ci=cos(inc);
668  s2i=si*si;
669  c2i=ci*ci;
670  P=FPlus;
671  Q=FCross;
672  p3 = pow (f0, 38./9.);
673  return(p3*((3254599.*c2b*e4*P)/(43776.) +(3254599.*c2b*c2i*e4*P)/(43776.) +(10815032477.*c2b*e6*P)/(1.9961856e7) +(10815032477.*c2b*c2i*e6*P)/(1.9961856e7) +(79934933490791.*c2b*e8*P)/(3.6410425344e10) +(79934933490791.*c2b*c2i*e8*P)/(3.6410425344e10)-(3254599.*ci*e4*Q*s2b)/(21888.) -(10815032477.*ci*e6*Q*s2b)/(9.980928e6) -(79934933490791.*ci*e8*Q*s2b)/(1.8205212672e10)-(3305.*e4*P*s2i)/(456.) -(10982515.*e6*P*s2i)/(207936.) -(81172812745.*e8*P*s2i)/(3.79275264e8)) );
674 }
675 
676 static REAL8 UNUSED
677 z9(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
678 {
679  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p4,s2i,c2i;
680  e2=e0*e0;
681  e3=e2*e0;
682  e4=e3*e0;
683  e5=e4*e0;
684  e6=e5*e0;
685  e7=e6*e0;
686  e8=e7*e0;
687  c2b=cos(2.*bet);
688  s2b=sin(2.*bet);
689  si=sin(inc);
690  ci=cos(inc);
691  s2i=si*si;
692  c2i=ci*ci;
693  P=FPlus;
694  Q=FCross;
695  p4 = pow (f0, 76./9.);
696  return(p4*((8340362777769439.*c2b*e8*P)/(2.18462552064e12) +(8340362777769439.*c2b*c2i*e8*P)/(2.18462552064e12)-(8340362777769439.*ci*e8*Q*s2b)/(1.09231276032e12)-(4442498396267.*e8*P*s2i)/(1.137825792e10)) );
697 }
698 
699  /* case 3 of zeta_real*/
700 
701 static REAL8 UNUSED
702 z10(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
703 {
704  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p5,s2i,c2i;
705  e2=e0*e0;
706  e3=e2*e0;
707  e4=e3*e0;
708  e5=e4*e0;
709  e6=e5*e0;
710  e7=e6*e0;
711  c2b=cos(2.*bet);
712  s2b=sin(2.*bet);
713  si=sin(inc);
714  ci=cos(inc);
715  s2i=si*si;
716  c2i=ci*ci;
717  P=FPlus;
718  Q=FCross;
719  p5 = pow (f0, 133./18.);
720  return(p5*((-74006050878931.*c2b*e7*P)/(4.045602816e10) -(74006050878931.*c2b*c2i*e7*P)/(4.045602816e10)+(74006050878931.*ci*e7*Q*s2b)/(2.022801408e10)+(2531702819.*e7*P*s2i)/(2.957312e7)) );
721 }
722 
723 static REAL8 UNUSED
724 z11(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
725 {
726  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p6,s2i,c2i;
727  e2=e0*e0;
728  e3=e2*e0;
729  e4=e3*e0;
730  e5=e4*e0;
731  e6=e5*e0;
732  e7=e6*e0;
733  c2b=cos(2.*bet);
734  s2b=sin(2.*bet);
735  si=sin(inc);
736  ci=cos(inc);
737  s2i=si*si;
738  c2i=ci*ci;
739  P=FPlus;
740  Q=FCross;
741  p6 = pow (f0, 95./18.);
742  return(p6*((1810486747.*c2b*e5*P)/(7.39328e6) +(1810486747.*c2b*c2i*e5*P)/(7.39328e6) +(6016247460281.*c2b*e7*P)/(2.697068544e9) +(6016247460281.*c2b*c2i*e7*P)/(2.697068544e9)-(1810486747.*ci*e5*Q*s2b)/(3.69664e6) -(6016247460281.*ci*e7*Q*s2b)/(1.348534272e9)-(50883.*e5*P*s2i)/(4864.) -(281807015.*e7*P*s2i)/(2.957312e6)) );
743 }
744 
745 static REAL8 UNUSED
746 z12(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
747 {
748  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p7,s2i,c2i;
749  e2=e0*e0;
750  e3=e2*e0;
751  e4=e3*e0;
752  e5=e4*e0;
753  e6=e5*e0;
754  e7=e6*e0;
755  c2b=cos(2.*bet);
756  s2b=sin(2.*bet);
757  si=sin(inc);
758  ci=cos(inc);
759  s2i=si*si;
760  c2i=ci*ci;
761  P=FPlus;
762  Q=FCross;
763  p7 = pow (f0, 19./6.);
764  return(p7*(-(40863.*c2b*e3*P)/(1216.) -(40863.*c2b*c2i*e3*P)/(1216.) -(135787749.*c2b*e5*P)/(739328.) -(135787749.*c2b*c2i*e5*P)/(739328.) -(518672547069.*c2b*e7*P)/(8.99022848e8) -(518672547069.*c2b*c2i*e7*P)/(8.99022848e8)+(40863.*ci*e3*Q*s2b)/(608.) +(135787749.*ci*e5*Q*s2b)/(369664.) +(518672547069.*ci*e7*Q*s2b)/(4.49511424e8)+(9.*e3*P*s2i)/(8.) +(29907.*e5*P*s2i)/(4864.) +(114236667.*e7*P*s2i)/(5.914624e6)) );
765 }
766 
767 static REAL8 UNUSED
768 z13(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
769 {
770  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p8,c2i;
771  e2=e0*e0;
772  e3=e2*e0;
773  e4=e3*e0;
774  e5=e4*e0;
775  e6=e5*e0;
776  e7=e6*e0;
777  c2b=cos(2.*bet);
778  s2b=sin(2.*bet);
779  ci=cos(inc);
780  c2i=ci*ci;
781  P=FPlus;
782  Q=FCross;
783  p8 = pow (f0, 19./18.);
784  return(p8*((9.*c2b*e0*P)/(2.) +(9.*c2b*c2i*e0*P)/(2.) +(9969.*c2b*e3*P)/(1216.) +(9969.*c2b*c2i*e3*P)/(1216.) +(15994231.*c2b*e5*P)/(1.478656e6) +(15994231.*c2b*c2i*e5*P)/(1.478656e6) +(105734339801.*c2b*e7*P)/(8.091205632e9) +(105734339801.*c2b*c2i*e7*P)/(8.091205632e9)-9.*ci*e0*Q*s2b -(9969.*ci*e3*Q*s2b)/(608.) -(15994231.*ci*e5*Q*s2b)/(739328.) -(105734339801.*ci*e7*Q*s2b)/(4.045602816e9)) );
785 }
786 
787  /*case 4 for z_real*/
788 
789 static REAL8 UNUSED
790 z14(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
791 {
792  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p1,c2i;
793  e2=e0*e0;
794  e3=e2*e0;
795  e4=e3*e0;
796  e5=e4*e0;
797  e6=e5*e0;
798  e7=e6*e0;
799  e8=e7*e0;
800  c2b=cos(2.*bet);
801  s2b=sin(2.*bet);
802  ci=cos(inc);
803  c2i=ci*ci;
804  P=FPlus;
805  Q=FCross;
806  p1 = pow (f0, 19./9.);
807  return(p1*(8.*c2b*e2*P +8.*c2b*c2i*e2*P +(3323.*c2b*e4*P)/(114.) +(3323.*c2b*c2i*e4*P)/(114.) +(1689785.*c2b*e6*P)/(25992.) +(1689785.*c2b*c2i*e6*P)/(25992.) +(8286900895.*c2b*e8*P)/(7.1114112e7) +(8286900895.*c2b*c2i*e8*P)/(7.1114112e7)-16.*ci*e2*Q*s2b -(3323.*ci*e4*Q*s2b)/(57.) -(1689785.*ci*e6*Q*s2b)/(12996.) -(8286900895.*ci*e8*Q*s2b)/(3.5557056e7) ) );
808 }
809 
810 static REAL8 UNUSED
811 z15(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
812 {
813  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p2,s2i,c2i;
814  e2=e0*e0;
815  e3=e2*e0;
816  e4=e3*e0;
817  e5=e4*e0;
818  e6=e5*e0;
819  e7=e6*e0;
820  e8=e7*e0;
821  c2b=cos(2.*bet);
822  s2b=sin(2.*bet);
823  si=sin(inc);
824  ci=cos(inc);
825  s2i=si*si;
826  c2i=ci*ci;
827  P=FPlus;
828  Q=FCross;
829  p2 = pow (f0, 19./3.);
830  return(p2*((71474367.*c2b*e6*P)/(115520.) +(71474367.*c2b*c2i*e6*P)/(115520.) +(237509321541.*c2b*e8*P)/(3.511808e7) +(237509321541.*c2b*c2i*e8*P)/(3.511808e7)-(71474367.*ci*e6*Q*s2b)/(57760.) -(237509321541.*ci*e8*Q*s2b)/(1.755904e7)-(51793.*e6*P*s2i)/(3420.) -(172108139.*e8*P*s2i)/(1.03968e6)) );
831 }
832 
833 static REAL8 UNUSED
834 z16(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
835 {
836  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p3,s2i,c2i;
837  e2=e0*e0;
838  e3=e2*e0;
839  e4=e3*e0;
840  e5=e4*e0;
841  e6=e5*e0;
842  e7=e6*e0;
843  e8=e7*e0;
844  c2b=cos(2.*bet);
845  s2b=sin(2.*bet);
846  si=sin(inc);
847  ci=cos(inc);
848  s2i=si*si;
849  c2i=ci*ci;
850  P=FPlus;
851  Q=FCross;
852  p3 = pow (f0, 38./9.);
853  return(p3*(-(1431.*c2b*e4*P)/(19.) -(1431.*c2b*c2i*e4*P)/(19.) -(1585071.*c2b*e6*P)/(2888.) -(1585071.*c2b*c2i*e6*P)/(2888.) -(3905136831.*c2b*e8*P)/(1.755904e6) -(3905136831.*c2b*c2i*e8*P)/(1.755904e6)+(2862.*ci*e4*Q*s2b)/(19.) +(1585071.*ci*e6*Q*s2b)/(1444.) +(3905136831.*ci*e8*Q*s2b)/(877952.)+(4.*e4*P*s2i)/(3.) +(3323.*e6*P*s2i)/(342.) +(24560609.*e8*P*s2i)/(623808.)) );
854 }
855 
856 static REAL8 UNUSED
857 z17(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
858 {
859  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p4,s2i,c2i;
860  e2=e0*e0;
861  e3=e2*e0;
862  e4=e3*e0;
863  e5=e4*e0;
864  e6=e5*e0;
865  e7=e6*e0;
866  e8=e7*e0;
867  c2b=cos(2.*bet);
868  s2b=sin(2.*bet);
869  si=sin(inc);
870  ci=cos(inc);
871  s2i=si*si;
872  c2i=ci*ci;
873  P=FPlus;
874  Q=FCross;
875  p4 = pow (f0, 76./9.);
876  return(p4*((-49436673769121.*c2b*e8*P)/(9.95597568e9) -(49436673769121.*c2b*c2i*e8*P)/(9.95597568e9)+(865873231.*e8*P*s2i)/(6.23808e6)+(49436673769121.*ci*e8*Q*s2b)/(4.97798784e9)) );
877 }
878 
879  /*case 5 for zeta_real*/
880 
881 static REAL8 UNUSED
882 z18(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
883 {
884  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p5,s2i,c2i;
885  e2=e0*e0;
886  e3=e2*e0;
887  e4=e3*e0;
888  e5=e4*e0;
889  e6=e5*e0;
890  e7=e6*e0;
891  c2b=cos(2.*bet);
892  s2b=sin(2.*bet);
893  si=sin(inc);
894  ci=cos(inc);
895  s2i=si*si;
896  c2i=ci*ci;
897  P=FPlus;
898  Q=FCross;
899  p5 = pow (f0, 133./18.);
900  return(p5*((3053741715625.*c2b*e7*P)/(2.235727872e9) +(3053741715625.*c2b*c2i*e7*P)/(2.235727872e9)-(3053741715625.*ci*e7*Q*s2b)/(1.117863936e9)-(15300625.*e7*P*s2i)/(700416.)) );
901 }
902 
903 static REAL8 UNUSED
904 z19(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
905 {
906  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p6,s2i,c2i;
907  e2=e0*e0;
908  e3=e2*e0;
909  e4=e3*e0;
910  e5=e4*e0;
911  e6=e5*e0;
912  e7=e6*e0;
913  c2b=cos(2.*bet);
914  s2b=sin(2.*bet);
915  si=sin(inc);
916  ci=cos(inc);
917  s2i=si*si;
918  c2i=ci*ci;
919  P=FPlus;
920  Q=FCross;
921  p6 = pow (f0, 95./18.);
922  return(p6*(-(13023125.*c2b*e5*P)/(87552.) -(13023125.*c2b*c2i*e5*P)/(87552.) -(216379221875.*c2b*e7*P)/(1.59694848e8) -(216379221875.*c2b*c2i*e7*P)/(1.59694848e8)+(13023125.*ci*e5*Q*s2b)/(43776.) +(216379221875.*ci*e7*Q*s2b)/(7.9847424e7)+(625.*e5*P*s2i)/(384.) + (10384375.*e7*P*s2i)/(700416.)) );
923 }
924 
925 static REAL8 UNUSED
926 z20(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
927 {
928  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p7,c2i;
929  e2=e0*e0;
930  e3=e2*e0;
931  e4=e3*e0;
932  e5=e4*e0;
933  e6=e5*e0;
934  e7=e6*e0;
935  c2b=cos(2.*bet);
936  s2b=sin(2.*bet);
937  ci=cos(inc);
938  c2i=ci*ci;
939  P=FPlus;
940  Q=FCross;
941  p7 = pow (f0, 19./6.);
942  return(p7*((625.*c2b*e3*P)/(48.) +(625.*c2b*c2i*e3*P)/(48.) +(2076875.*c2b*e5*P)/(29184.) +(2076875.*c2b*c2i*e5*P)/(29184.) +(7933101875.*c2b*e7*P)/(3.5487744e7) +(7933101875.*c2b*c2i*e7*P)/(3.5487744e7)-(625.*ci*e3*Q*s2b)/(24.) -(2076875.*ci*e5*Q*s2b)/(14592.) -(7933101875.*ci*e7*Q*s2b)/(1.7743872e7)) );
943 }
944 
945 /* case 6 for zeta_real*/
946 
947 static REAL8 UNUSED
948 z21(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
949 {
950  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p2,s2i,c2i;
951  e2=e0*e0;
952  e3=e2*e0;
953  e4=e3*e0;
954  e5=e4*e0;
955  e6=e5*e0;
956  e7=e6*e0;
957  e8=e7*e0;
958  c2b=cos(2.*bet);
959  s2b=sin(2.*bet);
960  si=sin(inc);
961  ci=cos(inc);
962  s2i=si*si;
963  c2i=ci*ci;
964  P=FPlus;
965  Q=FCross;
966  p2 = pow (f0, 19./3.);
967  return(p2*(-(1656963.*c2b*e6*P)/(6080.) -(1656963.*c2b*c2i*e6*P)/(6080.) -(5506088049.*c2b*e8*P)/(1.84832e6) -(5506088049.*c2b*c2i*e8*P)/(1.84832e6)+(1656963.*ci*e6*Q*s2b)/(3040.) +(5506088049.*ci*e8*Q*s2b)/(924160.)+(81.*e6*P*s2i)/(40.) +(269163.*e8*P*s2i)/(12160.)) );
968 }
969 
970 static REAL8 UNUSED
971 z22(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
972 {
973  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p3,c2i;
974  e2=e0*e0;
975  e3=e2*e0;
976  e4=e3*e0;
977  e5=e4*e0;
978  e6=e5*e0;
979  e7=e6*e0;
980  e8=e7*e0;
981  c2b=cos(2.*bet);
982  s2b=sin(2.*bet);
983  ci=cos(inc);
984  c2i=ci*ci;
985  P=FPlus;
986  Q=FCross;
987  p3 = pow (f0, 38./9.);
988  return(p3*((81.*c2b*e4*P)/(4.) +(81.*c2b*c2i*e4*P)/(4.) +(89721.*c2b*e6*P)/(608.) +(89721.*c2b*c2i*e6*P)/(608.) +(221045481.*c2b*e8*P)/(369664.) +(221045481.*c2b*c2i*e8*P)/(369664.)-(81.*ci*e4*Q*s2b)/(2.) -(89721.*ci*e6*Q*s2b)/(304.) -(221045481.*ci*e8*Q*s2b)/(184832.)) );
989 }
990 
991 static REAL8 UNUSED
992 z23(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
993 {
994  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p4,s2i,c2i;
995  e2=e0*e0;
996  e3=e2*e0;
997  e4=e3*e0;
998  e5=e4*e0;
999  e6=e5*e0;
1000  e7=e6*e0;
1001  e8=e7*e0;
1002  c2b=cos(2.*bet);
1003  s2b=sin(2.*bet);
1004  si=sin(inc);
1005  ci=cos(inc);
1006  s2i=si*si;
1007  c2i=ci*ci;
1008  P=FPlus;
1009  Q=FCross;
1010  p4 = pow (f0, 76./9.);
1011  return(p4*( (71730555921.*c2b*e8*P)/(2.587648e7) +(71730555921.*c2b*c2i*e8*P)/(2.587648e7)-(333693.*e8*P*s2i)/(10640.)-(71730555921.*ci*e8*Q*s2b)/(1.293824e7)) );
1012 }
1013 
1014  /* case 7 for zeta_real*/
1015 
1016 static REAL8 UNUSED
1017 z24(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1018 {
1019  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,si,ci,P,Q,p5,s2i,c2i;
1020  e2=e0*e0;
1021  e3=e2*e0;
1022  e4=e3*e0;
1023  e5=e4*e0;
1024  e6=e5*e0;
1025  e7=e6*e0;
1026  c2b=cos(2.*bet);
1027  s2b=sin(2.*bet);
1028  si=sin(inc);
1029  ci=cos(inc);
1030  s2i=si*si;
1031  c2i=ci*ci;
1032  P=FPlus;
1033  Q=FCross;
1034  p5 = pow (f0, 133./18.);
1035  return(p5*((-1109077123.*c2b*e7*P)/(2.33472e6) -(1109077123.*c2b*c2i*e7*P)/(2.33472e6)+(1109077123.*ci*e7*Q*s2b)/(1.16736e6)+(117649.*e7*P*s2i)/(46080.) ) );
1036 }
1037 
1038 static REAL8 UNUSED
1039 z25(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1040 {
1041  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p6,c2i;
1042  e2=e0*e0;
1043  e3=e2*e0;
1044  e4=e3*e0;
1045  e5=e4*e0;
1046  e6=e5*e0;
1047  e7=e6*e0;
1048  c2b=cos(2.*bet);
1049  s2b=sin(2.*bet);
1050  ci=cos(inc);
1051  c2i=ci*ci;
1052  P=FPlus;
1053  Q=FCross;
1054  p6 = pow (f0, 95./18.);
1055  return(p6*((117649.*c2b*e5*P)/(3840.) +(117649.*c2b*c2i*e5*P)/(3840.) +(390947627.*c2b*e7*P)/(1.400832e6) +(390947627.*c2b*c2i*e7*P)/(1.400832e6)-(117649.*ci*e5*Q*s2b)/(1920.) -(390947627.*ci*e7*Q*s2b)/(700416.)) );
1056 }
1057 
1058  /* case 8 for zeta_real*/
1059 
1060 static REAL8 UNUSED
1061 z26(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1062 {
1063  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p2,c2i;
1064  e2=e0*e0;
1065  e3=e2*e0;
1066  e4=e3*e0;
1067  e5=e4*e0;
1068  e6=e5*e0;
1069  e7=e6*e0;
1070  e8=e7*e0;
1071  c2b=cos(2.*bet);
1072  s2b=sin(2.*bet);
1073  ci=cos(inc);
1074  c2i=ci*ci;
1075  P=FPlus;
1076  Q=FCross;
1077  p2 = pow (f0, 19./3.);
1078  return(p2*((2048.*c2b*e6*P)/(45.) +(2048.*c2b*c2i*e6*P)/(45.) +(425344.*c2b*e8*P)/(855.) +(425344.*c2b*c2i*e8*P)/(855.)-(4096.*ci*e6*Q*s2b)/(45.) -(850688.*ci*e8*Q*s2b)/(855.)) );
1079 }
1080 
1081 
1082 static REAL8 UNUSED
1083 z27(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1084 {
1085  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,si,ci,P,Q,p4,s2i,c2i;
1086  e2=e0*e0;
1087  e3=e2*e0;
1088  e4=e3*e0;
1089  e5=e4*e0;
1090  e6=e5*e0;
1091  e7=e6*e0;
1092  e8=e7*e0;
1093  c2b=cos(2.*bet);
1094  s2b=sin(2.*bet);
1095  si=sin(inc);
1096  ci=cos(inc);
1097  s2i=si*si;
1098  c2i=ci*ci;
1099  P=FPlus;
1100  Q=FCross;
1101  p4 = pow (f0, 76./9.);
1102  return(p4*((-14348288.*c2b*e8*P)/(17955.) -(14348288.*c2b*c2i*e8*P)/(17955.)+(28696576.*ci*e8*Q*s2b)/(17955.)+(1024.*e8*P*s2i)/(315.)) );
1103 }
1104 
1105  /* case 9 for zeta_real*/
1106 
1107 static REAL8 UNUSED
1108 z28(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1109 {
1110  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p5,c2i;
1111  e2=e0*e0;
1112  e3=e2*e0;
1113  e4=e3*e0;
1114  e5=e4*e0;
1115  e6=e5*e0;
1116  e7=e6*e0;
1117  c2b=cos(2.*bet);
1118  s2b=sin(2.*bet);
1119  ci=cos(inc);
1120  c2i=ci*ci;
1121  P=FPlus;
1122  Q=FCross;
1123  p5 = pow (f0, 133./18.);
1124  return(p5*((4782969.*c2b*e7*P)/(71680.) +(4782969.*c2b*c2i*e7*P)/(71680.) -(4782969.*ci*e7*Q*s2b)/(35840.)) );
1125 }
1126 
1127  /*case 10 for zeta_real*/
1128 
1129 static REAL8 UNUSED
1130 z29(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1131 {
1132  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p4,c2i;
1133  e2=e0*e0;
1134  e3=e2*e0;
1135  e4=e3*e0;
1136  e5=e4*e0;
1137  e6=e5*e0;
1138  e7=e6*e0;
1139  e8=e7*e0;
1140  c2b=cos(2.*bet);
1141  s2b=sin(2.*bet);
1142  ci=cos(inc);
1143  c2i=ci*ci;
1144  P=FPlus;
1145  Q=FCross;
1146  p4 = pow (f0, 76./9.);
1147  return((390625.*e8*p4*(c2b*P + c2b*c2i*P - 2.*ci*Q*s2b))/(4032.) );
1148 }
1149 
1150  /*The following are the coefficients for the imaginary part of the zeta_l's*/
1151 
1152  /*case 1 for zeta_im*/
1153 
1154 static REAL8 UNUSED
1155 q1(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1156 {
1157  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p5,c2i;
1158  e2=e0*e0;
1159  e3=e2*e0;
1160  e4=e3*e0;
1161  e5=e4*e0;
1162  e6=e5*e0;
1163  e7=e6*e0;
1164  c2b=cos(2.*bet);
1165  s2b=sin(2.*bet);
1166  ci=cos(inc);
1167  c2i=ci*ci;
1168  P=FPlus;
1169  Q=FCross;
1170  p5 = pow (f0, 133./18.);
1171  return(p5*((16099508821573.*c2b*ci*e7*Q)/(2.022801408e10)+(16099508821573.*e7*P*s2b)/(4.045602816e10) +(16099508821573.*c2i*e7*P*s2b)/(4.045602816e10) ) );
1172 }
1173 
1174 static REAL8 UNUSED
1175 q2(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1176 {
1177  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p6,c2i;
1178  e2=e0*e0;
1179  e3=e2*e0;
1180  e4=e3*e0;
1181  e5=e4*e0;
1182  e6=e5*e0;
1183  e7=e6*e0;
1184  c2b=cos(2.*bet);
1185  s2b=sin(2.*bet);
1186  ci=cos(inc);
1187  c2i=ci*ci;
1188  P=FPlus;
1189  Q=FCross;
1190  p6 = pow (f0, 95./18.);
1191  return(p6*(-(749695861.*c2b*ci*e5*Q)/(6.653952e6) -(12456196730515.*c2b*ci*e7*Q)/(1.2136808448e10)-(749695861.*e5*P*s2b)/(1.3307904e7) -(749695861.*c2i*e5*P*s2b)/(1.3307904e7) -(12456196730515.*e7*P*s2b)/(2.4273616896e10) -(12456196730515.*c2i*e7*P*s2b)/(2.4273616896e10) ) );
1192 }
1193 
1194 static REAL8 UNUSED
1195 q3(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1196 {
1197  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p7,c2i;
1198  e2=e0*e0;
1199  e3=e2*e0;
1200  e4=e3*e0;
1201  e5=e4*e0;
1202  e6=e5*e0;
1203  e7=e6*e0;
1204  c2b=cos(2.*bet);
1205  s2b=sin(2.*bet);
1206  ci=cos(inc);
1207  c2i=ci*ci;
1208  P=FPlus;
1209  Q=FCross;
1210  p7 = pow (f0, 19./6.);
1211  return(p7*((31363.*c2b*ci*e3*Q)/(1824.) +(104219249.*c2b*ci*e5*Q)/(1.108992e6) +(398089398569.*c2b*ci*e7*Q)/(1.348534272e9)+(31363.*e3*P*s2b)/(3648.) +(31363.*c2i*e3*P*s2b)/(3648.) +(104219249.*e5*P*s2b)/(2.217984e6) +(104219249.*c2i*e5*P*s2b)/(2.217984e6) +(398089398569.*e7*P*s2b)/(2.697068544e9) +(398089398569.*c2i*e7*P*s2b)/(2.697068544e9) ) );
1212 }
1213 
1214 static REAL8 UNUSED
1215 q4(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1216 {
1217  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p8,c2i;
1218  e2=e0*e0;
1219  e3=e2*e0;
1220  e4=e3*e0;
1221  e5=e4*e0;
1222  e6=e5*e0;
1223  e7=e6*e0;
1224  c2b=cos(2.*bet);
1225  s2b=sin(2.*bet);
1226  ci=cos(inc);
1227  c2i=ci*ci;
1228  P=FPlus;
1229  Q=FCross;
1230  p8 = pow (f0, 19./18.);
1231  return(p8*(-3.*c2b*ci*e0*Q -(3323.*c2b*ci*e3*Q)/(608.) -(15994231.*c2b*ci*e5*Q)/(2.217984e6) -(105734339801.*c2b*ci*e7*Q)/(1.2136808448e10)-(3.*e0*P*s2b)/2. -(3.*c2i*e0*P*s2b)/2. -(3323.*e3*P*s2b)/1216. -(3323.*c2i*e3*P*s2b)/(1216.) -(15994231.*e5*P*s2b)/(4.435968e6) -(15994231.*c2i*e5*P*s2b)/(4.435968e6) -(105734339801.*e7*P*s2b)/(2.4273616896e10) -(105734339801.*c2i*e7*P*s2b)/(2.4273616896e10)) );
1232 }
1233 
1234  /* case 2 for zeta_im */
1235 
1236 static REAL8 UNUSED
1237 q5(REAL8 UNUSED f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1238 {
1239  REAL8 c2b,s2b,ci,c2i,P,Q;
1240  c2b=cos(2.*bet);
1241  s2b=sin(2.*bet);
1242  ci=cos(inc);
1243  c2i=ci*ci;
1244  P=FPlus;
1245  Q=FCross;
1246  return(4.*c2b*ci*Q + 2.*P*s2b + 2.*c2i*P*s2b);
1247 }
1248 
1249 static REAL8 UNUSED
1250 q6(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1251 {
1252  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p1,c2i;
1253  e2=e0*e0;
1254  e3=e2*e0;
1255  e4=e3*e0;
1256  e5=e4*e0;
1257  e6=e5*e0;
1258  e7=e6*e0;
1259  e8=e7*e0;
1260  c2b=cos(2.*bet);
1261  s2b=sin(2.*bet);
1262  ci=cos(inc);
1263  c2i=ci*ci;
1264  P=FPlus;
1265  Q=FCross;
1266  p1 = pow (f0, 19./9.);
1267  return(p1*(-(277.*c2b*ci*e2*Q)/(12.) -(920471.*c2b*ci*e4*Q)/(10944.) -(468070445.*c2b*ci*e6*Q)/(2.495232e6) -(2295471547915.*c2b*ci*e8*Q)/(6.826954752e9)-(277.*e2*P*s2b)/(24.) -(277.*c2i*e2*P*s2b)/(24.) -(920471.*e4*P*s2b)/(21888.) -(920471.*c2i*e4*P*s2b)/(21888.) -(468070445.*e6*P*s2b)/(4.990464e6) -(468070445.*c2i*e6*P*s2b)/(4.990464e6) -(2295471547915.*e8*P*s2b)/(1.3653909504e10) -(2295471547915.*c2i*e8*P*s2b)/(1.3653909504e10) ) );
1268 }
1269 
1270 static REAL8 UNUSED
1271 q7(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1272 {
1273  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p2,c2i;
1274  e2=e0*e0;
1275  e3=e2*e0;
1276  e4=e3*e0;
1277  e5=e4*e0;
1278  e6=e5*e0;
1279  e7=e6*e0;
1280  e8=e7*e0;
1281  c2b=cos(2.*bet);
1282  s2b=sin(2.*bet);
1283  ci=cos(inc);
1284  c2i=ci*ci;
1285  P=FPlus;
1286  Q=FCross;
1287  p2 = pow (f0, 19./3.);
1288  return(p2*(-(104238504751.*c2b*ci*e6*Q)/(9.980928e7) -(346384551287573.*c2b*ci*e8*Q)/(3.034202112e10)-(104238504751.*e6*P*s2b)/(1.9961856e8) -(104238504751.*c2i*e6*P*s2b)/(1.9961856e8) -(346384551287573.*e8*P*s2b)/(6.068404224e10) -(346384551287573.*c2i*e8*P*s2b)/(6.068404224e10)) );
1289 }
1290 
1291 static REAL8 UNUSED
1292 q8(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1293 {
1294  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p3,c2i;
1295  e2=e0*e0;
1296  e3=e2*e0;
1297  e4=e3*e0;
1298  e5=e4*e0;
1299  e6=e5*e0;
1300  e7=e6*e0;
1301  e8=e7*e0;
1302  c2b=cos(2.*bet);
1303  s2b=sin(2.*bet);
1304  ci=cos(inc);
1305  c2i=ci*ci;
1306  P=FPlus;
1307  Q=FCross;
1308  p3 = pow (f0, 38./9.);
1309  return(p3*((3265543.*c2b*ci*e4*Q)/(21888.) +(10851399389.*c2b*ci*e6*Q)/(9.980928e6) +(80203724795687.*c2b*ci*e8*Q)/(1.8205212672e10)+(3265543.*e4*P*s2b)/(43776.) +(3265543.*c2i*e4*P*s2b)/(43776.) +(10851399389.*e6*P*s2b)/(1.9961856e7) +(10851399389.*c2i*e6*P*s2b)/(1.9961856e7) +(80203724795687.*e8*P*s2b)/(3.6410425344e10) +(80203724795687.*c2i*e8*P*s2b)/(3.6410425344e10)) );
1310 }
1311 
1312 static REAL8 UNUSED
1313 q9(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1314 {
1315  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p4,c2i;
1316  e2=e0*e0;
1317  e3=e2*e0;
1318  e4=e3*e0;
1319  e5=e4*e0;
1320  e6=e5*e0;
1321  e7=e6*e0;
1322  e8=e7*e0;
1323  c2b=cos(2.*bet);
1324  s2b=sin(2.*bet);
1325  ci=cos(inc);
1326  c2i=ci*ci;
1327  P=FPlus;
1328  Q=FCross;
1329  p4 = pow (f0, 76./9.);
1330  return(p4*((8388155956587871.*e8*P*s2b)/(2.18462552064e12) +(8388155956587871.*c2i*e8*P*s2b)/(2.18462552064e12)+ (8388155956587871.*c2b*ci*e8*Q)/(1.09231276032e12)) );
1331 }
1332 
1333  /*case 3 for zeta_im */
1334 
1335 static REAL8 UNUSED
1336 q10(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1337 {
1338  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p5,c2i;
1339  e2=e0*e0;
1340  e3=e2*e0;
1341  e4=e3*e0;
1342  e5=e4*e0;
1343  e6=e5*e0;
1344  e7=e6*e0;
1345  c2b=cos(2.*bet);
1346  s2b=sin(2.*bet);
1347  ci=cos(inc);
1348  c2i=ci*ci;
1349  P=FPlus;
1350  Q=FCross;
1351  p5 = pow (f0, 133./18.);
1352  return(p5*((-74123294656819.*c2b*ci*e7*Q)/(2.022801408e10)-(74123294656819.*e7*P*s2b)/(4.045602816e10) -(74123294656819.*c2i*e7*P*s2b)/(4.045602816e10)) );
1353 }
1354 
1355 static REAL8 UNUSED
1356 q11(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1357 {
1358  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p6,c2i;
1359  e2=e0*e0;
1360  e3=e2*e0;
1361  e4=e3*e0;
1362  e5=e4*e0;
1363  e6=e5*e0;
1364  e7=e6*e0;
1365  c2b=cos(2.*bet);
1366  s2b=sin(2.*bet);
1367  ci=cos(inc);
1368  c2i=ci*ci;
1369  P=FPlus;
1370  Q=FCross;
1371  p6 = pow (f0, 95./18.);
1372  return(p6*( (1812254203.*c2b*ci*e5*Q)/(3.69664e6) +(6022120716569.*c2b*ci*e7*Q)/(1.348534272e9)+(1812254203.*e5*P*s2b)/(7.39328e6) +(1812254203.*c2i*e5*P*s2b)/(7.39328e6) +(6022120716569.*e7*P*s2b)/(2.697068544e9) +(6022120716569.*c2i*e7*P*s2b)/(2.697068544e9)) );
1373 }
1374 
1375 static REAL8 UNUSED
1376 q12(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1377 {
1378  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p7,c2i;
1379  e2=e0*e0;
1380  e3=e2*e0;
1381  e4=e3*e0;
1382  e5=e4*e0;
1383  e6=e5*e0;
1384  e7=e6*e0;
1385  c2b=cos(2.*bet);
1386  s2b=sin(2.*bet);
1387  ci=cos(inc);
1388  c2i=ci*ci;
1389  P=FPlus;
1390  Q=FCross;
1391  p7 = pow (f0, 19./6.);
1392  return(p7*( -(40863.*c2b*ci*e3*Q)/(608.) -(135787749.*c2b*ci*e5*Q)/(369664.) -(518672547069.*c2b*ci*e7*Q)/(4.49511424e8)-(40863.*e3*P*s2b)/(1216.) -(40863.*c2i*e3*P*s2b)/(1216.) -(135787749.*e5*P*s2b)/(739328.) -(135787749.*c2i*e5*P*s2b)/(739328.) -(518672547069.*e7*P*s2b)/(8.99022848e8) -(518672547069.*c2i*e7*P*s2b)/(8.99022848e8)) );
1393 }
1394 
1395 static REAL8 UNUSED
1396 q13(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1397 {
1398  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p8,c2i;
1399  e2=e0*e0;
1400  e3=e2*e0;
1401  e4=e3*e0;
1402  e5=e4*e0;
1403  e6=e5*e0;
1404  e7=e6*e0;
1405  c2b=cos(2.*bet);
1406  s2b=sin(2.*bet);
1407  ci=cos(inc);
1408  c2i=ci*ci;
1409  P=FPlus;
1410  Q=FCross;
1411  p8 = pow (f0, 19./18.);
1412  return(p8*(9.*c2b*ci*e0*Q +(9969.*c2b*ci*e3*Q)/(608.) +(15994231.*c2b*ci*e5*Q)/(739328.) +(105734339801.*c2b*ci*e7*Q)/(4.045602816e9)+(9.*e0*P*s2b)/(2.) +(9.*c2i*e0*P*s2b)/(2.) +(9969.*e3*P*s2b)/(1216.) +(9969.*c2i*e3*P*s2b)/(1216.) +(15994231.*e5*P*s2b)/(1.478656e6) +(15994231.*c2i*e5*P*s2b)/(1.478656e6) +(105734339801.*e7*P*s2b)/(8.091205632e9) +(105734339801.*c2i*e7*P*s2b)/(8.091205632e9) ) );
1413 }
1414 
1415  /* case 4 for zeta_im*/
1416 
1417 
1418 static REAL8 UNUSED
1419 q14(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1420 {
1421  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p1,c2i;
1422  e2=e0*e0;
1423  e3=e2*e0;
1424  e4=e3*e0;
1425  e5=e4*e0;
1426  e6=e5*e0;
1427  e7=e6*e0;
1428  e8=e7*e0;
1429  c2b=cos(2.*bet);
1430  s2b=sin(2.*bet);
1431  ci=cos(inc);
1432  c2i=ci*ci;
1433  P=FPlus;
1434  Q=FCross;
1435  p1 = pow (f0, 19./9.);
1436  return(p1*( 16.*c2b*ci*e2*Q +(3323.*c2b*ci*e4*Q)/(57.) +(1689785.*c2b*ci*e6*Q)/(12996.) +(8286900895.*c2b*ci*e8*Q)/(3.5557056e7)+8.*e2*P*s2b +8.*c2i*e2*P*s2b +(3323.*e4*P*s2b)/(114.) +(3323.*c2i*e4*P*s2b)/(114.) +(1689785.*e6*P*s2b)/(25992.) +(1689785.*c2i*e6*P*s2b)/(25992.) +(8286900895.*e8*P*s2b)/(7.1114112e7) +(8286900895.*c2i*e8*P*s2b)/(7.1114112e7)) );
1437 }
1438 
1439 static REAL8 UNUSED
1440 q15(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1441 {
1442  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p2,c2i;
1443  e2=e0*e0;
1444  e3=e2*e0;
1445  e4=e3*e0;
1446  e5=e4*e0;
1447  e6=e5*e0;
1448  e7=e6*e0;
1449  e8=e7*e0;
1450  c2b=cos(2.*bet);
1451  s2b=sin(2.*bet);
1452  ci=cos(inc);
1453  c2i=ci*ci;
1454  P=FPlus;
1455  Q=FCross;
1456  p2 = pow (f0, 19./3.);
1457  return(p2*((643523447.*c2b*ci*e6*Q)/(519840.) +(2138428414381.*c2b*ci*e8*Q)/(1.5803136e8)+(643523447.*e6*P*s2b)/(1.03968e6) +(643523447.*c2i*e6*P*s2b)/(1.03968e6) +(2138428414381.*e8*P*s2b)/(3.1606272e8) +(2138428414381.*c2i*e8*P*s2b)/(3.1606272e8) ) );
1458 }
1459 
1460 static REAL8 UNUSED
1461 q16(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1462 {
1463  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p3,c2i;
1464  e2=e0*e0;
1465  e3=e2*e0;
1466  e4=e3*e0;
1467  e5=e4*e0;
1468  e6=e5*e0;
1469  e7=e6*e0;
1470  e8=e7*e0;
1471  c2b=cos(2.*bet);
1472  s2b=sin(2.*bet);
1473  ci=cos(inc);
1474  c2i=ci*ci;
1475  P=FPlus;
1476  Q=FCross;
1477  p3 = pow (f0, 38./9.);
1478  return(p3*(-(2862.*c2b*ci*e4*Q)/(19.) -(1585071.*c2b*ci*e6*Q)/(1444.) -(3905136831.*c2b*ci*e8*Q)/(877952.)-(1431.*e4*P*s2b)/(19.) -(1431.*c2i*e4*P*s2b)/(19.) -(1585071.*e6*P*s2b)/(2888.) -(1585071.*c2i*e6*P*s2b)/(2888.) -(3905136831.*e8*P*s2b)/(1.755904e6) -(3905136831.*c2i*e8*P*s2b)/(1.755904e6) ) );
1479 }
1480 
1481 static REAL8 UNUSED
1482 q17(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1483 {
1484  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p4,c2i;
1485  e2=e0*e0;
1486  e3=e2*e0;
1487  e4=e3*e0;
1488  e5=e4*e0;
1489  e6=e5*e0;
1490  e7=e6*e0;
1491  e8=e7*e0;
1492  c2b=cos(2.*bet);
1493  s2b=sin(2.*bet);
1494  ci=cos(inc);
1495  c2i=ci*ci;
1496  P=FPlus;
1497  Q=FCross;
1498  p4 = pow (f0, 76./9.);
1499  return(p4*( (-49470967683233.*c2b*ci*e8*Q)/(4.97798784e9)-(49470967683233.*e8*P*s2b)/(9.95597568e9) -(49470967683233.*c2i*e8*P*s2b)/(9.95597568e9)) );
1500 }
1501 
1502  /* case 5 for zeta_im*/
1503 
1504 static REAL8 UNUSED
1505 q18(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1506 {
1507  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p5,c2i;
1508  e2=e0*e0;
1509  e3=e2*e0;
1510  e4=e3*e0;
1511  e5=e4*e0;
1512  e6=e5*e0;
1513  e7=e6*e0;
1514  c2b=cos(2.*bet);
1515  s2b=sin(2.*bet);
1516  ci=cos(inc);
1517  c2i=ci*ci;
1518  P=FPlus;
1519  Q=FCross;
1520  p5 = pow (f0, 133./18.);
1521  return(p5*( (3054326535625.*c2b*ci*e7*Q)/(1.117863936e9)+(3054326535625.*e7*P*s2b)/(2.235727872e9) +(3054326535625.*c2i*e7*P*s2b)/(2.235727872e9)) );
1522 }
1523 
1524 static REAL8 UNUSED
1525 q19(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1526 {
1527  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p6,c2i;
1528  e2=e0*e0;
1529  e3=e2*e0;
1530  e4=e3*e0;
1531  e5=e4*e0;
1532  e6=e5*e0;
1533  e7=e6*e0;
1534  c2b=cos(2.*bet);
1535  s2b=sin(2.*bet);
1536  ci=cos(inc);
1537  c2i=ci*ci;
1538  P=FPlus;
1539  Q=FCross;
1540  p6 = pow (f0, 95./18.);
1541  return(p6*( -(13023125.*c2b*ci*e5*Q)/(43776.) -(216379221875.*c2b*ci*e7*Q)/(7.9847424e7)-(13023125.*c2i*e5*P*s2b)/(87552.) -(216379221875.*e7*P*s2b)/(1.59694848e8) -(216379221875.*c2i*e7*P*s2b)/(1.59694848e8)-(13023125.*e5*P*s2b)/(87552.)) );
1542 }
1543 
1544 static REAL8 UNUSED
1545 q20(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1546 {
1547  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p7,c2i;
1548  e2=e0*e0;
1549  e3=e2*e0;
1550  e4=e3*e0;
1551  e5=e4*e0;
1552  e6=e5*e0;
1553  e7=e6*e0;
1554  c2b=cos(2.*bet);
1555  s2b=sin(2.*bet);
1556  ci=cos(inc);
1557  c2i=ci*ci;
1558  P=FPlus;
1559  Q=FCross;
1560  p7 = pow (f0, 19./6.);
1561  return(p7*( (625.*c2b*ci*e3*Q)/(24.) +(2076875.*c2b*ci*e5*Q)/(14592.) +(7933101875.*c2b*ci*e7*Q)/(1.7743872e7)+(625.*e3*P*s2b)/(48.) +(625.*c2i*e3*P*s2b)/(48.) +(2076875.*e5*P*s2b)/(29184.) +(2076875.*c2i*e5*P*s2b)/(29184.) +(7933101875.*e7*P*s2b)/(3.5487744e7) +(7933101875.*c2i*e7*P*s2b)/(3.5487744e7)) );
1562 }
1563 
1564 
1565  /*case 6 for zeta_im */
1566 
1567 static REAL8 UNUSED
1568 q21(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1569 {
1570  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p2,c2i;
1571  e2=e0*e0;
1572  e3=e2*e0;
1573  e4=e3*e0;
1574  e5=e4*e0;
1575  e6=e5*e0;
1576  e7=e6*e0;
1577  e8=e7*e0;
1578  c2b=cos(2.*bet);
1579  s2b=sin(2.*bet);
1580  ci=cos(inc);
1581  c2i=ci*ci;
1582  P=FPlus;
1583  Q=FCross;
1584  p2 = pow (f0, 19./3.);
1585  return(p2*( -(1656963.*c2b*ci*e6*Q)/(3040.) -(5506088049.*c2b*ci*e8*Q)/(924160.)-(1656963.*e6*P*s2b)/(6080.) -(1656963.*c2i*e6*P*s2b)/(6080.) -(5506088049.*e8*P*s2b)/(1.84832e6) -(5506088049.*c2i*e8*P*s2b)/(1.84832e6)) );
1586 }
1587 
1588 static REAL8 UNUSED
1589 q22(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1590 {
1591  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p3,c2i;
1592  e2=e0*e0;
1593  e3=e2*e0;
1594  e4=e3*e0;
1595  e5=e4*e0;
1596  e6=e5*e0;
1597  e7=e6*e0;
1598  e8=e7*e0;
1599  c2b=cos(2.*bet);
1600  s2b=sin(2.*bet);
1601  ci=cos(inc);
1602  c2i=ci*ci;
1603  P=FPlus;
1604  Q=FCross;
1605  p3 = pow (f0, 38./9.);
1606  return(p3*( (81.*c2b*ci*e4*Q)/(2.) +(89721.*c2b*ci*e6*Q)/(304.) +(221045481.*c2b*ci*e8*Q)/(184832.)+(81.*e4*P*s2b)/(4.) +(81.*c2i*e4*P*s2b)/(4.) +(89721.*e6*P*s2b)/(608.) +(89721.*c2i*e6*P*s2b)/(608.) +(221045481.*e8*P*s2b)/(369664.) +(221045481.*c2i*e8*P*s2b)/(369664.)) );
1607 }
1608 
1609 static REAL8 UNUSED
1610 q23(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1611 {
1612  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p4,c2i;
1613  e2=e0*e0;
1614  e3=e2*e0;
1615  e4=e3*e0;
1616  e5=e4*e0;
1617  e6=e5*e0;
1618  e7=e6*e0;
1619  e8=e7*e0;
1620  c2b=cos(2.*bet);
1621  s2b=sin(2.*bet);
1622  ci=cos(inc);
1623  c2i=ci*ci;
1624  P=FPlus;
1625  Q=FCross;
1626  p4 = pow (f0, 76./9.);
1627  return(p4*((71738041617.*c2b*ci*e8*Q)/(1.293824e7)+(71738041617.*e8*P*s2b)/(2.587648e7) +(71738041617.*c2i*e8*P*s2b)/(2.587648e7) ) );
1628 }
1629 
1630  /*case 7 for zeta_im*/
1631 
1632 
1633 static REAL8 UNUSED
1634 q24(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1635 {
1636  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p5,c2i;
1637  e2=e0*e0;
1638  e3=e2*e0;
1639  e4=e3*e0;
1640  e5=e4*e0;
1641  e6=e5*e0;
1642  e7=e6*e0;
1643  c2b=cos(2.*bet);
1644  s2b=sin(2.*bet);
1645  ci=cos(inc);
1646  c2i=ci*ci;
1647  P=FPlus;
1648  Q=FCross;
1649  p5 = pow (f0, 133./18.);
1650  return(p5*((-1109077123.*c2b*ci*e7*Q)/(1.16736e6)-(1109077123.*e7*P*s2b)/(2.33472e6) -(1109077123.*c2i*e7*P*s2b)/(2.33472e6) ) );
1651 }
1652 
1653 static REAL8 UNUSED
1654 q25(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1655 {
1656  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p6,c2i;
1657  e2=e0*e0;
1658  e3=e2*e0;
1659  e4=e3*e0;
1660  e5=e4*e0;
1661  e6=e5*e0;
1662  e7=e6*e0;
1663  c2b=cos(2.*bet);
1664  s2b=sin(2.*bet);
1665  ci=cos(inc);
1666  c2i=ci*ci;
1667  P=FPlus;
1668  Q=FCross;
1669  p6 = pow (f0, 95./18.);
1670  return(p6*((117649.*c2b*ci*e5*Q)/(1920.) +(390947627.*c2b*ci*e7*Q)/(700416.)+(117649.*e5*P*s2b)/(3840.) +(117649.*c2i*e5*P*s2b)/(3840.) +(390947627.*e7*P*s2b)/(1.400832e6) +(390947627.*c2i*e7*P*s2b)/(1.400832e6) ) );
1671 }
1672 
1673 
1674  /* case 8 for zeta_im*/
1675 
1676 static REAL8 UNUSED
1677 q26(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1678 {
1679  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p2,c2i;
1680  e2=e0*e0;
1681  e3=e2*e0;
1682  e4=e3*e0;
1683  e5=e4*e0;
1684  e6=e5*e0;
1685  e7=e6*e0;
1686  e8=e7*e0;
1687  c2b=cos(2.*bet);
1688  s2b=sin(2.*bet);
1689  ci=cos(inc);
1690  c2i=ci*ci;
1691  P=FPlus;
1692  Q=FCross;
1693  p2 = pow (f0, 19./3.);
1694  return(p2*((4096.*c2b*ci*e6*Q)/(45.) +(850688.*c2b*ci*e8*Q)/(855.)+(2048.*e6*P*s2b)/(45.) +(2048.*c2i*e6*P*s2b)/(45.) +(425344.*e8*P*s2b)/(855.) +(425344.*c2i*e8*P*s2b)/(855.)) );
1695 }
1696 
1697 static REAL8 UNUSED
1698 q27(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1699 {
1700  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p4,c2i;
1701  e2=e0*e0;
1702  e3=e2*e0;
1703  e4=e3*e0;
1704  e5=e4*e0;
1705  e6=e5*e0;
1706  e7=e6*e0;
1707  e8=e7*e0;
1708  c2b=cos(2.*bet);
1709  s2b=sin(2.*bet);
1710  ci=cos(inc);
1711  c2i=ci*ci;
1712  P=FPlus;
1713  Q=FCross;
1714  p4 = pow (f0, 76./9.);
1715  return(p4*( (-28696576.*c2b*ci*e8*Q)/(17955.)-(14348288.*e8*P*s2b)/(17955.) -(14348288.*c2i*e8*P*s2b)/(17955.)) );
1716 }
1717 
1718 
1719  /*case 9 of zeta_im*/
1720 
1721 static REAL8 UNUSED
1722 q28(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1723 {
1724  REAL8 e2,e3,e4,e5,e6,e7,c2b,s2b,ci,P,Q,p5,c2i;
1725  e2=e0*e0;
1726  e3=e2*e0;
1727  e4=e3*e0;
1728  e5=e4*e0;
1729  e6=e5*e0;
1730  e7=e6*e0;
1731  c2b=cos(2.*bet);
1732  s2b=sin(2.*bet);
1733  ci=cos(inc);
1734  c2i=ci*ci;
1735  P=FPlus;
1736  Q=FCross;
1737  p5 = pow (f0, 133./18.);
1738  return(p5*((4782969.*c2b*ci*e7*Q)/(35840.) +(4782969.*e7*P*s2b)/(71680.) +(4782969.*c2i*e7*P*s2b)/(71680.)));
1739 }
1740 
1741 
1742  /*case 10 for zeta_im*/
1743 
1744 static REAL8 UNUSED
1745 q29(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
1746 {
1747  REAL8 e2,e3,e4,e5,e6,e7,e8,c2b,s2b,ci,P,Q,p4;
1748  e2=e0*e0;
1749  e3=e2*e0;
1750  e4=e3*e0;
1751  e5=e4*e0;
1752  e6=e5*e0;
1753  e7=e6*e0;
1754  e8=e7*e0;
1755  c2b=cos(2.*bet);
1756  s2b=sin(2.*bet);
1757  ci=cos(inc);
1758  P=FPlus;
1759  Q=FCross;
1760  p4 = pow (f0, 76./9.);
1761  return(p4*(390625.*e8*(2.*c2b*ci*Q + P*s2b + ci*ci*P*s2b))/(4032.));
1762 }
z8
static REAL8 UNUSED z8(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:654
C25
static REAL8 UNUSED C25(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:318
C32
static REAL8 UNUSED C32(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:396
q24
static REAL8 UNUSED q24(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1634
q16
static REAL8 UNUSED q16(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1461
C22
static REAL8 UNUSED C22(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:282
C24
static REAL8 UNUSED C24(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:306
z19
static REAL8 UNUSED z19(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:904
q29
static REAL8 UNUSED q29(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1745
z10
static REAL8 UNUSED z10(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:702
C11
static REAL8 UNUSED C11(REAL8 UNUSED eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:161
z20
static REAL8 UNUSED z20(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:926
q8
static REAL8 UNUSED q8(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1292
z23
static REAL8 UNUSED z23(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:992
q23
static REAL8 UNUSED q23(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1610
C13
static REAL8 UNUSED C13(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:181
C28
static REAL8 UNUSED C28(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:347
C27
static REAL8 UNUSED C27(REAL8 UNUSED eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:341
q5
static REAL8 UNUSED q5(REAL8 UNUSED f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1237
q3
static REAL8 UNUSED q3(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1195
z16
static REAL8 UNUSED z16(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:834
C39
static REAL8 UNUSED C39(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:478
z18
static REAL8 UNUSED z18(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:882
C15
static REAL8 UNUSED C15(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:205
C19
static REAL8 UNUSED C19(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:249
C36
static REAL8 UNUSED C36(REAL8 eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:446
C31
static REAL8 UNUSED C31(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:384
C37
static REAL8 UNUSED C37(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:452
C8
static REAL8 UNUSED C8(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:120
q21
static REAL8 UNUSED q21(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1568
C35
static REAL8 UNUSED C35(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:432
C23
static REAL8 UNUSED C23(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:294
C17
static REAL8 UNUSED C17(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:225
q6
static REAL8 UNUSED q6(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1250
C7
static REAL8 UNUSED C7(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:107
C2
static REAL8 UNUSED C2(REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:50
z4
static REAL8 UNUSED z4(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:571
q12
static REAL8 UNUSED q12(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1376
z15
static REAL8 UNUSED z15(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:811
q20
static REAL8 UNUSED q20(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1545
q13
static REAL8 UNUSED q13(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1396
C1
static REAL8 UNUSED C1(REAL8 Mtotal)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:42
C4
static REAL8 UNUSED C4(REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:74
C10
static REAL8 UNUSED C10(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:146
q14
static REAL8 UNUSED q14(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1419
q22
static REAL8 UNUSED q22(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1589
q2
static REAL8 UNUSED q2(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1175
q10
static REAL8 UNUSED q10(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1336
z24
static REAL8 UNUSED z24(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1017
q9
static REAL8 UNUSED q9(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1313
z21
static REAL8 UNUSED z21(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:948
C38
static REAL8 UNUSED C38(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:466
z6
static REAL8 UNUSED z6(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:608
q1
static REAL8 UNUSED q1(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1155
C9
static REAL8 UNUSED C9(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:133
C16
static REAL8 UNUSED C16(REAL8 eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:219
C5
static REAL8 UNUSED C5(REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:87
C20
static REAL8 UNUSED C20(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:262
z26
static REAL8 UNUSED z26(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1061
C26
static REAL8 UNUSED C26(REAL8 eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:332
q11
static REAL8 UNUSED q11(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1356
q4
static REAL8 UNUSED q4(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1215
z12
static REAL8 UNUSED z12(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:746
q18
static REAL8 UNUSED q18(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1505
z29
static REAL8 UNUSED z29(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1130
q17
static REAL8 UNUSED q17(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1482
z14
static REAL8 UNUSED z14(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:790
z17
static REAL8 UNUSED z17(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:857
q15
static REAL8 UNUSED q15(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1440
C3
static REAL8 UNUSED C3(REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:62
q25
static REAL8 UNUSED q25(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1654
z7
static REAL8 UNUSED z7(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:631
C34
static REAL8 UNUSED C34(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:420
C14
static REAL8 UNUSED C14(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:193
z13
static REAL8 UNUSED z13(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:768
z25
static REAL8 UNUSED z25(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1039
z11
static REAL8 UNUSED z11(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:724
C21
static REAL8 UNUSED C21(REAL8 eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:276
q19
static REAL8 UNUSED q19(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1525
C0
static REAL8 UNUSED C0(REAL8 eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:36
q7
static REAL8 UNUSED q7(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1271
z9
static REAL8 UNUSED z9(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:677
z1
static REAL8 UNUSED z1(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:505
q27
static REAL8 UNUSED q27(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1698
C30
static REAL8 UNUSED C30(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:372
z5
static REAL8 UNUSED z5(REAL8 UNUSED f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:595
z22
static REAL8 UNUSED z22(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:971
C40
static REAL8 UNUSED C40(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:490
z3
static REAL8 UNUSED z3(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:549
C12
static REAL8 UNUSED C12(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:168
z27
static REAL8 UNUSED z27(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1083
z28
static REAL8 UNUSED z28(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1108
C18
static REAL8 UNUSED C18(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:237
q28
static REAL8 UNUSED q28(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1722
C29
static REAL8 UNUSED C29(REAL8 UNUSED eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:360
q26
static REAL8 UNUSED q26(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:1677
z2
static REAL8 UNUSED z2(REAL8 e0, REAL8 f0, REAL8 inc, REAL8 bet, REAL8 FPlus, REAL8 FCross)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:527
C6
static REAL8 UNUSED C6(REAL8 eta)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:101
C33
static REAL8 UNUSED C33(REAL8 eta, REAL8 e0, REAL8 f0)
Definition: LALSimInspiralOptimizedCoefficientsEccentricityFD.c:408
Q
const double Q
Definition: exact_derivatives-Hreal.c:33
LAL_PI
#define LAL_PI
REAL8
double REAL8