~ubuntu-branches/ubuntu/saucy/python-scipy/saucy

« back to all changes in this revision

Viewing changes to Lib/sandbox/fdfpack/src/pfdf.c

  • Committer: Bazaar Package Importer
  • Author(s): Ondrej Certik
  • Date: 2008-06-16 22:58:01 UTC
  • mfrom: (2.1.24 intrepid)
  • Revision ID: james.westby@ubuntu.com-20080616225801-irdhrpcwiocfbcmt
Tags: 0.6.0-12
* The description updated to match the current SciPy (Closes: #489149).
* Standards-Version bumped to 3.8.0 (no action needed)
* Build-Depends: netcdf-dev changed to libnetcdf-dev

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
/*****************************************************/
2
 
/******** Periodic Finite Difference Formulae ********/
3
 
/*****************************************************/
4
 
/* Author: Pearu Peterson, September 2002            */
5
 
/*                                                   */
6
 
/* This file is automatically generated using        */
7
 
/* genpfdf.cc utility from fdf package and it        */
8
 
/* contains the following C function                 */
9
 
/*                                                   */
10
 
/*  void periodic_finite_difference(int n,           */
11
 
/*    double *x, double *y, double h, int k, int m); */
12
 
/* where                                             */
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    */
21
 
/*       accuracy.                                   */
22
 
/*****************************************************/
23
 
/* Reference:                                        */
24
 
/*   http://epubs.siam.org/sam-bin/dbq/article/32250 */
25
 
/*****************************************************/
26
 
/* Parameters to genpfdf n=25 m=10 */
27
 
#include <stdio.h>
28
 
#include <math.h>
29
 
static double fd_coeffs_1_1[] = {
30
 
 0.5
31
 
};
32
 
static double fd_coeffs_2_1[] = {
33
 
 0.5
34
 
};
35
 
static double fd_coeffs_2_2[] = {
36
 
 -2.0,
37
 
 1.0
38
 
};
39
 
static double fd_coeffs_3_1[] = {
40
 
 0.5625,
41
 
 -0.020833333333333333334
42
 
};
43
 
static double fd_coeffs_3_2[] = {
44
 
 -0.125,
45
 
 0.125
46
 
};
47
 
static double fd_coeffs_3_3[] = {
48
 
 -0.375,
49
 
 0.125
50
 
};
51
 
static double fd_coeffs_4_1[] = {
52
 
 0.6666666666666666667,
53
 
 -0.083333333333333333336
54
 
};
55
 
static double fd_coeffs_4_2[] = {
56
 
 -2.5,
57
 
 1.3333333333333333334,
58
 
 -0.083333333333333333336
59
 
};
60
 
static double fd_coeffs_4_3[] = {
61
 
 -1.0,
62
 
 0.5
63
 
};
64
 
static double fd_coeffs_4_4[] = {
65
 
 6.0,
66
 
 -4.0,
67
 
 1.0
68
 
};
69
 
static double fd_coeffs_5_1[] = {
70
 
 0.5859375,
71
 
 -0.032552083333333333332,
72
 
 0.00234375
73
 
};
74
 
static double fd_coeffs_5_2[] = {
75
 
 -0.17708333333333333333,
76
 
 0.203125,
77
 
 -0.026041666666666666666
78
 
};
79
 
static double fd_coeffs_5_3[] = {
80
 
 -0.53125,
81
 
 0.203125,
82
 
 -0.015625
83
 
};
84
 
static double fd_coeffs_5_4[] = {
85
 
 0.0625,
86
 
 -0.09375,
87
 
 0.03125
88
 
};
89
 
static double fd_coeffs_5_5[] = {
90
 
 0.3125,
91
 
 -0.15625,
92
 
 0.03125
93
 
};
94
 
static double fd_coeffs_6_1[] = {
95
 
 0.75,
96
 
 -0.15,
97
 
 0.016666666666666666667
98
 
};
99
 
static double fd_coeffs_6_2[] = {
100
 
 -2.7222222222222222222,
101
 
 1.5,
102
 
 -0.15,
103
 
 0.011111111111111111111
104
 
};
105
 
static double fd_coeffs_6_3[] = {
106
 
 -1.625,
107
 
 1.0,
108
 
 -0.125
109
 
};
110
 
static double fd_coeffs_6_4[] = {
111
 
 9.333333333333333333,
112
 
 -6.5,
113
 
 2.0,
114
 
 -0.16666666666666666667
115
 
};
116
 
static double fd_coeffs_6_5[] = {
117
 
 2.5,
118
 
 -2.0,
119
 
 0.5
120
 
};
121
 
static double fd_coeffs_6_6[] = {
122
 
 -20.0,
123
 
 15.0,
124
 
 -6.0,
125
 
 1.0
126
 
};
127
 
static double fd_coeffs_7_1[] = {
128
 
 0.59814453125,
129
 
 -0.039876302083333333332,
130
 
 0.00478515625,
131
 
 -3.4877232142857142856E-4
132
 
};
133
 
static double fd_coeffs_7_2[] = {
134
 
 -0.20518663194444444444,
135
 
 0.2537109375,
136
 
 -0.054144965277777777777,
137
 
 0.0056206597222222222224
138
 
};
139
 
static double fd_coeffs_7_3[] = {
140
 
 -0.6155598958333333333,
141
 
 0.2537109375,
142
 
 -0.032486979166666666666,
143
 
 0.0024088541666666666667
144
 
};
145
 
static double fd_coeffs_7_4[] = {
146
 
 0.108072916666666666664,
147
 
 -0.17578125,
148
 
 0.076822916666666666664,
149
 
 -0.009114583333333333333
150
 
};
151
 
static double fd_coeffs_7_5[] = {
152
 
 0.5403645833333333333,
153
 
 -0.29296875,
154
 
 0.076822916666666666664,
155
 
 -0.0065104166666666666665
156
 
};
157
 
static double fd_coeffs_7_6[] = {
158
 
 -0.0390625,
159
 
 0.0703125,
160
 
 -0.0390625,
161
 
 0.0078125
162
 
};
163
 
static double fd_coeffs_7_7[] = {
164
 
 -0.2734375,
165
 
 0.1640625,
166
 
 -0.0546875,
167
 
 0.0078125
168
 
};
169
 
static double fd_coeffs_8_1[] = {
170
 
 0.8,
171
 
 -0.2,
172
 
 0.03809523809523809524,
173
 
 -0.0035714285714285714287
174
 
};
175
 
static double fd_coeffs_8_2[] = {
176
 
 -2.8472222222222222222,
177
 
 1.6,
178
 
 -0.2,
179
 
 0.025396825396825396825,
180
 
 -0.0017857142857142857143
181
 
};
182
 
static double fd_coeffs_8_3[] = {
183
 
 -2.0333333333333333333,
184
 
 1.4083333333333333333,
185
 
 -0.3,
186
 
 0.029166666666666666667
187
 
};
188
 
static double fd_coeffs_8_4[] = {
189
 
 11.375,
190
 
 -8.133333333333333333,
191
 
 2.8166666666666666666,
192
 
 -0.4,
193
 
 0.029166666666666666667
194
 
};
195
 
static double fd_coeffs_8_5[] = {
196
 
 4.8333333333333333335,
197
 
 -4.3333333333333333335,
198
 
 1.5,
199
 
 -0.16666666666666666667
200
 
};
201
 
static double fd_coeffs_8_6[] = {
202
 
 -37.5,
203
 
 29.0,
204
 
 -13.0,
205
 
 3.0,
206
 
 -0.25
207
 
};
208
 
static double fd_coeffs_8_7[] = {
209
 
 -7.0,
210
 
 7.0,
211
 
 -3.0,
212
 
 0.5
213
 
};
214
 
static double fd_coeffs_8_8[] = {
215
 
 70.0,
216
 
 -56.0,
217
 
 28.0,
218
 
 -8.0,
219
 
 1.0
220
 
};
221
 
static double fd_coeffs_9_1[] = {
222
 
 0.605621337890625,
223
 
 -0.04486083984375,
224
 
 0.00692138671875,
225
 
 -8.828299386160714286E-4,
226
 
 5.9339735243055555557E-5
227
 
};
228
 
static double fd_coeffs_9_2[] = {
229
 
 -0.222705078125,
230
 
 0.28874782986111111112,
231
 
 -0.07917131696428571429,
232
 
 0.0143798828125,
233
 
 -0.0012513175843253968254
234
 
};
235
 
static double fd_coeffs_9_3[] = {
236
 
 -0.668115234375,
237
 
 0.28874782986111111112,
238
 
 -0.04750279017857142857,
239
 
 0.0061628069196428571427,
240
 
 -4.171058614417989418E-4
241
 
};
242
 
static double fd_coeffs_9_4[] = {
243
 
 0.14020182291666666667,
244
 
 -0.2400390625,
245
 
 0.122721354166666666664,
246
 
 -0.025179036458333333333,
247
 
 0.002294921875
248
 
};
249
 
static double fd_coeffs_9_5[] = {
250
 
 0.7010091145833333333,
251
 
 -0.40006510416666666666,
252
 
 0.122721354166666666664,
253
 
 -0.017985026041666666666,
254
 
 0.0012749565972222222222
255
 
};
256
 
static double fd_coeffs_9_6[] = {
257
 
 -0.080078125,
258
 
 0.15234375,
259
 
 -0.09765625,
260
 
 0.0283203125,
261
 
 -0.0029296875
262
 
};
263
 
static double fd_coeffs_9_7[] = {
264
 
 -0.560546875,
265
 
 0.35546875,
266
 
 -0.13671875,
267
 
 0.0283203125,
268
 
 -0.0022786458333333333333
269
 
};
270
 
static double fd_coeffs_9_8[] = {
271
 
 0.02734375,
272
 
 -0.0546875,
273
 
 0.0390625,
274
 
 -0.013671875,
275
 
 0.001953125
276
 
};
277
 
static double fd_coeffs_9_9[] = {
278
 
 0.24609375,
279
 
 -0.1640625,
280
 
 0.0703125,
281
 
 -0.017578125,
282
 
 0.001953125
283
 
};
284
 
static double fd_coeffs_10_1[] = {
285
 
 0.8333333333333333333,
286
 
 -0.23809523809523809523,
287
 
 0.05952380952380952381,
288
 
 -0.009920634920634920635,
289
 
 7.936507936507936508E-4
290
 
};
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
298
 
};
299
 
