1
/*****************************************************/
2
/******** Periodic Finite Difference Formulae ********/
3
/*****************************************************/
4
/* Author: Pearu Peterson, September 2002 */
6
/* This file is automatically generated using */
7
/* genpfdf.cc utility from fdf package and it */
8
/* contains the following C function */
10
/* void periodic_finite_difference(int n, */
11
/* double *x, double *y, double h, int k, int m); */
13
/* n -- length of arrays x and y */
14
/* x -- input array */
15
/* y -- output array */
16
/* h -- discretization step of arrays x and y */
17
/* k -- the order of derivative */
18
/* m -- number of grid intervals used to evaluate */
19
/* finite differences at the center point. The */
20
/* error is O(h^(2*(m-1))) within numerical */
22
/*****************************************************/
24
/* http://epubs.siam.org/sam-bin/dbq/article/32250 */
25
/*****************************************************/
26
/* Parameters to genpfdf n=25 m=10 */
29
static double fd_coeffs_1_1[] = {
32
static double fd_coeffs_2_1[] = {
35
static double fd_coeffs_2_2[] = {
39
static double fd_coeffs_3_1[] = {
41
-0.020833333333333333334
43
static double fd_coeffs_3_2[] = {
47
static double fd_coeffs_3_3[] = {
51
static double fd_coeffs_4_1[] = {
52
0.6666666666666666667,
53
-0.083333333333333333336
55
static double fd_coeffs_4_2[] = {
57
1.3333333333333333334,
58
-0.083333333333333333336
60
static double fd_coeffs_4_3[] = {
64
static double fd_coeffs_4_4[] = {
69
static double fd_coeffs_5_1[] = {
71
-0.032552083333333333332,
74
static double fd_coeffs_5_2[] = {
75
-0.17708333333333333333,
77
-0.026041666666666666666
79
static double fd_coeffs_5_3[] = {
84
static double fd_coeffs_5_4[] = {
89
static double fd_coeffs_5_5[] = {
94
static double fd_coeffs_6_1[] = {
97
0.016666666666666666667
99
static double fd_coeffs_6_2[] = {
100
-2.7222222222222222222,
103
0.011111111111111111111
105
static double fd_coeffs_6_3[] = {
110
static double fd_coeffs_6_4[] = {
111
9.333333333333333333,
114
-0.16666666666666666667
116
static double fd_coeffs_6_5[] = {
121
static double fd_coeffs_6_6[] = {
127
static double fd_coeffs_7_1[] = {
129
-0.039876302083333333332,
131
-3.4877232142857142856E-4
133
static double fd_coeffs_7_2[] = {
134
-0.20518663194444444444,
136
-0.054144965277777777777,
137
0.0056206597222222222224
139
static double fd_coeffs_7_3[] = {
140
-0.6155598958333333333,
142
-0.032486979166666666666,
143
0.0024088541666666666667
145
static double fd_coeffs_7_4[] = {
146
0.108072916666666666664,
148
0.076822916666666666664,
149
-0.009114583333333333333
151
static double fd_coeffs_7_5[] = {
152
0.5403645833333333333,
154
0.076822916666666666664,
155
-0.0065104166666666666665
157
static double fd_coeffs_7_6[] = {
163
static double fd_coeffs_7_7[] = {
169
static double fd_coeffs_8_1[] = {
172
0.03809523809523809524,
173
-0.0035714285714285714287
175
static double fd_coeffs_8_2[] = {
176
-2.8472222222222222222,
179
0.025396825396825396825,
180
-0.0017857142857142857143
182
static double fd_coeffs_8_3[] = {
183
-2.0333333333333333333,
184
1.4083333333333333333,
186
0.029166666666666666667
188
static double fd_coeffs_8_4[] = {
190
-8.133333333333333333,
191
2.8166666666666666666,
193
0.029166666666666666667
195
static double fd_coeffs_8_5[] = {
196
4.8333333333333333335,
197
-4.3333333333333333335,
199
-0.16666666666666666667
201
static double fd_coeffs_8_6[] = {
208
static double fd_coeffs_8_7[] = {
214
static double fd_coeffs_8_8[] = {
221
static double fd_coeffs_9_1[] = {
225
-8.828299386160714286E-4,
226
5.9339735243055555557E-5
228
static double fd_coeffs_9_2[] = {
230
0.28874782986111111112,
231
-0.07917131696428571429,
233
-0.0012513175843253968254
235
static double fd_coeffs_9_3[] = {
237
0.28874782986111111112,
238
-0.04750279017857142857,
239
0.0061628069196428571427,
240
-4.171058614417989418E-4
242
static double fd_coeffs_9_4[] = {
243
0.14020182291666666667,
245
0.122721354166666666664,
246
-0.025179036458333333333,
249
static double fd_coeffs_9_5[] = {
250
0.7010091145833333333,
251
-0.40006510416666666666,
252
0.122721354166666666664,
253
-0.017985026041666666666,
254
0.0012749565972222222222
256
static double fd_coeffs_9_6[] = {
263
static double fd_coeffs_9_7[] = {
268
-0.0022786458333333333333
270
static double fd_coeffs_9_8[] = {
277
static double fd_coeffs_9_9[] = {
284
static double fd_coeffs_10_1[] = {
285
0.8333333333333333333,
286
-0.23809523809523809523,
287
0.05952380952380952381,
288
-0.009920634920634920635,
289
7.936507936507936508E-4
291
static double fd_coeffs_10_2[] = {
292
-2.9272222222222222222,
293
1.6666666666666666666,
294
-0.23809523809523809523,
295
0.03968253968253968254,
296
-0.0049603174603174603174,
297
3.174603174603174603E-4
299
static double fd_coeffs_10_3[] = {
300
-2.3180555555555555556,
301
1.7337301587301587302,
302
-0.4830357142857142857,
303
0.08339947089947089947,
304
-0.0067791005291005291007
306
static double fd_coeffs_10_4[] = {
307
12.741666666666666667,
308
-9.272222222222222222,
309
3.4674603174603174603,
310
-0.6440476190476190476,
311
0.08339947089947089947,
312
-0.0054232804232804232802
314
static double fd_coeffs_10_5[] = {
315
6.7291666666666666665,
318
-0.52777777777777777775,
319
0.04513888888888888889
321
static double fd_coeffs_10_6[] = {
326
-0.7916666666666666667,
327
0.054166666666666666667
329
static double fd_coeffs_10_7[] = {
333
2.1666666666666666667,
334
-0.20833333333333333333
336
static double fd_coeffs_10_8[] = {
341
4.3333333333333333335,
342
-0.33333333333333333334
344
static double fd_coeffs_10_9[] = {
351
static double fd_coeffs_10_10[] = {
359
static double fd_coeffs_11_1[] = {
360
0.610668182373046875,
361
-0.048465728759765625,
362
0.0087238311767578125,
363
-0.0014836447579520089286,
364
1.7950269911024305556E-4,
365
-1.0923905806107954545E-5
367
static double fd_coeffs_11_2[] = {
368
-0.23465464274088541667,
369
0.31435403975229414683,
370
-0.10050982520693824405,
371
0.024337853325737847222,
372
-0.0038119385734437003968,
373
2.845134432353670635E-4
375
static double fd_coeffs_11_3[] = {
376
-0.70396392822265625,
377
0.31435403975229414683,
378
-0.060305895124162946428,
379
0.0104305085681733630955,
380
-0.0012706461911479001323,
381
7.7594575427827380955E-5
383
static double fd_coeffs_11_4[] = {
384
0.16364067925347222222,
385
-0.29026518322172619049,
386
0.16457645476810515873,
387
-0.04471141673900462963,
388
0.0073175339471726190477,
389
-5.580680080191798942E-4
391
static double fd_coeffs_11_5[] = {
392
0.8182033962673611111,
393
-0.48377530536954365078,
394
0.16457645476810515873,
395
-0.031936726242146164022,
396
0.004065296637318121693,
397
-2.5366727637235449736E-4
399
static double fd_coeffs_11_6[] = {
403
0.058085123697916666668,
405
8.5042317708333333334E-4
407
static double fd_coeffs_11_7[] = {
411
0.058085123697916666668,
412
-0.008231608072916666667,
413
5.411783854166666667E-4
415
static double fd_coeffs_11_8[] = {
419
-0.045003255208333333332,
421
-8.951822916666666667E-4
423
static double fd_coeffs_11_9[] = {
431
static double fd_coeffs_11_10[] = {
439
static double fd_coeffs_12_1[] = {
440
0.85714285714285714287,
441
-0.26785714285714285715,
442
0.07936507936507936508,
443
-0.017857142857142857143,
444
0.0025974025974025974027,
445
-1.8037518037518037518E-4
447
static double fd_coeffs_12_2[] = {
448
-2.9827777777777777777,
449
1.7142857142857142857,
450
-0.26785714285714285715,
451
0.052910052910052910053,
452
-0.008928571428571428571,
453
0.0010389610389610389611,
454
-6.012506012506012506E-5
456
static double fd_coeffs_12_3[] = {
457
-2.5271428571428571428,
458
1.9950892857142857143,
459
-0.65727513227513227514,
460
0.15309523809523809524,
461
-0.022619047619047619048,
462
0.0015839947089947089947
464
static double fd_coeffs_12_4[] = {
465
13.717407407407407408,
466
-10.108571428571428571,
467
3.9901785714285714286,
468
-0.8763668430335097002,
469
0.15309523809523809524,
470
-0.018095238095238095238,
471
0.0010559964726631393298
473
static double fd_coeffs_12_5[] = {
474
8.246031746031746032,
475
-8.396081349206349206,
476
3.9828042328042328043,
477
-1.0333994708994708995,
478
0.16005291005291005292,
479
-0.011491402116402116402
481
static double fd_coeffs_12_6[] = {
482
-61.768055555555555554,
483
49.476190476190476192,
484
-25.18824404761904762,
485
7.9656084656084656086,
486
-1.5500992063492063492,
487
0.19206349206349206349,
488
-0.011491402116402116402
490
static double fd_coeffs_12_7[] = {
493
-15.729166666666666667,
494
5.0083333333333333332,
495
-0.8541666666666666667,
496
0.06458333333333333333
498
static double fd_coeffs_12_8[] = {
499
233.56666666666666667,
502
-41.944444444444444444,
503
10.0166666666666666664,
504
-1.3666666666666666667,
505
0.08611111111111111111
507
static double fd_coeffs_12_9[] = {
515
static double fd_coeffs_12_10[] = {
519
136.66666666666666667,
522
-0.41666666666666666666
524
static double fd_coeffs_13_1[] = {
525
0.6143031120300292969,
526
-0.051191926002502441406,
527
0.010238385200500488281,
528
-0.002089466367449079241,
529
3.4472677442762586804E-4,
530
-3.8461251692338423294E-5,
531
2.1182573758638822115E-6
533
static double fd_coeffs_13_2[] = {
534
-0.24332125588068886408,
535
0.33385391931685190353,
536
-0.118565269248195426174,
537
0.034448901988841869214,
538
-0.0073573712215451107053,
539
0.0010067312048856543485,
540
-6.5656160150026116824E-5
542
static double fd_coeffs_13_3[] = {
543
-0.72996376764206659225,
544
0.33385391931685190353,
545
-0.0711391615489172557,
546
0.0147638151380750868054,
547
-0.0024524570738483702352,
548
2.7456305587790573143E-4,
549
-1.5151421573082950036E-5
551
static double fd_coeffs_13_4[] = {
552
0.18137577682575851521,
553
-0.3301691527593703497,
554
0.20152457471036826912,
555
-0.06540236390667197145,
556
0.014572801135835193452,
557
-0.00203599280570970431,
558
1.3435679979004767416E-4
560
static double fd_coeffs_13_5[] = {
561
0.90687888412879257607,
562
-0.5502819212656172495,
563
0.20152457471036826912,
564
-0.04671597421905140818,
565
0.008096000631019551918,
566
-9.2545127532259286815E-4,
567
5.167569222694141314E-5
569
static double fd_coeffs_13_6[] = {
570
-0.14570670960441468254,
571
0.2961811610630580357,
572
-0.22375217821232225529,
573
0.09298107005931712963,
574
-0.022819737025669642856,
575
0.003342990774326223545,
576
-2.2659705429480820106E-4
578
static double fd_coeffs_13_7[] = {
579
-1.0199469672309027778,
580
0.69108937581380208335,
581
-0.3132530494972511574,
582
0.09298107005931712963,
583
-0.017748684353298611112,
584
0.002127357765480324074,
585
-1.2201379846643518518E-4
587
static double fd_coeffs_13_8[] = {
588
0.10392659505208333333,
590
0.18742031521267361111,
591
-0.09048597547743055556,
593
-0.004143948025173611111,
594
2.9534233940972222222E-4
596
static double fd_coeffs_13_9[] = {
605
static double fd_coeffs_13_10[] = {
608
-0.109354654947916666664,
609
0.057820638020833333332,
611
0.0033976236979166666667,
612
-2.6448567708333333334E-4
614
static double fd_coeffs_14_1[] = {
616
-0.29166666666666666666,
617
0.09722222222222222222,
618
-0.026515151515151515152,
619
0.0053030303030303030304,
620
-6.7987567987567987567E-4,
621
4.1625041625041625042E-5
623
static double fd_coeffs_14_2[] = {
624
-3.0235941043083900226,
626
-0.29166666666666666666,
627
0.064814814814814814815,
628
-0.013257575757575757576,
629
0.0021212121212121212122,
630
-2.2662522662522662522E-4,
631
1.1892869035726178583E-5
633
static double fd_coeffs_14_3[] = {
634
-2.6869345238095238095,
636
-0.8170667989417989418,
637
0.23056998556998556999,
638
-0.046829906204906204906,
639
0.006053691678691678692,
640
-3.7247474747474747474E-4
642
static double fd_coeffs_14_4[] = {
643
14.447883597883597884,
644
-10.747738095238095238,
646
-1.0894223985890652557,
647
0.23056998556998556999,
648
-0.037463924963924963926,
649
0.0040357944524611191277,
650
-2.1284271284271284271E-4
652
static double fd_coeffs_14_5[] = {
653
9.4708002645502645505,
654
-10.029106040564373898,
655
5.207572751322751323,
656
-1.62722663139329806,
657
0.34562389770723104057,
658
-0.04575066137566137566,
659
0.0028549382716049382716
661
static double fd_coeffs_14_6[] = {
662
-70.16646825396825397,
663
56.824801587301587303,
664
-30.087318121693121692,
665
10.415145502645502646,
666
-2.44083994708994709,
667
0.4147486772486772487,
668
-0.04575066137566137566,
669
0.00244708994708994709
671
static double fd_coeffs_14_7[] = {
672
-31.996006944444444445,
673
37.950925925925925924,
674
-23.450173611111111111,
675
8.751851851851851852,
676
-2.0240162037037037037,
677
0.28055555555555555556,
678
-0.017997685185185185186
680
static double fd_coeffs_14_8[] = {
681
304.15873015873015872,
682
-255.96805555555555556,
683
151.8037037037037037,
684
-62.533796296296296297,
685
17.503703703703703704,
686
-3.238425925925925926,
687
0.37407407407407407407,
688
-0.020568783068783068783
690
static double fd_coeffs_14_9[] = {
699
static double fd_coeffs_14_10[] = {
703
261.79166666666666666,
706
-2.1666666666666666667,
709
static double fd_coeffs_15_1[] = {
710
0.6170455366373062134,
711
-0.053324922919273376465,
712
0.011518183350563049316,
713
-0.0026711927992956978933,
714
5.386355850431654188E-4,
715
-8.320943875746293501E-5,
716
8.510855528024526743E-6,
717
-4.2617321014404296876E-7
719
static double fd_coeffs_15_2[] = {
720
-0.24989236070185291524,
721
0.3491864972329013562,
722
-0.13389784716424487886,
723
0.04420599702632788456,
724
-0.011538983380467688711,
725
0.0021861602753509968628,
726
-2.6478054867014887898E-4,
727
1.5317260655394058626E-5
729
static double fd_coeffs_15_3[] = {
730
-0.7496770821055587457,
731
0.3491864972329013562,
732
-0.080338708298546927314,
733
0.01894542729699766481,
734
-0.0038463277934892295702,
735
5.962255296411809626E-4,
736
-6.1103203539265125915E-5,
737
3.0634521310788117254E-6
739
static double fd_coeffs_15_4[] = {
740
0.19522055739448184059,
741
-0.36247364075305810893,
742
0.23382906270405602833,
743
-0.085959765357200545494,
744
0.023383116043204582327,
745
-0.0045209534206087627107,
746
5.538956049028757159E-4,
747
-3.227221577790984936E-5
749
static double fd_coeffs_15_5[] = {
750
0.97610278697240920295,
751
-0.6041227345884301815,
752
0.23382906270405602833,
753
-0.06139983239800038964,
754
0.0129906200240025457374,
755
-0.0020549788275494375958,
756
2.1303677111649065995E-4,
757
-1.075740525930328312E-5
759
static double fd_coeffs_15_6[] = {
760
-0.17064858466859847781,
761
0.354378869546153558,
762
-0.28194988669541777758,
763
0.13001597545765064381,
764
-0.038691839339241148935,
765
0.0078197375807181867976,
766
-9.824114501791656333E-4,
767
5.813956891418134094E-5
769
static double fd_coeffs_15_7[] = {
770
-1.1945400926801893446,
771
0.82688402894102496865,
772
-0.3947298413735848886,
773
0.13001597545765064381,
774
-0.030093652819409782504,
775
0.0049761966422752097802,
776
-5.289907808657045718E-4,
777
2.7131798826617959105E-5
779
static double fd_coeffs_15_8[] = {
780
0.140817230466812376,
781
-0.30904873035572193288,
782
0.27349846451370804398,
783
-0.14526297957808883102,
784
0.049434110853407118054,
785
-0.0107653441252531828705,
786
0.0014132403822802992725,
787
-8.5992157143890542326E-5
789
static double fd_coeffs_15_9[] = {
790
1.267355074201311384,
791
-0.92714619106716579863,
792
0.49229723612467447916,
793
-0.18676668802897135416,
794
0.049434110853407118054,
795
-0.008808008829752604166,
796
9.78397187732514881E-4,
797
-5.1595294286334325397E-5
799
static double fd_coeffs_15_10[] = {
800
-0.097423553466796875,
801
0.22050603230794270833,
802
-0.207363128662109375,
803
0.120189666748046875,
804
-0.045406341552734375,
805
0.010936737060546875,
806
-0.001537322998046875,
807
9.7910563151041666664E-5
809
static double fd_coeffs_16_1[] = {
810
0.8888888888888888889,
811
-0.3111111111111111111,
812
0.11313131313131313131,
813
-0.035353535353535353535,
814
0.008702408702408702409,
815
-0.0015540015540015540016,
816
1.7760017760017760018E-4,
817
-9.7125097125097125096E-6
819
static double fd_coeffs_16_2[] = {
820
-3.0548441043083900226,
821
1.7777777777777777778,
822
-0.3111111111111111111,
823
0.07542087542087542087,
824
-0.017676767676767676768,
825
0.0034809634809634809634,
826
-5.1800051800051800053E-4,
827
5.074290788576502862E-5,
828
-2.4281274281274281274E-6
830
static double fd_coeffs_16_3[] = {
831
-2.8129176114890400604,
832
2.3845211640211640212,
833
-0.9613746993746993747,
834
0.31074104136604136604,
835
-0.077664927664927664926,
836
0.013982697196982911268,
837
-0.0016058756058756058756,
838
8.810006131434702863E-5
840
static double fd_coeffs_16_4[] = {
841
15.014807492441421013,
842
-11.251670445956160242,
843
4.7690423280423280424,
844
-1.2818329324995991663,
845
0.31074104136604136604,
846
-0.06213194213194213194,
847
0.0093217981313219408455,
848
-9.176432033574890718E-4,
849
4.4050030657173514316E-5
851
static double fd_coeffs_16_5[] = {
852
10.474125514403292181,
853
-11.433761390358612581,
854
6.3568362193362193363,
855
-2.265706335845224734,
856
0.59119301480412591524,
857
-0.108897005772005772006,
858
0.012677702955480733258,
859
-7.016260488482710705E-4
861
static double fd_coeffs_16_6[] = {
862
-76.93891369047619048,
863
62.844753086419753087,
864
-34.301284171075837742,
865
12.713672438672438673,
866
-3.3985595037678371013,
867
0.70943161776495109827,
868
-0.108897005772005772006,
869
0.010866602533269199936,
870
-5.2621953663620330285E-4
872
static double fd_coeffs_16_7[] = {
873
-38.817746913580246912,
874
47.50136188271604938,
875
-31.264166666666666666,
876
13.092959104938271605,
877
-3.6936728395061728396,
878
0.70989583333333333335,
879
-0.08478395061728395062,
880
0.004770447530864197531
882
static double fd_coeffs_16_8[] = {
883
365.55438988095238095,
884
-310.5419753086419753,
885
190.00544753086419753,
886
-83.371111111111111114,
887
26.18591820987654321,
888
-5.909876543209876543,
889
0.9465277777777777778,
890
-0.096895943562610229275,
891
0.004770447530864197531
893
static double fd_coeffs_16_9[] = {
894
129.55734126984126985,
895
-168.52777777777777778,
897
-57.394444444444444445,
898
17.993055555555555556,
899
-3.6928571428571428572,
900
0.45972222222222222223,
901
-0.026587301587301587302
903
static double fd_coeffs_16_10[] = {
904
-1493.7232142857142857,
905
1295.5734126984126984,
906
-842.6388888888888889,
907
406.95833333333333334,
908
-143.4861111111111111,
909
35.986111111111111112,
910
-6.154761904761904762,
911
0.65674603174603174606,
912
-0.03323412698412698413
914
static double fd_coeffs_17_1[] = {
915
0.6191880558617413044,
916
-0.05503893829882144928,
917
0.012608920410275459289,
918
-0.00321656132915190288,
919
7.483927119109365675E-4,
920
-1.431400464339689775E-4,
921
2.0496977063325735239E-5,
922
-1.9244384020566940308E-6,
923
8.8133246583097121296E-8
925
static double fd_coeffs_17_2[] = {
926
-0.25504508095871564213,
927
0.3615530258493719007,
928
-0.14701386236353182001,
929
0.053387207665828743367,
930
-0.01607914248791316834,
931
0.0037716126620779897496,
932
-6.3952384007834719766E-4,
933
6.936677383926881613E-5,
934
-3.6033008789249838336E-6
936
static double fd_coeffs_17_3[] = {
937
-0.76513524287614692636,
938
0.3615530258493719007,
939
-0.08820831741811909201,
940
0.022880231856783747156,
941
-0.00535971416263772278,
942
0.0010286216351121790226,
943
-1.4758242463346473793E-4,
944
1.3873354767853763226E-5,
945
-6.358766256926442059E-7
947
static double fd_coeffs_17_4[] = {
948
0.20631058824801210733,
949
-0.38908971480153074912,
950
0.26205823214940579825,
951
-0.10572018396894538443,
952
0.03315475162044104114,
953
-0.007933270606310383248,
954
0.0013604433033414405703,
955
-1.4860121074501054951E-4,
956
7.755266331140046677E-6
958
static double fd_coeffs_17_5[] = {
959
1.0315529412400605367,
960
-0.6484828580025512485,
961
0.26205823214940579825,
962
-0.07551441712067527459,
963
0.018419306455800578411,
964
-0.0036060320937774469312,
965
5.2324742436209252704E-4,
966
-4.9533736915003516506E-5,
967
2.2809606856294254933E-6
969
static double fd_coeffs_17_6[] = {
970
-0.19157658920506741425,
971
0.40460608043367900545,
972
-0.33522117097006597943,
973
0.16730587444990438511,
974
-0.057131899280465526498,
975
0.014259123591939398011,
976
-0.0025044481437405428291,
977
2.7766409202399535956E-4,
978
-1.4634968207320934575E-5
980
static double fd_coeffs_17_7[] = {
981
-1.3410361244354718997,
982
0.94408085434525101275,
983
-0.4693096393580923712,
984
0.16730587444990438511,
985
-0.04443592166258429839,
986
0.009073987740325071461,
987
-0.0013485490004756769079,
988
1.2957657627786450113E-4,
989
-6.0261633794850907076E-6
991
static double fd_coeffs_17_8[] = {
992
0.17448784842482740076,
993
-0.3898582134549579923,
994
0.3592054920432008343,
995
-0.20525789884873378424,
996
0.079101928075154622395,
997
-0.021125534266180882799,
998
0.0038620125974086647101,
999
-4.3918045740278940353E-4,
1000
2.3545886683926590746E-5
1002
static double fd_coeffs_17_9[] = {
1003
1.5703906358234466068,
1004
-1.169574640364873977,
1005
0.6465698856777615017,
1006
-0.26390301280551486544,
1007
0.079101928075154622395,
1008
-0.017284528035966176835,
1009
0.0026737010289752294147,
1010
-2.635082744416736421E-4,
1011
1.2465469420902312748E-5
1013
static double fd_coeffs_17_10[] = {
1014
-0.14176721421499100942,
1015
0.32693081810360863094,
1016
-0.32023790147569444444,
1017
0.19920200771755642361,
1018
-0.0844783782958984375,
1019
0.024580940367683531746,
1020
-0.0047623165070064484126,
1021
5.630538577125186012E-4,
1022
-3.1009552970765128967E-5
1024
static double fd_coeffs_18_1[] = {
1026
-0.32727272727272727271,
1027
0.12727272727272727273,
1028
-0.044055944055944055943,
1029
0.012587412587412587413,
1030
-0.0027972027972027972027,
1031
4.4955044955044955046E-4,
1032
-4.627725215960510078E-5,
1033
2.285296402943461767E-6
1035
static double fd_coeffs_18_2[] = {
1036
-3.0795354623330813807,
1038
-0.32727272727272727271,
1039
0.084848484848484848486,
1040
-0.022027972027972027972,
1041
0.005034965034965034965,
1042
-9.324009324009324009E-4,
1043
1.284429855858427287E-4,
1044
-1.1569313039901275195E-5,
1045
5.0784364509854705933E-7
1047
static double fd_coeffs_18_3[] = {
1048
-2.9147457482993197279,
1049
2.532634817563388992,
1050
-1.0909741462241462241,
1051
0.39049454711954711955,
1052
-0.1132691713048855906,
1053
0.025376055161769447483,
1054
-0.0040981726606726606727,
1055
4.2319882498453927025E-4,
1056
-2.0943672729387015101E-5
1058
static double fd_coeffs_18_4[] = {
1059
15.467376989375997313,
1060
-11.6589829931972789115,
1061
5.065269635126777984,
1062
-1.4546321949655282988,
1063
0.39049454711954711955,
1064
-0.09061533704390847248,
1065
0.016917370107846298323,
1066
-0.0023418129489558060987,
1067
2.1159941249226963512E-4,
1068
-9.308298990838673378E-6
1070
static double fd_coeffs_18_5[] = {
1071
11.308281486205593348,
1072
-12.647079167525596097,
1073
7.4184892743573299127,
1074
-2.9190312927812927812,
1075
0.88285594200772772204,
1076
-0.20222914247715835017,
1077
0.033094107859732859734,
1078
-0.003446688893117464546,
1079
1.7156642776682459222E-4
1081
static double fd_coeffs_18_6[] = {
1082
-82.499953502491531594,
1083
67.84968891723356009,
1084
-37.94123750257678829,
1085
14.836978548714659825,
1086
-4.378546939171939172,
1087
1.0594271304092732664,
1088
-0.20222914247715835017,
1089
0.028366378165485308343,
1090
-0.0025850166698380984095,
1091
1.1437761851121639481E-4
1093
static double fd_coeffs_18_7[] = {
1094
-44.801884645061728397,
1095
56.205562219416386083,
1096
-38.88034196127946128,
1097
17.779836209315375981,
1098
-5.786028689674523008,
1099
1.3794497053872053872,
1100
-0.2312488601290684624,
1101
0.024463208473625140291,
1102
-0.0012307975589225589225
1104
static double fd_coeffs_18_8[] = {
1105
418.74672527189888302,
1106
-358.41507716049382717,
1107
224.82224887766554433,
1108
-103.68091189674523008,
1109
35.559672418630751962,
1110
-9.2576459034792368125,
1111
1.8392662738496071829,
1112
-0.26428441157607824276,
1113
0.024463208473625140291,
1114
-0.0010940422745978301534
1116
static double fd_coeffs_18_9[] = {
1117
166.11278025793650794,
1118
-221.69932539682539682,
1119
168.61260416666666667,
1120
-86.02527777777777778,
1121
30.774677579365079366,
1122
-7.782976190476190476,
1123
1.3544357638888888889,
1124
-0.14688492063492063492,
1125
0.007518601190476190476
1127
static double fd_coeffs_18_10[] = {
1128
-1899.8947585978835979,
1129
1661.1278025793650793,
1130
-1108.4966269841269841,
1131
562.0420138888888889,
1132
-215.06319444444444444,
1133
61.549355158730158732,
1134
-12.971626984126984127,
1135
1.934908234126984127,
1136
-0.18360615079365079365,
1137
0.008354001322751322751
1139
static double fd_coeffs_19_1[] = {
1140
0.6209080226835794747,
1141
-0.05644618388032540679,
1142
0.013547084131278097629,
1143
-0.003721726409691785063,
1144
9.6489203214231464597E-4,
1145
-2.1530648651109500365E-4,
1146
3.8538587082607241776E-5,
1147
-5.1082519348710775376E-6,
1148
4.4189030578469528872E-7,
1149
-1.8618792589557798286E-8
1151
static double fd_coeffs_19_2[] = {
1152
-0.25919347316061063725,
1153
0.37173544307220507063,
1154
-0.15832765927779089768,
1155
0.06191606995503943268,
1156
-0.020778719667682323676,
1157
0.005686255216798756738,
1158
-0.0012052136857913010806,
1159
1.8455248902968929003E-4,
1160
-1.8108168717718673139E-5,
1161
8.532275199290405474E-7
1163
static double fd_coeffs_19_3[] = {
1164
-0.7775804194818319118,
1165
0.37173544307220507063,
1166
-0.094996595566674538606,
1167
0.026535458552159756863,
1168
-0.0069262398892274412254,
1169
0.0015507968773087518376,
1170
-2.7812623518260794168E-4,
1171
3.6910497805937858005E-5,
1172
-3.1955591854797658481E-6,
1173
1.3472013472563798117E-7
1175
static double fd_coeffs_19_4[] = {
1176
0.21538517591399155125,
1177
-0.41136370270893483873,
1178
0.28680710760207700895,
1179
-0.12437702854095906634,
1180
0.04343505373155062097,
1181
-0.012121541836762434291,
1182
0.0025978870759750011053,
1183
-4.0056939974279436888E-4,
1184
3.9484593834564675782E-5,
1185
-1.8664310296132134769E-6
1187
static double fd_coeffs_19_5[] = {
1188
1.0769258795699577563,
1189
-0.68560617118155806456,
1190
0.28680710760207700895,
1191
-0.088840734672113618816,
1192
0.024130585406417011651,
1193
-0.0055097917439829246775,
1194
9.991873369134619636E-4,
1195
-1.3352313324759812296E-4,
1196
1.1613115833695492877E-5,
1197
-4.9116606042452986233E-7
1199
static double fd_coeffs_19_6[] = {
1200
-0.20930129541797138823,
1201
0.4481121775017160325,
1202
-0.38356127882344045392,
1203
0.20374687883167898897,
1204
-0.07721163638879030822,
1205
0.02243975722866430908,
1206
-0.004921453536409266554,
1207
7.6981451134568118596E-4,
1208
-7.6609465455236927526E-5,
1209
3.6455586616421172325E-6
1211
static double fd_coeffs_19_7[] = {
1212
-1.4651090679257997176,
1213
1.0455950808373374092,
1214
-0.5369857903528166355,
1215
0.20374687883167898897,
1216
-0.06005349496905912861,
1217
0.014279845509150014869,
1218
-0.00265001344268191276,
1219
3.5924677196131788679E-4,
1220
-3.1545074010979911334E-5,
1221
1.3431005595523589804E-6
1223
static double fd_coeffs_19_8[] = {
1224
0.20477333965790685242,
1225
-0.46419532829978937366,
1226
0.44180228631523570248,
1227
-0.26752317453072930025,
1228
0.11341136569584602918,
1229
-0.03510345329683293741,
1230
0.007991852311010408119,
1231
-0.0012800980461452258441,
1232
1.2943921267371488327E-4,
1233
-6.22901917586989956E-6
1235
static double fd_coeffs_19_9[] = {
1236
1.8429600569211616718,
1237
-1.392585984899368121,
1238
0.79524411536742426447,
1239
-0.3439583672537948146,
1240
0.11341136569584602918,
1241
-0.028721007242863312428,
1242
0.005532820830699513313,
1243
-7.680588276871355065E-4,
1244
6.852664200373140879E-5,
1245
-2.9505880306752155812E-6
1247
static double fd_coeffs_19_10[] = {
1248
-0.18578297413846172354,
1249
0.43496950155212765648,
1250
-0.44028088308516002836,
1251
0.28969594770007663304,
1252
-0.13434238604136875698,
1253
0.044895906486208476717,
1254
-0.010764465587479727609,
1255
0.0017852109102975754511,
1256
-1.849108114444389545E-4,
1257
9.053015204333754444E-6
1259
static double fd_coeffs_20_1[] = {
1260
0.90909090909090909093,
1261
-0.3409090909090909091,
1262
0.13986013986013986014,
1263
-0.05244755244755244755,
1264
0.016783216783216783217,
1265
-0.0043706293706293706294,
1266
8.814714697067638244E-4,
1267
-1.2854792266556972439E-4,
1268
1.2027875804965588247E-5,
1269
-5.412544112234514711E-7
1271
static double fd_coeffs_20_2[] = {
1272
-3.0995354623330813807,
1273
1.8181818181818181819,
1274
-0.3409090909090909091,
1275
0.093240093240093240094,
1276
-0.026223776223776223776,
1277
0.006713286713286713287,
1278
-0.0014568764568764568765,
1279
2.5184899134478966413E-4,
1280
-3.2136980666392431097E-5,
1281
2.672861289992352944E-6,
1282
-1.0825088224469029423E-7
1284
static double fd_coeffs_20_3[] = {
1285
-2.9987330790902219475,
1286
2.6586158137497423212,
1287
-1.2072642965500108357,
1288
0.4680213140034568606,
1289
-0.15203255474684046113,
1290
0.03991232395250252393,
1291
-0.0080885209561680149916,
1292
0.0011832651669836543786,
1293
-1.10951529018755909514E-4,
1294
5.000436460520494134E-6
1296
static double fd_coeffs_20_4[] = {
1297
15.836921244855967078,
1298
-11.99493231636088779,
1299
5.3172316274994846424,
1300
-1.609685728733347781,
1301
0.4680213140034568606,
1302
-0.121626043797472368903,
1303
0.026608215968335015954,
1304
-0.004622011974953151424,
1305
5.916325834918271893E-4,
1306
-4.9311790675002626448E-5,
1307
2.0001745842081976535E-6
1309
static double fd_coeffs_20_5[] = {
1310
12.011344076631775044,
1311
-13.701673053164868641,
1312
8.391960553408966107,
1313
-3.5680121454823835777,
1314
1.2073463683582731201,
1315
-0.32391305235861287448,
1316
0.06649753410170076837,
1317
-0.0098092462725399233335,
1318
9.250271700668526065E-4,
1319
-4.1858930127779334128E-5
1321
static double fd_coeffs_20_6[] = {
1322
-87.140166599304330786,
1323
72.06806445979065027,
1324
-41.105019159494605923,
1325
16.783921106817932214,
1326
-5.3520182182235753665,
1327
1.4488156420299277442,
1328
-0.32391305235861287448,
1329
0.056997886372886372885,
1330
-0.0073569347044049425004,
1331
6.1668478004456840436E-4,
1332
-2.5115358076667600477E-5
1334
static double fd_coeffs_20_7[] = {
1335
-50.051881686129371316,
1336
64.08055778101785046,
1337
-46.149568633526966862,
1338
22.62598732414704637,
1339
-8.209104247090358202,
1340
2.2881030394181435849,
1341
-0.4806831086865809088,
1342
0.07197449391315132056,
1343
-0.0068571339925506592175,
1344
3.125742463126722386E-4
1346
static double fd_coeffs_20_8[] = {
1347
464.9466992332941407,
1348
-400.41505348903497052,
1349
256.32223112407140184,
1350
-123.06551635607191163,
1351
45.25197464829409274,
1352
-13.134566795344573122,
1353
3.0508040525575247798,
1354
-0.54935212421323532434,
1355
0.07197449391315132056,
1356
-0.006095230215600585971,
1357
2.5005939705013779087E-4
1359
static double fd_coeffs_20_9[] = {
1360
200.373875962000962,
1361
-273.09096895292207793,
1362
216.0510443722943723,
1363
-117.65090458152958153,
1364
46.58749098124098124,
1365
-13.7127812161796536795,
1366
2.982225378787878788,
1367
-0.45694008537758537758,
1368
0.04423566017316017316,
1369
-0.002039836610149110149
1371
static double fd_coeffs_20_10[] = {
1372
-2276.7668113425925926,
1373
2003.73875962000962,
1374
-1365.4548447646103896,
1375
720.17014790764790766,
1376
-294.12726145382395382,
1377
93.17498196248196248,
1378
-22.8546353602994228,
1379
4.260321969696969697,
1380
-0.571175106721981722,
1381
0.049150733525733525733,
1382
-0.002039836610149110149
1384
static double fd_coeffs_21_1[] = {
1385
0.62231917728058760986,
1386
-0.057622146044498852763,
1387
0.014361211783398175612,
1388
-0.0041869422109032581958,
1389
0.0011819927393743354413,
1390
-2.9671925172310280191E-4,
1391
6.248351802731541773E-5,
1392
-1.0429347700361783306E-5,
1393
1.2820633213884909363E-6,
1394
-1.0263609414993736305E-7,
1395
4.000823883827598322E-9
1397
static double fd_coeffs_21_2[] = {
1398
-0.26260485842999193177,
1399
0.3802639062456583069,
1400
-0.168168193708698478,
1401
0.06978849749976549695,
1402
-0.025502176194517962235,
1403
0.007851172791598424411,
1404
-0.0019577251422724689884,
1405
3.7750414453768106126E-4,
1406
-5.2636359703359305883E-5,
1407
4.7122606300888759717E-6,
1408
-2.0310700579788607497E-7
1410
static double fd_coeffs_21_3[] = {
1411
-0.78781457528997579534,
1412
0.3802639062456583069,
1413
-0.1009009162252190868,
1414
0.029909356071328070118,
1415
-0.008500725398172654078,
1416
0.0021412289431632066575,
1417
-4.5178272513980053578E-4,
1418
7.550082890753621225E-5,
1419
-9.288769359416348097E-6,
1420
7.4404115211929620605E-7,
1421
-2.9015286542555153567E-8
1423
static double fd_coeffs_21_4[] = {
1424
0.22294360058181273126,
1425
-0.43025976437848778873,
1426
0.30861025568233041278,
1427
-0.14181954700516178941,
1428
0.053900564810072254814,
1429
-0.016918234414418183135,
1430
0.004265186635053202575,
1431
-8.280821071987434637E-4,
1432
1.1598686780036609275E-4,
1433
-1.0416685178732195373E-5,
1434
4.500133762694200998E-7
1436
static double fd_coeffs_21_5[] = {
1437
1.1147180029090636563,
1438
-0.7170996072974796479,
1439
0.30861025568233041278,
1440
-0.10129967643225842101,
1441
0.02994475822781791934,
1442
-0.0076901065520082650613,
1443
0.0016404563980973856059,
1444
-2.7602736906624782124E-4,
1445
3.4113784647166497868E-5,
1446
-2.7412329417716303614E-6,
1447
1.07146041968909547575E-7
1449
static double fd_coeffs_21_6[] = {
1450
-0.22446233308351165624,
1451
0.48601477166556670255,
1452
-0.42729504132019122706,
1453
0.23873388882907960746,
1454
-0.09820384238723067931,
1455
0.032061184977949479167,
1456
-0.008265800080278443322,
1457
0.0016273392661839316392,
1458
-2.3006126368945016653E-4,
1459
2.0796053758407126298E-5,
1460
-9.026576366718425824E-7
1462
static double fd_coeffs_21_7[] = {
1463
-1.5712363315845815937,
1464
1.1340344672196556393,
1465
-0.5982130578482677179,
1466
0.23873388882907960746,
1467
-0.076380766301179417245,
1468
0.020402572258695123105,
1469
-0.0044508154278422387117,
1470
7.59424990885834765E-4,
1471
-9.4731108578008892104E-5,
1472
7.661704016255257057E-6,
1473
-3.0088587889061419412E-7
1475
static double fd_coeffs_21_8[] = {
1476
0.2318770803012348174,
1477
-0.5319546799081092861,
1478
0.5199861535556048322,
1479
-0.33007026832302460405,
1480
0.15093962197122321146,
1481
-0.05230390408971414596,
1482
0.013970618629391576865,
1483
-0.0028131150508583460356,
1484
4.037685714118521807E-4,
1485
-3.688935927013230339E-5,
1486
1.6137021102243370436E-6
1488
static double fd_coeffs_21_9[] = {
1489
2.0868937227111133566,
1490
-1.5958640397243278584,
1491
0.935975076400088698,
1492
-0.42437605927246020522,
1493
0.15093962197122321146,
1494
-0.042794103346129755784,
1495
0.00967196674342493783,
1496
-0.0016878690305150076214,
1497
2.1375983192392174271E-4,
1498
-1.7473907022694248974E-5,
1499
6.915866186675730187E-7
1501
static double fd_coeffs_21_10[] = {
1502
-0.2280906640097119018,
1503
0.5407387262302531022,
1504
-0.5623222961753047734,
1505
0.38732907817219242906,
1506
-0.19292226432463823459,
1507
0.071745017366040320626,
1508
-0.020097044235549619876,
1509
0.0041781797944180606476,
1510
-6.1312629597126262125E-4,
1511
5.6912392886743458377E-5,
1512
-2.518914614863668628E-6
1514
static double fd_coeffs_22_1[] = {
1515
0.9166666666666666667,
1516
-0.3525641025641025641,
1517
0.1510989010989010989,
1518
-0.06043956043956043956,
1519
0.021153846153846153846,
1520
-0.006221719457013574661,
1521
0.0014813617754794225382,
1522
-2.7288243232515678337E-4,
1523
3.6384324310020904448E-5,
1524
-3.1186563694303632384E-6,
1525
1.2887009791034558836E-7
1527
static double fd_coeffs_22_2[] = {
1528
-3.1160643879529160915,
1529
1.8333333333333333334,
1530
-0.3525641025641025641,
1531
0.1007326007326007326,
1532
-0.03021978021978021978,
1533
0.008461538461538461539,
1534
-0.0020739064856711915535,
1535
4.2324622156554929663E-4,
1536
-6.8220608081289195844E-5,
1537
8.085405402226867655E-6,
1538
-6.2373127388607264766E-7,
1539
2.3430926892790106976E-8
1541
static double fd_coeffs_22_3[] = {
1542
-3.0691770668705192515,
1543
2.7669911795655843276,
1544
-1.3117691135867156275,
1545
0.54233585056289137924,
1546
-0.1926733169277812135,
1547
0.05712488205266566611,
1548
-0.0136666647849245888466,
1549
0.0025253749603686947046,
1550
-3.3743255665248146452E-4,
1551
2.8966682770967642812E-5,
1552
-1.1983123155223574339E-6
1554
static double fd_coeffs_22_4[] = {
1555
16.144313191533628041,
1556
-12.276708267482077006,
1557
5.533982359131168655,
1558
-1.7490254847822875034,
1559
0.54233585056289137924,
1560
-0.1541386535422249708,
1561
0.038083254701777110742,
1562
-0.007809522734242622198,
1563
0.0012626874801843473523,
1564
-1.4997002517888065089E-4,
1565
1.1586673108387057124E-5,
1566
-4.3574993291722088505E-7
1568
static double fd_coeffs_22_5[] = {
1569
12.611227457118743494,
1570
-14.624570561606358564,
1571
9.2818974365489742475,
1572
-4.2008561512708338104,
1573
1.5534329340238318413,
1574
-0.4704908919346142152,
1575
0.113999611742071573226,
1576
-0.021238317584508838788,
1577
0.0028536829539616070896,
1578
-2.4594948927008139583E-4,
1579
1.0204527957115103084E-5
1581
static double fd_coeffs_22_6[] = {
1582
-91.066675998855397,
1583
75.66736474271246097,
1584
-43.87371168481907569,
1585
18.563794873097948495,
1586
-6.3012842269062507156,
1587
1.8641195208285982096,
1588
-0.4704908919346142152,
1589
0.09771395292177563419,
1590
-0.01592873818838162909,
1591
0.001902455302641071393,
1592
-1.475696935620488375E-4,
1593
5.566106158426419864E-6
1595
static double fd_coeffs_22_7[] = {
1596
-54.672951127001570674,
1597
71.18989538235969563,
1598
-53.005001320535174703,
1599
27.500961679352883056,
1600
-10.87510584759355014,
1601
3.4172331290430248763,
1602
-0.84660489699094058664,
1603
0.16001582643750853628,
1604
-0.02171410885603593937,
1605
0.0018847408985333368049,
1606
-7.8608332611033228316E-5
1608
static double fd_coeffs_22_8[] = {
1609
505.27603253545151693,
1610
-437.3836090160125654,
1611
284.75958152943878252,
1612
-141.34667018809379921,
1613
55.001923358705766112,
1614
-17.400169356149680223,
1615
4.5563108387240331683,
1616
-0.967548453703932099,
1617
0.16001582643750853628,
1618
-0.019301430094254168328,
1619
0.001507792718826669444,
1620
-5.7169696444387802412E-5
1622
static double fd_coeffs_22_9[] = {
1623
232.07478727336192614,
1624
-321.86160173963125353,
1625
263.0798688451925059,
1626
-151.09362420670158765,
1627
64.87647827625692209,
1628
-21.45870524700993451,
1629
5.4924785369273216497,
1630
-1.0609107700427144872,
1631
0.1461557132104007104,
1632
-0.012825027407740701392,
1633
5.3925953987957956213E-4
1635
static double fd_coeffs_22_10[] = {
1636
-2622.5949347392576558,
1637
2320.7478727336192614,
1638
-1609.3080086981562676,
1639
876.93289615064168635,
1640
-377.73406051675396913,
1641
129.75295655251384418,
1642
-35.764508745016557516,
1643
7.846397909896173785,
1644
-1.326138462553393109,
1645
0.16239523690044523378,
1646
-0.012825027407740701392,
1647
4.902359453450723292E-4
1649
static double fd_coeffs_23_1[] = {
1650
0.6234978120860432682,
1651
-0.058619452418345948292,
1652
0.015073573479003243847,
1653
-0.0046143592282662991366,
1654
0.0013957012480558559117,
1655
-3.8471687294490534856E-4,
1656
9.181605843458293328E-5,
1657
-1.8148437281221655819E-5,
1658
2.8258812375604654388E-6,
1659
-3.2318151074593372054E-7,
1660
2.4050407210736358095E-8,
1661
-8.7172101421342433797E-10
1663
static double fd_coeffs_23_2[] = {
1664
-0.26545948572732368774,
1665
0.38751026784657737975,
1666
-0.17679481466217356472,
1667
0.07703485910068456979,
1668
-0.030160551509394509063,
1669
0.010195583897778124449,
1670
-0.0028812810325856841546,
1671
6.579360950376451806E-4,
1672
-1.162009351500178396E-4,
1673
1.4861226457706625052E-5,
1674
-1.2228595530704828726E-6,
1675
4.8559645108218895124E-8
1677
static double fd_coeffs_23_3[] = {
1678
-0.7963784571819710632,
1679
0.38751026784657737975,
1680
-0.10607688879730413883,
1681
0.033014939614579101338,
1682
-0.010053517169798169687,
1683
0.0027806137903031248497,
1684
-6.6491100751977326646E-4,
1685
1.3158721900752903612E-4,
1686
-2.050604737941491287E-5,
1687
2.3465094406905197452E-6,
1688
-1.7469422186721183894E-7,
1689
6.3338667532459428423E-9
1691
static double fd_coeffs_23_4[] = {
1692
0.22933413448662658715,
1693
-0.44648188890609219215,
1694
0.327922308691383274,
1695
-0.15804167153276619284,
1696
0.06432907343496079987,
1697
-0.022166568820407843066,
1698
0.006332712310140038306,
1699
-0.0014558733040874587666,
1700
2.5828620576180822806E-4,
1701
-3.313674754232379681E-5,
1702
2.732890455864748474E-6,
1703
-1.0870843236168230352E-7
1705
static double fd_coeffs_23_5[] = {
1706
1.1466706724331329357,
1707
-0.7441364815101536536,
1708
0.327922308691383274,
1709
-0.11288690823769013774,
1710
0.035738374130533777708,
1711
-0.010075713100185383212,
1712
0.0024356585808230916562,
1713
-4.852911013624862555E-4,
1714
7.596653110641418473E-5,
1715
-8.720196721664157055E-6,
1716
6.506882037773210652E-7,
1717
-2.3632267904713544245E-8
1719
static double fd_coeffs_23_6[] = {
1720
-0.23755470117165160623,
1721
0.51924924450469119096,
1722
-0.46685988993819657042,
1723
0.2719683616682040959,
1724
-0.119568860640953564734,
1725
0.04281351442590151954,
1726
-0.012501566226441368318,
1727
0.0029135030757475728325,
1728
-5.215917271905421703E-4,
1729
6.7342934485472236145E-5,
1730
-5.5796169441759923757E-6,
1731
2.227123479763880854E-7
1733
static double fd_coeffs_23_7[] = {
1734
-1.6628829082015612436,
1735
1.2115815705109461123,
1736
-0.65360384591347519856,
1737
0.2719683616682040959,
1738
-0.092998002720741661456,
1739
0.027244963725573694253,
1740
-0.0067316125834684290943,
1741
0.0013596347686822006551,
1742
-2.1477306413728207013E-4,
1743
2.4810554810437139633E-5,
1744
-1.859872314725330792E-6,
1745
6.778201894933550425E-8
1747
static double fd_coeffs_23_8[] = {
1748
0.2561228525227838469,
1749
-0.59350164016281066866,
1750
0.59325634433501124,
1751
-0.3916172285777259866,
1752
0.19050552499210267167,
1753
-0.07221615593682341678,
1754
0.021814839054010380521,
1755
-0.005194963357928832976,
1756
9.436541876811625538E-4,
1757
-1.2308958371649278313E-4,
1758
1.0274968681389744099E-5,
1759
-4.1244126529359081215E-7
1761
static double fd_coeffs_23_9[] = {
1762
2.305105672705054622,
1763
-1.780504920488432006,
1764
1.067861419803020232,
1765
-0.50350786531421912564,
1766
0.19050552499210267167,
1767
-0.059085945766491886457,
1768
0.015102580883545648053,
1769
-0.0031169780147572997854,
1770
4.9958162877238017556E-4,
1771
-5.8305592286759739378E-5,
1772
4.403558006309890328E-6,
1773
-1.6139006033227466562E-7
1775
static double fd_coeffs_23_10[] = {
1776
-0.2680471562277438061,
1777
0.64216674493756485925,
1778
-0.6830699374935330556,
1779
0.4887570968795041861,
1780
-0.258125990636481507,
1781
0.10455996459487647732,
1782
-0.03302414465903053009,
1783
0.008103412716527648769,
1784
-0.0015028457583161025955,
1785
1.9896844149642379039E-4,
1786
-1.6792488877080352706E-5,
1787
6.796940124865087656E-7
1789
static double fd_coeffs_24_1[] = {
1790
0.92307692307692307694,
1791
-0.36263736263736263737,
1792
0.16117216117216117216,
1793
-0.06799450549450549451,
1794
0.025597931480284421461,
1795
-0.008295625942684766214,
1796
0.0022454325859898615317,
1797
-4.9118837818528221004E-4,
1798
8.3164169851476353024E-5,
1799
-1.02065117544993705984E-5,
1800
8.068388738734680315E-7,
1801
-3.081676254377829287E-8
1803
static double fd_coeffs_24_2[] = {
1804
-3.1299532768418049803,
1805
1.8461538461538461539,
1806
-0.36263736263736263737,
1807
0.10744810744810744811,
1808
-0.033997252747252747254,
1809
0.010239172592113768585,
1810
-0.0027652086475615887382,
1811
6.4155216742567472335E-4,
1812
-1.2279709454632055251E-4,
1813
1.8480926633661411784E-5,
1814
-2.0413023508998741198E-6,
1815
1.4669797706790327845E-7,
1816
-5.1361270906297154785E-9
1818
static double fd_coeffs_24_3[] = {
1819
-3.1291013820234599456,
1820
2.8611579605202054181,
1821
-1.4059358945413367182,
1822
0.6129609362788571972,
1823
-0.23421748499599640056,
1824
0.076512160484499420075,
1825
-0.020809346312442287675,
1826
0.004566141111088037227,
1827
-7.747395889494834336E-4,
1828
9.522532402808915328E-5,
1829
-7.5360953922905019134E-6,
1830
2.8808104894400656725E-7
1832
static double fd_coeffs_24_4[] = {
1833
16.403985223863037715,
1834
-12.516405528093839782,
1835
5.7223159210404108362,
1836
-1.8745811927217822909,
1837
0.6129609362788571972,
1838
-0.18737398799679712044,
1839
0.05100810698966628005,
1840
-0.011891055035681307242,
1841
0.0022830705555440186135,
1842
-3.4432870619977041494E-4,
1843
3.809012961123566131E-5,
1844
-2.740398324469273423E-6,
1845
9.602701631466885575E-8
1847
static double fd_coeffs_24_5[] = {
1848
13.128673392744821316,
1849
-15.43769988901876657,
1850
10.0950267639613822535,
1851
-4.810703146830139815,
1852
1.9121664608234236085,
1853
-0.63789987110775703993,
1854
0.1756766040690189297,
1855
-0.038860315392208083495,
1856
0.0066298253413257309555,
1857
-8.180922752343425876E-4,
1858
6.4931229223261825776E-5,
1859
-2.487577330279396486E-6
1861
static double fd_coeffs_24_6[] = {
1862
-94.43007458042490285,
1863
78.7720403564689279,
1864
-46.31309966705629971,
1865
20.190053527922764507,
1866
-7.2160547202452097223,
1867
2.2945997529881083303,
1868
-0.63789987110775703993,
1869
0.15057994634487336832,
1870
-0.02914523654415606262,
1871
0.004419883560883820637,
1872
-4.9085536514060555256E-4,
1873
3.5417034121779177698E-5,
1874
-1.243788665139698243E-6
1876
static double fd_coeffs_24_7[] = {
1877
-58.75927633207841541,
1878
77.61126356176616593,
1879
-59.426369499941645004,
1880
32.316987813907735783,
1881
-13.708062397331698803,
1882
4.7392795189208275857,
1883
-1.3336746195774994794,
1884
0.29917860431938250566,
1885
-0.051534704116437504237,
1886
0.0064030129076850890578,
1887
-5.1079087261685300903E-4,
1888
1.9644660909355444577E-5
1890
static double fd_coeffs_24_8[] = {
1891
540.690850979450838,
1892
-470.0742106566273233,
1893
310.44505424706466373,
1894
-158.47031866651105335,
1895
64.633975627815471565,
1896
-21.932899835730718083,
1897
6.3190393585611034474,
1898
-1.5241995652314279765,
1899
0.29917860431938250566,
1900
-0.045808625881277781544,
1901
0.005122410326148071246,
1902
-3.7148427099407491565E-4,
1903
1.3096440606236963051E-5
1905
static double fd_coeffs_24_9[] = {
1906
261.2253276119456675,
1907
-367.6695937002628471,
1908
308.88786080582409949,
1909
-185.44961817717528282,
1910
85.08588649418262514,
1911
-30.889762415375262598,
1912
8.967078546325074103,
1913
-2.0536536298706437596,
1914
0.35888632603067126877,
1915
-0.045056938441115028417,
1916
0.0036223118995936456255,
1917
-1.4013874362336663924E-4
1919
static double fd_coeffs_24_10[] = {
1920
-2938.392455073914854,
1921
2612.253276119456675,
1922
-1838.3479685013142354,
1923
1029.6262026860803316,
1924
-463.62404544293820707,
1925
170.17177298836525028,
1926
-51.482937358958770997,
1927
12.810112209035820147,
1928
-2.5670670373383046995,
1929
0.39876258447852363197,
1930
-0.045056938441115028417,
1931
0.0032930108178124051139,
1932
-1.167822863528055327E-4
1934
static double fd_coeffs_25_1[] = {
1935
0.6244970072977196196,
1936
-0.059475905456925678056,
1937
0.015701639040628379007,
1938
-0.0050069002042820086118,
1939
0.0016035170588877021045,
1940
-4.7707945553683698978E-4,
1941
1.2584437833687353795E-4,
1942
-2.8356933251908837216E-5,
1943
5.2564755162955086288E-6,
1944
-7.6510774324321430053E-7,
1945
8.169295927559904244E-8,
1946
-5.6752670196186480336E-9,
1947
1.9214184021620894783E-10
1949
static double fd_coeffs_25_2[] = {
1950
-0.2678832919931881272,
1951
0.3937429125302287955,
1952
-0.18441249149774751732,
1953
0.083700326331811778306,
1954
-0.034697550212934877888,
1955
0.012660126403404991465,
1956
-0.003954359348432788932,
1957
0.0010293862812924122189,
1958
-2.1643352509178037374E-4,
1959
3.5229346338920402042E-5,
1960
-4.1592246618038191433E-6,
1961
3.1656122249261069763E-7,
1962
-1.165224249497355663E-8
1964
static double fd_coeffs_25_3[] = {
1965
-0.8036498759795643816,
1966
0.3937429125302287955,
1967
-0.110647494898648510394,
1968
0.03587156842791933356,
1969
-0.0115658500709782926296,
1970
0.0034527617463831794902,
1971
-9.12544465022951292E-4,
1972
2.0587725625848244376E-4,
1973
-3.8194151486784771836E-5,
1974
5.562528369303221375E-6,
1975
-5.9417495168625987764E-7,
1976
4.1290594238166612735E-8,
1977
-1.3982690993968267957E-9
1979
static double fd_coeffs_25_4[] = {
1980
0.23480664724818830324,
1981
-0.46055406457867946213,
1982
0.34512163451343438173,
1983
-0.17309108162706091208,
1984
0.07457278954956476844,
1985
-0.027731056586365554385,
1986
0.008755527557518754717,
1987
-0.0022945401204877836782,
1988
4.8459312447300701375E-4,
1989
-7.912424974032141101E-5,
1990
9.362667660610399949E-6,
1991
-7.138071454932298588E-7,
1992
2.6308639701371632839E-8
1994
static double fd_coeffs_25_5[] = {
1995
1.1740332362409415162,
1996
-0.76759010763113243685,
1997
0.34512163451343438173,
1998
-0.12363648687647208006,
1999
0.041429327527535982465,
2000
-0.012605025721075251993,
2001
0.0033675105990456748912,
2002
-7.648467068292612261E-4,
2003
1.4252738955088441581E-4,
2004
-2.082217098429510816E-5,
2005
2.2292065858596190354E-6,
2006
-1.5517546641157170843E-7,
2007
5.2617279402743265677E-9
2009
static double fd_coeffs_25_6[] = {
2010
-0.2489610773507709159,
2011
0.54857992610814084437,
2012
-0.50270850078685725787,
2013
0.30333589616078219741,
2014
-0.14091987151405294477,
2015
0.054411594406350565486,
2016
-0.017551447940416790238,
2017
0.004661539053662141958,
2018
-9.932839751992354264E-4,
2019
1.6319483515034038446E-4,
2020
-1.9398083172749661354E-5,
2021
1.4839215672509689842E-6,
2022
-5.4835183446720908645E-8
2024
static double fd_coeffs_25_7[] = {
2025
-1.7427275414553964112,
2026
1.2800198275856619701,
2027
-0.70379190110160016103,
2028
0.30333589616078219741,
2029
-0.10960434451093006815,
2030
0.034625560076768541672,
2031
-0.009450779660224425512,
2032
0.0021753848917089995806,
2033
-4.0899928390556752854E-4,
2034
6.01244129501254048E-5,
2035
-6.466027724249887118E-6,
2036
4.5162830307638186475E-7,
2037
-1.535385136508185442E-8
2039
static double fd_coeffs_25_8[] = {
2040
0.27785231745569293952,
2041
-0.64937740713314833543,
2042
0.66154894840986838825,
2043
-0.4513732571432259913,
2044
0.23117979653668670848,
2045
-0.094310821961041905917,
2046
0.031435004643502641343,
2047
-0.008525020677368461721,
2048
0.0018422410834029671359,
2049
-3.0569012937119104046E-4,
2050
3.659953247142633497E-5,
2051
-2.8150800239080415665E-6,
2052
1.044625547223674241E-7
2054
static double fd_coeffs_25_9[] = {
2055
2.5006708571012364556,
2056
-1.9481322213994450062,
2057
1.1907881071377630989,
2058
-0.5803370448984334174,
2059
0.23117979653668670848,
2060
-0.07716339978630701393,
2061
0.021762695522424905545,
2062
-0.0051150124064210770325,
2063
9.753041029780414249E-4,
2064
-1.4480058759687996654E-4,
2065
1.568551391632557213E-5,
2066
-1.1015530528335814825E-6,
2067
3.7606519700052272676E-8
2069
static double fd_coeffs_25_10[] = {
2070
-0.30541751501504875807,
2071
0.73826195324777759283,
2072
-0.8005196365393486189,
2073
0.591525583544592804,
2074
-0.32807764962700401162,
2075
0.14255839663911092427,
2076
-0.049568916506041700774,
2077
0.013830449125108438621,
2078
-0.0030482365352347284287,
2079
5.130050699611712859E-4,
2080
-6.2065503964675529406E-5,
2081
4.811754913338449496E-6,
2082
-1.796548217761713361E-7
2084
static void apply_periodic_fd(int n,double *x,double *y,double coeff,int m,double *c) {
2087
for(i=m;i<n-m;++i) {
2088
double d = 0.0, *c_ptr = c + l;
2090
d += (*(--c_ptr)) * (x[i+k]-x[i-k]);
2094
double d = 0.0, *c_ptr = c + l;
2096
d += (*(--c_ptr)) * (x[i+k]-x[(i+n-k)%n]);
2099
for(i=n-m;i<n;++i) {
2100
double d = 0.0, *c_ptr = c + l;
2102
d += (*(--c_ptr)) * (x[(i+k)%n]-x[i-k]);
2106
for(i=l;i<n-l;++i) {
2107
double d = 0.0, *c_ptr = c + l;
2109
d += (*(--c_ptr)) * (x[i+k]-x[i-k]);
2113
double d = 0.0, *c_ptr = c + l;
2115
d += (*(--c_ptr)) * (x[i+k]-x[(i+n-k)%n]);
2118
for(i=n-l;i<n;++i) {
2119
double d = 0.0, *c_ptr = c + l;
2121
d += (*(--c_ptr)) * (x[(i+k)%n]-x[i-k]);
2126
static void apply_periodic_fd2(int n,double *x,double *y,double coeff,int m,double *c) {
2129
for(i=m;i<n-m;++i) {
2130
double d = 0.0, *c_ptr = c + l;
2132
d += (*(--c_ptr)) * (x[i+k]+x[i-k]);
2136
double d = 0.0, *c_ptr = c + l;
2138
d += (*(--c_ptr)) * (x[i+k]+x[(i+n-k)%n]);
2141
for(i=n-m;i<n;++i) {
2142
double d = 0.0, *c_ptr = c + l;
2144
d += (*(--c_ptr)) * (x[(i+k)%n]+x[i-k]);
2148
for(i=l;i<n-l;++i) {
2149
double d = 0.0, *c_ptr = c + l + 1;
2151
d += (*(--c_ptr)) * (x[i+k]+x[i-k]);
2152
y[i] = coeff * (d+(*c)*x[i]);
2155
double d = 0.0, *c_ptr = c + l + 1;
2157
d += (*(--c_ptr)) * (x[i+k]+x[(i+n-k)%n]);
2158
y[i] = coeff * (d+(*c)*x[i]);
2160
for(i=n-l;i<n;++i) {
2161
double d = 0.0, *c_ptr = c + l + 1;
2163
d += (*(--c_ptr)) * (x[(i+k)%n]+x[i-k]);
2164
y[i] = coeff * (d+(*c)*x[i]);
2168
extern void periodic_finite_difference(
2169
int n, double *x, double *y, double h, int k, int m) {
2170
double coeff = pow(1.0/h,k), *c=NULL;
2172
case 1001: c=fd_coeffs_1_1; break;
2173
case 1002: c=fd_coeffs_2_1; break;
2174
case 2002: c=fd_coeffs_2_2; break;
2175
case 1003: c=fd_coeffs_3_1; break;
2176
case 2003: c=fd_coeffs_3_2; break;
2177
case 3003: c=fd_coeffs_3_3; break;
2178
case 1004: c=fd_coeffs_4_1; break;
2179
case 2004: c=fd_coeffs_4_2; break;
2180
case 3004: c=fd_coeffs_4_3; break;
2181
case 4004: c=fd_coeffs_4_4; break;
2182
case 1005: c=fd_coeffs_5_1; break;
2183
case 2005: c=fd_coeffs_5_2; break;
2184
case 3005: c=fd_coeffs_5_3; break;
2185
case 4005: c=fd_coeffs_5_4; break;
2186
case 5005: c=fd_coeffs_5_5; break;
2187
case 1006: c=fd_coeffs_6_1; break;
2188
case 2006: c=fd_coeffs_6_2; break;
2189
case 3006: c=fd_coeffs_6_3; break;
2190
case 4006: c=fd_coeffs_6_4; break;
2191
case 5006: c=fd_coeffs_6_5; break;
2192
case 6006: c=fd_coeffs_6_6; break;
2193
case 1007: c=fd_coeffs_7_1; break;
2194
case 2007: c=fd_coeffs_7_2; break;
2195
case 3007: c=fd_coeffs_7_3; break;
2196
case 4007: c=fd_coeffs_7_4; break;
2197
case 5007: c=fd_coeffs_7_5; break;
2198
case 6007: c=fd_coeffs_7_6; break;
2199
case 7007: c=fd_coeffs_7_7; break;
2200
case 1008: c=fd_coeffs_8_1; break;
2201
case 2008: c=fd_coeffs_8_2; break;
2202
case 3008: c=fd_coeffs_8_3; break;
2203
case 4008: c=fd_coeffs_8_4; break;
2204
case 5008: c=fd_coeffs_8_5; break;
2205
case 6008: c=fd_coeffs_8_6; break;
2206
case 7008: c=fd_coeffs_8_7; break;
2207
case 8008: c=fd_coeffs_8_8; break;
2208
case 1009: c=fd_coeffs_9_1; break;
2209
case 2009: c=fd_coeffs_9_2; break;
2210
case 3009: c=fd_coeffs_9_3; break;
2211
case 4009: c=fd_coeffs_9_4; break;
2212
case 5009: c=fd_coeffs_9_5; break;
2213
case 6009: c=fd_coeffs_9_6; break;
2214
case 7009: c=fd_coeffs_9_7; break;
2215
case 8009: c=fd_coeffs_9_8; break;
2216
case 9009: c=fd_coeffs_9_9; break;
2217
case 1010: c=fd_coeffs_10_1; break;
2218
case 2010: c=fd_coeffs_10_2; break;
2219
case 3010: c=fd_coeffs_10_3; break;
2220
case 4010: c=fd_coeffs_10_4; break;
2221
case 5010: c=fd_coeffs_10_5; break;
2222
case 6010: c=fd_coeffs_10_6; break;
2223
case 7010: c=fd_coeffs_10_7; break;
2224
case 8010: c=fd_coeffs_10_8; break;
2225
case 9010: c=fd_coeffs_10_9; break;
2226
case 10010: c=fd_coeffs_10_10; break;
2227
case 1011: c=fd_coeffs_11_1; break;
2228
case 2011: c=fd_coeffs_11_2; break;
2229
case 3011: c=fd_coeffs_11_3; break;
2230
case 4011: c=fd_coeffs_11_4; break;
2231
case 5011: c=fd_coeffs_11_5; break;
2232
case 6011: c=fd_coeffs_11_6; break;
2233
case 7011: c=fd_coeffs_11_7; break;
2234
case 8011: c=fd_coeffs_11_8; break;
2235
case 9011: c=fd_coeffs_11_9; break;
2236
case 10011: c=fd_coeffs_11_10; break;
2237
case 1012: c=fd_coeffs_12_1; break;
2238
case 2012: c=fd_coeffs_12_2; break;
2239
case 3012: c=fd_coeffs_12_3; break;
2240
case 4012: c=fd_coeffs_12_4; break;
2241
case 5012: c=fd_coeffs_12_5; break;
2242
case 6012: c=fd_coeffs_12_6; break;
2243
case 7012: c=fd_coeffs_12_7; break;
2244
case 8012: c=fd_coeffs_12_8; break;
2245
case 9012: c=fd_coeffs_12_9; break;
2246
case 10012: c=fd_coeffs_12_10; break;
2247
case 1013: c=fd_coeffs_13_1; break;
2248
case 2013: c=fd_coeffs_13_2; break;
2249
case 3013: c=fd_coeffs_13_3; break;
2250
case 4013: c=fd_coeffs_13_4; break;
2251
case 5013: c=fd_coeffs_13_5; break;
2252
case 6013: c=fd_coeffs_13_6; break;
2253
case 7013: c=fd_coeffs_13_7; break;
2254
case 8013: c=fd_coeffs_13_8; break;
2255
case 9013: c=fd_coeffs_13_9; break;
2256
case 10013: c=fd_coeffs_13_10; break;
2257
case 1014: c=fd_coeffs_14_1; break;
2258
case 2014: c=fd_coeffs_14_2; break;
2259
case 3014: c=fd_coeffs_14_3; break;
2260
case 4014: c=fd_coeffs_14_4; break;
2261
case 5014: c=fd_coeffs_14_5; break;
2262
case 6014: c=fd_coeffs_14_6; break;
2263
case 7014: c=fd_coeffs_14_7; break;
2264
case 8014: c=fd_coeffs_14_8; break;
2265
case 9014: c=fd_coeffs_14_9; break;
2266
case 10014: c=fd_coeffs_14_10; break;
2267
case 1015: c=fd_coeffs_15_1; break;
2268
case 2015: c=fd_coeffs_15_2; break;
2269
case 3015: c=fd_coeffs_15_3; break;
2270
case 4015: c=fd_coeffs_15_4; break;
2271
case 5015: c=fd_coeffs_15_5; break;
2272
case 6015: c=fd_coeffs_15_6; break;
2273
case 7015: c=fd_coeffs_15_7; break;
2274
case 8015: c=fd_coeffs_15_8; break;
2275
case 9015: c=fd_coeffs_15_9; break;
2276
case 10015: c=fd_coeffs_15_10; break;
2277
case 1016: c=fd_coeffs_16_1; break;
2278
case 2016: c=fd_coeffs_16_2; break;
2279
case 3016: c=fd_coeffs_16_3; break;
2280
case 4016: c=fd_coeffs_16_4; break;
2281
case 5016: c=fd_coeffs_16_5; break;
2282
case 6016: c=fd_coeffs_16_6; break;
2283
case 7016: c=fd_coeffs_16_7; break;
2284
case 8016: c=fd_coeffs_16_8; break;
2285
case 9016: c=fd_coeffs_16_9; break;
2286
case 10016: c=fd_coeffs_16_10; break;
2287
case 1017: c=fd_coeffs_17_1; break;
2288
case 2017: c=fd_coeffs_17_2; break;
2289
case 3017: c=fd_coeffs_17_3; break;
2290
case 4017: c=fd_coeffs_17_4; break;
2291
case 5017: c=fd_coeffs_17_5; break;
2292
case 6017: c=fd_coeffs_17_6; break;
2293
case 7017: c=fd_coeffs_17_7; break;
2294
case 8017: c=fd_coeffs_17_8; break;
2295
case 9017: c=fd_coeffs_17_9; break;
2296
case 10017: c=fd_coeffs_17_10; break;
2297
case 1018: c=fd_coeffs_18_1; break;
2298
case 2018: c=fd_coeffs_18_2; break;
2299
case 3018: c=fd_coeffs_18_3; break;
2300
case 4018: c=fd_coeffs_18_4; break;
2301
case 5018: c=fd_coeffs_18_5; break;
2302
case 6018: c=fd_coeffs_18_6; break;
2303
case 7018: c=fd_coeffs_18_7; break;
2304
case 8018: c=fd_coeffs_18_8; break;
2305
case 9018: c=fd_coeffs_18_9; break;
2306
case 10018: c=fd_coeffs_18_10; break;
2307
case 1019: c=fd_coeffs_19_1; break;
2308
case 2019: c=fd_coeffs_19_2; break;
2309
case 3019: c=fd_coeffs_19_3; break;
2310
case 4019: c=fd_coeffs_19_4; break;
2311
case 5019: c=fd_coeffs_19_5; break;
2312
case 6019: c=fd_coeffs_19_6; break;
2313
case 7019: c=fd_coeffs_19_7; break;
2314
case 8019: c=fd_coeffs_19_8; break;
2315
case 9019: c=fd_coeffs_19_9; break;
2316
case 10019: c=fd_coeffs_19_10; break;
2317
case 1020: c=fd_coeffs_20_1; break;
2318
case 2020: c=fd_coeffs_20_2; break;
2319
case 3020: c=fd_coeffs_20_3; break;
2320
case 4020: c=fd_coeffs_20_4; break;
2321
case 5020: c=fd_coeffs_20_5; break;
2322
case 6020: c=fd_coeffs_20_6; break;
2323
case 7020: c=fd_coeffs_20_7; break;
2324
case 8020: c=fd_coeffs_20_8; break;
2325
case 9020: c=fd_coeffs_20_9; break;
2326
case 10020: c=fd_coeffs_20_10; break;
2327
case 1021: c=fd_coeffs_21_1; break;
2328
case 2021: c=fd_coeffs_21_2; break;
2329
case 3021: c=fd_coeffs_21_3; break;
2330
case 4021: c=fd_coeffs_21_4; break;
2331
case 5021: c=fd_coeffs_21_5; break;
2332
case 6021: c=fd_coeffs_21_6; break;
2333
case 7021: c=fd_coeffs_21_7; break;
2334
case 8021: c=fd_coeffs_21_8; break;
2335
case 9021: c=fd_coeffs_21_9; break;
2336
case 10021: c=fd_coeffs_21_10; break;
2337
case 1022: c=fd_coeffs_22_1; break;
2338
case 2022: c=fd_coeffs_22_2; break;
2339
case 3022: c=fd_coeffs_22_3; break;
2340
case 4022: c=fd_coeffs_22_4; break;
2341
case 5022: c=fd_coeffs_22_5; break;
2342
case 6022: c=fd_coeffs_22_6; break;
2343
case 7022: c=fd_coeffs_22_7; break;
2344
case 8022: c=fd_coeffs_22_8; break;
2345
case 9022: c=fd_coeffs_22_9; break;
2346
case 10022: c=fd_coeffs_22_10; break;
2347
case 1023: c=fd_coeffs_23_1; break;
2348
case 2023: c=fd_coeffs_23_2; break;
2349
case 3023: c=fd_coeffs_23_3; break;
2350
case 4023: c=fd_coeffs_23_4; break;
2351
case 5023: c=fd_coeffs_23_5; break;
2352
case 6023: c=fd_coeffs_23_6; break;
2353
case 7023: c=fd_coeffs_23_7; break;
2354
case 8023: c=fd_coeffs_23_8; break;
2355
case 9023: c=fd_coeffs_23_9; break;
2356
case 10023: c=fd_coeffs_23_10; break;
2357
case 1024: c=fd_coeffs_24_1; break;
2358
case 2024: c=fd_coeffs_24_2; break;
2359
case 3024: c=fd_coeffs_24_3; break;
2360
case 4024: c=fd_coeffs_24_4; break;
2361
case 5024: c=fd_coeffs_24_5; break;
2362
case 6024: c=fd_coeffs_24_6; break;
2363
case 7024: c=fd_coeffs_24_7; break;
2364
case 8024: c=fd_coeffs_24_8; break;
2365
case 9024: c=fd_coeffs_24_9; break;
2366
case 10024: c=fd_coeffs_24_10; break;
2367
case 1025: c=fd_coeffs_25_1; break;
2368
case 2025: c=fd_coeffs_25_2; break;
2369
case 3025: c=fd_coeffs_25_3; break;
2370
case 4025: c=fd_coeffs_25_4; break;
2371
case 5025: c=fd_coeffs_25_5; break;
2372
case 6025: c=fd_coeffs_25_6; break;
2373
case 7025: c=fd_coeffs_25_7; break;
2374
case 8025: c=fd_coeffs_25_8; break;
2375
case 9025: c=fd_coeffs_25_9; break;
2376
case 10025: c=fd_coeffs_25_10; break;
2379
"Unsupported derivative/order combination: k,m=%d,%d\n",
2384
apply_periodic_fd(n,x,y,coeff,m,c);
2386
apply_periodic_fd2(n,x,y,coeff,m,c);