/
YB-slip.cpp
9156 lines (9135 loc) · 783 KB
/
YB-slip.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
/* twobody solutions for slip particles
* Copyright (C) 2007 Kengo Ichiki <kichiki@users.sourceforge.net>
* $Id: YB-slip.c,v 1.4 2008/10/22 06:13:09 kichiki Exp $
*
* This program is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License
* as published by the Free Software Foundation; either version 2
* of the License, or (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
*/
#include <cstdio> // fprintf()
#include <cstdlib> // malloc(), free()
#include "memory-check.h" // CHECK_MALLOC
static void twobody_YB_slip_15 (int n, double l, double *f, double *fk);
static void twobody_YB_slip_16 (int n, double l, double *f, double *fk);
static void twobody_YB_slip_17 (int n, double l, double *f, double *fk);
static void twobody_YB_slip_18 (int n, double l, double *f, double *fk);
//static void twobody_YB_slip_19 (int n, double l, double *f, double *fk);
//static void twobody_YB_slip_20 (int n, double l, double *f, double *fk);
double SL_G1 (int m, int n);
double SL_G2 (int m, int n);
/* YB1 */
void twobody_YB_slip (int n, double l, double *f)
{
int i;
for (i = 0; i <= n; i ++) f[i] = 0.0;
int q;
double *fk = (double *)malloc (sizeof (double) * (n+1));
CHECK_MALLOC (fk, "twobody_YB_slip");
/* k = 2 */
for (q = 0; q <= 2; q++) fk[q] = 0.0;
fk[1] =
- 6.0 * SL_G1(0,3) * SL_G2(2,3)
;
f[2] = 0.0;
for (q = 0; q <= 2; q++) f[2] = f[2] * l + fk[2-q];
if (n == 2)
{
free (fk);
return;
}
/* k = 3 */
for (q = 0; q <= 3; q++) fk[q] = 0.0;
fk[1] =
- 9.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3)
;
f[3] = 0.0;
for (q = 0; q <= 3; q++) f[3] = f[3] * l + fk[3-q];
if (n == 3)
{
free (fk);
return;
}
/* k = 4 */
for (q = 0; q <= 4; q++) fk[q] = 0.0;
fk[2] =
- 27.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) / 2.0
;
f[4] = 0.0;
for (q = 0; q <= 4; q++) f[4] = f[4] * l + fk[4-q];
if (n == 4)
{
free (fk);
return;
}
/* k = 5 */
for (q = 0; q <= 5; q++) fk[q] = 0.0;
fk[1] =
- 12.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3)
;
fk[2] =
- 81.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) / 4.0
;
fk[3] =
- 36.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3)
;
f[5] = 0.0;
for (q = 0; q <= 5; q++) f[5] = f[5] * l + fk[5-q];
if (n == 5)
{
free (fk);
return;
}
/* k = 6 */
for (q = 0; q <= 6; q++) fk[q] = 0.0;
fk[2] =
- 108.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3)
;
fk[3] =
- 243.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) / 8.0
;
fk[4] =
- 72.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,3)
;
f[6] = 0.0;
for (q = 0; q <= 6; q++) f[6] = f[6] * l + fk[6-q];
if (n == 6)
{
free (fk);
return;
}
/* k = 7 */
for (q = 0; q <= 7; q++) fk[q] = 0.0;
fk[2] =
- 189.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
;
fk[3] =
- 3.0 * SL_G1(0,3) * ( 2560.0 * SL_G1(0,3) * SL_G2(2,5)
+ 243.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) ) / 16.0
;
fk[4] =
- 243.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
;
fk[5] =
48.0 * SL_G1(0,3) * SL_G1(2,3) * ( 7.0 * SL_G2(2,7)
- 6.0 * SL_G2(0,5)
- 4.0 * SL_G2(-1,4) )
;
f[7] = 0.0;
for (q = 0; q <= 7; q++) f[7] = f[7] * l + fk[7-q];
if (n == 7)
{
free (fk);
return;
}
/* k = 8 */
for (q = 0; q <= 8; q++) fk[q] = 0.0;
fk[2] =
- 24.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 21.0 * SL_G1(2,7)
+ 5.0 * SL_G1(0,2) * SL_G1(0,3)
+ 60.0 * SL_G1(-1,4)
+ 4.0 * SL_G1(-3,2) ) / 5.0
;
fk[3] =
- 486.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
;
fk[4] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * ( 23040.0 * SL_G1(0,3) * SL_G2(2,5)
+ 729.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 2048.0 * SL_G1(0,3) * SL_G2(0,3) ) / 32.0
;
fk[5] =
- 405.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
;
fk[6] =
72.0 * SL_G1(0,3) * SL_G1(2,3) * ( 7.0 * SL_G2(2,3) * SL_G2(2,7)
- 6.0 * SL_G2(0,5) * SL_G2(2,3)
- 4.0 * SL_G2(-1,4) * SL_G2(2,3)
- SL_G2(0,3) * SL_G2(0,3) )
;
f[8] = 0.0;
for (q = 0; q <= 8; q++) f[8] = f[8] * l + fk[8-q];
if (n == 8)
{
free (fk);
return;
}
/* k = 9 */
for (q = 0; q <= 9; q++) fk[q] = 0.0;
fk[2] =
- 36.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 21.0 * SL_G1(2,3) * SL_G1(2,7)
+ 60.0 * SL_G1(-1,4) * SL_G1(2,3)
+ 35.0 * SL_G1(0,3) * SL_G1(0,3)
+ 4.0 * SL_G1(-3,3) ) / 5.0
;
fk[3] =
- 3159.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) / 4.0
;
fk[4] =
- 27.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * ( 7680.0 * SL_G1(0,3) * SL_G2(2,5)
+ 243.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 2560.0 * SL_G1(0,3) * SL_G2(0,3) ) / 64.0
;
fk[5] =
- 15.0 * SL_G1(0,3) * ( 1792.0 * SL_G1(0,3) * SL_G2(2,7)
+ 243.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) ) / 4.0
;
fk[6] =
36.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * ( 84.0 * SL_G2(2,3) * SL_G2(2,7)
- 90.0 * SL_G2(0,5) * SL_G2(2,3)
- 120.0 * SL_G2(-1,4) * SL_G2(2,3)
- 95.0 * SL_G2(0,3) * SL_G2(0,3)
- 4.0 * SL_G2(-3,3) ) / 5.0
;
fk[7] =
576.0 * SL_G1(0,3) * SL_G1(2,3) * ( 6.0 * SL_G2(2,9)
- 5.0 * SL_G2(0,7)
- 2.0 * SL_G2(-2,5) )
;
f[9] = 0.0;
for (q = 0; q <= 9; q++) f[9] = f[9] * l + fk[9-q];
if (n == 9)
{
free (fk);
return;
}
/* k = 10 */
for (q = 0; q <= 10; q++) fk[q] = 0.0;
fk[2] =
- 288.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 180.0 * SL_G1(2,9)
- 105.0 * SL_G1(0,2) * SL_G1(2,7)
- 105.0 * SL_G1(0,7)
+ 84.0 * SL_G1(0,2) * SL_G1(0,5)
+ 140.0 * SL_G1(-2,5)
+ 16.0 * SL_G1(-5,2) ) / 35.0
;
fk[3] =
- 108.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 21.0 * SL_G1(2,3) * SL_G1(2,7)
+ 60.0 * SL_G1(-1,4) * SL_G1(2,3)
+ 95.0 * SL_G1(0,3) * SL_G1(0,3)
+ 4.0 * SL_G1(-3,3) ) / 5.0
;
fk[4] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * ( 14336.0 * SL_G2(2,5) * SL_G1(2,7)
- 7168.0 * SL_G2(0,3) * SL_G1(2,7)
- 7680.0 * SL_G1(0,5) * SL_G2(2,5)
+ 1280.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G2(2,5)
+ 10240.0 * SL_G1(-1,4) * SL_G2(2,5)
+ 1024.0 * SL_G1(-3,2) * SL_G2(2,5)
+ 2187.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 4608.0 * SL_G2(0,3) * SL_G1(0,5)
+ 1024.0 * SL_G1(-1,4) * SL_G2(0,3) ) / 4.0
;
fk[5] =
- 27.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 23040.0 * SL_G1(0,3) * SL_G2(2,5)
+ 729.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 17920.0 * SL_G1(0,3) * SL_G2(0,3) ) / 128.0
;
fk[6] =
- 6.0 * SL_G1(0,3) * ( 3808.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
- 1600.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
+ 160.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 640.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 243.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 96.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 320.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3) )
;
fk[7] =
108.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * ( 42.0 * SL_G2(2,3) * SL_G2(2,7)
- 45.0 * SL_G2(0,5) * SL_G2(2,3)
- 60.0 * SL_G2(-1,4) * SL_G2(2,3)
- 70.0 * SL_G2(0,3) * SL_G2(0,3)
- 2.0 * SL_G2(-3,3) ) / 5.0
;
fk[8] =
96.0 * SL_G1(0,3) * SL_G1(2,3) * ( 54.0 * SL_G2(2,3) * SL_G2(2,9)
+ 7.0 * SL_G2(0,3) * SL_G2(2,7)
- 45.0 * SL_G2(0,7) * SL_G2(2,3)
- 18.0 * SL_G2(-2,5) * SL_G2(2,3)
- 6.0 * SL_G2(0,3) * SL_G2(0,5)
- 4.0 * SL_G2(-1,4) * SL_G2(0,3) )
;
f[10] = 0.0;
for (q = 0; q <= 10; q++) f[10] = f[10] * l + fk[10-q];
if (n == 10)
{
free (fk);
return;
}
/* k = 11 */
for (q = 0; q <= 11; q++) fk[q] = 0.0;
fk[2] =
- 48.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 1620.0 * SL_G1(2,3) * SL_G1(2,9)
- 798.0 * SL_G1(0,3) * SL_G1(2,7)
- 945.0 * SL_G1(0,7) * SL_G1(2,3)
+ 1260.0 * SL_G1(-2,5) * SL_G1(2,3)
+ 756.0 * SL_G1(0,3) * SL_G1(0,5)
+ 35.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,3)
+ 420.0 * SL_G1(-1,4) * SL_G1(0,3)
+ 28.0 * SL_G1(-3,2) * SL_G1(0,3)
+ 144.0 * SL_G1(-5,3) ) / 35.0
;
fk[3] =
- 162.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 21.0 * SL_G1(2,3) * SL_G1(2,7)
+ 60.0 * SL_G1(-1,4) * SL_G1(2,3)
+ 125.0 * SL_G1(0,3) * SL_G1(0,3)
+ 4.0 * SL_G1(-3,3) ) / 5.0
;
fk[4] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * ( 430080.0 * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,7)
- 317184.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,7)
- 230400.0 * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,5)
+ 307200.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G2(2,5)
+ 384000.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5)
+ 30720.0 * SL_G1(-3,3) * SL_G2(2,5)
+ 138240.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 69255.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 92160.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3)
+ 107520.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G1(2,3)
+ 11520.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3)
+ 1024.0 * SL_G1(-3,3) * SL_G2(0,3) ) / 80.0
;
fk[5] =
- 81.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 23040.0 * SL_G1(0,3) * SL_G2(2,5)
+ 729.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 32768.0 * SL_G1(0,3) * SL_G2(0,3) ) / 256.0
;
fk[6] =
- 3.0 * SL_G1(0,3) * SL_G1(2,3) * ( 623616.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
- 1152000.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
+ 506880.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 460800.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 76545.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 184320.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 168960.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3)
+ 24320.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3)
+ 1024.0 * SL_G2(-3,3) * SL_G1(0,3) ) / 80.0
;
fk[7] =
- 9.0 * SL_G1(0,3) * ( 35840.0 * SL_G1(0,3) * SL_G2(2,9)
- 567.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 810.0 * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1620.0 * SL_G2(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 3330.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
+ 72.0 * SL_G2(-3,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) ) / 5.0
;
fk[8] =
144.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * ( 1350.0 * SL_G2(2,3) * SL_G2(2,9)
+ 1638.0 * SL_G2(0,3) * SL_G2(2,7)
- 1260.0 * SL_G2(0,7) * SL_G2(2,3)
- 1050.0 * SL_G2(-2,5) * SL_G2(2,3)
- 1512.0 * SL_G2(0,3) * SL_G2(0,5)
- 35.0 * SL_G2(0,2) * SL_G2(0,3) * SL_G2(0,3)
- 1260.0 * SL_G2(-1,4) * SL_G2(0,3)
- 28.0 * SL_G2(-3,2) * SL_G2(0,3)
- 48.0 * SL_G2(-5,3) ) / 35.0
;
fk[9] =
768.0 * SL_G1(0,3) * SL_G1(2,3) * ( 33.0 * SL_G2(2,11)
- 28.0 * SL_G2(0,9)
- 8.0 * SL_G2(-3,6) )
;
f[11] = 0.0;
for (q = 0; q <= 11; q++) f[11] = f[11] * l + fk[11-q];
if (n == 11)
{
free (fk);
return;
}
/* k = 12 */
for (q = 0; q <= 12; q++) fk[q] = 0.0;
fk[2] =
- 1152.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 77.0 * SL_G1(2,11)
- 56.0 * SL_G1(0,2) * SL_G1(2,9)
- 56.0 * SL_G1(0,9)
+ 45.0 * SL_G1(0,2) * SL_G1(0,7)
+ 28.0 * SL_G1(-3,6)
+ 4.0 * SL_G1(-7,2) ) / 7.0
;
fk[3] =
- 432.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 540.0 * SL_G1(2,3) * SL_G1(2,9)
- 21.0 * SL_G1(0,3) * SL_G1(2,7)
- 315.0 * SL_G1(0,7) * SL_G1(2,3)
+ 420.0 * SL_G1(-2,5) * SL_G1(2,3)
+ 252.0 * SL_G1(0,3) * SL_G1(0,5)
+ 70.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,3)
+ 840.0 * SL_G1(-1,4) * SL_G1(0,3)
+ 56.0 * SL_G1(-3,2) * SL_G1(0,3)
+ 48.0 * SL_G1(-5,3) ) / 35.0
;
fk[4] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * ( 3456000.0 * SL_G2(2,5) * SL_G1(2,9)
- 1612800.0 * SL_G2(0,3) * SL_G1(2,9)
- 1344000.0 * SL_G1(0,2) * SL_G2(2,5) * SL_G1(2,7)
+ 35721.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,7)
+ 403200.0 * SL_G1(0,2) * SL_G2(0,3) * SL_G1(2,7)
- 2150400.0 * SL_G1(0,7) * SL_G2(2,5)
+ 860160.0 * SL_G1(0,2) * SL_G1(0,5) * SL_G2(2,5)
+ 896000.0 * SL_G1(-2,5) * SL_G2(2,5)
+ 163840.0 * SL_G1(-5,2) * SL_G2(2,5)
+ 102060.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 314685.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 6804.0 * SL_G1(-3,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 672000.0 * SL_G2(0,3) * SL_G1(0,7)
+ 107520.0 * SL_G1(-2,5) * SL_G2(0,3) ) / 70.0
;
fk[5] =
- 27.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 17920.0 * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,7)
- 23296.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,7)
- 9600.0 * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,5)
+ 12800.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G2(2,5)
+ 59200.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5)
+ 1280.0 * SL_G1(-3,3) * SL_G2(2,5)
+ 7680.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 3645.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 9600.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3)
+ 10240.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G1(2,3)
+ 5440.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3)
+ 256.0 * SL_G1(-3,3) * SL_G2(0,3) ) / 20.0
;
fk[6] =
- SL_G1(0,3) * ( 97943552.0 * SL_G2(2,3) * SL_G1(2,7) * SL_G2(2,7)
- 44040192.0 * SL_G1(0,5) * SL_G2(2,3) * SL_G2(2,7)
+ 6881280.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
+ 45875200.0 * SL_G1(-1,4) * SL_G2(2,3) * SL_G2(2,7)
+ 5505024.0 * SL_G1(-3,2) * SL_G2(2,3) * SL_G2(2,7)
- 110100480.0 * SL_G2(0,3) * SL_G2(2,5) * SL_G1(2,7)
- 41287680.0 * SL_G2(0,5) * SL_G2(2,3) * SL_G1(2,7)
- 9175040.0 * SL_G2(-1,4) * SL_G2(2,3) * SL_G1(2,7)
+ 34406400.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,7)
+ 5598720.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
+ 31457280.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G2(2,5)
+ 177147.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 12690432.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 18874368.0 * SL_G1(0,5) * SL_G2(0,5) * SL_G2(2,3)
+ 6291456.0 * SL_G2(-1,4) * SL_G1(0,5) * SL_G2(2,3)
+ 8388608.0 * SL_G1(-1,4) * SL_G2(-1,4) * SL_G2(2,3) ) / 512.0
;
fk[7] =
- 9.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * ( 193536.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
- 576000.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
+ 357120.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 345600.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 40095.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 149760.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 161280.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3)
+ 66560.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3)
+ 1024.0 * SL_G2(-3,3) * SL_G1(0,3) ) / 80.0
;
fk[8] =
- SL_G1(0,3) * ( 1866240.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,9)
- 1433600.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,7)
- 15309.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 134400.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,7)
+ 614400.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G1(2,5)
+ 204800.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G1(2,5)
+ 21870.0 * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 43740.0 * SL_G2(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 108135.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1944.0 * SL_G2(-3,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 57600.0 * SL_G1(0,3) * SL_G2(0,7) * SL_G2(2,3)
+ 161280.0 * SL_G2(-2,5) * SL_G1(0,3) * SL_G2(2,3) ) / 10.0
;
fk[9] =
144.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * ( 2025.0 * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,9)
+ 2751.0 * SL_G2(0,3) * SL_G2(2,3) * SL_G2(2,7)
- 1890.0 * SL_G2(0,7) * SL_G2(2,3) * SL_G2(2,3)
- 1575.0 * SL_G2(-2,5) * SL_G2(2,3) * SL_G2(2,3)
- 2583.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G2(2,3)
- 2310.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G2(2,3)
- 72.0 * SL_G2(-5,3) * SL_G2(2,3)
- 385.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,3)
- 56.0 * SL_G2(-3,3) * SL_G2(0,3) ) / 35.0
;
fk[10] =
1152.0 * SL_G1(0,3) * SL_G1(2,3) * ( 33.0 * SL_G2(2,3) * SL_G2(2,11)
+ 6.0 * SL_G2(0,3) * SL_G2(2,9)
- 28.0 * SL_G2(0,9) * SL_G2(2,3)
- 8.0 * SL_G2(-3,6) * SL_G2(2,3)
- 5.0 * SL_G2(0,3) * SL_G2(0,7)
- 2.0 * SL_G2(-2,5) * SL_G2(0,3) )
;
f[12] = 0.0;
for (q = 0; q <= 12; q++) f[12] = f[12] * l + fk[12-q];
if (n == 12)
{
free (fk);
return;
}
/* k = 13 */
for (q = 0; q <= 13; q++) fk[q] = 0.0;
fk[2] =
- 576.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 1155.0 * SL_G1(2,3) * SL_G1(2,11)
- 660.0 * SL_G1(0,3) * SL_G1(2,9)
- 105.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(2,7)
- 840.0 * SL_G1(0,9) * SL_G1(2,3)
+ 420.0 * SL_G1(-3,6) * SL_G1(2,3)
+ 570.0 * SL_G1(0,3) * SL_G1(0,7)
+ 84.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,5)
+ 140.0 * SL_G1(-2,5) * SL_G1(0,3)
+ 16.0 * SL_G1(-5,2) * SL_G1(0,3)
+ 60.0 * SL_G1(-7,3) ) / 35.0
;
fk[3] =
- 216.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 1620.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,9)
+ 84.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,7)
- 945.0 * SL_G1(0,7) * SL_G1(2,3) * SL_G1(2,3)
+ 1260.0 * SL_G1(-2,5) * SL_G1(2,3) * SL_G1(2,3)
+ 756.0 * SL_G1(0,3) * SL_G1(0,5) * SL_G1(2,3)
+ 2940.0 * SL_G1(-1,4) * SL_G1(0,3) * SL_G1(2,3)
+ 144.0 * SL_G1(-5,3) * SL_G1(2,3)
+ 875.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(0,3)
+ 196.0 * SL_G1(-3,3) * SL_G1(0,3) ) / 35.0
;
fk[4] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * ( 10368000.0 * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,9)
- 6635520.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,9)
- 1021440.0 * SL_G1(0,3) * SL_G2(2,5) * SL_G1(2,7)
+ 107163.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,7)
+ 2240000.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,7)
- 6451200.0 * SL_G1(0,7) * SL_G1(2,3) * SL_G2(2,5)
+ 2688000.0 * SL_G1(-2,5) * SL_G1(2,3) * SL_G2(2,5)
+ 967680.0 * SL_G1(0,3) * SL_G1(0,5) * SL_G2(2,5)
+ 268800.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5)
+ 2150400.0 * SL_G1(-1,4) * SL_G1(0,3) * SL_G2(2,5)
+ 215040.0 * SL_G1(-3,2) * SL_G1(0,3) * SL_G2(2,5)
+ 491520.0 * SL_G1(-5,3) * SL_G2(2,5)
+ 306180.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1097145.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 20412.0 * SL_G1(-3,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1935360.0 * SL_G2(0,3) * SL_G1(0,7) * SL_G1(2,3)
+ 1075200.0 * SL_G1(-2,5) * SL_G2(0,3) * SL_G1(2,3)
+ 387072.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(0,5)
+ 71680.0 * SL_G1(-1,4) * SL_G1(0,3) * SL_G2(0,3)
+ 12288.0 * SL_G1(-5,3) * SL_G2(0,3) ) / 140.0
;
fk[5] =
- 27.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 430080.0 * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,7)
- 634368.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,7)
- 230400.0 * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,5)
+ 307200.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G2(2,5)
+ 1766400.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5)
+ 30720.0 * SL_G1(-3,3) * SL_G2(2,5)
+ 276480.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 91125.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 230400.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3)
+ 399360.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G1(2,3)
+ 366080.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3)
+ 12288.0 * SL_G1(-3,3) * SL_G2(0,3) ) / 320.0
;
fk[6] =
- 3.0 * SL_G1(0,3) * ( 391774208.0 * SL_G1(2,3) * SL_G2(2,3) * SL_G1(2,7) * SL_G2(2,7)
+ 44040192.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,7)
- 220200960.0 * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 183500800.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 339476480.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
+ 22020096.0 * SL_G1(-3,3) * SL_G2(2,3) * SL_G2(2,7)
- 1247805440.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,7)
- 162398208.0 * SL_G2(0,5) * SL_G1(2,3) * SL_G2(2,3) * SL_G1(2,7)
- 91750400.0 * SL_G2(-1,4) * SL_G1(2,3) * SL_G2(2,3) * SL_G1(2,7)
+ 383057920.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,7)
- 5505024.0 * SL_G2(-3,3) * SL_G1(2,3) * SL_G1(2,7)
+ 327680000.0 * SL_G1(0,3) * SL_G1(2,5) * SL_G2(2,5) * SL_G2(2,5)
+ 196608000.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5) * SL_G2(2,5)
+ 308019200.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
+ 27993600.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
- 19660800.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,5)
+ 131072000.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,5)
+ 43253760.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 2621440.0 * SL_G1(-3,3) * SL_G2(0,3) * SL_G2(2,5)
- 23592960.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3) * SL_G1(2,5)
+ 31457280.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3) * SL_G1(2,5)
+ 17039360.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 3145728.0 * SL_G2(-3,3) * SL_G1(0,3) * SL_G1(2,5)
+ 885735.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 92689920.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 94371840.0 * SL_G1(0,5) * SL_G2(0,5) * SL_G1(2,3) * SL_G2(2,3)
+ 31457280.0 * SL_G1(-1,4) * SL_G2(0,5) * SL_G1(2,3) * SL_G2(2,3)
+ 31457280.0 * SL_G2(-1,4) * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,3)
+ 83886080.0 * SL_G1(-1,4) * SL_G2(-1,4) * SL_G1(2,3) * SL_G2(2,3)
+ 5898240.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 3145728.0 * SL_G1(-3,3) * SL_G2(0,5) * SL_G2(2,3)
+ 6553600.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3)
+ 7864320.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3)
+ 13107200.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) ) / 5120.0
;
fk[7] =
- 27.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * ( 26880.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
- 230400.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
+ 267264.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 230400.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 19683.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 92160.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 104448.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3)
+ 94720.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3)
+ 2048.0 * SL_G2(-3,3) * SL_G1(0,3) ) / 64.0
;
fk[8] =
- 3.0 * SL_G1(0,3) * SL_G1(2,3) * ( 6912000.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,9)
- 10035200.0 * SL_G1(2,5) * SL_G2(2,5) * SL_G2(2,7)
- 3386880.0 * SL_G1(0,3) * SL_G2(2,5) * SL_G2(2,7)
- 26091520.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,7)
- 71442.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 9031680.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,7)
+ 5376000.0 * SL_G2(0,5) * SL_G1(2,5) * SL_G2(2,5)
- 7168000.0 * SL_G2(-1,4) * SL_G1(2,5) * SL_G2(2,5)
+ 3225600.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,5)
+ 3225600.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,5)
+ 9856000.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G1(2,5)
+ 358400.0 * SL_G2(0,2) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 7168000.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G1(2,5)
+ 286720.0 * SL_G2(-3,2) * SL_G2(0,3) * SL_G1(2,5)
- 716800.0 * SL_G2(-3,5) * SL_G1(2,5)
+ 153090.0 * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 408240.0 * SL_G2(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1403325.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 20412.0 * SL_G2(-3,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 3548160.0 * SL_G1(0,3) * SL_G2(0,7) * SL_G2(2,3)
+ 2042880.0 * SL_G2(-2,5) * SL_G1(0,3) * SL_G2(2,3)
+ 1121792.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,5)
+ 645120.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(0,3)
+ 71680.0 * SL_G2(-3,5) * SL_G1(0,3)
+ 12288.0 * SL_G2(-5,3) * SL_G1(0,3) ) / 140.0
;
fk[9] =
- 36.0 * SL_G1(0,3) * ( 492800.0 * SL_G1(0,3) * SL_G2(2,11)
- 7290.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,9)
- 24192.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 8505.0 * SL_G2(0,7) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 13230.0 * SL_G2(-2,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 27216.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
+ 37800.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
+ 864.0 * SL_G2(-5,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
+ 13230.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3)
+ 1512.0 * SL_G2(-3,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) ) / 35.0
;
fk[10] =
192.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * ( 6930.0 * SL_G2(2,3) * SL_G2(2,11)
+ 12240.0 * SL_G2(0,3) * SL_G2(2,9)
+ 1029.0 * SL_G2(2,7) * SL_G2(2,7)
- 882.0 * SL_G2(0,5) * SL_G2(2,7)
+ 1190.0 * SL_G2(0,2) * SL_G2(0,3) * SL_G2(2,7)
+ 2352.0 * SL_G2(-1,4) * SL_G2(2,7)
+ 196.0 * SL_G2(-3,2) * SL_G2(2,7)
- 6300.0 * SL_G2(0,9) * SL_G2(2,3)
- 3780.0 * SL_G2(-3,6) * SL_G2(2,3)
- 10530.0 * SL_G2(0,3) * SL_G2(0,7)
- 966.0 * SL_G2(0,2) * SL_G2(0,3) * SL_G2(0,5)
- 2520.0 * SL_G2(-1,4) * SL_G2(0,5)
- 168.0 * SL_G2(-3,2) * SL_G2(0,5)
- 140.0 * SL_G2(-1,4) * SL_G2(0,2) * SL_G2(0,3)
- 5040.0 * SL_G2(-2,5) * SL_G2(0,3)
- 144.0 * SL_G2(-5,2) * SL_G2(0,3)
- 1680.0 * SL_G2(-1,4) * SL_G2(-1,4)
- 112.0 * SL_G2(-3,2) * SL_G2(-1,4)
- 180.0 * SL_G2(-7,3) ) / 35.0
;
fk[11] =
3072.0 * SL_G1(0,3) * SL_G1(2,3) * ( 52.0 * SL_G2(2,13)
- 45.0 * SL_G2(0,11)
- 10.0 * SL_G2(-4,7) )
;
f[13] = 0.0;
for (q = 0; q <= 13; q++) f[13] = f[13] * l + fk[13-q];
if (n == 13)
{
free (fk);
return;
}
/* k = 14 */
for (q = 0; q <= 14; q++) fk[q] = 0.0;
fk[2] =
- 512.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 1872.0 * SL_G1(2,13)
- 1485.0 * SL_G1(0,2) * SL_G1(2,11)
- 1485.0 * SL_G1(0,11)
+ 1232.0 * SL_G1(0,2) * SL_G1(0,9)
+ 396.0 * SL_G1(-4,7)
+ 64.0 * SL_G1(-9,2) ) / 11.0
;
fk[3] =
- 96.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 103950.0 * SL_G1(2,3) * SL_G1(2,11)
+ 21600.0 * SL_G1(0,3) * SL_G1(2,9)
+ 3087.0 * SL_G1(2,7) * SL_G1(2,7)
- 55230.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(2,7)
+ 17640.0 * SL_G1(-1,4) * SL_G1(2,7)
+ 1176.0 * SL_G1(-3,2) * SL_G1(2,7)
- 75600.0 * SL_G1(0,9) * SL_G1(2,3)
+ 37800.0 * SL_G1(-3,6) * SL_G1(2,3)
+ 4050.0 * SL_G1(0,3) * SL_G1(0,7)
+ 45360.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,5)
+ 175.0 * SL_G1(0,2) * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,3)
+ 4200.0 * SL_G1(-1,4) * SL_G1(0,2) * SL_G1(0,3)
+ 280.0 * SL_G1(-3,2) * SL_G1(0,2) * SL_G1(0,3)
+ 75600.0 * SL_G1(-2,5) * SL_G1(0,3)
+ 8640.0 * SL_G1(-5,2) * SL_G1(0,3)
+ 25200.0 * SL_G1(-1,4) * SL_G1(-1,4)
+ 3360.0 * SL_G1(-3,2) * SL_G1(-1,4)
+ 112.0 * SL_G1(-3,2) * SL_G1(-3,2)
+ 5400.0 * SL_G1(-7,3) ) / 175.0
;
fk[4] =
- 6.0 * SL_G1(0,3) * SL_G2(2,3) * ( 7884800.0 * SL_G2(2,5) * SL_G1(2,11)
- 3548160.0 * SL_G2(0,3) * SL_G1(2,11)
- 4300800.0 * SL_G1(0,2) * SL_G2(2,5) * SL_G1(2,9)
+ 131220.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,9)
+ 1003520.0 * SL_G1(0,2) * SL_G2(0,3) * SL_G1(2,9)
+ 66339.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,7)
- 5376000.0 * SL_G1(0,9) * SL_G2(2,5)
+ 2880000.0 * SL_G1(0,2) * SL_G1(0,7) * SL_G2(2,5)
+ 1075200.0 * SL_G1(-3,6) * SL_G2(2,5)
+ 256000.0 * SL_G1(-7,2) * SL_G2(2,5)
- 76545.0 * SL_G1(0,7) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 102060.0 * SL_G1(-2,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 61236.0 * SL_G1(0,3) * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 408240.0 * SL_G1(-1,4) * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 11664.0 * SL_G1(-5,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 238140.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 27216.0 * SL_G1(-3,3) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1505280.0 * SL_G2(0,3) * SL_G1(0,9)
+ 143360.0 * SL_G1(-3,6) * SL_G2(0,3) ) / 35.0
;
fk[5] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 7776000.0 * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,9)
- 8812800.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,9)
+ 11370240.0 * SL_G1(0,3) * SL_G2(2,5) * SL_G1(2,7)
+ 107163.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,7)
- 2325120.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,7)
- 4838400.0 * SL_G1(0,7) * SL_G1(2,3) * SL_G2(2,5)
+ 2016000.0 * SL_G1(-2,5) * SL_G1(2,3) * SL_G2(2,5)
- 5322240.0 * SL_G1(0,3) * SL_G1(0,5) * SL_G2(2,5)
+ 1411200.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5)
+ 12096000.0 * SL_G1(-1,4) * SL_G1(0,3) * SL_G2(2,5)
+ 1128960.0 * SL_G1(-3,2) * SL_G1(0,3) * SL_G2(2,5)
+ 368640.0 * SL_G1(-5,3) * SL_G2(2,5)
+ 306180.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 1403325.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 20412.0 * SL_G1(-3,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 3265920.0 * SL_G2(0,3) * SL_G1(0,7) * SL_G1(2,3)
+ 1693440.0 * SL_G1(-2,5) * SL_G2(0,3) * SL_G1(2,3)
+ 3677184.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(0,5)
+ 22400.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3)
+ 1236480.0 * SL_G1(-1,4) * SL_G1(0,3) * SL_G2(0,3)
+ 17920.0 * SL_G1(-3,2) * SL_G1(0,3) * SL_G2(0,3)
+ 55296.0 * SL_G1(-5,3) * SL_G2(0,3) ) / 70.0
;
fk[6] =
- 3.0 * SL_G1(0,3) * ( 333250560.0 * SL_G2(2,3) * SL_G2(2,7) * SL_G1(2,9)
- 294912000.0 * SL_G2(0,3) * SL_G2(2,5) * SL_G1(2,9)
- 141557760.0 * SL_G2(0,5) * SL_G2(2,3) * SL_G1(2,9)
- 23592960.0 * SL_G2(-1,4) * SL_G2(2,3) * SL_G1(2,9)
+ 82575360.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,9)
- 129024000.0 * SL_G1(0,2) * SL_G2(2,3) * SL_G1(2,7) * SL_G2(2,7)
- 174899200.0 * SL_G1(0,7) * SL_G2(2,3) * SL_G2(2,7)
+ 68812800.0 * SL_G1(0,2) * SL_G1(0,5) * SL_G2(2,3) * SL_G2(2,7)
+ 57344000.0 * SL_G1(-2,5) * SL_G2(2,3) * SL_G2(2,7)
+ 13107200.0 * SL_G1(-5,2) * SL_G2(2,3) * SL_G2(2,7)
+ 5806080.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5) * SL_G1(2,7)
+ 122880000.0 * SL_G1(0,2) * SL_G2(0,3) * SL_G2(2,5) * SL_G1(2,7)
- 11321856.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,7)
- 3110400.0 * SL_G1(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
+ 4147200.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
+ 37843200.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
+ 414720.0 * SL_G1(-3,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
+ 4354560.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,5)
+ 1476225.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 4976640.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 8709120.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 14100480.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 331776.0 * SL_G1(-3,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 73728000.0 * SL_G2(0,5) * SL_G1(0,7) * SL_G2(2,3)
+ 16384000.0 * SL_G2(-1,4) * SL_G1(0,7) * SL_G2(2,3)
+ 13107200.0 * SL_G1(-2,5) * SL_G2(-1,4) * SL_G2(2,3) ) / 320.0
;
fk[7] =
- 3.0 * SL_G1(0,3) * SL_G2(2,3) * ( 852590592.0 * SL_G1(2,3) * SL_G2(2,3) * SL_G1(2,7) * SL_G2(2,7)
+ 132120576.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,7)
- 528482304.0 * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 330301440.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 3049783296.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,7)
+ 44040192.0 * SL_G1(-3,3) * SL_G2(2,3) * SL_G2(2,7)
- 3501195264.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,5) * SL_G1(2,7)
- 462422016.0 * SL_G2(0,5) * SL_G1(2,3) * SL_G2(2,3) * SL_G1(2,7)
- 396361728.0 * SL_G2(-1,4) * SL_G1(2,3) * SL_G2(2,3) * SL_G1(2,7)
+ 1546911744.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,7)
- 33030144.0 * SL_G2(-3,3) * SL_G1(2,3) * SL_G1(2,7)
+ 2949120000.0 * SL_G1(0,3) * SL_G1(2,5) * SL_G2(2,5) * SL_G2(2,5)
+ 1769472000.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,5) * SL_G2(2,5)
- 1533542400.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
+ 83980800.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
- 188743680.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,5)
+ 566231040.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G1(2,3) * SL_G2(2,5)
+ 774635520.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 25165824.0 * SL_G1(-3,3) * SL_G2(0,3) * SL_G2(2,5)
- 70778880.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3) * SL_G1(2,5)
+ 94371840.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3) * SL_G1(2,5)
+ 802160640.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 9437184.0 * SL_G2(-3,3) * SL_G1(0,3) * SL_G1(2,5)
+ 2657205.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 382579200.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 283115520.0 * SL_G1(0,5) * SL_G2(0,5) * SL_G1(2,3) * SL_G2(2,3)
+ 165150720.0 * SL_G1(-1,4) * SL_G2(0,5) * SL_G1(2,3) * SL_G2(2,3)
+ 188743680.0 * SL_G2(-1,4) * SL_G1(0,5) * SL_G1(2,3) * SL_G2(2,3)
+ 424673280.0 * SL_G1(-1,4) * SL_G2(-1,4) * SL_G1(2,3) * SL_G2(2,3)
+ 88473600.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 14155776.0 * SL_G1(-3,3) * SL_G2(0,5) * SL_G2(2,3)
+ 220200960.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G1(0,3) * SL_G2(2,3)
+ 3145728.0 * SL_G1(-3,3) * SL_G2(-1,4) * SL_G2(2,3)
+ 165150720.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(0,5) * SL_G1(2,3)
+ 9437184.0 * SL_G2(-3,3) * SL_G1(0,5) * SL_G1(2,3)
+ 125829120.0 * SL_G1(-1,4) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3)
+ 14417920.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3)
+ 1048576.0 * SL_G1(-3,3) * SL_G2(0,3) * SL_G2(0,3)
+ 262144.0 * SL_G2(-3,3) * SL_G1(0,3) * SL_G1(0,3) ) / 10240.0
;
fk[8] =
- SL_G1(0,3) * ( 349913088.0 * SL_G2(2,3) * SL_G1(2,7) * SL_G2(2,9)
- 132710400.0 * SL_G1(0,5) * SL_G2(2,3) * SL_G2(2,9)
+ 20643840.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,9)
+ 123863040.0 * SL_G1(-1,4) * SL_G2(2,3) * SL_G2(2,9)
+ 16515072.0 * SL_G1(-3,2) * SL_G2(2,3) * SL_G2(2,9)
- 612147200.0 * SL_G2(0,3) * SL_G1(2,7) * SL_G2(2,7)
- 1524096.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 172032000.0 * SL_G2(0,3) * SL_G1(0,5) * SL_G2(2,7)
- 193536000.0 * SL_G2(0,7) * SL_G2(2,3) * SL_G1(2,7)
- 25804800.0 * SL_G2(-2,5) * SL_G2(2,3) * SL_G1(2,7)
+ 258048000.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G1(2,7)
+ 57344000.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G1(2,7)
- 23328000.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,5)
+ 31259520.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,5)
+ 27993600.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G1(2,5)
+ 2066715.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 11664000.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 15240960.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 19595520.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 331776.0 * SL_G2(-3,3) * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 73728000.0 * SL_G1(0,5) * SL_G2(0,7) * SL_G2(2,3)
+ 14745600.0 * SL_G2(-2,5) * SL_G1(0,5) * SL_G2(2,3)
+ 29491200.0 * SL_G2(-2,5) * SL_G1(-1,4) * SL_G2(2,3) ) / 160.0
;
fk[9] =
- 3.0 * SL_G1(0,3) * SL_G1(2,3) * ( 4043520.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,9)
- 15052800.0 * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,5) * SL_G2(2,7)
- 5080320.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,5) * SL_G2(2,7)
- 36126720.0 * SL_G2(0,3) * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,7)
- 107163.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 15262464.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,3) * SL_G2(2,7)
+ 8064000.0 * SL_G2(0,5) * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,5)
- 10752000.0 * SL_G2(-1,4) * SL_G2(2,3) * SL_G1(2,5) * SL_G2(2,5)
- 6988800.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,5)
- 1075200.0 * SL_G2(-3,3) * SL_G1(2,5) * SL_G2(2,5)
+ 4838400.0 * SL_G1(0,3) * SL_G2(0,5) * SL_G2(2,3) * SL_G2(2,5)
+ 4838400.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,5)
+ 3360000.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G2(2,5)
+ 107520.0 * SL_G2(-3,3) * SL_G1(0,3) * SL_G2(2,5)
+ 16128000.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G2(2,3) * SL_G1(2,5)
+ 12902400.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G2(2,3) * SL_G1(2,5)
+ 4838400.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,5)
+ 645120.0 * SL_G2(-3,3) * SL_G2(0,3) * SL_G1(2,5)
+ 229635.0 * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 612360.0 * SL_G2(-1,4) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 2372895.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 30618.0 * SL_G2(-3,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 8225280.0 * SL_G1(0,3) * SL_G2(0,7) * SL_G2(2,3) * SL_G2(2,3)
+ 5402880.0 * SL_G2(-2,5) * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3)
+ 3257856.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,5) * SL_G2(2,3)
+ 3333120.0 * SL_G2(-1,4) * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,3)
+ 36864.0 * SL_G2(-5,3) * SL_G1(0,3) * SL_G2(2,3)
+ 98560.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,3)
+ 14336.0 * SL_G2(-3,3) * SL_G1(0,3) * SL_G2(0,3) ) / 140.0
;
fk[10] =
- 6.0 * SL_G1(0,3) * ( 7687680.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,11)
- 8064000.0 * SL_G2(0,3) * SL_G1(2,5) * SL_G2(2,9)
- 65610.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,9)
+ 752640.0 * SL_G1(0,3) * SL_G2(0,3) * SL_G2(2,9)
- 229635.0 * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,7)
+ 4480000.0 * SL_G2(0,3) * SL_G2(0,7) * SL_G1(2,5)
+ 896000.0 * SL_G2(-2,5) * SL_G2(0,3) * SL_G1(2,5)
+ 76545.0 * SL_G2(0,7) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 119070.0 * SL_G2(-2,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 261954.0 * SL_G2(0,3) * SL_G2(0,5) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 374220.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 7776.0 * SL_G2(-5,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3) * SL_G2(2,3)
+ 189000.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
+ 15120.0 * SL_G2(-3,3) * SL_G2(0,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,3) * SL_G2(2,3)
+ 250880.0 * SL_G1(0,3) * SL_G2(0,9) * SL_G2(2,3)
+ 645120.0 * SL_G2(-3,6) * SL_G1(0,3) * SL_G2(2,3) ) / 35.0
;
fk[11] =
288.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,3) * ( 6930.0 * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,11)
+ 13590.0 * SL_G2(0,3) * SL_G2(2,3) * SL_G2(2,9)
+ 1029.0 * SL_G2(2,3) * SL_G2(2,7) * SL_G2(2,7)
- 882.0 * SL_G2(0,5) * SL_G2(2,3) * SL_G2(2,7)
+ 2352.0 * SL_G2(-1,4) * SL_G2(2,3) * SL_G2(2,7)
+ 2828.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G2(2,7)
+ 196.0 * SL_G2(-3,3) * SL_G2(2,7)
- 6300.0 * SL_G2(0,9) * SL_G2(2,3) * SL_G2(2,3)
- 3780.0 * SL_G2(-3,6) * SL_G2(2,3) * SL_G2(2,3)
- 11790.0 * SL_G2(0,3) * SL_G2(0,7) * SL_G2(2,3)
- 2520.0 * SL_G2(-1,4) * SL_G2(0,5) * SL_G2(2,3)
- 6090.0 * SL_G2(-2,5) * SL_G2(0,3) * SL_G2(2,3)
- 1680.0 * SL_G2(-1,4) * SL_G2(-1,4) * SL_G2(2,3)
- 180.0 * SL_G2(-7,3) * SL_G2(2,3)
- 2478.0 * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,5)
- 168.0 * SL_G2(-3,3) * SL_G2(0,5)
- 35.0 * SL_G2(0,2) * SL_G2(0,3) * SL_G2(0,3) * SL_G2(0,3)
- 1400.0 * SL_G2(-1,4) * SL_G2(0,3) * SL_G2(0,3)
- 28.0 * SL_G2(-3,2) * SL_G2(0,3) * SL_G2(0,3)
- 192.0 * SL_G2(-5,3) * SL_G2(0,3)
- 112.0 * SL_G2(-3,3) * SL_G2(-1,4) ) / 35.0
;
fk[12] =
1536.0 * SL_G1(0,3) * SL_G1(2,3) * ( 156.0 * SL_G2(2,3) * SL_G2(2,13)
+ 33.0 * SL_G2(0,3) * SL_G2(2,11)
- 135.0 * SL_G2(0,11) * SL_G2(2,3)
- 30.0 * SL_G2(-4,7) * SL_G2(2,3)
- 28.0 * SL_G2(0,3) * SL_G2(0,9)
- 8.0 * SL_G2(-3,6) * SL_G2(0,3) )
;
f[14] = 0.0;
for (q = 0; q <= 14; q++) f[14] = f[14] * l + fk[14-q];
if (n == 14)
{
free (fk);
return;
}
/* k = 15 */
twobody_YB_slip_15 (n, l, f, fk);
if (n == 15)
{
free (fk);
return;
}
/* k = 16 */
twobody_YB_slip_16 (n, l, f, fk);
if (n == 16)
{
free (fk);
return;
}
/* k = 17 */
twobody_YB_slip_17 (n, l, f, fk);
if (n == 17)
{
free (fk);
return;
}
/* k = 18 */
twobody_YB_slip_18 (n, l, f, fk);
if (n == 18)
{
free (fk);
return;
}
/* k = 19
twobody_YB_slip_19 (n, l, f, fk);
if (n == 19)
{
free (fk);
return;
}
*/
/* k = 20
twobody_YB_slip_20 (n, l, f, fk);
if (n == 20)
{
free (fk);
return;
}
*/
//fprintf (stderr, "twobody_YB_slip is implemented only up to k=20\n");
fprintf (stderr, "twobody_YB_slip is implemented only up to k=18\n");
free (fk);
return;
}
static void twobody_YB_slip_15 (int n, double l, double *f, double *fk)
{
int q;
for (q = 0; q <= 15; q++) fk[q] = 0.0;
fk[2] =
- 768.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * ( 13104.0 * SL_G1(2,3) * SL_G1(2,13)
- 7854.0 * SL_G1(0,3) * SL_G1(2,11)
- 1848.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(2,9)
- 10395.0 * SL_G1(0,11) * SL_G1(2,3)
+ 2772.0 * SL_G1(-4,7) * SL_G1(2,3)
+ 6776.0 * SL_G1(0,3) * SL_G1(0,9)
+ 1485.0 * SL_G1(0,2) * SL_G1(0,3) * SL_G1(0,7)
+ 924.0 * SL_G1(-3,6) * SL_G1(0,3)
+ 132.0 * SL_G1(-7,2) * SL_G1(0,3)
+ 448.0 * SL_G1(-9,3) ) / 77.0
;
fk[3] =
- 144.0 * SL_G1(0,3) * SL_G2(2,3) * SL_G2(2,3) * SL_G2(2,3) * ( 103950.0 * SL_G1(2,3) * SL_G1(2,3) * SL_G1(2,11)
+ 37800.0 * SL_G1(0,3) * SL_G1(2,3) * SL_G1(2,9)
+ 3087.0 * SL_G1(2,3) * SL_G1(2,7) * SL_G1(2,7)
+ 17640.0 * SL_G1(-1,4) * SL_G1(2,3) * SL_G1(2,7)
- 55860.0 * SL_G1(0,3) * SL_G1(0,3) * SL_G1(2,7)
+ 1176.0 * SL_G1(-3,3) * SL_G1(2,7)
- 75600.0 * SL_G1(0,9) * SL_G1(2,3) * SL_G1(2,3)
+ 37800.0 * SL_G1(-3,6) * SL_G1(2,3) * SL_G1(2,3)
- 5400.0 * SL_G1(0,3) * SL_G1(0,7) * SL_G1(2,3)