static double fd_coeffs_10_3[] = {
300
 
 -2.3180555555555555556,
301
 
 1.7337301587301587302,
302
 
 -0.4830357142857142857,
303
 
 0.08339947089947089947,
304
 
 -0.0067791005291005291007
305
 
};
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
313
 
};
314
 
static double fd_coeffs_10_5[] = {
315
 
 6.7291666666666666665,
316
 
 -6.5,
317
 
 2.71875,
318
 
 -0.52777777777777777775,
319
 
 0.04513888888888888889
320
 
};
321
 
static double fd_coeffs_10_6[] = {
322
 
 -51.15,
323
 
 40.375,
324
 
 -19.5,
325
 
 5.4375,
326
 
 -0.7916666666666666667,
327
 
 0.054166666666666666667
328
 
};
329
 
static double fd_coeffs_10_7[] = {
330
 
 -15.75,
331
 
 17.0,
332
 
 -8.625,
333
 
 2.1666666666666666667,
334
 
 -0.20833333333333333333
335
 
};
336
 
static double fd_coeffs_10_8[] = {
337
 
 154.0,
338
 
 -126.0,
339
 
 68.0,
340
 
 -23.0,
341
 
 4.3333333333333333335,
342
 
 -0.33333333333333333334
343
 
};
344
 
