Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALSimulation 6.2.0.1-b246709
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
LALSimInspiralPNCoefficients.c
Go to the documentation of this file.
1/*
2* Copyright (C) 2011 Drew Keppel, 2012 Riccardo Sturani
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 * Computes the PN Coefficients for using in the PN energy equation.
33 *
34 * Terms given in equation 3.1 of: Alessandra Buonanno, Bala R Iyer, Evan
35 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
36 * templates for compact binary inspiral signals in gravitational-wave
37 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
38 * For the spin terms a good reference are (3.15) and (3.16) of 1303.7412
39 *
40 * In the latest version coefficients of the terms n.S and L.S are reported
41 * "Averaged" spin coefficients refer to the ones obtained by orbital averaging,
42 * i.e. by using
43 * n_i n_j = 1/2 (\f$\delta_{ij} - \hat LN_i \hat LN_j\f$)
44 * However such orbital averaging at 2PN would introduce corrections
45 * at 3PN, as LNh is not constant.
46 */
47
48static REAL8 UNUSED
50 REAL8 eta)
51{
52 return -eta / 2.0;
53}
54
55static REAL8 UNUSED
57 REAL8 eta)
58{
59 return -(0.75 + eta/12.0);
60}
61
62static REAL8 UNUSED
64 REAL8 eta)
65{
66 return -(27.0/8.0 - 19.0/8.0 * eta + 1./24.0 * eta*eta);
67}
68
69static REAL8 UNUSED
71 REAL8 eta)
72{
73 return -(67.5/6.4 - (344.45/5.76 - 20.5/9.6 * LAL_PI*LAL_PI) * eta + 15.5/9.6 * eta*eta + 3.5/518.4 * eta*eta*eta);
74}
75
76static REAL8 UNUSED
78 REAL8 eta)
79{
80 return (-39.69/1.28 + (-123.671/5.76 + 9.037/1.536 *LAL_PI*LAL_PI+ 1792./15.*log(2)+89.6/1.5*LAL_GAMMA)* eta + (-498.449/3.456 +31.57/5.76*LAL_PI*LAL_PI)*eta*eta + 3.01/17.28 *eta*eta*eta + 7.7/3110.4*eta*eta*eta*eta);
81 /*see arXiv:1305.4884, or eq.(26) of arXiv:1309.3474
82 note that in getting a4 from PRD 62, 084011 (2000),
83 the first reference is using the fact that \omega_{static} = 0
84 (see arXiv:gr-qc/0105038) */
85}
86
87static REAL8 UNUSED
89 REAL8 eta)
90{
91 return 896./15.*eta;
92 /* arXiv:1305.4884 has a misprint, it should have a factor of nu
93 See for instance arXiv:1002.0726
94 Also note that this term is usually given as 448*log(x)/15
95 since x=v^2 the log(v) term is twice this */
96}
97
98/* Eq. (4.6) of arXiv:1212.5520
99 */
100static REAL8 UNUSED
102 REAL8 mByM)
103{
104 return 2. / 3. + 2. / mByM;
105}
106
107/* Eq. (6) of arXiv:astro-ph/0504538v2
108 */
109static REAL8 UNUSED
111 REAL8 eta)
112{
113 return 1./eta;
114}
115
116static REAL8 UNUSED
118 REAL8 eta)
119{
120 return -2./eta;
121}
122
123/* Eq. (6) of arXiv:astro-ph/0504538v2
124 */
125static REAL8 UNUSED
127 REAL8 eta)
128{
129 return -3./eta;
130}
131
132static REAL8 UNUSED
134 REAL8 UNUSED eta)
135{
136 return 6./eta;
137}
138
139/* Eq. (6) of arXiv:astro-ph/0504538v2
140 */
141static REAL8 UNUSED
143 REAL8 mByM)
144{
145 return .5/mByM/mByM;
146}
147
148/* Eq. (6) of arXiv:astro-ph/0504538v2
149 */
150static REAL8 UNUSED
152 REAL8 mByM)
153{
154 return -1.5/mByM/mByM;
155}
156
157static REAL8 UNUSED
159 REAL8 mByM)
160{
161 return -1./mByM/mByM;
162}
163
164static REAL8 UNUSED
166 REAL8 mByM)
167{
168 return 3./mByM/mByM;
169}
170
171/* Eq. 4.6 of arXiv:1212.5520
172 */
173static REAL8 UNUSED
175 REAL8 mByM)
176{
177 return 5./3. + 3./mByM + 29.*mByM/9. + mByM*mByM/9.;
178}
179
180/* From (3.30) of arXiv:1501.01529
181 */
182static REAL8 UNUSED
184 REAL8 eta)
185{
186 return -7./eta -1./3.;
187}
188
189/* From (3.30) of arXiv:1501.01529
190 */
191static REAL8 UNUSED
193 REAL8 eta)
194{
195 return 16./3./eta - 2./9.;
196}
197
198static REAL8 UNUSED
200 REAL8 eta)
201{
202 return 13./eta - 3.;
203}
204
205static REAL8 UNUSED
207 REAL8 eta)
208{
209 return 5./eta;
210}
211
212static REAL8 UNUSED
214 REAL8 eta)
215{
216 return 2./eta -11./6.;
217}
218
219static REAL8 UNUSED
221 REAL8 eta)
222{
223 return -11./3./eta + 2.3/1.8;
224}
225
226static REAL8 UNUSED
228 REAL8 mByM)
229{
230 return 2./(mByM*mByM) - 1./mByM -1.;
231}
232
233static REAL8 UNUSED
235 REAL8 mByM)
236{
237 return 3./(mByM*mByM) -2./(3.*mByM) -1./9.;
238}
239
240static REAL8 UNUSED
242 REAL8 mByM)
243{
244 return -8./(mByM*mByM)+11./3./mByM+1.;
245}
246
247static REAL8 UNUSED
249 REAL8 mByM)
250{
251 return 2./(mByM*mByM)-2./mByM;
252}
253
254static REAL8 UNUSED
256 REAL8 mByM)
257{
258 return -1./(mByM*mByM) - 1./6./mByM -0.5;
259}
260
261static REAL8 UNUSED
263 REAL8 mByM)
264{
265 return 6./(mByM*mByM) -1.5/mByM -1.1/1.8;
266}
267
268/* From (3.30) of arXiv:1501.01529
269 */
270static REAL8 UNUSED
272 REAL8 mByM)
273{
274 return -2.5/mByM/mByM - 2.5/mByM - 5./6.;
275}
276
277/* From (3.30) of arXiv:1501.01529
278 */
279static REAL8 UNUSED
281 REAL8 mByM)
282{
283 return 6.5/mByM/mByM + 8./5./mByM + 2.5;
284}
285
286static REAL8 UNUSED
288 REAL8 mByM)
289{
290 return 1./mByM/mByM -1./mByM;
291}
292
293static REAL8 UNUSED
295 REAL8 mByM)
296{
297 return 1.25/mByM/mByM + 1.25/mByM + 5./12.;
298}
299
300static REAL8 UNUSED
302 REAL8 mByM)
303{
304 return -3.75/mByM/mByM - 3.75/mByM - 1.25;
305}
306
307/* Eq. (4.6) of arXiv:1212.5520
308 * Symbol definitions right above eq. (3.1)
309 */
310static REAL8 UNUSED
312 REAL8 mByM)
313{
314 return -75./4. + 27./(4.*mByM) + 53.*mByM/2. + 67*mByM*mByM/6. + 17.*mByM*mByM*mByM/12. - mByM*mByM*mByM*mByM/12.;
315}
316
317/*
318 * Tidal correction coefficients to Energy
319 */
320
321static REAL8 UNUSED
323 REAL8 mByM)
324{
325 return -9.0 * mByM*mByM*mByM*mByM*(1.-mByM);
326}
327
328static REAL8 UNUSED
330 REAL8 mByM)
331{
332 return (-33./2. + 11./2.*mByM - 11./2.*mByM*mByM + 33./2.*mByM*mByM*mByM)*mByM*mByM*mByM*mByM;
333}
334
335/**
336 * Computes the flux PN Coefficients.
337 *
338 * Terms given in equation 3.2 of: Alessandra Buonanno, Bala R Iyer, Evan
339 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
340 * templates for compact binary inspiral signals in gravitational-wave
341 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
342 * For terms involving spins see eq.(3.13) of arXiv:1303.7412
343 */
344
345static REAL8 UNUSED
347 REAL8 eta)
348{
349 return 32.0 * eta*eta / 5.0;
350}
351
352
353static REAL8 UNUSED
355 REAL8 eta)
356{
357 return -(12.47/3.36 + 3.5/1.2 * eta);
358}
359
360static REAL8 UNUSED
362 REAL8 UNUSED eta)
363{
364 return 4.0 * LAL_PI;
365}
366
367static REAL8 UNUSED
369 REAL8 mByM)
370{
371 return -3./2. - 5./4./mByM;
372}
373
374static REAL8 UNUSED
376 REAL8 eta)
377{
378 return -(44.711/9.072 - 92.71/5.04 * eta - 6.5/1.8 * eta*eta);
379}
380
381static REAL8 UNUSED
383 REAL8 eta)
384{
385 return 31./8./eta;
386}
387
388static REAL8 UNUSED
390 REAL8 eta)
391{
392 return -15./eta;
393}
394
395static REAL8 UNUSED
397 REAL8 eta)
398{
399 return 71./24./eta;
400}
401
402static REAL8 UNUSED
404 REAL8 mByM)
405{
406 return 1./(16.*mByM*mByM);
407}
408
409static REAL8 UNUSED
411 REAL8 mByM)
412{
413 return 1./48./(mByM*mByM);
414}
415
416static REAL8 UNUSED
418 REAL8 eta)
419{
420 return -103./48./eta;
421}
422
423static REAL8 UNUSED
425 REAL8 eta)
426{
427 return 289./48./eta;
428}
429
430static REAL8 UNUSED
432 REAL8 mByM)
433{
434 return 7./96./(mByM*mByM);
435}
436
437static REAL8 UNUSED
439 REAL8 mByM)
440{
441 return -1./96./(mByM*mByM);
442}
443
444static REAL8 UNUSED
446 REAL8 mByM)
447{
448 return 2./mByM/mByM;
449}
450
451static REAL8 UNUSED
453 REAL8 mByM)
454{
455 return -7.5/mByM/mByM;
456}
457
458static REAL8 UNUSED
460 REAL8 mByM)
461{
462 return 1.5/mByM/mByM;
463}
464
465static REAL8 UNUSED
467 REAL8 mByM)
468{
469 return -1./mByM/mByM;
470}
471
472static REAL8 UNUSED
474 REAL8 mByM)
475{
476 return 3./mByM/mByM;
477}
478
479static REAL8 UNUSED
481 REAL8 eta)
482{
483 return -(81.91/6.72 + 58.3/2.4 * eta) * LAL_PI;
484}
485
486/* Eq. (4.9) of arXiv:1307.6793
487 */
488static REAL8 UNUSED
490 REAL8 mByM)
491{
492 return 63./8. - 13./(16.*mByM) - (73.*mByM)/36. - (157.*mByM*mByM)/18.;
493}
494
495static REAL8 UNUSED
497 REAL8 eta)
498{
499 return (664.3739519/6.9854400 + 16.0/3.0 * LAL_PI*LAL_PI - 17.12/1.05 * LAL_GAMMA - 17.12/1.05*log(4.) + (4.1/4.8 * LAL_PI*LAL_PI - 134.543/7.776) * eta - 94.403/3.024 * eta*eta - 7.75/3.24 * eta*eta*eta);
500}
501
502/* Note that this coefficient multiplies log(v)*/
503static REAL8 UNUSED
505 REAL8 UNUSED eta)
506{
507 return -17.12/1.05;
508}
509
510/* Eq. (4.9) of arXiv:1307.6793
511 * (symbol definitions around eq. 3.1)
512 */
513
514static REAL8 UNUSED
516 REAL8 mByM)
517{
518 return LAL_PI*( -17./3. - 31./(6.*mByM) );
519}
520
521/* From (4.12) of 1501.01529
522 */
523static REAL8 UNUSED
525 REAL8 eta)
526{
527 return -2.9/16.8/eta - 25.9/1.2;
528}
529
530static REAL8 UNUSED
532 REAL8 eta)
533{
534 return -49./6./eta + 44./9.;
535}
536
537static REAL8 UNUSED
539 REAL8 eta)
540{
541 return 349.9/4.2/eta + 117.59/1.68;
542}
543
544static REAL8 UNUSED
546 REAL8 eta)
547{
548 return -38.9/1.2/eta - 20.27/5.04;
549}
550
551static REAL8 UNUSED
553 REAL8 eta)
554{
555 return 212.3/8.4/eta + 82.1/7.2;
556}
557
558static REAL8 UNUSED
560 REAL8 eta)
561{
562 return -564.7/16.8/eta - 202.3/7.2;
563}
564
565static REAL8 UNUSED
567 REAL8 mByM)
568{
569 return -21./(8.*mByM*mByM) + 21.5/2.4/mByM - 1./24.;
570}
571
572static REAL8 UNUSED
574 REAL8 mByM)
575{
576 return -.5/(mByM*mByM) - 43./(6.*mByM) + 22./9.;
577}
578
579static REAL8 UNUSED
581 REAL8 mByM)
582{
583 return 23.53/(1.12*mByM*mByM) - 6.47/1.68/mByM + 2.27/3.36;
584}
585
586static REAL8 UNUSED
588 REAL8 mByM)
589{
590 return 8.81/(1.12*mByM*mByM) - 36.67/2.52/mByM + 6.1/100.8;
591}
592
593
594static REAL8 UNUSED
596 REAL8 mByM)
597{
598 return 18.9/(1.6*mByM*mByM) - 3.5/14.4/mByM +4.7/14.4;
599}
600
601
602static REAL8 UNUSED
604 REAL8 mByM)
605{
606 return -23.9/(1.6*mByM*mByM) + 2.93/1.44/mByM + 2.99/1.44;
607}
608
609static REAL8 UNUSED
611 REAL8 mByM)
612{
613 return -27.9/(5.6*mByM*mByM) - 45./(8.*mByM) + 43./4.;
614}
615
616static REAL8 UNUSED
618 REAL8 mByM)
619{
620 return 154.1/(4.2*mByM*mByM) + 15.17/(1.68*mByM) - 96.1/2.8;
621}
622
623static REAL8 UNUSED
625 REAL8 mByM)
626{
627 return -36.53/(1.68*mByM*mByM) + 65.9/(8.4*mByM) + 2.9/1.4;
628}
629
630
631static REAL8 UNUSED
633 REAL8 mByM)
634{
635 return 27.9/(11.2*mByM*mByM) + 4.5/(1.6*mByM) - 43./8.;
636}
637
638static REAL8 UNUSED
640 REAL8 mByM)
641{
642 return -8.37/(1.12*mByM*mByM) - 13.5/(1.6*mByM) + 12.9/8.;
643}
644
645/*
646 */
647static REAL8 UNUSED
649 REAL8 eta)
650{
651 return -(162.85/5.04 - 214.745/1.728 * eta - 193.385/3.024 * eta*eta) * LAL_PI;
652}
653
654/* Eq. (4.9) of arXiv:1307.6793
655 */
656static REAL8 UNUSED
658 REAL8 mByM)
659{
660 return (380.647/13.608) + 95.35/(3.36*mByM) - 401.15*mByM/7.56 + 3742.*mByM*mByM/63. - 35.*mByM*mByM*mByM/108. - 1117.*mByM*mByM*mByM*mByM/54.;
661}
662
663/* Eq. (4.9) of arXiv:1307.6793
664 */
665static REAL8 UNUSED
667 REAL8 mByM)
668{
669 return LAL_PI * (125.47/2.52 - 71.63/(6.72*mByM) -3.137*mByM/2.016 - 212.41*mByM*mByM/3.36);
670}
671
672/*
673 * Tidal correction coefficients to Flux
674 */
675
676static REAL8 UNUSED
678 REAL8 mByM)
679{
680 return 6. *(3. - 2.*mByM) * mByM*mByM*mByM*mByM;
681}
682
683static REAL8 UNUSED
685 REAL8 mByM)
686{
687 return (-176./7. - 1803./28.*mByM + 643./4.*mByM*mByM -155./2.*mByM*mByM*mByM) * mByM*mByM*mByM*mByM;
688}
689
690/* Non-spin phasing terms - see arXiv:0907.0700, Eq. 3.18 */
691static REAL8 UNUSED
693 REAL8 eta
694 )
695{
696 return 5.*(74.3/8.4 + 11.*eta)/9.;
697}
698
699static REAL8 UNUSED
701 REAL8 UNUSED eta)
702{
703 return -16.*LAL_PI;
704}
705
706static REAL8 UNUSED
708 REAL8 eta
709 )
710{
711 return 5.*(3058.673/7.056 + 5429./7.*eta+617.*eta*eta)/72.;
712}
713
714static REAL8 UNUSED
716 REAL8 eta
717 )
718{
719 return 5./9.*(772.9/8.4-13.*eta)*LAL_PI;
720}
721
722static REAL8 UNUSED
724 REAL8 eta
725 )
726{
727 return 5./3.*(772.9/8.4-13.*eta)*LAL_PI;
728}
729
731 REAL8 UNUSED eta
732 )
733{
734 return -684.8/2.1;
735}
736
737static REAL8 UNUSED
739 REAL8 eta
740 )
741{
742 return 11583.231236531/4.694215680 - 640./3.*LAL_PI*LAL_PI - 684.8/2.1*LAL_GAMMA + eta*(-15737.765635/3.048192 + 225.5/1.2*LAL_PI*LAL_PI) + eta*eta*76.055/1.728 - eta*eta*eta*127.825/1.296 + XLALSimInspiralTaylorF2Phasing_6PNLogCoeff(eta)*log(4.);
743}
744
745static REAL8 UNUSED
747 REAL8 eta
748 )
749{
750 return LAL_PI*(770.96675/2.54016 + 378.515/1.512*eta - 740.45/7.56*eta*eta);
751}
752
753/* Spin-orbit terms - can be derived from arXiv:1303.7412, Eq. 3.15-16 */
754
755static REAL8 UNUSED
757 REAL8 mByM
758 )
759{
760 return mByM*(25.+38./3.*mByM);
761}
762
763static REAL8 UNUSED
765 REAL8 mByM
766 )
767{
768 return -mByM*(1391.5/8.4-mByM*(1.-mByM)*10./3.+ mByM*(1276./8.1+mByM*(1.-mByM)*170./9.));
769
770}
771
772static REAL8 UNUSED
774 REAL8 mByM)
775{
776 return LAL_PI*mByM*(1490./3. + mByM*260.);
777}
778
779static REAL8 UNUSED
781 REAL8 mByM
782 )
783{
784 REAL8 eta=mByM*(1.-mByM);
785 return mByM*(-17097.8035/4.8384+eta*28764.25/6.72+eta*eta*47.35/1.44 + mByM*(-7189.233785/1.524096+eta*458.555/3.024-eta*eta*534.5/7.2));
786}
787
788/*
789 * Spin-squared corrections to TF2 phasing
790 * Compute 2.0PN SS, QM, and self-spin
791 * See Eq. (6.24) in arXiv:0810.5336
792 * 9b,c,d in arXiv:astro-ph/0504538
793 * Note that these terms are understood to multiply
794 * dimensionless spin magnitudes \chi_i=S_i/m_i^2
795 * differently from the convention adopted for SpinTaylorTX
796 * whose spinning coefficients multiply \chi_LAL=S_i/M^2
797 * where M=m_1+m_2.
798 * See also https://dcc.ligo.org/T1800298
799 */
800
801static REAL8 UNUSED
803 REAL8 eta
804 )
805{
806 return 247./4.8*eta;
807}
808
809static REAL8 UNUSED
811 REAL8 eta
812 )
813{
814 return -721./4.8*eta;
815}
816
817static REAL8 UNUSED
819 REAL8 mByM
820 )
821{
822 return -720./9.6*mByM*mByM;
823}
824
825static REAL8 UNUSED
827 REAL8 mByM
828 )
829{
830 return 1./9.6*mByM*mByM;
831}
832
833static REAL8 UNUSED
835 REAL8 mByM
836 )
837{
838 return 240./9.6*mByM*mByM;
839}
840
841static REAL8 UNUSED
843 REAL8 mByM
844 )
845{
846 return -7./9.6*mByM*mByM;
847}
848
849static REAL8 UNUSED
851 REAL8 eta
852 )
853{
854 return (326.75/1.12L + 557.5/1.8*eta)*eta;
855}
856
857static REAL8 UNUSED
859 REAL8 mByM
860 )
861{
862 return (-4108.25/6.72-108.5/1.2*mByM+125.5/3.6*mByM*mByM)*mByM*mByM;
863}
864
865static REAL8 UNUSED
867 REAL8 mByM
868 )
869{
870 return (4703.5/8.4+2935./6.*mByM-120.*mByM*mByM)*mByM*mByM;
871}
872
873/*
874 * Tidal corrections to F2 phasing
875 * See arXiv:1101.1673
876 */
877
878static REAL8 UNUSED
880 REAL8 mByM /**< ratio of object mass to total mass */
881 )
882{
883 return (-288. + 264.*mByM)*mByM*mByM*mByM*mByM;
884
885}
886
887static REAL8 UNUSED
889 REAL8 mByM /**< ratio of object mass to total mass */
890 )
891{
892 return (-15895./28. + 4595./28.*mByM + 5715./14.*mByM*mByM - 325./7.*mByM*mByM*mByM)*mByM*mByM*mByM*mByM;
893}
894
895static REAL8 UNUSED
897 REAL8 mByM /**< ratio of object mass to total mass */
898 )
899/* literature: Agathos et al (arxiv 1503.0545) eq (5)
900 * the coefficient mByM4 conversion & transformation (6.5PN, 7PN, 7.5PN):
901 * mByM=mA/M: mA= mass star A, M is total mass (mA+mB)
902 * Lambda (unitless) = lambda(m) / mA^5
903 * to call the function:
904 * Lambda * XLALSimInspiralTaylorF2Phasing_13PNTidalCoeff
905 * lambda(m)*mByM^4/mA^5= lambda(m)*(mA/M)^4/(mA)^5= lambda/(M^4*mA)
906 * =lambda/(mByM*M^5) eq (5)
907 */
908{
909 return mByM*mByM*mByM*mByM * 24.L*(12.L - 11.L*mByM)*LAL_PI;
910}
911
912static REAL8 UNUSED
914 REAL8 mByM /**< ratio of object mass to total mass */
915 )
916/* literature: Agathos et al (arxiv 1503.0545) eq (5)
917 * caveat: these are incomplete terms
918 * conversion see XLALSimInspiralTaylorF2Phasing_13PNTidalCoeff above
919 * --> completed by the terms given in equation (4) of :
920 * Tatsuya Narikawa, Nami Uchikata, Takahiro Tanaka,
921 * "Gravitational-wave constraints on the GWTC-2 events by measuring
922 * the tidal deformability and the spin-induced quadrupole moment",
923 * Phys. Rev. D 104, 084056 (2021), arXiv:2106.09193
924 */
925{
926 REAL8 mByM3 = mByM*mByM*mByM;
927 REAL8 mByM4 = mByM3 * mByM;
928 return - mByM4 * 5.L*(193986935.L/571536.L - 14415613.L/381024.L*mByM - 57859.L/378.L*mByM*mByM - 209495.L/1512.L*mByM3 + 965.L/54.L*mByM4 - 4.L*mByM4*mByM);
929}
930
931static REAL8 UNUSED
933 REAL8 mByM /**< ratio of object mass to total mass */
934 )
935/* literature: Agathos et al (arxiv 1503.0545) eq (5)
936 * conversion see XLALSimInspiralTaylorF2Phasing_13PNTidalCoeff above
937 * --> corrected by the terms given in equation (4) of :
938 * Tatsuya Narikawa, Nami Uchikata, Takahiro Tanaka,
939 * "Gravitational-wave constraints on the GWTC-2 events by measuring
940 * the tidal deformability and the spin-induced quadrupole moment",
941 * Phys. Rev. D 104, 084056 (2021), arXiv:2106.09193
942 */
943{
944 REAL8 mByM2 = mByM*mByM;
945 REAL8 mByM3 = mByM2*mByM;
946 REAL8 mByM4 = mByM3*mByM;
947 return mByM4 * 1.L/28.L*LAL_PI*(27719.L - 22415.L*mByM + 7598.L*mByM2 - 10520.L*mByM3) ;
948}
949
950/* The phasing function for TaylorF2 frequency-domain waveform.
951 * This function is tested in ../test/PNCoefficients.c for consistency
952 * with the energy and flux in this file.
953 */
954static void UNUSED
956 PNPhasingSeries *pfa, /**< \todo UNDOCUMENTED */
957 const REAL8 m1, /**< Mass of body 1, in Msol */
958 const REAL8 m2, /**< Mass of body 2, in Msol */
959 const REAL8 chi1L, /**< Component of dimensionless spin 1 along Lhat */
960 const REAL8 chi2L, /**< Component of dimensionless spin 2 along Lhat */
961 const REAL8 chi1sq,/**< Magnitude of dimensionless spin 1 */
962 const REAL8 chi2sq, /**< Magnitude of dimensionless spin 2 */
963 const REAL8 chi1dotchi2, /**< Dot product of dimensionles spin 1 and spin 2 */
964 LALDict *p /**< LAL dictionary containing accessory parameters */
965 )
966{
967 const REAL8 mtot = m1 + m2;
968 const REAL8 eta = m1*m2/mtot/mtot;
969 const REAL8 m1M = m1/mtot;
970 const REAL8 m2M = m2/mtot;
971
972 const REAL8 pfaN = 3.L/(128.L * eta);
973
974 memset(pfa, 0, sizeof(PNPhasingSeries));
975
976 pfa->v[0] = 1.L;
977 pfa->v[1] = 0.L;
986
987 /* modify the PN coefficients if a non null LALSimInspiralTestGRParam structure is passed */
988 /* BEWARE: this is for the non-spinning case only!*/
999
1002
1004 {
1008#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1009 __attribute__ ((fallthrough));
1010#endif
1012 pfa->v[6] += XLALSimInspiralTaylorF2Phasing_6PNSOCoeff(m1M)*chi1L
1017#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1018 __attribute__ ((fallthrough));
1019#endif
1021 pfa->v[5] += XLALSimInspiralTaylorF2Phasing_5PNSOCoeff(m1M)*chi1L
1023 pfa->vlogv[5] += 3.*(XLALSimInspiralTaylorF2Phasing_5PNSOCoeff(m1M)*chi1L
1025#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1026 __attribute__ ((fallthrough));
1027#endif
1029 /* 2PN SS, QM, and self-spin */
1035#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1036 __attribute__ ((fallthrough));
1037#endif
1040#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1041 __attribute__ ((fallthrough));
1042#endif
1046 break;
1047 default:
1048 XLALPrintError("XLAL Error - %s: Invalid spin PN order %i\n",
1051 break;
1052 }
1053
1057 {
1060#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1061 __attribute__ ((fallthrough));
1062#endif
1064#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1065 __attribute__ ((fallthrough));
1066#endif
1069#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1070 __attribute__ ((fallthrough));
1071#endif
1074#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1075 __attribute__ ((fallthrough));
1076#endif
1079#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1080 __attribute__ ((fallthrough));
1081#endif
1084#if __GNUC__ >= 7 && !defined __INTEL_COMPILER
1085 __attribute__ ((fallthrough));
1086#endif
1088 break;
1089 default:
1090 XLALPrintError("XLAL Error - %s: Invalid tidal PN order %i\n",
1093 }
1094
1095
1096 /* At the very end, multiply everything in the series by pfaN */
1097 for(int ii = 0; ii <= PN_PHASING_SERIES_MAX_ORDER; ii++)
1098 {
1099 pfa->v[ii] *= pfaN;
1100 pfa->vlogv[ii] *= pfaN;
1101 pfa->vlogvsq[ii] *= pfaN;
1102 }
1103}
1104
1105/**
1106 * Computes the PN Coefficients for using in the TaylorT2 phasing equation.
1107 *
1108 * Terms given in equation 3.8a of: Alessandra Buonanno, Bala R Iyer, Evan
1109 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
1110 * templates for compact binary inspiral signals in gravitational-wave
1111 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
1112 */
1113
1114static REAL8 UNUSED
1116 REAL8 eta)
1117{
1118 return -1./(32.*eta);
1119}
1120
1121static REAL8 UNUSED
1123 REAL8 eta)
1124{
1125 return 3.715/1.008 + 5.5/1.2 * eta;
1126}
1127
1128static REAL8 UNUSED
1130 REAL8 UNUSED eta)
1131{
1132 return -10. * LAL_PI;
1133}
1134
1135static REAL8 UNUSED
1137 REAL8 eta)
1138{
1139 return 15.293365/1.016064 + 27.145/1.008 * eta + 30.85/1.44 * eta*eta;
1140}
1141
1142static REAL8 UNUSED
1144 REAL8 eta)
1145{
1146 return (386.45/6.72 - 65./8. * eta) * LAL_PI;
1147}
1148
1149static REAL8 UNUSED
1151 REAL8 eta)
1152{
1153 return 1234.8611926451/1.8776862720 - 160./3. * LAL_PI*LAL_PI - 171.2/2.1 * LAL_GAMMA
1154 + (225.5/4.8 * LAL_PI*LAL_PI - 1573.7765635/1.2192768) * eta
1155 + 76.055/6.912 * eta*eta - 127.825/5.184 * eta*eta*eta;
1156}
1157
1158static REAL8 UNUSED
1160 REAL8 UNUSED eta)
1161{
1162 return -85.6/2.1;
1163}
1164
1165static REAL8 UNUSED
1167 REAL8 eta)
1168{
1169 return (77.096675/2.032128 + 37.8515/1.2096 * eta - 74.045/6.048 * eta*eta) * LAL_PI;
1170}
1171
1172/*
1173 * TaylorT2 derivatives dt/dv
1174 */
1175
1176/* The expression for dt/dv has an extra factor of M not implemented here */
1177static REAL8 UNUSED
1179 REAL8 eta)
1180{
1181 return 5./(32.*eta);
1182}
1183
1184static REAL8 UNUSED
1186 REAL8 eta)
1187{
1188 return 743./336. + 11.*eta/4.;
1189}
1190
1191static REAL8 UNUSED
1193 REAL8 UNUSED eta)
1194{
1195 return -4.*LAL_PI;
1196}
1197
1198static REAL8 UNUSED
1200 REAL8 eta)
1201{
1202 return 3058673./1016064. + 5429.*eta/1008. + 617.*eta*eta/144.;
1203}
1204
1205static REAL8 UNUSED
1207 REAL8 eta)
1208{
1209 return (-7729./672.+13.*eta/8.)*LAL_PI;
1210}
1211
1212static REAL8 UNUSED
1214 REAL8 eta)
1215{
1216 return -10817850546611./93884313600. + 32.*LAL_PI*LAL_PI/3.
1217 + 1712.*LAL_GAMMA/105.
1218 + (3147553127./12192768. - 451.*LAL_PI*LAL_PI/48.)*eta
1219 - 15211.*eta*eta/6912. + 25565.*eta*eta*eta/5184.
1220 + 856.*log(16.)/105.;
1221}
1222
1223/* The convention here is that this is the coefficient in front of v^6 log(v)
1224 * in the dt/dv expansion, NOT the one proportional to v^6 log(16 v^2).
1225 * Hence the term above containing log(16).
1226 */
1227static REAL8 UNUSED
1229 REAL8 UNUSED eta)
1230{
1231 return 1712./105.;
1232}
1233
1234static REAL8 UNUSED
1236 REAL8 eta)
1237{
1238 return LAL_PI*(-15419335./1016064. -75703.*eta/6048. + 14809.*eta*eta/3024);
1239}
1240
1241static REAL8 UNUSED
1243 REAL8 mByM)
1244{
1245 return 19./6. + 25./mByM/4.;
1246}
1247
1248static REAL8 UNUSED
1250 REAL8 mByM)
1251{
1252 return -17.*mByM*mByM/4. + 5.*mByM + 1249./36. + 8349./mByM/224.;
1253}
1254
1255static REAL8 UNUSED
1257 REAL8 mByM)
1258{
1259 return LAL_PI*( -13. - 149./mByM/6.);
1260}
1261
1262static REAL8 UNUSED
1264 REAL8 mByM)
1265{
1266 const REAL8 mByMsq = mByM*mByM;
1267 return 1069.*mByMsq*mByMsq/288. - 1741.*mByMsq*mByM/192. + 176383.*mByMsq/12096. + 707767.*mByM/3456. + 133100377./6096384. + 34195607./mByM/193536.;
1268}
1269
1270/* At 2 PN there are several spin^2 terms; see arXiv:astro-ph/0504538
1271 * The dt/dv spin^2 term at 2 PN is just -sigma (Eq. 9b-9d)
1272 * The terms 4PNSS and 4PNSSL are spin1-spin1 terms.
1273 */
1274static REAL8 UNUSED
1276 REAL8 eta)
1277{
1278 return -79./8./eta;
1279}
1280
1281static REAL8 UNUSED
1283 REAL8 eta)
1284{
1285 return -71./24./eta;
1286}
1287
1288static REAL8 UNUSED
1290 REAL8 eta)
1291{
1292 return 33./eta;
1293}
1294
1295static REAL8 UNUSED
1297 REAL8 eta)
1298{
1299 return 247./48./eta;
1300}
1301
1302static REAL8 UNUSED
1304 REAL8 eta)
1305{
1306 return -721./48./eta;
1307}
1308
1309static REAL8 UNUSED
1311 REAL8 mByM)
1312{
1313 return -1./16./mByM/mByM;
1314}
1315
1316static REAL8 UNUSED
1318 REAL8 mByM)
1319{
1320 return -1./48./mByM/mByM;
1321}
1322
1323static REAL8 UNUSED
1325 REAL8 mByM)
1326{
1327 return - 7./96./mByM/mByM;
1328}
1329
1330static REAL8 UNUSED
1332 REAL8 mByM)
1333{
1334 return 1./96./mByM/mByM;
1335}
1336
1337static REAL8 UNUSED
1339 REAL8 mByM)
1340{
1341 return -5./mByM/mByM;
1342}
1343
1344static REAL8 UNUSED
1346 REAL8 mByM)
1347{
1348 return 16.5/mByM/mByM;
1349}
1350
1351static REAL8 UNUSED
1353 REAL8 mByM)
1354{
1355 return -1.5/mByM/mByM;
1356}
1357
1358static REAL8 UNUSED
1360 REAL8 mByM)
1361{
1362 return 2.5/mByM/mByM;
1363}
1364
1365static REAL8 UNUSED
1367 REAL8 mByM)
1368{
1369 return -7.5/mByM/mByM;
1370}
1371
1372static REAL8 UNUSED
1374 REAL8 eta)
1375{
1376 return -98.173/(1.344*eta) - 46.1/2.4;
1377}
1378
1379static REAL8 UNUSED
1381 REAL8 eta)
1382{
1383 return 140.3/(2.4*eta) + 6.7/1.8;
1384}
1385
1386static REAL8 UNUSED
1388 REAL8 eta)
1389{
1390 return 373./(3.*eta) + 93.25/1.68;
1391}
1392
1393static REAL8 UNUSED
1395 REAL8 eta)
1396{
1397 return 140.699/(4.032*eta) - 32.11/2.52;
1398}
1399
1400static REAL8 UNUSED
1402 REAL8 mByM)
1403{
1404 return 27.565/(2.688*mByM*mByM) - 21.3/(1.6*mByM) - 17.3/4.8;
1405}
1406
1407static REAL8 UNUSED
1409 REAL8 mByM)
1410{
1411 return 32.5/(1.6*mByM*mByM) + 107./6./mByM + 67./36.;
1412}
1413
1414static REAL8 UNUSED
1416 REAL8 mByM)
1417{
1418 return -59.37/(1.12*mByM*mByM) + 103.7/5.6/mByM + 11.17/3.36;
1419}
1420
1421static REAL8 UNUSED
1423 REAL8 mByM)
1424{
1425 return 8.5/(806.4*mByM*mByM) + 6.485/1.008/mByM +2.9/50.4;
1426}
1427
1428static REAL8 UNUSED
1430 REAL8 mByM)
1431{
1432 return -94.07/(3.36*mByM*mByM) -58.7/2.4/mByM + 6.;
1433}
1434
1435static REAL8 UNUSED
1437 REAL8 mByM)
1438{
1439 return 563.9/(8.4*mByM*mByM) + 157.45/1.68/mByM - 171./7;
1440}
1441
1442static REAL8 UNUSED
1444 REAL8 mByM)
1445{
1446 return 56.65/(3.36*mByM*mByM) - 170.9/8.4/mByM + 45./7.;
1447}
1448
1449static REAL8 UNUSED
1451 REAL8 eta)
1452{
1453 return 52.973/(8.064*eta) + 3.13/1.44;
1454}
1455
1456static REAL8 UNUSED
1458 REAL8 eta)
1459{
1460 return -170.603/(8.064*eta) - 25.43/1.44;
1461}
1462
1463static REAL8 UNUSED
1465 REAL8 mByM)
1466{
1467 return -37.427/(2.304*mByM*mByM) - 2.41/(2.88*mByM) - 5.51/2.88;
1468}
1469
1470static REAL8 UNUSED
1472 REAL8 mByM)
1473{
1474 return 75.4979/(1.6128*mByM*mByM) + 15.43/2.88/mByM + 4.9/28.8;
1475}
1476
1477static REAL8 UNUSED
1479 REAL8 mByM)
1480{
1481 return 94.07/(6.72*mByM*mByM) + 58.7/4.8/mByM -3.;
1482}
1483
1484static REAL8 UNUSED
1486 REAL8 mByM)
1487{
1488 return -94.07/(2.24*mByM*mByM) -58.7/1.6/mByM + 9.;
1489}
1490
1491static REAL8 UNUSED
1493 REAL8 mByM)
1494{
1495 return 6.*mByM*mByM*mByM*mByM * (-12.+11.*mByM);
1496}
1497
1498static REAL8 UNUSED
1500 REAL8 mByM)
1501{
1502 return mByM*mByM*mByM*mByM * (-3179/8. + 919/8.*mByM + 1143/4.*mByM*mByM - 65./2.*mByM*mByM*mByM);
1503}
1504
1505/*
1506 * Tidal correction coefficients to Phasing
1507 */
1508
1509static REAL8 UNUSED
1511 REAL8 chi,
1512 REAL8 lambda)
1513{
1514 return lambda * (-66.*chi + 72.) * chi*chi*chi*chi;
1515}
1516
1517static REAL8 UNUSED
1519 REAL8 eta,
1520 REAL8 chi,
1521 REAL8 lambda)
1522{
1523 return lambda * chi*chi*chi*chi * ( -1497.5*chi/5.6 - 225.5*eta*chi/1.4 + 1589.5/5.6
1524 + 259.5*eta/1.4 + 398.5*chi*chi/2.8 - 965.*chi*chi*chi/7.);
1525}
1526
1527/**
1528 * Computes the PN Coefficients for using in the TaylorT2 timing equation.
1529 *
1530 * Terms given in equation 3.8b of: Alessandra Buonanno, Bala R Iyer, Evan
1531 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
1532 * templates for compact binary inspiral signals in gravitational-wave
1533 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
1534 */
1535
1536static REAL8 UNUSED
1538 REAL8 totalmass,
1539 REAL8 eta)
1540{
1541 totalmass *= LAL_G_SI / pow(LAL_C_SI, 3.0); /* convert totalmass from kilograms to seconds */
1542 return -5.*totalmass/(256.*eta);
1543}
1544
1545static REAL8 UNUSED
1547 REAL8 eta)
1548{
1549 return 7.43/2.52 + 11./3. * eta;
1550}
1551
1552static REAL8 UNUSED
1554 REAL8 UNUSED eta)
1555{
1556 return -32./5. * LAL_PI;;
1557}
1558
1559static REAL8 UNUSED
1561 REAL8 eta)
1562{
1563 return 30.58673/5.08032 + 54.29/5.04*eta + 61.7/7.2*eta*eta;
1564}
1565
1566static REAL8 UNUSED
1568 REAL8 eta)
1569{
1570 return -(77.29/2.52 -13./3.*eta) * LAL_PI;
1571}
1572
1573static REAL8 UNUSED
1575 REAL8 eta)
1576{
1577 return -1005.2469856691/2.3471078400 + 128./3. * LAL_PI*LAL_PI + 68.48/1.05 * LAL_GAMMA
1578 + (3147.553127/3.048192 - 45.1/1.2 * LAL_PI*LAL_PI) * eta
1579 - 15.211/1.728 * eta*eta + 25.565/1.296 * eta*eta*eta;
1580}
1581
1582static REAL8 UNUSED
1584 REAL8 UNUSED eta)
1585{
1586 return 34.24/1.05;
1587}
1588
1589static REAL8 UNUSED
1591 REAL8 eta)
1592{
1593 return (-154.19335/1.27008 - 757.03/7.56 * eta + 148.09/3.78 * eta*eta) * LAL_PI;
1594}
1595
1596/*
1597 * Tidal correction coefficients to Timing
1598 */
1599
1600static REAL8 UNUSED
1602 REAL8 chi,
1603 REAL8 lambda)
1604{
1605 return lambda * (-264.*chi + 288.) * chi*chi*chi*chi;
1606}
1607
1608static REAL8 UNUSED
1610 REAL8 eta,
1611 REAL8 chi,
1612 REAL8 lambda)
1613{
1614 return lambda * chi*chi*chi*chi * (-2995.*chi/4. - 451.*eta*chi + 3179./4. + 519.*eta
1615 + (797.*chi*chi)/2. - 386.*chi*chi*chi);
1616}
1617
1618
1619
1620/**
1621 * Computes the PN Coefficients for using in the TaylorT3 phasing equation.
1622 *
1623 * Terms given in equation 3.10a of: Alessandra Buonanno, Bala R Iyer, Evan
1624 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
1625 * templates for compact binary inspiral signals in gravitational-wave
1626 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
1627 */
1628
1629
1630static REAL8 UNUSED
1632 REAL8 eta)
1633{
1634 return -1./eta;
1635}
1636
1637static REAL8 UNUSED
1639 REAL8 eta)
1640{
1641 return 3.715/8.064 + 5.5/9.6 * eta;
1642}
1643
1644static REAL8 UNUSED
1646 REAL8 UNUSED eta)
1647{
1648 return -3./4. * LAL_PI;
1649}
1650
1651static REAL8 UNUSED
1653 REAL8 eta)
1654{
1655 return 9.275495/14.450688 + 2.84875/2.58048 * eta + 1.855/2.048 * eta*eta;
1656}
1657
1658static REAL8 UNUSED
1660 REAL8 eta)
1661{
1662 return (3.8645/2.1504 - 6.5/25.6 * eta) * LAL_PI;
1663}
1664
1665static REAL8 UNUSED
1667 REAL8 eta)
1668{
1669 return 83.1032450749357/5.7682522275840 - 5.3/4.0 * LAL_PI*LAL_PI - 10.7/5.6 * LAL_GAMMA
1670 + (-126.510089885/4.161798144 + 2.255/2.048 * LAL_PI*LAL_PI) * eta
1671 + 1.54565/18.35008 * eta*eta - 1.179625/1.769472 * eta*eta*eta;
1672}
1673
1674static REAL8 UNUSED
1676 REAL8 UNUSED eta)
1677{
1678 return -10.7/5.6;
1679}
1680
1681static REAL8 UNUSED
1683 REAL8 eta)
1684{
1685 return (1.88516689/1.73408256 + 4.88825/5.16096 * eta - 1.41769/5.16096 * eta*eta) * LAL_PI;
1686}
1687
1688/*
1689 * Tidal correction coefficients to Phasing
1690 */
1691
1692static REAL8 UNUSED
1694 REAL8 chi,
1695 REAL8 lambda)
1696{
1697 return lambda * (-3.3*chi/51.2 + 9./128.) * chi*chi*chi*chi;
1698}
1699
1700static REAL8 UNUSED
1702 REAL8 eta,
1703 REAL8 chi,
1704 REAL8 lambda)
1705{
1706 return lambda * chi*chi*chi*chi * (-1.30715*chi/13.76256
1707 - 8.745*eta*chi/114.688 + 2.3325/22.9376 + 4.905*eta/57.344
1708 + 3.985*chi*chi/114.688 - 9.65*chi*chi*chi/286.72);
1709}
1710
1711
1712
1713/**
1714 * Computes the PN Coefficients for using in the TaylorT3 frequency equation.
1715 *
1716 * Terms given in equation 3.10b of: Alessandra Buonanno, Bala R Iyer, Evan
1717 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
1718 * templates for compact binary inspiral signals in gravitational-wave
1719 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
1720 */
1721
1722static REAL8 UNUSED
1724 REAL8 totalmass)
1725{
1726 totalmass *= LAL_G_SI / pow(LAL_C_SI, 3.0); /* convert totalmass from kilograms to seconds */
1727 return 1. / (8. * LAL_PI * totalmass);
1728}
1729
1730static REAL8 UNUSED
1732 REAL8 eta)
1733{
1734 return 7.43/26.88 + 1.1/3.2 * eta;
1735}
1736
1737static REAL8 UNUSED
1739 REAL8 UNUSED eta)
1740{
1741 return -3./10. * LAL_PI;
1742}
1743
1744static REAL8 UNUSED
1746 REAL8 eta)
1747{
1748 return 1.855099/14.450688 + 5.6975/25.8048 * eta + 3.71/20.48 * eta*eta;
1749}
1750
1751static REAL8 UNUSED
1753 REAL8 eta)
1754{
1755 return (-7.729/21.504 + 1.3/25.6 * eta) * LAL_PI;
1756}
1757
1758static REAL8 UNUSED
1760 REAL8 eta)
1761{
1762 return -7.20817631400877/2.88412611379200 + 5.3/20.0 * LAL_PI*LAL_PI + 1.07/2.80 * LAL_GAMMA
1763 + (25.302017977/4.161798144 - 4.51/20.48 * LAL_PI*LAL_PI) * eta
1764 - 3.0913/183.5008 * eta*eta + 2.35925/17.69472 * eta*eta*eta;
1765}
1766
1767static REAL8 UNUSED
1769 REAL8 UNUSED eta)
1770{
1771 return 1.07/2.80;
1772}
1773
1774static REAL8 UNUSED
1776 REAL8 eta)
1777{
1778 return (-1.88516689/4.33520640 - 9.7765/25.8048 * eta + 1.41769/12.90240 * eta*eta) * LAL_PI;
1779}
1780
1781/*
1782 * Tidal correction coefficients to Frequency
1783 */
1784
1785static REAL8 UNUSED
1787 REAL8 chi,
1788 REAL8 lambda)
1789{
1790 return lambda * (-9.9*chi/102.4 + 2.7/25.6) * chi*chi*chi*chi;
1791}
1792
1793static REAL8 UNUSED
1795 REAL8 eta,
1796 REAL8 chi,
1797 REAL8 lambda)
1798{
1799 return lambda * chi*chi*chi*chi * (-8.579*chi/65.536 - 1.947*eta*chi/16.384
1800 + 1.8453/13.1072 + 4.329*eta/32.768 + 2.391*chi*chi/65.536
1801 - 5.79*chi*chi*chi/163.84);
1802}
1803
1804/**
1805 * Computes the PN Coefficients for using in the TaylorT4 frequency equation.
1806 *
1807 * Spin-less terms given in equation 3.6 of: Alessandra Buonanno, Bala R Iyer,
1808 * Evan Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
1809 * templates for compact binary inspiral signals in gravitational-wave
1810 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
1811 * ( with (domega/dt)/omega = 3 (dv/dt)/v ),
1812 * spin-orbit terms in eq. 4.10a of S.Marsat, A.Bohe, L.Blanchet, A.Buonanno
1813 * Class.Quant.Grav. 31 (2014) 025023, arXiv:1307.6793,
1814 * spin-spin terms in eqs. 5.15-17 of E.Racine, A.Buonanno, L.Kidder PRD80(2009)
1815 * 044010, arXiv:0812.4413
1816 *** UNREVIEWED ****
1817 */
1818
1819static REAL8 UNUSED
1821 REAL8 eta)
1822{
1823 return 96. / 5. * eta;
1824}
1825
1826static REAL8 UNUSED
1828 REAL8 eta)
1829{
1830 return ( -(1.0/336.0) * (743.0 + 924.0*eta) );
1831}
1832
1833static REAL8 UNUSED
1835 REAL8 UNUSED eta)
1836{
1837 return 4.0 * LAL_PI;
1838}
1839
1840static REAL8 UNUSED
1842 REAL8 mByM)
1843{
1844 return - 19./6. - 25./4./mByM;
1845}
1846
1847static REAL8 UNUSED
1849 REAL8 eta)
1850{
1851 return (34103. + 122949.*eta + 59472.*eta*eta)/18144.;
1852}
1853
1854static REAL8 UNUSED
1856 REAL8 eta)
1857{
1858 return 79. / (8. * eta);
1859}
1860
1861static REAL8 UNUSED
1863 REAL8 eta)
1864{
1865 return -33. / eta;
1866}
1867
1868static REAL8 UNUSED
1870 REAL8 eta)
1871{
1872 return 7.1 / 2.4 / eta;
1873}
1874
1875static REAL8 UNUSED
1877 REAL8 eta)
1878{
1879 return - 247. / 48. / eta;
1880}
1881
1882static REAL8 UNUSED
1884 REAL8 eta)
1885{
1886 return 721. / 48. / eta;
1887}
1888
1889static REAL8 UNUSED
1891 REAL8 mByM)
1892{
1893 return 1./(16.*mByM*mByM);
1894}
1895
1896static REAL8 UNUSED
1898 REAL8 mByM)
1899{
1900 return 1./48./mByM/mByM;
1901}
1902
1903static REAL8 UNUSED
1905 REAL8 mByM)
1906{
1907 return 7./96./mByM/mByM;
1908}
1909
1910static REAL8 UNUSED
1912 REAL8 mByM)
1913{
1914 return -1./96./mByM/mByM;
1915}
1916
1917static REAL8 UNUSED
1919 REAL8 mByM)
1920{
1921 return 5./mByM/mByM;
1922}
1923
1924static REAL8 UNUSED
1926 REAL8 mByM)
1927{
1928 return -16.5/mByM/mByM;
1929}
1930
1931static REAL8 UNUSED
1933 REAL8 mByM)
1934{
1935 return 1.5/mByM/mByM;
1936}
1937
1938static REAL8 UNUSED
1940 REAL8 mByM)
1941{
1942 return -2.5/mByM/mByM;
1943}
1944
1945static REAL8 UNUSED
1947 REAL8 mByM)
1948{
1949 return 7.5/mByM/mByM;
1950}
1951
1952static REAL8 UNUSED
1954 REAL8 eta)
1955{
1956 return ( -(1.0/672.0) * LAL_PI * (4159.0 + 15876.0*eta) );
1957 /* coefficient 15876 corrected (from 14532) according
1958 to 2005 erratas for L. Blanchet, Phys. Rev. D 54, 1417 (1996)
1959 (see Phys. Rev. D 71 129904 (E) (2005)) and L. Blanchet,
1960 B. R. Iyer, and B. Joguet, Phys. Rev. D 65, 064005 (2002)
1961 (see Phys. Rev. D 71 129903 (E) (2005)).
1962 See errata for Arun et al., Phys. Rev. D 71, 084008
1963 (2005) (see Phys. Rev. D 72 069903 (E) (2005))
1964 for corrected coefficients
1965 */
1966}
1967
1968static REAL8 UNUSED
1970 REAL8 mByM)
1971{
1972 return -809./(84.*mByM) + 13.795/1.008 - 527.*mByM/24. - 79.*mByM*mByM/6.;
1973}
1974
1975static REAL8 UNUSED
1977 REAL8 eta)
1978{
1979 return ( 16447.322263/139.7088 - 1712./105.
1980 * LAL_GAMMA - 561.98689/2.17728 * eta + LAL_PI * LAL_PI
1981 * (16./3. + 451./48. * eta) + 541./896. * eta * eta
1982 - 5605./2592. * eta * eta * eta - 856./105. * log(16.) );
1983}
1984
1985/* The coefficient of the log is normalized for the argument of the log to be v=(M omega)^(1/3) */
1986static REAL8 UNUSED
1988 REAL8 UNUSED eta)
1989{
1990 return -(1712.0/105.0);
1991}
1992
1993static REAL8 UNUSED
1995 REAL8 mByM)
1996{
1997 return LAL_PI * ( -37./3. - 151./6./mByM );
1998}
1999
2000static REAL8 UNUSED
2002 REAL8 eta)
2003{
2004 return 98.69/(3.36*eta) - 168.5/4.8;
2005}
2006
2007static REAL8 UNUSED
2009 REAL8 eta)
2010{
2011 return 237./(4.*eta) + 49./3.;
2012}
2013
2014static REAL8 UNUSED
2016 REAL8 eta)
2017{
2018 return 36.31/(1.68*eta) + 211.67/1.68;
2019}
2020
2021static REAL8 UNUSED
2023 REAL8 eta)
2024{
2025 return -230.3/(4.8*eta) - 3.557/1.008;
2026}
2027
2028static REAL8 UNUSED
2030 REAL8 eta)
2031{
2032 return 108.79/(6.72*eta) + 75.25/2.88;
2033}
2034
2035static REAL8 UNUSED
2037 REAL8 eta)
2038{
2039 return 162.25/(2.24*eta) - 129.31/2.88;
2040}
2041
2042static REAL8 UNUSED
2044 REAL8 mByM)
2045{
2046 return -33.7/(3.2*mByM*mByM) + 41.5/(3.2*mByM) + 37.9/9.6;
2047}
2048
2049static REAL8 UNUSED
2051 REAL8 mByM)
2052{
2053 return 75./(4.*mByM*mByM) + 87./(4.*mByM) + 49./6.;
2054}
2055
2056static REAL8 UNUSED
2058 REAL8 mByM)
2059{
2060 return 59.37/(1.12*mByM*mByM) - 103.7/(5.6*mByM) - 11.17/3.36;
2061}
2062
2063static REAL8 UNUSED
2065 REAL8 mByM)
2066{
2067 return -2.3/(22.4*mByM*mByM) -13.201/(2.016*mByM) + 1.15/20.16;
2068}
2069
2070static REAL8 UNUSED
2072 REAL8 mByM)
2073{
2074 return 101.9/(6.4*mByM*mByM) + 2.51/(5.76*mByM) + 13.33/5.76;
2075}
2076
2077static REAL8 UNUSED
2079 REAL8 mByM)
2080{
2081 return -49.3/(6.4*mByM*mByM) + 197.47/(5.76*mByM) + 56.45/5.76;
2082}
2083
2084static REAL8 UNUSED
2086 REAL8 mByM)
2087{
2088 return 6.59/(1.12*mByM*mByM) - 7.3/(2.4*mByM) + 21.5;
2089}
2090
2091static REAL8 UNUSED
2093 REAL8 mByM)
2094{
2095 return 19.63/(3.36*mByM*mByM) - 4.99/(1.68*mByM) - 185.7/2.8;
2096}
2097
2098static REAL8 UNUSED
2100 REAL8 mByM)
2101{
2102 return -39.47/(1.68*mByM*mByM) + 25.4/(2.1*mByM) + 5.1/2.8;
2103}
2104
2105static REAL8 UNUSED
2107 REAL8 mByM)
2108{
2109 return -6.59/(2.24*mByM*mByM) + 7.3/(4.8*mByM) - 43./4.;
2110}
2111
2112static REAL8 UNUSED
2114 REAL8 mByM)
2115{
2116 return 19.77/(2.24*mByM*mByM) - 7.3/(1.6*mByM) + 129./4.;
2117}
2118
2119static REAL8 UNUSED
2121 REAL8 eta)
2122{
2123 return (LAL_PI/12096.0) * (-13245.0 + 717350.0*eta + 731960.0*eta*eta);
2124 /* coefficients 717350 and 731960 corrected (from 661775 and 599156) according
2125 to 2005 erratas for L. Blanchet, Phys. Rev. D 54, 1417 (1996)
2126 (see Phys. Rev. D 71 129904 (E) (2005)) and L. Blanchet,
2127 B. R. Iyer, and B. Joguet, Phys. Rev. D 65, 064005 (2002)
2128 (see Phys. Rev. D 71 129903 (E) (2005)).
2129 See errata for Arun et al., Phys. Rev. D 71, 084008
2130 (2005) (see Phys. Rev. D 72 069903 (E) (2005))
2131 for corrected coefficients
2132 */
2133}
2134
2135static REAL8 UNUSED
2137 REAL8 mByM)
2138{
2139 return -1195.759 / 18.144 / mByM + 2694.373 / 18.144 - 432.2 / 2.1 * mByM + 1425.7 / 86.4 *mByM*mByM - 351.05 / 8.64 *mByM*mByM*mByM - 108.19 / 4.32 *mByM*mByM*mByM*mByM;
2140}
2141
2142static REAL8 UNUSED
2144 REAL8 mByM)
2145{
2146 return LAL_PI*(266.519/2.016 - 166.5/(2.8*mByM) - 828.43*mByM/6.72 -343.03*mByM*mByM/3.36);
2147 // see eq.(4.10a) of arXiv:1307.6793
2148}
2149
2150/*
2151 * Tidal correction coefficients to DOmega/dt
2152 */
2153
2154static REAL8 UNUSED
2156 REAL8 mByM)
2157{
2158 return 6.*mByM*mByM*mByM*mByM * (12.-11.*mByM);
2159}
2160
2161static REAL8 UNUSED
2163 REAL8 mByM)
2164{
2165 return mByM*mByM*mByM*mByM * (4421./56. - 12263./56.*mByM + 1893./4.*mByM*mByM - 661./2.*mByM*mByM*mByM);
2166}
2167
2168/*
2169 * For L see eq. 2.9 of arxiv:gr-qc/9506022
2170 */
2171
2172static REAL8 UNUSED
2174 REAL8 eta,
2175 REAL8 v)
2176{
2177 return M*M*eta/v;
2178}
2179/* eq. 4.7 of http://arxiv.org/pdf/1212.5520.pdf */
2180static REAL8 UNUSED
2182 REAL8 eta)
2183{
2184 return 1.5 + eta/6.;
2185}
2186
2187/* Orbital averaged from eq.(2.9c) of
2188 * \cite Kidder:1995zr, see also eq.(4.7) of \cite Bohe:2012mr.
2189 * Explicit formula can be found in eq.(9) of https://dcc.ligo.org/T1500554/public.
2190 */
2191
2192static REAL8 UNUSED
2194 REAL8 mByM)
2195{
2196 return -(1.+1./mByM);
2197}
2198
2199static REAL8 UNUSED
2201 REAL8 mByM)
2202{
2203 return 0.5*(1.+3./mByM);
2204}
2205
2206static REAL8 UNUSED
2208 REAL8 mByM)
2209{
2210 return 0.5*(1./3.-3./mByM);
2211}
2212
2213static REAL8 UNUSED
2215 REAL8 mByM)
2216{
2217 return -0.75-0.25/mByM;
2218}
2219
2220static REAL8 UNUSED
2222 REAL8 mByM)
2223{
2224 return -(1./3.+9./mByM)/4.;
2225}
2226
2227static REAL8 UNUSED
2229 REAL8 mByM)
2230{
2231 REAL8 eta=mByM*(1.-mByM);
2232 return -0.5*(5.+1./mByM)+eta/3.*(1.+4./mByM);
2233}
2234
2235static REAL8 UNUSED
2237 REAL8 mByM)
2238{
2239 REAL8 eta=mByM*(1.-mByM);
2240 return 3./8.*(3.+5./mByM)-7.*eta/4.*(1./6.+1./mByM);
2241}
2242
2243static REAL8 UNUSED
2245 REAL8 mByM)
2246{
2247 REAL8 eta=mByM*(1.-mByM);
2248 return -1./8.*(15.+17./mByM)+eta/4.*(-1.7/1.8+19./3./mByM);
2249}
2250
2251/* eq. 4.7 of http://arxiv.org/pdf/1212.5520.pdf */
2252static REAL8 UNUSED
2254 REAL8 eta)
2255{
2256 return 27./8. - 19./8.*eta + eta*eta/24.;
2257}
2258
2259static REAL8 UNUSED
2261 REAL8 eta)
2262{
2263 return 13.5/1.6 + (-68.89/1.44 + 4.1/2.4 * LAL_PI*LAL_PI)*eta + 3.1/2.4*eta*eta + 7./1296.*eta*eta*eta;
2264}
2265
2266/*
2267 * dLh
2268 *
2269 * \f$d \hat{L_N}/d \hat{t} = M * d\hat{L_N} / dt = \Omega_L x \hat{L_N}\f$
2270 * This is Eq. (10) of gr-qc/0405090 ( times M b/c we use \f$\hat{t}\f$)
2271 */
2272
2273static REAL8 UNUSED
2275 REAL8 mByM)
2276{
2277 return 0.5+1.5/mByM;
2278}
2279
2280/* Using spin-self^2 derivatives at v^6 from
2281 * eq. A.2 of Blanchet et al. 1501.01529
2282 * and relating to LNh derivative through (4.25)
2283 * of arXiv:0812.4413.
2284 */
2285static REAL8 UNUSED
2287 REAL8 eta)
2288{
2289 return -1.5/eta;
2290}
2291
2292static REAL8 UNUSED
2294 REAL8 mByM)
2295{
2296 return -1.5/(mByM*mByM);
2297}
2298
2299/* Using spin derivatives at v^7 from
2300 * eq. 7.8 of Blanchet et al. gr-qc/0605140
2301 * and relating to LNh derivative through (4.25)
2302 * of arXiv:0812.4413.
2303 */
2304static REAL8 UNUSED
2306 REAL8 mByM)
2307{
2308 return ( 9./8./mByM + 5./8 + 29./24.*mByM +mByM*mByM/24.);
2309}
2310
2311// See (3.4) of arXiv:1212.5520
2312static REAL8 UNUSED
2314 REAL8 mByM)
2315{
2316 return -7.5/1.6 + 2.7/(1.6*mByM) + 53.*mByM/8. + 6.7*mByM*mByM/2.4 + 1.7*mByM*mByM*mByM/4.8 - mByM*mByM*mByM*mByM/48.;
2317}
2318
2319/*
2320 * dS1
2321 * d S_1 / d \hat{t} = M * d S_1 / dt = \Omega_{S1,S2,LN,v} x S_1
2322 * However, that paper uses spin variables which are M^2 times our spins
2323 */
2324
2325/* dS1, 1.5PN: eq. (8) of gr-qc/0405090.
2326 */
2327static REAL8 UNUSED
2329 REAL8 mByM)
2330{
2331 return 3./2. -mByM - mByM*mByM/2.;
2332}
2333
2334/* S1S2 contribution
2335 * see. eq. A.2 of arXiv:1501.01529
2336 */
2337static const REAL8 UNUSED
2339
2340static const REAL8 UNUSED
2342
2343static const REAL8 UNUSED
2345
2346static const REAL8 UNUSED
2348
2349/* S1S1 contribution
2350 * again eq. A.2 of arXiv:1501.01529
2351 */
2352static REAL8 UNUSED
2354 REAL8 mByM)
2355{
2356 return 1.5 * (1. - 1./mByM);
2357}
2358
2359static REAL8 UNUSED
2361 REAL8 mByM)
2362{
2363 return 3. * (1./mByM - 1.);
2364}
2365
2366/* dS1, 2.5PN
2367 * eq. 7.8 of Blanchet et al. gr-qc/0605140
2368 */
2369static REAL8 UNUSED
2371 REAL8 mByM)
2372{
2373 return 9./8. - mByM/2. + 7.*mByM*mByM/12. - 7.*mByM*mByM*mByM/6. - mByM*mByM*mByM*mByM/24.;
2374}
2375
2376/* S1S2 contribution
2377 * again eq. A.2 of arXiv:1501.01529
2378 */
2379
2380static REAL8 UNUSED
2382 REAL8 mByM)
2383{
2384 return -1.5 -mByM;
2385}
2386
2387static REAL8 UNUSED
2389 REAL8 mByM)
2390{
2391 return 3.5-3./mByM-.5*mByM*mByM;
2392}
2393
2394static REAL8 UNUSED
2396 REAL8 mByM)
2397{
2398 return 1.5 +2.*mByM+mByM*mByM;
2399}
2400
2401static REAL8 UNUSED
2403 REAL8 mByM)
2404{
2405 return 3. -1.5*mByM-1.5/mByM;
2406}
2407
2408static REAL8 UNUSED
2410 REAL8 mByM)
2411{
2412 return 1.5 +mByM;
2413}
2414
2415static REAL8 UNUSED
2417 REAL8 mByM)
2418{
2420}
2421
2422static REAL8 UNUSED
2424 REAL8 mByM)
2425{
2427}
2428
2429static REAL8 UNUSED
2431 REAL8 mByM)
2432{
2434}
2435
2436/* S1S1 contribution
2437 * again eq. A.2 of arXiv:1501.01529
2438 */
2439
2440static REAL8 UNUSED
2442 REAL8 mByM)
2443{
2444 return 3. * (.5/mByM + 1. - mByM - .5*mByM*mByM);
2445}
2446
2447static REAL8 UNUSED
2449 REAL8 mByM)
2450{
2451 return 3. * (1./mByM -1.);
2452}
2453
2454static REAL8 UNUSED
2456 REAL8 mByM)
2457{
2459}
2460
2461/* dS1, 3.5PN
2462 * eq. 3.4 of Bohe' et al. arXiv:1212.5520
2463 */
2464
2465static REAL8 UNUSED
2467 REAL8 mByM)
2468{
2469 return (mByM*mByM*mByM*mByM*mByM*mByM/48. - 3./8.*mByM*mByM*mByM*mByM*mByM - 3.9/1.6*mByM*mByM*mByM*mByM - 23./6.*mByM*mByM*mByM +18.1/1.6*mByM*mByM -51./8.*mByM + 2.7/1.6);
2470}
2471
2472/**
2473 * Computes the PN Coefficients for using in the TaylorEt v(zeta) equation,
2474 * which is the square root of the x(zeta) equation.
2475 *
2476 * Terms given in equation 3.11 of: Alessandra Buonanno, Bala R Iyer, Evan
2477 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
2478 * templates for compact binary inspiral signals in gravitational-wave
2479 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
2480 */
2481
2482static REAL8 UNUSED
2484 REAL8 eta)
2485{
2486 return (3.0/4.0 + 1.0/12.0 * eta);
2487}
2488
2489static REAL8 UNUSED
2491 REAL8 eta)
2492{
2493 return (9.0/2.0 - 17.0/8.0 * eta + 1.0/18.0 * eta*eta);
2494}
2495
2496static REAL8 UNUSED
2498 REAL8 eta)
2499{
2500 return (40.5/1.6 + (20.5/9.6 * LAL_PI*LAL_PI - 479.5/7.2) * eta
2501 + 5.5/6.4 * eta*eta + 3.5/129.6 * eta*eta*eta);
2502}
2503
2504
2505/**
2506 * Computes the PN Coefficients for using in the TaylorEt dPhase/dt equation.
2507 *
2508 * Terms given in equation 3.13a of: Alessandra Buonanno, Bala R Iyer, Evan
2509 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
2510 * templates for compact binary inspiral signals in gravitational-wave
2511 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
2512 */
2513
2514static REAL8 UNUSED
2516 REAL8 m)
2517{
2518 return 1.0/m;
2519}
2520
2521static REAL8 UNUSED
2523 REAL8 eta)
2524{
2525 return (9.0/8.0 + 1.0/8.0 * eta);
2526}
2527
2528static REAL8 UNUSED
2530 REAL8 eta)
2531{
2532 return (8.91/1.28 - 20.1/6.4 * eta + 1.1/12.8 * eta*eta);
2533}
2534
2535static REAL8 UNUSED
2537 REAL8 eta)
2538{
2539 return (41.445/1.024 - (309.715/3.072 - 20.5/6.4 * LAL_PI*LAL_PI) * eta
2540 + 1.215/1.024 * eta*eta + 4.5/102.4 * eta*eta*eta);
2541}
2542
2543
2544/**
2545 * Computes the PN Coefficients for using in the TaylorEt dZeta/dt equation.
2546 *
2547 * Terms given in equation 3.13b of: Alessandra Buonanno, Bala R Iyer, Evan
2548 * Ochsner, Yi Pan, and B S Sathyaprakash, "Comparison of post-Newtonian
2549 * templates for compact binary inspiral signals in gravitational-wave
2550 * detectors", Phys. Rev. D 80, 084043 (2009), arXiv:0907.0700v1
2551 */
2552
2553static REAL8 UNUSED
2555 REAL8 m,
2556 REAL8 eta)
2557{
2558 return 64.0 * eta / (5.0 * m);
2559}
2560
2561static REAL8 UNUSED
2563 REAL8 eta)
2564{
2565 return (1.3/33.6 - 5.0/2.0 * eta);
2566}
2567
2568static REAL8 UNUSED
2570 REAL8 UNUSED eta)
2571{
2572 return 4.0 * LAL_PI;
2573}
2574
2575static REAL8 UNUSED
2577 REAL8 eta)
2578{
2579 return (11.7857/1.8144 - 12.017/2.016 * eta + 5.0/2.0 * eta*eta);
2580}
2581
2582static REAL8 UNUSED
2584 REAL8 eta)
2585{
2586 return (49.13/6.72 - 177.0/8.0 * eta) * LAL_PI;
2587}
2588
2589static REAL8 UNUSED
2591 REAL8 eta)
2592{
2593 return (379.99588601/2.79417600 + 16.0/3.0 * LAL_PI*LAL_PI - 17.12/1.05 * LAL_GAMMA
2594 + (36.9/3.2 * LAL_PI*LAL_PI - 2486.1497/7.2576) * eta
2595 + 48.8849/1.6128 * eta*eta - 8.5/6.4 * eta*eta*eta);
2596}
2597
2598static REAL8 UNUSED
2600 REAL8 UNUSED eta)
2601{
2602 return -8.56/1.05;
2603}
2604
2605static REAL8 UNUSED
2607 REAL8 eta)
2608{
2609 return (129.817/2.304 - 320.7739/4.8384 * eta + 61.3373/1.2096 * eta*eta) * LAL_PI;
2610}
2611
2612
2613/**
2614 * Computes the PN Coefficients for using in the TaylorF2Ecc equation.
2615 *
2616 * 3-dimensional REAL8 array eccPNCoeffs[ORDER][v_power][v0_power] are calculated,
2617 * where ORDER is relative PN order, v_power is power of v, and v0_power is power of v0.
2618 * Note that ORDER = v_power + v0_power.
2619 *
2620 * Terms given in equation 6.26 of: Blake Moore, Marc Favata,
2621 * K.G.Arun, and Chandra Kant Mishra, "Gravitational-wave phasing
2622 * for low-eccentricity inspiralling compact binaries to 3PN order",
2623 * Phys. Rev. D 93, 124061 (2016), arXiv:1605.00304
2624 */
2625
2626static INT4 UNUSED
2628{
2629 INT4 ret = 0;
2630 memset(eccPNCoeffs, 0x00, (LAL_MAX_ECC_PN_ORDER+1)*(LAL_MAX_ECC_PN_ORDER+1)*(LAL_MAX_ECC_PN_ORDER+1)*sizeof(REAL8));
2631 eccPNCoeffs[0][0][0] = 1.0; // lowest order constant term
2632
2633 eccPNCoeffs[2][2][0] = 29.9076223/8.1976608 + 18.766963/2.927736*eta; //v^2 term
2634 eccPNCoeffs[2][1][1] = 0.0; //v*v0 term
2635 eccPNCoeffs[2][0][2] = 2.833/1.008 - 19.7/3.6*eta; //v0^2 term
2636
2637 eccPNCoeffs[3][3][0] = -28.19123/2.82600*LAL_PI; //v^3 term
2638 eccPNCoeffs[3][0][3] = 37.7/7.2*LAL_PI; //v0^3 term
2639
2640 eccPNCoeffs[4][4][0] = 16.237683263/3.330429696 + 241.33060753/9.71375328*eta+156.2608261/6.9383952*eta*eta; //v^4 term
2641 eccPNCoeffs[4][2][2] = 84.7282939759/8.2632420864-7.18901219/3.68894736*eta-36.97091711/1.05398496*eta*eta; //v^2*v0^2 term
2642 eccPNCoeffs[4][0][4] = -1.193251/3.048192 - 66.317/9.072*eta +18.155/1.296*eta*eta; //v0^4 term
2643
2644 eccPNCoeffs[5][5][0] = -28.31492681/1.18395270*LAL_PI - 115.52066831/2.70617760*LAL_PI*eta; //v^5 term
2645 eccPNCoeffs[5][3][2] = -79.86575459/2.84860800*LAL_PI + 55.5367231/1.0173600*LAL_PI*eta; //v^3*v0^2 term
2646 eccPNCoeffs[5][2][3] = 112.751736071/5.902315776*LAL_PI + 70.75145051/2.10796992*LAL_PI*eta; //v^2*v0^3 term
2647 eccPNCoeffs[5][0][5] = 76.4881/9.0720*LAL_PI - 94.9457/2.2680*LAL_PI*eta; //v0^5 term
2648
2649 eccPNCoeffs[6][6][0] = -436.03153867072577087/1.32658535116800000 + 53.6803271/1.9782000*LAL_GAMMA + 157.22503703/3.25555200*LAL_PI*LAL_PI
2650 +(2991.72861614477/6.89135247360 - 15.075413/1.446912*LAL_PI*LAL_PI)*eta
2651 +345.5209264991/4.1019955200*eta*eta + 506.12671711/8.78999040*eta*eta*eta
2652 + 384.3505163/5.9346000*log(2.0) - 112.1397129/1.7584000*log(3.0); //v^6 term except log(16*v^2) term
2653 eccPNCoeffs[6][4][2] = 46.001356684079/3.357073133568 + 253.471410141755/5.874877983744*eta
2654 - 169.3852244423/2.3313007872*eta*eta - 307.833827417/2.497822272*eta*eta*eta; //v^4*v0^2 term
2655 eccPNCoeffs[6][3][3] = -106.2809371/2.0347200*LAL_PI*LAL_PI; //v^3*v0^3 term
2656 eccPNCoeffs[6][2][4] = -3.56873002170973/2.49880440692736 - 260.399751935005/8.924301453312*eta
2657 + 15.0484695827/3.5413894656*eta*eta + 340.714213265/3.794345856*eta*eta*eta; //v^2*v0^4 term
2658 eccPNCoeffs[6][0][6] = 265.31900578691/1.68991764480 - 33.17/1.26*LAL_GAMMA + 12.2833/1.0368*LAL_PI*LAL_PI
2659 + (91.55185261/5.48674560 - 3.977/1.152*LAL_PI*LAL_PI)*eta - 5.732473/1.306368*eta*eta
2660 - 30.90307/1.39968*eta*eta*eta + 87.419/1.890*log(2.0) - 260.01/5.60*log(3.0); //v0^6 term except log(16*v0^2) term
2661 //printPNCoeffs_F2(eccPNCoeffs);
2662 return ret;
2663}
2664
2665/**
2666 * Compute eccentric phase correction term using eccPNCeoffs[k][i][j]
2667 *
2668 */
2669
2670static REAL8 UNUSED
2672{
2673 static REAL8 v0_power[LAL_MAX_ECC_PN_ORDER+1];
2674 /* following code is not efficient in memory usage, need to be improved later */
2675 static REAL8 eccPNCoeffs[LAL_MAX_ECC_PN_ORDER+1][LAL_MAX_ECC_PN_ORDER+1][LAL_MAX_ECC_PN_ORDER+1]; // we want to calculate just one time
2676 REAL8 v_power[LAL_MAX_ECC_PN_ORDER+1];
2677 REAL8 phasing = 0.0;
2678 REAL8 global_factor;
2679 v0_power[0] = 1.0;
2680 for(int i=1; i<=LAL_MAX_ECC_PN_ORDER; i++)
2681 {
2682 v0_power[i] = v0_power[i-1]*v0;
2683 }
2684 eccentricityPNCoeffs_F2(eta, eccPNCoeffs);
2685 //printPNCoeffs_F2(eccPNCoeffs);
2686 v_power[0] = 1.0;
2687 for(int i=1; i<=LAL_MAX_ECC_PN_ORDER; i++)
2688 {
2689 v_power[i] = v_power[i-1]*v;
2690 }
2691
2692 global_factor = -2.355/1.462*ecc*ecc*pow(v0/v, 19.0/3.0);
2693 global_factor *= (3.0/128.0/eta); // overall factor except v^-5 in phase term, this is Newtonian phase term
2694 if(ecc_order == -1) {
2695 ecc_order = LAL_MAX_ECC_PN_ORDER;
2696 }
2697 if(ecc_order > LAL_MAX_ECC_PN_ORDER) {
2698 return XLAL_REAL8_FAIL_NAN;
2699 }
2700
2701 REAL8 phaseOrder = 0;
2702 for(int i=0; i<=ecc_order; i++)
2703 {
2704 phaseOrder = 0;
2705 INT4 k = 0;
2706 for(int j=i; j>=0; j--)
2707 {
2708 k = i - j;
2709 if( j==6 )
2710 {
2711 phaseOrder += (eccPNCoeffs[i][j][k]+53.6803271/3.9564000*log(16.0*v_power[2]))*v_power[j]*v0_power[k];
2712 //phasing += (eccPNCoeffs[i][j][k]+53.6803271/3.9564000*log(16.0*v_power[2]))*v_power[j]*v0_power[k];
2713 }
2714 else if( k == 6 )
2715 {
2716 phaseOrder += (eccPNCoeffs[i][j][k] - 33.17/2.52*log(16.0*v0_power[2]))*v_power[j]*v0_power[k];
2717 //phasing += (eccPNCoeffs[i][j][k] - 33.17/2.52*log(16.0*v0_power[2]))*v_power[j]*v0_power[k];
2718 }
2719 else
2720 {
2721 phaseOrder += eccPNCoeffs[i][j][k]*v_power[j]*v0_power[k];
2722 //phasing += eccPNCoeffs[i][j][k]*v_power[j]*v0_power[k];
2723 }
2724 }
2725 phasing += phaseOrder;
2726 //ecc_phase_order[i] = phaseOrder*global_factor;
2727 }
2728 //fprintf(stdout, "======== DEBUG for eccentricity ================\n");
2729 //fprintf(stdout, "eccentricityPhasing_F2 phasing = %g, global_factor = %g, ecc_order = %d, ecc = %g\n", phasing, global_factor, ecc_order, ecc);
2730 return phasing*global_factor;
2731}
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_12PNTidalCoeff(REAL8 eta, REAL8 chi, REAL8 lambda)
static const REAL8 UNUSED XLALSimInspiralSpinDot_4PNS2OCoeffAvg
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_12PNTidalCoeff(REAL8 eta, REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_5PNSicoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorEtVOfZeta_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtVOfZeta_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_12PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_3PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_12PNTidalCoeff(REAL8 eta, REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralPNEnergy_5PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorEtPhasing_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_5PNLogCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_12PNTidalCoeff(REAL8 eta, REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_0PNCoeff(REAL8 eta)
Computes the flux PN Coefficients.
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_10PNTidalCoeff(REAL8 mByM)
static INT4 UNUSED eccentricityPNCoeffs_F2(REAL8 eta, REAL8 eccPNCoeffs[LAL_MAX_ECC_PN_ORDER+1][LAL_MAX_ECC_PN_ORDER+1][LAL_MAX_ECC_PN_ORDER+1])
Computes the PN Coefficients for using in the TaylorF2Ecc equation.
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_12PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_5PNCoeff(REAL8 eta)
static const REAL8 UNUSED XLALSimInspiralSpinDot_4PNS2CoeffAvg
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_6PNS1S2OCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNSelf2SCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_14PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_0PNCoeff(REAL8 eta)
Computes the PN Coefficients for using in the TaylorT4 frequency equation.
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_0PNCoeff(REAL8 eta)
Computes the PN Coefficients for using in the TaylorT2 phasing equation.
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS2vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_7PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_4PNQMSOCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS2Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNQM2SOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_3PNSicoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorEtPhasing_0PNCoeff(REAL8 m)
Computes the PN Coefficients for using in the TaylorEt dPhase/dt equation.
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_12PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_4PN(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralL_5PNSincoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1OS2OCoeff(REAL8 eta)
static const REAL8 UNUSED XLALSimInspiralSpinDot_4PNS2Coeff
static REAL8 UNUSED XLALSimInspiralPNEnergy_8PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS2OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNSelf2SOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_6PNCoeff(REAL8 eta)
static REAL8 UNUSED eccentricityPhasing_F2(REAL8 v, REAL8 v0, REAL8 ecc, REAL8 eta, INT4 ecc_order)
Compute eccentric phase correction term using eccPNCeoffs[k][i][j].
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_7PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_0PNCoeff(REAL8 eta)
Computes the PN Coefficients for using in the PN energy equation.
static REAL8 UNUSED XLALSimInspiralPNFlux_10PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1OS1OCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralLDot_5PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNQMSOCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralLN(REAL8 M, REAL8 eta, REAL8 v)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_0PNCoeff(REAL8 eta)
Computes the PN Coefficients for using in the TaylorT3 phasing equation.
static REAL8 UNUSED XLALSimInspiralPNFlux_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralLDot_7PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_6PN(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_0PNCoeff(REAL8 m, REAL8 eta)
Computes the PN Coefficients for using in the TaylorEt dZeta/dt equation.
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_6PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_5PNSiLcoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_7PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_3PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_5PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_3PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralLDot_4PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_10PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_13PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralL_3PNSincoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_12PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1OS1OCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS2nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1OS2OCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralL_3PNSiLcoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_7PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_3PNSicoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_7PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNS1nS2nCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralSpinDot_3PNCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_8PNLogCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtPhasing_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNS1S2OCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1OS2OCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralSpinDot_7PNCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_15PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtZeta_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1nS2nCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorEtPhasing_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNQMS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_0PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_0PNCoeff(REAL8 totalmass)
Computes the PN Coefficients for using in the TaylorT3 frequency equation.
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_5PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNQMS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_6PNQM2SCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Phasing_10PNTidalCoeff(REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralSpinDot_5PNCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNQMS1S1Coeff(REAL8 mByM)
static void UNUSED XLALSimInspiralPNPhasing_F2(PNPhasingSeries *pfa, const REAL8 m1, const REAL8 m2, const REAL8 chi1L, const REAL8 chi2L, const REAL8 chi1sq, const REAL8 chi2sq, const REAL8 chi1dotchi2, LALDict *p)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1OS1OCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_3PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_4PNQMSnCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_10PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_3PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_10PNTidalCoeff(REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralL_2PN(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNQMSvCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_3PNCoeff(REAL8 UNUSED eta)
static const REAL8 UNUSED XLALSimInspiralSpinDot_4PNS2nCoeff
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_6PNSelf2SCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_12PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_0PNCoeff(REAL8 totalmass, REAL8 eta)
Computes the PN Coefficients for using in the TaylorT2 timing equation.
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNQMS1S1Coeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralL_3PNSiLcoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralLDot_4PNQMSSCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNQMS1OS1OCoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_4PNQM2SCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorEtVOfZeta_2PNCoeff(REAL8 eta)
Computes the PN Coefficients for using in the TaylorEt v(zeta) equation, which is the square root of ...
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1OS2OCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_8PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNS1OS1OCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNQMS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_5PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_5PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1nS1nCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_10PNTidalCoeff(REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_6PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNQMSnCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNS1S2CoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT3Frequency_10PNTidalCoeff(REAL8 chi, REAL8 lambda)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1vS2vCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1S2Coeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralLDot_3PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2dtdv_4PNS1vS1vCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralSpinDot_6PNS2CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralPNFlux_8PNSOCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT2Timing_6PNLogCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_4PNS1S1CoeffAvg(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorF2Phasing_10PNTidalCoeff(REAL8 mByM)
static REAL8 UNUSED XLALSimInspiralTaylorT4wdot_6PNS1OS2OCoeffAvg(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralTaylorT2Phasing_7PNCoeff(REAL8 eta)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi5(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi6(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi6L(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupdQuadMon1(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi0(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupTidalLambda2(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi1(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupTidalLambda1(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi4(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi3(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi7(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi5L(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupdQuadMon2(LALDict *params)
INT4 XLALSimInspiralWaveformParamsLookupPNTidalOrder(LALDict *params)
INT4 XLALSimInspiralWaveformParamsLookupPNSpinOrder(LALDict *params)
REAL8 XLALSimInspiralWaveformParamsLookupNonGRDChi2(LALDict *params)
REAL8 M
Definition: bh_qnmode.c:133
double i
Definition: bh_ringdown.c:118
#define __attribute__(x)
#define LAL_C_SI
#define LAL_PI
#define LAL_GAMMA
#define LAL_G_SI
double REAL8
int32_t INT4
#define LAL_SIM_INSPIRAL_TIDAL_ORDER_DEFAULT
#define LAL_MAX_ECC_PN_ORDER
#define PN_PHASING_SERIES_MAX_ORDER
Structure for passing around PN phasing coefficients.
@ LAL_SIM_INSPIRAL_SPIN_ORDER_0PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_25PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_35PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_ALL
@ LAL_SIM_INSPIRAL_SPIN_ORDER_2PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_15PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_1PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_05PN
@ LAL_SIM_INSPIRAL_SPIN_ORDER_3PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_5PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_6PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_75PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_7PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_0PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_65PN
static const INT4 m
#define XLAL_ERROR_VOID(...)
#define XLAL_REAL8_FAIL_NAN
int XLALPrintError(const char *fmt,...) _LAL_GCC_PRINTF_FORMAT_(1
XLAL_EINVAL
p
REAL8 vlogv[PN_PHASING_SERIES_MAX_ORDER+1]
REAL8 vlogvsq[PN_PHASING_SERIES_MAX_ORDER+1]
REAL8 v[PN_PHASING_SERIES_MAX_ORDER+1]