Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALSimulation 6.2.0.1-6c6b863
•All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
LALSimInspiralEccentricityFD.c
Go to the documentation of this file.
1/* Copyright (C) 2014 E. A. Huerta, Prayush Kumar
2 *
3 * This program is free software; you can redistribute it and/or modify
4 * it under the terms of the GNU General Public License as published by
5 * the Free Software Foundation; either version 2 of the License, or
6 * (at your option) any later version.
7 *
8 * This program is distributed in the hope that it will be useful,
9 * but WITHOUT ANY WARRANTY; without even the implied warranty of
10 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
11 * GNU General Public License for more details.
12 *
13 * You should have received a copy of the GNU General Public License
14 * along with with program; see the file COPYING. If not, write to the
15 * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
16 * MA 02110-1301 USA
17 *
18 * This code describes the gravitational radiation emitted by compact binaries with low to moderate values of eccentricity [0 < e < 0.5]. The waveform phase includes post-Newtonian corrections up to 3.5 order and eccentricity corrections up to order e^8 at each post-Newtonian order. The waveform amplitude is modeled within the restricted post-Newtonian approach.
19 */
20
21
22#include <stdio.h>
23#include <stdlib.h>
24#include <math.h>
25#include <complex.h>
26#include <lal/Date.h>
27#include <lal/FrequencySeries.h>
28#include <lal/LALConstants.h>
29#include <lal/LALDatatypes.h>
30#include <lal/LALSimInspiral.h>
31#include <lal/Units.h>
32#include <lal/XLALError.h>
33
34
35#include "gsl/gsl_matrix.h"
36#include "gsl/gsl_vector.h"
37#include "gsl/gsl_linalg.h"
38#include "gsl/gsl_eigen.h"
39#include "gsl/gsl_permutation.h"
40#include <gsl/gsl_errno.h>
41#include <gsl/gsl_math.h>
42#include <gsl/gsl_roots.h>
43#include <gsl/gsl_sf_bessel.h>
44#include <gsl/gsl_sf_gamma.h>
45#include <gsl/gsl_complex.h>
46#include <gsl/gsl_complex_math.h>
47#include <gsl/gsl_errno.h>
48#include <gsl/gsl_spline.h>
49#include <gsl/gsl_errno.h>
50#include <gsl/gsl_multiroots.h>
51#include <gsl/gsl_rng.h>
52
53#include "LALSimInspiralOptimizedCoefficientsEccentricityFD.c"
54
55#define gamma (0.577215664901532860606512090)
56
57#ifdef __GNUC__
58#define UNUSED __attribute__ ((unused))
59#else
60#define UNUSED
61#endif
62
63
64typedef struct
65tagexpnCoeffsEPC {
66
67 /*Input params*/
68 REAL8 f0, e0;
69
70 /* Coefficients to reconstruct the phase*/
71 REAL8 a0, a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11, a12,a13,a14, a15, a16, a17, a18, a19, a20, a21, a22, a23, a24, a25, a26, a27, a28, a29, a30, a31, a32, a33, a34, a35, a36, a37, a38, a39, a40;
72
73 /*Coefficients to reconstruct z_real*/
74 /*plus component*/
75 REAL8 b1p, b2p, b3p, b4p, b5p, b6p, b7p, b8p, b9p, b10p, b11p, b12p, b13p, b14p, b15p, b16p, b17p, b18p, b19p, b20p, b21p, b22p, b23p, b24p, b25p, b26p, b27p, b28p, b29p;
76 /*cross component*/
77 REAL8 b1c, b2c, b3c, b4c, b5c, b6c, b7c, b8c, b9c, b10c, b11c, b12c, b13c, b14c, b15c, b16c, b17c, b18c, b19c, b20c, b21c, b22c, b23c, b24c, b25c, b26c, b27c, b28c, b29c;
78
79 /*Coefficients to reconstruct z_imaginary*/
80 /*plus component*/
81 REAL8 d1p, d2p, d3p, d4p, d5p, d6p, d7p, d8p, d9p, d10p, d11p, d12p, d13p, d14p, d15p, d16p, d17p, d18p, d19p, d20p, d21p, d22p, d23p, d24p, d25p, d26p, d27p, d28p, d29p;
82 /*cross component*/
83 REAL8 d1c, d2c, d3c, d4c, d5c, d6c, d7c, d8c, d9c, d10c, d11c, d12c, d13c, d14c, d15c, d16c, d17c, d18c, d19c, d20c, d21c, d22c, d23c, d24c, d25c, d26c, d27c, d28c, d29c;
84 /* symmetric mass ratio, total mass, component masses*/
85 //REAL8 nu,m,m1,m2,mu,Mtotal;
86
87}expnCoeffsEPC;
88
89static REAL8 PhaseAccTaylorF2(int k, REAL8 f, expnCoeffsEPC *ak){
90
91 REAL8 k4, k5, k19, k19ov3, k19ov9, k19ov18, f1, f2, f3, f4, f8, fact, lo1, lo2, lo3, lo4,lfact;
92 k4=f*f*f*f;
93 k5=f*f*f*f*f;
94 k19=k5*k5*k5*k4;
95 k19ov3=cbrt(k19);
96 k19ov9=cbrt(k19ov3);
97 k19ov18=sqrt(k19ov9);
98 f8=1./(k19ov18);
99 f1=f8*f8;
100 f2=f1*f1*f1;
101 f3=f1*f1;
102 f4=f3*f3;
103 fact=ak->a1*f;
104 lo1=fact;
105 lo4=cbrt(lo1);
106 lo3=lo4*lo4;
107 lo2=lo3*lo3*lo4;
108 lfact=6.*sqrt(6.);
109
110 switch( k )
111 {
112 case 1: return( ak->a0*(1./lo2)*( 1.+ ak->a2*f1 + ak->a3*f2 + ak->a4*f3 + ak->a5*f4) );
113 break;
114
115 case 2: return( ak->a0*(1./lo1)*(ak->a6 + ak->a7*f1 + ak->a8*f2 + ak->a9*f3 + ak->a10*f4) );
116 break;
117
118 case 3: return( ak->a0*(1./lo3)*(ak->a11 + ak->a12*f1 + ak->a13*f2 + ak->a14*f3 + ak->a15*f4 ) );
119 break;
120
121 case 4: return(ak->a0*(1./lo4)*(ak->a16 + ak->a17*f1 + ak->a18*f2 + ak->a19*f3 + ak->a20*f4) );
122 break;
123
124 case 5: return( ak->a0*(ak->a21 + ak->a21*log(lfact*lo1) +ak->a22*f1+ak->a23*f2 + ak->a24*f3+ak->a25*f4 ) );
125 break;
126
127 case 6: return(ak->a0*(lo4)*(ak->a26 + log(4.*lo4)*(ak->a27 + ak->a28*f1 + ak->a29*f2 + ak->a30*f3 + ak->a31*f4 ) + ak->a32*f1 + ak->a33*f2 + ak->a34*f3 + ak->a35*f4 ) );
128 break;
129
130 case 7: return(ak->a0*(lo3)*( ak->a36 + ak->a37*f1+ ak->a38*f2+ ak->a39*f3+ak->a40*f4) );
131 break;
132
133 default: return 0;
134 break;
135
136 }
137}
138
139
140static REAL8 zeta_generic_re_plus(int k, REAL8 f, expnCoeffsEPC *ak) {
141 REAL8 k4, k5, k19, k19ov3, k19ov9, k19ov18, f1,f2,f3,f4,f5,f6,f7,f8;
142 k4=f*f*f*f;
143 k5=f*f*f*f*f;
144 k19=k5*k5*k5*k4;
145 k19ov3=cbrt(k19);
146 k19ov9=cbrt(k19ov3);
147 k19ov18=sqrt(k19ov9);
148 f8=1./(k19ov18);
149 f1=f8*f8;
150 f2=f1*f1*f1;
151 f3=f1*f1;
152 f4=f3*f3;
153 f5=f2*f8;
154 f6=f3*f8;
155 f7=f1*f8;
156
157 switch( k )
158 {
159 case 1: return(ak->b1p*f5 + ak->b2p*f6 + ak->b3p*f7 + ak->b4p*f8);
160 break;
161
162 case 2: return(ak->b5p + ak->b6p*f1 + ak->b7p*f2 + ak->b8p*f3 + ak->b9p*f4);
163 break;
164
165 case 3: return(ak->b10p*f5 + ak->b11p*f6 + ak->b12p*f7 + ak->b13p*f8);
166 break;
167
168 case 4: return(ak->b14p*f1 + ak->b15p*f2 + ak->b16p*f3 + ak->b17p*f4);
169 break;
170
171 case 5: return(ak->b18p*f5 + ak->b19p*f6 + ak->b20p*f7);
172 break;
173
174 case 6: return(ak->b21p*f2 + ak->b22p*f3 + ak->b23p*f4);
175 break;
176
177 case 7: return(ak->b24p*f5 + ak->b25p*f6);
178 break;
179
180 case 8: return(ak->b26p*f2 + ak->b27p*f4);
181 break;
182
183 case 9: return(ak->b28p*f5);
184 break;
185
186 case 10: return(ak->b29p*f4);
187 break;
188
189 default: return 0;
190 break;
191
192 }
193}
194
195static REAL8 zeta_generic_re_cross(int k, REAL8 f, expnCoeffsEPC *ak) {
196 REAL8 k4, k5, k19, k19ov3, k19ov9, k19ov18, f1,f2,f3,f4,f5,f6,f7,f8;
197 k4=f*f*f*f;
198 k5=f*f*f*f*f;
199 k19=k5*k5*k5*k4;
200 k19ov3=cbrt(k19);
201 k19ov9=cbrt(k19ov3);
202 k19ov18=sqrt(k19ov9);
203 f8=1./(k19ov18);
204 f1=f8*f8;
205 f2=f1*f1*f1;
206 f3=f1*f1;
207 f4=f3*f3;
208 f5=f2*f8;
209 f6=f3*f8;
210 f7=f1*f8;
211
212 switch( k )
213 {
214 case 1: return(ak->b1c*f5 + ak->b2c*f6 + ak->b3c*f7 + ak->b4c*f8);
215 break;
216
217 case 2: return(ak->b5c + ak->b6c*f1 + ak->b7c*f2 + ak->b8c*f3 + ak->b9c*f4);
218 break;
219
220 case 3: return(ak->b10c*f5 + ak->b11c*f6 + ak->b12c*f7 + ak->b13c*f8);
221 break;
222
223 case 4: return(ak->b14c*f1 + ak->b15c*f2 + ak->b16c*f3 + ak->b17c*f4);
224 break;
225
226 case 5: return(ak->b18c*f5 + ak->b19c*f6 + ak->b20c*f7);
227 break;
228
229 case 6: return(ak->b21c*f2 + ak->b22c*f3 + ak->b23c*f4);
230 break;
231
232 case 7: return(ak->b24c*f5 + ak->b25c*f6);
233 break;
234
235 case 8: return(ak->b26c*f2 + ak->b27c*f4);
236 break;
237
238 case 9: return(ak->b28c*f5);
239 break;
240
241 case 10: return(ak->b29c*f4);
242 break;
243
244 default: return 0;
245 break;
246
247 }
248}
249
250
251static REAL8 zeta_generic_im_plus(int k, REAL8 f, expnCoeffsEPC *ak) {
252 REAL8 k4, k5, k19, k19ov3, k19ov9, k19ov18, f1,f2,f3,f4,f5,f6,f7,f8;
253 k4=f*f*f*f;
254 k5=f*f*f*f*f;
255 k19=k5*k5*k5*k4;
256 k19ov3=cbrt(k19);
257 k19ov9=cbrt(k19ov3);
258 k19ov18=sqrt(k19ov9);
259 f8=1./(k19ov18);
260 f1=f8*f8;
261 f2=f1*f1*f1;
262 f3=f1*f1;
263 f4=f3*f3;
264 f5=f2*f8;
265 f6=f3*f8;
266 f7=f1*f8;
267
268 switch( k )
269 {
270 case 1: return(ak->d1p*f5 + ak->d2p*f6 + ak->d3p*f7 + ak->d4p*f8);
271 break;
272
273 case 2: return(ak->d5p + ak->d6p*f1 + ak->d7p*f2 + ak->d8p*f3 + ak->d9p*f4);
274 break;
275
276 case 3: return(ak->d10p*f5 + ak->d11p*f6 + ak->d12p*f7 + ak->d13p*f8);
277 break;
278
279 case 4: return(ak->d14p*f1 + ak->d15p*f2 + ak->d16p*f3 + ak->d17p*f4);
280 break;
281
282 case 5: return(ak->d18p*f5 + ak->d19p*f6 + ak->d20p*f7);
283 break;
284
285 case 6: return(ak->d21p*f2 + ak->d22p*f3 + ak->d23p*f4);
286 break;
287
288 case 7: return(ak->d24p*f5 + ak->d25p*f6);
289 break;
290
291 case 8: return(ak->d26p*f2 + ak->d27p*f4);
292 break;
293
294 case 9: return(ak->d28p*f5);
295 break;
296
297 case 10: return(ak->d29p*f4);
298 break;
299
300 default: return 0;
301 break;
302
303 }
304}
305
306static REAL8 zeta_generic_im_cross(int k, REAL8 f, expnCoeffsEPC *ak) {
307 REAL8 k4, k5, k19, k19ov3, k19ov9, k19ov18, f1,f2,f3,f4,f5,f6,f7,f8;
308 k4=f*f*f*f;
309 k5=f*f*f*f*f;
310 k19=k5*k5*k5*k4;
311 k19ov3=cbrt(k19);
312 k19ov9=cbrt(k19ov3);
313 k19ov18=sqrt(k19ov9);
314 f8=1./(k19ov18);
315 f1=f8*f8;
316 f2=f1*f1*f1;
317 f3=f1*f1;
318 f4=f3*f3;
319 f5=f2*f8;
320 f6=f3*f8;
321 f7=f1*f8;
322
323 switch( k )
324 {
325 case 1: return(ak->d1c*f5 + ak->d2c*f6 + ak->d3c*f7 + ak->d4c*f8);
326 break;
327
328 case 2: return(ak->d5c + ak->d6c*f1 + ak->d7c*f2 + ak->d8c*f3 + ak->d9c*f4);
329 break;
330
331 case 3: return(ak->d10c*f5 + ak->d11c*f6 + ak->d12c*f7 + ak->d13c*f8);
332 break;
333
334 case 4: return(ak->d14c*f1 + ak->d15c*f2 + ak->d16c*f3 + ak->d17c*f4);
335 break;
336
337 case 5: return(ak->d18c*f5 + ak->d19c*f6 + ak->d20c*f7);
338 break;
339
340 case 6: return(ak->d21c*f2 + ak->d22c*f3 + ak->d23c*f4);
341 break;
342
343 case 7: return(ak->d24c*f5 + ak->d25c*f6);
344 break;
345
346 case 8: return(ak->d26c*f2 + ak->d27c*f4);
347 break;
348
349 case 9: return(ak->d28c*f5);
350 break;
351
352 case 10: return(ak->d29c*f4);
353 break;
354
355 default: return 0;
356 break;
357
358 }
359}
360
361
362static int
363EPCSetup(
364 expnCoeffsEPC *ak, /**< coefficients for EPC evolution [modified] */
365 const REAL8 m1, /**< Mass of companion 1 (kg) */
366 const REAL8 m2, /**< Mass of companion 2 (kg) */
367 const REAL8 fStart, /**< Start GW frequency (Hz) */
368 const REAL8 i, /**< Polar inclination of source (rad) */
369 const REAL8 inclination_azimuth, /**< Azimuthal component of inclination angles */
370 const REAL8 e_min /**< Initial eccentricity at f_min: range [0, 0.5] */
371 )
372{
373
374 const REAL8 m = m1 + m2;
375 const REAL8 Mtotal = m * LAL_MTSUN_SI; /* total mass in seconds */
376 const REAL8 eta = m1 * m2 / (m * m);
377 const REAL8 e0=e_min;
378 const REAL8 inc=i;
379 const REAL8 bet=inclination_azimuth;
380 const REAL8 f0=fStart;
381
382 /* Coefficients to reconstruct the phase*/
383
384 ak->a0 = C0(eta); ak->b1p = z1(e0, f0, inc, bet, 1., 0.);
385 ak->a1 = C1(Mtotal); ak->b2p = z2(e0, f0, inc, bet, 1., 0.);
386 ak->a2 = C2(e0, f0); ak->b3p = z3(e0, f0, inc, bet, 1., 0.);
387 ak->a3 = C3(e0, f0); ak->b4p = z4(e0, f0, inc, bet, 1., 0.);
388 ak->a4 = C4(e0, f0); ak->b5p = z5(f0, inc, bet, 1., 0.);
389 ak->a5 = C5(e0, f0); ak->b6p = z6(e0, f0, inc, bet, 1., 0.);
390 ak->a6 = C6(eta); ak->b7p = z7(e0, f0, inc, bet, 1., 0.);
391 ak->a7 = C7(eta, e0, f0); ak->b8p = z8(e0, f0, inc, bet, 1., 0.);
392 ak->a8 = C8(eta, e0, f0); ak->b9p = z9(e0, f0, inc, bet, 1., 0.);
393 ak->a9 = C9(eta, e0, f0); ak->b10p = z10(e0, f0, inc, bet, 1., 0.);
394 ak->a10 = C10(eta, e0, f0); ak->b11p = z11(e0, f0, inc, bet, 1., 0.);
395 ak->a11 = C11(eta); ak->b12p = z12(e0, f0, inc, bet, 1., 0.);
396 ak->a12 = C12(eta, e0, f0); ak->b13p = z13(e0, f0, inc, bet, 1., 0.);
397 ak->a13 = C13(eta, e0, f0); ak->b14p = z14(e0, f0, inc, bet, 1., 0.);
398 ak->a14 = C14(eta, e0, f0); ak->b15p = z15(e0, f0, inc, bet, 1., 0.);
399 ak->a15 = C15(eta, e0, f0); ak->b16p = z16(e0, f0, inc, bet, 1., 0.);
400 ak->a16 = C16(eta); ak->b17p = z17(e0, f0, inc, bet, 1., 0.);
401 ak->a17 = C17(eta, e0, f0); ak->b18p = z18(e0, f0, inc, bet, 1., 0.);
402 ak->a18 = C18(eta, e0, f0); ak->b19p = z19(e0, f0, inc, bet, 1., 0.);
403 ak->a19 = C19(eta, e0, f0); ak->b20p = z20(e0, f0, inc, bet, 1., 0.);
404 ak->a20 = C20(eta, e0, f0); ak->b21p = z21(e0, f0, inc, bet, 1., 0.);
405 ak->a21 = C21(eta); ak->b22p = z22(e0, f0, inc, bet, 1., 0.);
406 ak->a22 = C22(eta, e0, f0); ak->b23p = z23(e0, f0, inc, bet, 1., 0.);
407 ak->a23 = C23(eta, e0, f0); ak->b24p = z24(e0, f0, inc, bet, 1., 0.);
408 ak->a24 = C24(eta, e0, f0); ak->b25p = z25(e0, f0, inc, bet, 1., 0.);
409 ak->a25 = C25(eta, e0, f0); ak->b26p = z26(e0, f0, inc, bet, 1., 0.);
410 ak->a26 = C26(eta); ak->b27p = z27(e0, f0, inc, bet, 1., 0.);
411 ak->a27 = C27(eta); ak->b28p = z28(e0, f0, inc, bet, 1., 0.);
412 ak->a28 = C28(eta, e0, f0); ak->b29p = z29(e0, f0, inc, bet, 1., 0.);
413 ak->a29 = C29(eta, e0, f0); ak->d1p = q1(e0, f0, inc, bet, 1., 0.);
414 ak->a30 = C30(eta, e0, f0); ak->d2p = q2(e0, f0, inc, bet,1., 0.);
415 ak->a31 = C31(eta, e0, f0); ak->d3p = q3(e0, f0, inc, bet, 1., 0.);
416 ak->a32 = C32(eta, e0, f0); ak->d4p = q4(e0, f0, inc, bet, 1., 0.);
417 ak->a33 = C33(eta, e0, f0); ak->d5p = q5(f0, inc, bet, 1., 0.);
418 ak->a34 = C34(eta, e0, f0); ak->d6p = q6(e0, f0, inc, bet, 1., 0.);
419 ak->a35 = C35(eta, e0, f0); ak->d7p = q7(e0, f0, inc, bet, 1., 0.);
420 ak->a36 = C36(eta); ak->d8p = q8(e0, f0, inc, bet, 1., 0.);
421 ak->a37 = C37(eta, e0, f0); ak->d9p = q9(e0, f0, inc, bet, 1., 0.);
422 ak->a38 = C38(eta, e0, f0); ak->d10p = q10(e0, f0, inc, bet, 1., 0.);
423 ak->a39 = C39(eta, e0, f0); ak->d11p = q11(e0, f0, inc, bet, 1., 0.);
424 ak->a40 = C40(eta, e0, f0); ak->d12p = q12(e0, f0, inc, bet, 1., 0.);
425 ak->d13p = q13(e0, f0, inc, bet, 1., 0.); ak->d14p = q14(e0, f0, inc, bet, 1., 0.);
426 ak->d15p = q15(e0, f0, inc, bet, 1., 0.); ak->d16p = q16(e0, f0, inc, bet, 1., 0.);
427 ak->d17p = q17(e0, f0, inc, bet, 1., 0.); ak->d18p = q18(e0, f0, inc, bet, 1., 0.);
428 ak->d19p = q19(e0, f0, inc, bet, 1., 0.); ak->d20p = q20(e0, f0, inc, bet, 1., 0.);
429 ak->d21p = q21(e0, f0, inc, bet, 1., 0.); ak->d22p = q22(e0, f0, inc, bet, 1., 0.);
430 ak->d23p = q23(e0, f0, inc, bet, 1., 0.); ak->d24p = q24(e0, f0, inc, bet, 1., 0.);
431 ak->d25p = q25(e0, f0, inc, bet, 1., 0.); ak->d26p = q26(e0, f0, inc, bet, 1., 0.);
432 ak->d27p = q27(e0, f0, inc, bet, 1., 0.); ak->d28p = q28(e0, f0, inc, bet, 1., 0.);
433 ak->d29p = q29(e0, f0, inc, bet, 1., 0.);
434
435 ak->b1c = z1(e0, f0, inc, bet, 0., 1.); ak->b2c = z2(e0, f0, inc, bet, 0., 1.);
436 ak->b3c = z3(e0, f0, inc, bet, 0., 1.); ak->b4c = z4(e0, f0, inc, bet, 0., 1.);
437 ak->b5c = z5(f0, inc, bet, 0., 1.); ak->b6c = z6(e0, f0, inc, bet, 0., 1.);
438 ak->b7c = z7(e0, f0, inc, bet, 0., 1.); ak->b8c = z8(e0, f0, inc, bet, 0., 1.);
439 ak->b9c = z9(e0, f0, inc, bet, 0., 1.); ak->b10c = z10(e0, f0, inc, bet, 0., 1.);
440 ak->b11c = z11(e0, f0, inc, bet, 0., 1.); ak->b12c = z12(e0, f0, inc, bet, 0., 1.);
441 ak->b13c = z13(e0, f0, inc, bet, 0., 1.); ak->b14c = z14(e0, f0, inc, bet, 0., 1.);
442 ak->b15c = z15(e0, f0, inc, bet, 0., 1.); ak->b16c = z16(e0, f0, inc, bet, 0., 1.);
443 ak->b17c = z17(e0, f0, inc, bet, 0., 1.); ak->b18c = z18(e0, f0, inc, bet, 0., 1.);
444 ak->b19c = z19(e0, f0, inc, bet, 0., 1.); ak->b20c = z20(e0, f0, inc, bet, 0., 1.);
445 ak->b21c = z21(e0, f0, inc, bet, 0., 1.); ak->b22c = z22(e0, f0, inc, bet, 0., 1.);
446 ak->b23c = z23(e0, f0, inc, bet, 0., 1.); ak->b24c = z24(e0, f0, inc, bet, 0., 1.);
447 ak->b25c = z25(e0, f0, inc, bet, 0., 1.); ak->b26c = z26(e0, f0, inc, bet, 0., 1.);
448 ak->b27c = z27(e0, f0, inc, bet, 0., 1.); ak->b28c = z28(e0, f0, inc, bet, 0., 1.);
449 ak->b29c = z29(e0, f0, inc, bet, 0., 1.);
450
451 ak->d1c = q1(e0, f0, inc, bet, 0., 1.); ak->d2c = q2(e0, f0, inc, bet, 0., 1.);
452 ak->d3c = q3(e0, f0, inc, bet, 0., 1.); ak->d4c = q4(e0, f0, inc, bet, 0., 1.);
453 ak->d5c = q5(f0, inc, bet, 0., 1.); ak->d6c = q6(e0, f0, inc, bet, 0., 1.);
454 ak->d7c = q7(e0, f0, inc, bet, 0., 1.); ak->d8c = q8(e0, f0, inc, bet, 0., 1.);
455 ak->d9c = q9(e0, f0, inc, bet, 0., 1.); ak->d10c = q10(e0, f0, inc, bet, 0., 1.);
456 ak->d11c =q11(e0, f0, inc, bet, 0., 1.); ak->d12c = q12(e0, f0, inc, bet, 0., 1.);
457 ak->d13c = q13(e0, f0, inc, bet, 0., 1.); ak->d14c = q14(e0, f0, inc, bet, 0., 1.);
458 ak->d15c = q15(e0, f0, inc, bet, 0., 1.); ak->d16c = q16(e0, f0, inc, bet, 0., 1.);
459 ak->d17c = q17(e0, f0, inc, bet, 0., 1.); ak->d18c = q18(e0, f0, inc, bet, 0., 1.);
460 ak->d19c = q19(e0, f0, inc, bet, 0., 1.); ak->d20c = q20(e0, f0, inc, bet, 0., 1.);
461 ak->d21c = q21(e0, f0, inc, bet, 0., 1.); ak->d22c = q22(e0, f0, inc, bet, 0., 1.);
462 ak->d23c = q23(e0, f0, inc, bet, 0., 1.); ak->d24c = q24(e0, f0, inc, bet, 0., 1.);
463 ak->d25c = q25(e0, f0, inc, bet, 0., 1.); ak->d26c = q26(e0, f0, inc, bet, 0., 1.);
464 ak->d27c = q27(e0, f0, inc, bet, 0., 1.); ak->d28c = q28(e0, f0, inc, bet, 0., 1.);
465 ak->d29c = q29(e0, f0, inc, bet, 0., 1.);
466
467 return 0;
468
469}
470
471
472static REAL8 Heaviside( REAL8 x){
473if( x >= 0 )
474 return 1;
475
476return 0;
477}
478
479
480
481/**
482 * @addtogroup LALSimInspiralEccentricityFD_c
483 * @brief Routines to generate frequency-domain eccentric inspiral waveforms.
484 * @{
485 */
486
487
488int XLALSimInspiralEFD(
489 COMPLEX16FrequencySeries **hptilde, /**< FD plus polarization */
490 COMPLEX16FrequencySeries **hctilde, /**< FD cross polarization */
491 const REAL8 phiRef, /**< Orbital coalescence phase (rad) */
492 const REAL8 deltaF, /**< Frequency resolution */
493 const REAL8 m1_SI, /**< Mass of companion 1 (kg) */
494 const REAL8 m2_SI, /**< Mass of companion 2 (kg) */
495 const REAL8 fStart, /**< Start GW frequency (Hz) */
496 const REAL8 fEnd, /**< Highest GW frequency (Hz): end at Schwarzschild ISCO */
497 const REAL8 i, /**< Polar inclination of source (rad) */
498 const REAL8 r, /**< Distance of source (m) */
499 const REAL8 inclination_azimuth, /**< Azimuthal component of inclination angles [0, 2 M_PI]*/
500 const REAL8 e_min, /**< Initial eccentricity at frequency f_min: range [0, 0.4] */
501 const INT4 phaseO /**< Twice PN phase order */
502 )
503
504 {
505
506 gsl_complex cphase, exphase, czeta_FPlus,czeta_FCross, czeta_FPlus_times_exphase,czeta_FCross_times_exphase;
507
508
509 const REAL8 m1 = m1_SI / LAL_MSUN_SI;
510 const REAL8 m2 = m2_SI / LAL_MSUN_SI;
511 const REAL8 m = m1 + m2;
512 const REAL8 Mtotal = m * LAL_MTSUN_SI; /* total mass in seconds */
513 const REAL8 eta = m1 * m2 / (m * m);
514 const REAL8 piM = LAL_PI * Mtotal;
515 const REAL8 vISCO = 1. / sqrt(6.);
516 const REAL8 fISCO = vISCO * vISCO * vISCO / piM;
517 const REAL8 fupper = 2.*fISCO;
518 const REAL8 mchirp= pow(eta, 3./5.)*Mtotal;
519 REAL8 shft, f_max;
520 REAL8 f;
521 const REAL8 zr=0.;
522 REAL8 Amplitude, zim, Phaseorder, pc_re_hplus, pc_im_hplus, pc_re_hcross, pc_im_hcross, re_hplus, im_hplus, re_hcross, im_hcross;
523 size_t j, n, jStart;
524
525 expnCoeffsEPC ak;
526
527 expnCoeffsEPC* ak_ptr = &ak;
528 EPCSetup( ak_ptr, m1, m2, fStart, i, inclination_azimuth, e_min);
529
530
531 COMPLEX16 *data_p = NULL;
532 COMPLEX16 *data_c = NULL;
533 LIGOTimeGPS tC = {0, 0};
534
535 COMPLEX16FrequencySeries *htilde_p;
536 COMPLEX16FrequencySeries *htilde_c;
537
538 /* Perform some initial checks */
539 if (!hptilde) XLAL_ERROR(XLAL_EFAULT);
540 if (*hptilde) XLAL_ERROR(XLAL_EFAULT);
541 if (!hctilde) XLAL_ERROR(XLAL_EFAULT);
542 if (*hctilde) XLAL_ERROR(XLAL_EFAULT);
543 if (m1_SI <= 0) XLAL_ERROR(XLAL_EDOM);
544 if (m2_SI <= 0) XLAL_ERROR(XLAL_EDOM);
545 if (fStart <= 0) XLAL_ERROR(XLAL_EDOM);
546 if (r <= 0) XLAL_ERROR(XLAL_EDOM);
547
548
549 /* allocate htilde_p and htilde_c*/
550 if ( fEnd == 0. ) // End at ISCO
551 f_max = fISCO;
552 else // End at user-specified freq.
553 f_max = fEnd;
554 n = (size_t) (f_max / deltaF + 1);
555 XLALGPSAdd(&tC, -1 / deltaF); /* coalesce at t=0 */
556
557
558 htilde_p = XLALCreateCOMPLEX16FrequencySeries("htilde_p: FD waveform", &tC, 0.0, deltaF, &lalStrainUnit, n);
559 if (!htilde_p) XLAL_ERROR(XLAL_EFUNC);
560 memset(htilde_p->data->data, 0, n * sizeof(COMPLEX16));
561 XLALUnitDivide(&htilde_p->sampleUnits, &htilde_p->sampleUnits, &lalSecondUnit);
562
563 htilde_c = XLALCreateCOMPLEX16FrequencySeries("htilde_c: FD waveform", &tC, 0.0, deltaF, &lalStrainUnit, n);
564 if (!htilde_c) XLAL_ERROR(XLAL_EFUNC);
565 memset(htilde_c->data->data, 0, n * sizeof(COMPLEX16));
566 XLALUnitDivide(&htilde_c->sampleUnits, &htilde_c->sampleUnits, &lalSecondUnit);
567
568
569 /* extrinsic parameters*/
570 Amplitude = -sqrt(5./384.)*pow(M_PI, -2./3.)*(pow(mchirp,5./6.)/r)*LAL_MRSUN_SI/LAL_MTSUN_SI;
571 shft = LAL_TWOPI * (tC.gpsSeconds + 1e-9 * tC.gpsNanoSeconds);
572
573 jStart = (size_t) ceil(fStart / deltaF);
574 f=jStart*deltaF;
575 data_p = htilde_p->data->data;
576 data_c = htilde_c->data->data;
577
578 /* In order to decompose the waveform in the form h = F_+ h_+ + F_x h_x we decompose the amplitude function using a two basis decomposition. Note that the functions zeta_real and zeta_im depend on several paramenters, including F_+ and F_x. Since zeta_real and zeta_im are linear function in the antenna pattern functions, czeta_FPlus and czeta_FCross are used to extract F_+ and F_x from the waveform amplitude*/
579
580
581 for ( j=jStart;j<n;j++) {
582 switch (phaseO) {
583 case -1:
584 case 7:
585 pc_re_hplus=0.;
586 pc_im_hplus=0.;
587 pc_re_hcross=0.;
588 pc_im_hcross=0.;
589 Phaseorder=0.;
590
591 for(int k=1; k<8; k++){
592 Phaseorder+=-PhaseAccTaylorF2(k, f, &ak);
593 }
594
595 for(int lm=1;lm<11;lm++){
596
597 zim=M_PI/4. + pow(((REAL8)lm)/2.,8./3.)*Phaseorder - shft*f + ((REAL8)lm)*phiRef;
598
599 cphase = gsl_complex_rect (zr,zim);
600
601 czeta_FPlus = gsl_complex_rect (((double)zeta_generic_re_plus(lm, f, &ak)),((double)zeta_generic_im_plus(lm, f, &ak)));
602
603 czeta_FCross=gsl_complex_rect (((double)zeta_generic_re_cross(lm, f, &ak)),((double)zeta_generic_im_cross(lm, f, &ak)));
604
605 exphase=gsl_complex_exp (cphase);
606
607 czeta_FPlus_times_exphase = gsl_complex_mul(czeta_FPlus, exphase);
608 czeta_FCross_times_exphase = gsl_complex_mul(czeta_FCross, exphase);
609
610
611 re_hplus=Heaviside(((double)lm)*fupper-2.*f)*pow(((double)lm)/2., 2./3.)*GSL_REAL (czeta_FPlus_times_exphase);
612 im_hplus=Heaviside(((double)lm)*fupper-2.*f)*pow(((double)lm)/2., 2./3.)*GSL_IMAG (czeta_FPlus_times_exphase);
613
614 re_hcross=Heaviside(((double)lm)*fupper-2.*f)*pow(((double)lm)/2., 2./3.)*GSL_REAL (czeta_FCross_times_exphase);
615 im_hcross=Heaviside(((double)lm)*fupper-2.*f)*pow(((double)lm)/2., 2./3.)*GSL_IMAG (czeta_FCross_times_exphase);
616
617
618 pc_re_hplus+=re_hplus;
619 pc_im_hplus+=im_hplus;
620
621 pc_re_hcross+=re_hcross;
622 pc_im_hcross+=im_hcross;
623 }
624 break;
625 default:
626 XLAL_ERROR(XLAL_ETYPE, "Invalid phase PN order %d", phaseO);
627 }
628
629
630
631 /*Note that h(f)= FPlus*data_p + FCross*data_c, where the polarizations are given by:*/
632
633 data_p[j] = Amplitude*pow(f,-7./6.)*(pc_re_hplus + pc_im_hplus*1.0j);
634 data_c[j] = Amplitude*pow(f,-7./6.)*(pc_re_hcross + pc_im_hcross*1.0j);
635
636 f+=deltaF;
637}
638
639 *hptilde = htilde_p;
640 *hctilde = htilde_c;
641 return XLAL_SUCCESS;
642
643}
644
645/** @} */
Heaviside
static REAL8 Heaviside(REAL8 x)
Definition: LALSimInspiralEccentricityFD.c:472
EPCSetup
static int EPCSetup(expnCoeffsEPC *ak, const REAL8 m1, const REAL8 m2, const REAL8 fStart, const REAL8 i, const REAL8 inclination_azimuth, const REAL8 e_min)
Definition: LALSimInspiralEccentricityFD.c:363
zeta_generic_re_plus
static REAL8 zeta_generic_re_plus(int k, REAL8 f, expnCoeffsEPC *ak)
Definition: LALSimInspiralEccentricityFD.c:140
zeta_generic_re_cross
static REAL8 zeta_generic_re_cross(int k, REAL8 f, expnCoeffsEPC *ak)
Definition: LALSimInspiralEccentricityFD.c:195
PhaseAccTaylorF2
static REAL8 PhaseAccTaylorF2(int k, REAL8 f, expnCoeffsEPC *ak)
Definition: LALSimInspiralEccentricityFD.c:89
zeta_generic_im_cross
static REAL8 zeta_generic_im_cross(int k, REAL8 f, expnCoeffsEPC *ak)
Definition: LALSimInspiralEccentricityFD.c:306
zeta_generic_im_plus
static REAL8 zeta_generic_im_plus(int k, REAL8 f, expnCoeffsEPC *ak)
Definition: LALSimInspiralEccentricityFD.c:251
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
i
double i
Definition: bh_ringdown.c:118
e
double e
Definition: bh_ringdown.c:117
k5
coeffs k5
Definition: exact_derivatives-HrealHeader.c:43
a4
const double a4
Definition: exact_derivatives-HrealHeader.c:34
a2
const double a2
Definition: exact_derivatives-HrealHeader.c:33
k4
coeffs k4
Definition: exact_derivatives-HrealHeader.c:42
XLALCreateCOMPLEX16FrequencySeries
COMPLEX16FrequencySeries * XLALCreateCOMPLEX16FrequencySeries(const CHAR *name, const LIGOTimeGPS *epoch, REAL8 f0, REAL8 deltaF, const LALUnit *sampleUnits, size_t length)
LAL_MSUN_SI
#define LAL_MSUN_SI
LAL_PI
#define LAL_PI
LAL_TWOPI
#define LAL_TWOPI
LAL_MTSUN_SI
#define LAL_MTSUN_SI
LAL_MRSUN_SI
#define LAL_MRSUN_SI
COMPLEX16
double complex COMPLEX16
REAL8
double REAL8
INT4
int32_t INT4
XLALSimInspiralEFD
int XLALSimInspiralEFD(COMPLEX16FrequencySeries **hptilde, COMPLEX16FrequencySeries **hctilde, const REAL8 phiRef, const REAL8 deltaF, const REAL8 m1_SI, const REAL8 m2_SI, const REAL8 fStart, const REAL8 fEnd, const REAL8 i, const REAL8 r, const REAL8 inclination_azimuth, const REAL8 e_min, const INT4 phaseO)
Definition: LALSimInspiralEccentricityFD.c:488
r
static const INT4 r
m
static const INT4 m
lalStrainUnit
const LALUnit lalStrainUnit
lalSecondUnit
const LALUnit lalSecondUnit
XLALUnitDivide
LALUnit * XLALUnitDivide(LALUnit *output, const LALUnit *unit1, const LALUnit *unit2)
XLAL_ERROR
#define XLAL_ERROR(...)
XLAL_SUCCESS
XLAL_SUCCESS
XLAL_EFAULT
XLAL_EFAULT
XLAL_EFUNC
XLAL_EFUNC
XLAL_EDOM
XLAL_EDOM
XLAL_ETYPE
XLAL_ETYPE
XLALGPSAdd
LIGOTimeGPS * XLALGPSAdd(LIGOTimeGPS *epoch, REAL8 dt)
x
x
COMPLEX16FrequencySeries::data
COMPLEX16Sequence * data
COMPLEX16Vector::data
COMPLEX16 * data
LIGOTimeGPS::gpsNanoSeconds
INT4 gpsNanoSeconds
f_max
double f_max
Definition: unicorn.c:23