static double fd_coeffs_10_9[] = {
345
 
 21.0,
346
 
 -24.0,
347
 
 13.5,
348
 
 -4.0,
349
 
 0.5
350
 
};
351
 
static double fd_coeffs_10_10[] = {
352
 
 -252.0,
353
 
 210.0,
354
 
 -120.0,
355
 
 45.0,
356
 
 -10.0,
357
 
 1.0
358
 
};
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
366
 
};
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
374
 
};
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
382
 
};
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
390
 
};
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
398
 
};
399
 
static double fd_coeffs_11_6[] = {
400
 
 -0.1157958984375,
401
 
 0.2288818359375,
402
 
 -0.16143798828125,
403
 
 0.058085123697916666668,
404
 
 -0.01058349609375,
405
 
 8.5042317708333333334E-4
406
 
};
407
 
static double fd_coeffs_11_7[] = {
408
 
 -0.8105712890625,
409
 
 0.5340576171875,
410
 
 -0.22601318359375,
411
 
 0.058085123697916666668,
412
 
 -0.008231608072916666667,
413
 
 5.411783854166666667E-4
414
 
};
415
 
static double fd_coeffs_11_8[] = {
416
 
 0.06494140625,
417
 
 -0.13525390625,
418
 
 0.106201171875,
419
 
 -0.045003255208333333332,
420
 
 0.010009765625,
421
 
 -8.951822916666666667E-4
422
 
};
423
 
static double fd_coeffs_11_9[] = {
424
 
 0.58447265625,
425
 
 -0.40576171875,
426
 
 0.191162109375,
427
 
 -0.057861328125,
428
 
 0.010009765625,
429
 
 -7.32421875E-4
430
 
};
431
 
static double fd_coeffs_11_10[] = {
432
 
 -0.0205078125,
433
 
 0.0439453125,
434
 
 -0.03662109375,
435
 
 0.01708984375,
436
 
 -0.00439453125,
437
 
 4.8828125E-4
438
 
};
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
446
 
};
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
455
 
};
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
463
 
};
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
472
 
};
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
480
 
};
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
489
 
};
490
 
static double fd_coeffs_12_7[] = {
491
 
 -24.275,
492
 
 27.65625,
493
 
 -15.729166666666666667,
494
 
 5.0083333333333333332,
495
 
 -0.8541666666666666667,
496
 
 0.06458333333333333333
497
 
};
498
 
static double fd_coeffs_12_8[] = {
499
 
 233.56666666666666667,
500
 
 -194.2,
501
 
 110.625,
502
 
 -41.944444444444444444,
503
 
 10.0166666666666666664,
504
 
 -1.3666666666666666667,
505
 
 0.08611111111111111111
506
 
};
507
 
static double fd_coeffs_12_9[] = {
508
 
 54.0,
509
 
 -65.25,
510
 
 41.0,
511
 
 -15.0,
512
 
 3.0,
513
 
 -0.25
514
 
};
515
 
static double fd_coeffs_12_10[] = {
516
 
 -637.0,
517
 
 540.0,
518
 
 -326.25,
519
 
 136.66666666666666667,
520
 
 -37.5,
521
 
 6.0,
522
 
 -0.41666666666666666666
523
 
};
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
532
 
};
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
541
 
};
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
550
 
};
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
559
 
};
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
568
 
};
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
577
 
};
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
586
 
};
587
 
static double fd_coeffs_13_8[] = {
588
 
 0.10392659505208333333,
589
 
 -0.2229705810546875,
590
 
 0.18742031521267361111,
591
 
 -0.09048597547743055556,
592
 
 0.025958251953125,
593
 
 -0.004143948025173611111,
594
 
 2.9534233940972222222E-4
595
 
};
596
 
static double fd_coeffs_13_9[] = {
597
 
 0.93533935546875,
598
 
 -0.6689117431640625,
599
 
 0.3373565673828125,
600
 
 -0.116339111328125,
601
 
 0.025958251953125,
602
 
 -0.0033905029296875,
603
 
 2.044677734375E-4
604
 
};
605
 
static double fd_coeffs_13_10[] = {
606
 
 -0.055419921875,
607
 
 0.12249755859375,
608
 
 -0.109354654947916666664,
609
 
 0.057820638020833333332,
610
 
 -0.0186767578125,
611
 
 0.0033976236979166666667,
612
 
 -2.6448567708333333334E-4
613
 
};
614
 
static double fd_coeffs_14_1[] = {
615
 
 0.875,
616
 
 -0.29166666666666666666,
617
 
 0.09722222222222222222,
618
 
 -0.026515151515151515152,
619
 
 0.0053030303030303030304,
620
 
 -6.7987567987567987567E-4,
621
 
 4.1625041625041625042E-5
622
 
};
623
 
static double fd_coeffs_14_2[] = {
624
 
 -3.0235941043083900226,
625
 
 1.75,
626
 
 -0.29166666666666666666,
627
 
 0.064814814814814814815,
628
 
 -0.013257575757575757576,
629
 
 0.0021212121212121212122,
630
 
 -2.2662522662522662522E-4,
631
 
 1.1892869035726178583E-5
632
 
};
633
 
static double fd_coeffs_14_3[] = {
634
 
 -2.6869345238095238095,
635
 
 2.20814484126984127,
636
 
 -0.8170667989417989418,
637
 
 0.23056998556998556999,
638
 
 -0.046829906204906204906,
639
 
 0.006053691678691678692,
640
 
 -3.7247474747474747474E-4
641
 
};
642
 
static double fd_coeffs_14_4[] = {
643
 
 14.447883597883597884,
644
 
 -10.747738095238095238,
645
 
 4.41628968253968254,
646
 
 -1.0894223985890652557,
647
 
 0.23056998556998556999,
648
 
 -0.037463924963924963926,
649
 
 0.0040357944524611191277,
650
 
 -2.1284271284271284271E-4
651
 
};
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
660
 
};
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
670
 
};
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
679
 
};
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
689
 
};
690
 
static double fd_coeffs_14_9[] = {
691
 
 91.5375,
692
 
 -115.3,
693
 
 78.5375,
694
 
 -33.2,
695
 
 8.6875,
696
 
 -1.3,
697
 
 0.0875
698
 
};
699
 
static double fd_coeffs_14_10[] = {
700
 
 -1066.0,
701
 
 915.375,
702
 
 -576.5,
703
 
 261.79166666666666666,
704
 
 -83.0,
705
 
 17.375,
706
 
 -2.1666666666666666667,
707
 
 0.125
708
 
};
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
718
 
};
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
728
 
};
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
738
 
};
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
748
 
};
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
758
 
};
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
768
 
};
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
778
 
};
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
788
 
};
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
798
 
};
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
808
 
};
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
818
 
};
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
829
 
};
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
839
 
};
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
850
 
};
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
860
 
};
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
871
 
};
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
881
 
};
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
892
 
};
893
 
static double fd_coeffs_16_9[] = {
894
 
 129.55734126984126985,
895
 
 -168.52777777777777778,
896
 
 122.0875,
897
 
 -57.394444444444444445,
898
 
 17.993055555555555556,
899
 
 -3.6928571428571428572,
900
 
 0.45972222222222222223,
901
 
 -0.026587301587301587302
902
 
};
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
913
 
};
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
924
 
};
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
935
 
};
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
946
 
};
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
957
 
};
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
968
 
};
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
979
 
};
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
990
 
};
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
1001
 
};
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
1012
 
};
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
1023
 
};
1024
 
static double fd_coeffs_18_1[] = {
1025
 
 0.9,
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
1034
 
};
1035
 
static double fd_coeffs_18_2[] = {
1036
 
 -3.0795354623330813807,
1037
 
 1.8,
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
1046
 
};
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
1057
 
};
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
1069
 
};
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
1080
 
};
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
1092
 
};
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
1103
 
};
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
1115
 
};
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
1126
 
};
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
1138
 
};
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
1150
 
};
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
1162
 
};
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
1174
 
};
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
1186
 
};
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
1198
 
};
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
1210
 
};
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
1222
 
};
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
1234
 
};
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
1246
 
};
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
1258
 
};
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
1270
 
};
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
1283
 
};
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
1295
 
};
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
1308
 
};
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
1320
 
};
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
1333
 
};
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
1345
 
};
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
1358
 
};
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
1370
 
};
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
1383
 
};
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
1396
 
};
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
1409
 
};
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
1422
 
};
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
1435
 
};
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
1448
 
};
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
1461
 
};
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
1474
 
};
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
1487
 
};
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
1500
 
};
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
1513
 
};
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
1526
 
};
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
1540
 
};
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
1553
 
};
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
1567
 
};
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
1580
 
};
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
1594
 
};
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
1607
 
};
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
1621
 
};
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
1634
 
};
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
1648
 
};
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
1662
 
};
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
1676
 
};
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
1690
 
};
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
1704
 
};
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
1718
 
};
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
1732
 
};
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
1746
 
};
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
1760
 
};
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
1774
 
};
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
1788
 
};
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
1802
 
};
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
1817
 
};
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
1831
 
};
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
1846
 
};
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
1860
 
};
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
1875
 
};
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
1889
 
};
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
1904
 
};
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
1918
 
};
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
1933
 
};
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
1948
 
};
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
1963
 
};
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
1978
 
};
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
1993
 
};
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
2008
 
};
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
2023
 
};
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
2038
 
};
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
2053
 
};
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
2068
 
};
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
2083
 
};
2084
 
static void apply_periodic_fd(int n,double *x,double *y,double coeff,int m,double *c) {
2085
 
 int i,k,l=(m+1)/2;
2086
 
 if (m%2) {
2087
 
  for(i=m;i<n-m;++i) {
2088
 
   double d = 0.0, *c_ptr = c + l;
2089
 
   for(k=m;k>0;k-=2)
2090
 
    d += (*(--c_ptr)) * (x[i+k]-x[i-k]);
2091
 
   y[i] = coeff * d;
2092
 
  }
2093
 
  for(i=0;i<m;++i) {
2094
 
   double d = 0.0, *c_ptr = c + l;
2095
 
   for(k=m;k>0;k-=2)
2096
 
    d += (*(--c_ptr)) * (x[i+k]-x[(i+n-k)%n]);
2097
 
   y[i] = coeff * d;
2098
 
  }
2099
 
  for(i=n-m;i<n;++i) {
2100
 
   double d = 0.0, *c_ptr = c + l;
2101
 
   for(k=m;k>0;k-=2)
2102
 
    d += (*(--c_ptr)) * (x[(i+k)%n]-x[i-k]);
2103
 
   y[i] = coeff * d;
2104
 
  }
2105
 
 } else {
2106
 
  for(i=l;i<n-l;++i) {
2107
 
   double d = 0.0, *c_ptr = c + l;
2108
 
   for(k=l;k>0;--k)
2109
 
    d += (*(--c_ptr)) * (x[i+k]-x[i-k]);
2110
 
   y[i] = coeff * d;
2111
 
  }
2112
 
  for(i=0;i<l;++i) {
2113
 
   double d = 0.0, *c_ptr = c + l;
2114
 
   for(k=l;k>0;--k)
2115
 
    d += (*(--c_ptr)) * (x[i+k]-x[(i+n-k)%n]);
2116
 
   y[i] = coeff * d;
2117
 
  }
2118
 
  for(i=n-l;i<n;++i) {
2119
 
   double d = 0.0, *c_ptr = c + l;
2120
 
   for(k=l;k>0;--k)
2121
 
    d += (*(--c_ptr)) * (x[(i+k)%n]-x[i-k]);
2122
 
   y[i] = coeff * d;
2123
 
  }
2124
 
 }
2125
 
}
2126
 
static void apply_periodic_fd2(int n,double *x,double *y,double coeff,int m,double *c) {
2127
 
 int i,k,l=(m+1)/2;
2128
 
 if (m%2) {
2129
 
  for(i=m;i<n-m;++i) {
2130
 
   double d = 0.0, *c_ptr = c + l;
2131
 
   for(k=m;k>0;k-=2)
2132
 
    d += (*(--c_ptr)) * (x[i+k]+x[i-k]);
2133
 
   y[i] = coeff * d;
2134
 
  }
2135
 
  for(i=0;i<m;++i) {
2136
 
   double d = 0.0, *c_ptr = c + l;
2137
 
   for(k=m;k>0;k-=2)
2138
 
    d += (*(--c_ptr)) * (x[i+k]+x[(i+n-k)%n]);
2139
 
   y[i] = coeff * d;
2140
 
  }
2141
 
  for(i=n-m;i<n;++i) {
2142
 
   double d = 0.0, *c_ptr = c + l;
2143
 
   for(k=m;k>0;k-=2)
2144
 
    d += (*(--c_ptr)) * (x[(i+k)%n]+x[i-k]);
2145
 
   y[i] = coeff * d;
2146
 
  }
2147
 
 } else {
2148
 
  for(i=l;i<n-l;++i) {
2149
 
   double d = 0.0, *c_ptr = c + l + 1;
2150
 
   for(k=l;k>0;--k)
2151
 
    d += (*(--c_ptr)) * (x[i+k]+x[i-k]);
2152
 
   y[i] = coeff * (d+(*c)*x[i]);
2153
 
  }
2154
 
  for(i=0;i<l;++i) {
2155
 
   double d = 0.0, *c_ptr = c + l + 1;
2156
 
   for(k=l;k>0;--k)
2157
 
    d += (*(--c_ptr)) * (x[i+k]+x[(i+n-k)%n]);
2158
 
   y[i] = coeff * (d+(*c)*x[i]);
2159
 
  }
2160
 
  for(i=n-l;i<n;++i) {
2161
 
   double d = 0.0, *c_ptr = c + l + 1;
2162
 
   for(k=l;k>0;--k)
2163
 
    d += (*(--c_ptr)) * (x[(i+k)%n]+x[i-k]);
2164
 
   y[i] = coeff * (d+(*c)*x[i]);
2165
 
  }
2166
 
 }
2167
 
}
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;
2171
 
  switch (k*1000+m) {
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;
2377
 
  default:
2378
 
    fprintf(stderr,
2379
 
     "Unsupported derivative/order combination: k,m=%d,%d\n",
2380
 
      k,m);
2381
 
  }
2382
 
  if (c!=NULL) {
2383
 
   if (k%2)
2384
 
    apply_periodic_fd(n,x,y,coeff,m,c);
2385
 
   else
2386
 
    apply_periodic_fd2(n,x,y,coeff,m,c);
2387
 
  }
2388
 
}