|
Packit |
67cb25 |
#define hahn1_N 236
|
|
Packit |
67cb25 |
#define hahn1_P 7
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#define hahn1_NTRIES 1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* double hahn1_x0[7] = { 10, -1, 0.05, -0.00001, -0.05, 0.001, -0.000001 }; */
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double hahn1_x0[hahn1_P] = { 1, -0.1, 0.005, -0.000001, -0.005, 0.0001, -0.0000001};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double hahn1_epsrel = 1.0e-5;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double hahn1_sigma[hahn1_P] = {
|
|
Packit |
67cb25 |
1.7070154742E-01,
|
|
Packit |
67cb25 |
1.2000289189E-02,
|
|
Packit |
67cb25 |
2.2508314937E-04,
|
|
Packit |
67cb25 |
2.7578037666E-07,
|
|
Packit |
67cb25 |
2.4712888219E-04,
|
|
Packit |
67cb25 |
1.0449373768E-05,
|
|
Packit |
67cb25 |
1.3027335327E-08
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double hahn1_F1[hahn1_N] = {
|
|
Packit |
67cb25 |
.591E0,
|
|
Packit |
67cb25 |
1.547E0,
|
|
Packit |
67cb25 |
2.902E0,
|
|
Packit |
67cb25 |
2.894E0,
|
|
Packit |
67cb25 |
4.703E0,
|
|
Packit |
67cb25 |
6.307E0,
|
|
Packit |
67cb25 |
7.03E0 ,
|
|
Packit |
67cb25 |
7.898E0,
|
|
Packit |
67cb25 |
9.470E0,
|
|
Packit |
67cb25 |
9.484E0,
|
|
Packit |
67cb25 |
10.072E0,
|
|
Packit |
67cb25 |
10.163E0,
|
|
Packit |
67cb25 |
11.615E0,
|
|
Packit |
67cb25 |
12.005E0,
|
|
Packit |
67cb25 |
12.478E0,
|
|
Packit |
67cb25 |
12.982E0,
|
|
Packit |
67cb25 |
12.970E0,
|
|
Packit |
67cb25 |
13.926E0,
|
|
Packit |
67cb25 |
14.452E0,
|
|
Packit |
67cb25 |
14.404E0,
|
|
Packit |
67cb25 |
15.190E0,
|
|
Packit |
67cb25 |
15.550E0,
|
|
Packit |
67cb25 |
15.528E0,
|
|
Packit |
67cb25 |
15.499E0,
|
|
Packit |
67cb25 |
16.131E0,
|
|
Packit |
67cb25 |
16.438E0,
|
|
Packit |
67cb25 |
16.387E0,
|
|
Packit |
67cb25 |
16.549E0,
|
|
Packit |
67cb25 |
16.872E0,
|
|
Packit |
67cb25 |
16.830E0,
|
|
Packit |
67cb25 |
16.926E0,
|
|
Packit |
67cb25 |
16.907E0,
|
|
Packit |
67cb25 |
16.966E0,
|
|
Packit |
67cb25 |
17.060E0,
|
|
Packit |
67cb25 |
17.122E0,
|
|
Packit |
67cb25 |
17.311E0,
|
|
Packit |
67cb25 |
17.355E0,
|
|
Packit |
67cb25 |
17.668E0,
|
|
Packit |
67cb25 |
17.767E0,
|
|
Packit |
67cb25 |
17.803E0,
|
|
Packit |
67cb25 |
17.765E0,
|
|
Packit |
67cb25 |
17.768E0,
|
|
Packit |
67cb25 |
17.736E0,
|
|
Packit |
67cb25 |
17.858E0,
|
|
Packit |
67cb25 |
17.877E0,
|
|
Packit |
67cb25 |
17.912E0,
|
|
Packit |
67cb25 |
18.046E0,
|
|
Packit |
67cb25 |
18.085E0,
|
|
Packit |
67cb25 |
18.291E0,
|
|
Packit |
67cb25 |
18.357E0,
|
|
Packit |
67cb25 |
18.426E0,
|
|
Packit |
67cb25 |
18.584E0,
|
|
Packit |
67cb25 |
18.610E0,
|
|
Packit |
67cb25 |
18.870E0,
|
|
Packit |
67cb25 |
18.795E0,
|
|
Packit |
67cb25 |
19.111E0,
|
|
Packit |
67cb25 |
.367E0,
|
|
Packit |
67cb25 |
.796E0,
|
|
Packit |
67cb25 |
0.892E0,
|
|
Packit |
67cb25 |
1.903E0,
|
|
Packit |
67cb25 |
2.150E0,
|
|
Packit |
67cb25 |
3.697E0,
|
|
Packit |
67cb25 |
5.870E0,
|
|
Packit |
67cb25 |
6.421E0,
|
|
Packit |
67cb25 |
7.422E0,
|
|
Packit |
67cb25 |
9.944E0,
|
|
Packit |
67cb25 |
11.023E0,
|
|
Packit |
67cb25 |
11.87E0 ,
|
|
Packit |
67cb25 |
12.786E0,
|
|
Packit |
67cb25 |
14.067E0,
|
|
Packit |
67cb25 |
13.974E0,
|
|
Packit |
67cb25 |
14.462E0,
|
|
Packit |
67cb25 |
14.464E0,
|
|
Packit |
67cb25 |
15.381E0,
|
|
Packit |
67cb25 |
15.483E0,
|
|
Packit |
67cb25 |
15.59E0 ,
|
|
Packit |
67cb25 |
16.075E0,
|
|
Packit |
67cb25 |
16.347E0,
|
|
Packit |
67cb25 |
16.181E0,
|
|
Packit |
67cb25 |
16.915E0,
|
|
Packit |
67cb25 |
17.003E0,
|
|
Packit |
67cb25 |
16.978E0,
|
|
Packit |
67cb25 |
17.756E0,
|
|
Packit |
67cb25 |
17.808E0,
|
|
Packit |
67cb25 |
17.868E0,
|
|
Packit |
67cb25 |
18.481E0,
|
|
Packit |
67cb25 |
18.486E0,
|
|
Packit |
67cb25 |
19.090E0,
|
|
Packit |
67cb25 |
16.062E0,
|
|
Packit |
67cb25 |
16.337E0,
|
|
Packit |
67cb25 |
16.345E0,
|
|
Packit |
67cb25 |
16.388E0,
|
|
Packit |
67cb25 |
17.159E0,
|
|
Packit |
67cb25 |
17.116E0,
|
|
Packit |
67cb25 |
17.164E0,
|
|
Packit |
67cb25 |
17.123E0,
|
|
Packit |
67cb25 |
17.979E0,
|
|
Packit |
67cb25 |
17.974E0,
|
|
Packit |
67cb25 |
18.007E0,
|
|
Packit |
67cb25 |
17.993E0,
|
|
Packit |
67cb25 |
18.523E0,
|
|
Packit |
67cb25 |
18.669E0,
|
|
Packit |
67cb25 |
18.617E0,
|
|
Packit |
67cb25 |
19.371E0,
|
|
Packit |
67cb25 |
19.330E0,
|
|
Packit |
67cb25 |
0.080E0,
|
|
Packit |
67cb25 |
0.248E0,
|
|
Packit |
67cb25 |
1.089E0,
|
|
Packit |
67cb25 |
1.418E0,
|
|
Packit |
67cb25 |
2.278E0,
|
|
Packit |
67cb25 |
3.624E0,
|
|
Packit |
67cb25 |
4.574E0,
|
|
Packit |
67cb25 |
5.556E0,
|
|
Packit |
67cb25 |
7.267E0,
|
|
Packit |
67cb25 |
7.695E0,
|
|
Packit |
67cb25 |
9.136E0,
|
|
Packit |
67cb25 |
9.959E0,
|
|
Packit |
67cb25 |
9.957E0,
|
|
Packit |
67cb25 |
11.600E0,
|
|
Packit |
67cb25 |
13.138E0,
|
|
Packit |
67cb25 |
13.564E0,
|
|
Packit |
67cb25 |
13.871E0,
|
|
Packit |
67cb25 |
13.994E0,
|
|
Packit |
67cb25 |
14.947E0,
|
|
Packit |
67cb25 |
15.473E0,
|
|
Packit |
67cb25 |
15.379E0,
|
|
Packit |
67cb25 |
15.455E0,
|
|
Packit |
67cb25 |
15.908E0,
|
|
Packit |
67cb25 |
16.114E0,
|
|
Packit |
67cb25 |
17.071E0,
|
|
Packit |
67cb25 |
17.135E0,
|
|
Packit |
67cb25 |
17.282E0,
|
|
Packit |
67cb25 |
17.368E0,
|
|
Packit |
67cb25 |
17.483E0,
|
|
Packit |
67cb25 |
17.764E0,
|
|
Packit |
67cb25 |
18.185E0,
|
|
Packit |
67cb25 |
18.271E0,
|
|
Packit |
67cb25 |
18.236E0,
|
|
Packit |
67cb25 |
18.237E0,
|
|
Packit |
67cb25 |
18.523E0,
|
|
Packit |
67cb25 |
18.627E0,
|
|
Packit |
67cb25 |
18.665E0,
|
|
Packit |
67cb25 |
19.086E0,
|
|
Packit |
67cb25 |
0.214E0,
|
|
Packit |
67cb25 |
0.943E0,
|
|
Packit |
67cb25 |
1.429E0,
|
|
Packit |
67cb25 |
2.241E0,
|
|
Packit |
67cb25 |
2.951E0,
|
|
Packit |
67cb25 |
3.782E0,
|
|
Packit |
67cb25 |
4.757E0,
|
|
Packit |
67cb25 |
5.602E0,
|
|
Packit |
67cb25 |
7.169E0,
|
|
Packit |
67cb25 |
8.920E0,
|
|
Packit |
67cb25 |
10.055E0,
|
|
Packit |
67cb25 |
12.035E0,
|
|
Packit |
67cb25 |
12.861E0,
|
|
Packit |
67cb25 |
13.436E0,
|
|
Packit |
67cb25 |
14.167E0,
|
|
Packit |
67cb25 |
14.755E0,
|
|
Packit |
67cb25 |
15.168E0,
|
|
Packit |
67cb25 |
15.651E0,
|
|
Packit |
67cb25 |
15.746E0,
|
|
Packit |
67cb25 |
16.216E0,
|
|
Packit |
67cb25 |
16.445E0,
|
|
Packit |
67cb25 |
16.965E0,
|
|
Packit |
67cb25 |
17.121E0,
|
|
Packit |
67cb25 |
17.206E0,
|
|
Packit |
67cb25 |
17.250E0,
|
|
Packit |
67cb25 |
17.339E0,
|
|
Packit |
67cb25 |
17.793E0,
|
|
Packit |
67cb25 |
18.123E0,
|
|
Packit |
67cb25 |
18.49E0 ,
|
|
Packit |
67cb25 |
18.566E0,
|
|
Packit |
67cb25 |
18.645E0,
|
|
Packit |
67cb25 |
18.706E0,
|
|
Packit |
67cb25 |
18.924E0,
|
|
Packit |
67cb25 |
19.1E0 ,
|
|
Packit |
67cb25 |
0.375E0,
|
|
Packit |
67cb25 |
0.471E0,
|
|
Packit |
67cb25 |
1.504E0,
|
|
Packit |
67cb25 |
2.204E0,
|
|
Packit |
67cb25 |
2.813E0,
|
|
Packit |
67cb25 |
4.765E0,
|
|
Packit |
67cb25 |
9.835E0,
|
|
Packit |
67cb25 |
10.040E0,
|
|
Packit |
67cb25 |
11.946E0,
|
|
Packit |
67cb25 |
12.596E0,
|
|
Packit |
67cb25 |
13.303E0,
|
|
Packit |
67cb25 |
13.922E0,
|
|
Packit |
67cb25 |
14.440E0,
|
|
Packit |
67cb25 |
14.951E0,
|
|
Packit |
67cb25 |
15.627E0,
|
|
Packit |
67cb25 |
15.639E0,
|
|
Packit |
67cb25 |
15.814E0,
|
|
Packit |
67cb25 |
16.315E0,
|
|
Packit |
67cb25 |
16.334E0,
|
|
Packit |
67cb25 |
16.430E0,
|
|
Packit |
67cb25 |
16.423E0,
|
|
Packit |
67cb25 |
17.024E0,
|
|
Packit |
67cb25 |
17.009E0,
|
|
Packit |
67cb25 |
17.165E0,
|
|
Packit |
67cb25 |
17.134E0,
|
|
Packit |
67cb25 |
17.349E0,
|
|
Packit |
67cb25 |
17.576E0,
|
|
Packit |
67cb25 |
17.848E0,
|
|
Packit |
67cb25 |
18.090E0,
|
|
Packit |
67cb25 |
18.276E0,
|
|
Packit |
67cb25 |
18.404E0,
|
|
Packit |
67cb25 |
18.519E0,
|
|
Packit |
67cb25 |
19.133E0,
|
|
Packit |
67cb25 |
19.074E0,
|
|
Packit |
67cb25 |
19.239E0,
|
|
Packit |
67cb25 |
19.280E0,
|
|
Packit |
67cb25 |
19.101E0,
|
|
Packit |
67cb25 |
19.398E0,
|
|
Packit |
67cb25 |
19.252E0,
|
|
Packit |
67cb25 |
19.89E0 ,
|
|
Packit |
67cb25 |
20.007E0,
|
|
Packit |
67cb25 |
19.929E0,
|
|
Packit |
67cb25 |
19.268E0,
|
|
Packit |
67cb25 |
19.324E0,
|
|
Packit |
67cb25 |
20.049E0,
|
|
Packit |
67cb25 |
20.107E0,
|
|
Packit |
67cb25 |
20.062E0,
|
|
Packit |
67cb25 |
20.065E0,
|
|
Packit |
67cb25 |
19.286E0,
|
|
Packit |
67cb25 |
19.972E0,
|
|
Packit |
67cb25 |
20.088E0,
|
|
Packit |
67cb25 |
20.743E0,
|
|
Packit |
67cb25 |
20.83E0 ,
|
|
Packit |
67cb25 |
20.935E0,
|
|
Packit |
67cb25 |
21.035E0,
|
|
Packit |
67cb25 |
20.93E0 ,
|
|
Packit |
67cb25 |
21.074E0,
|
|
Packit |
67cb25 |
21.085E0,
|
|
Packit |
67cb25 |
20.935E0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double hahn1_F0[hahn1_N] = {
|
|
Packit |
67cb25 |
24.41E0,
|
|
Packit |
67cb25 |
34.82E0,
|
|
Packit |
67cb25 |
44.09E0,
|
|
Packit |
67cb25 |
45.07E0,
|
|
Packit |
67cb25 |
54.98E0,
|
|
Packit |
67cb25 |
65.51E0,
|
|
Packit |
67cb25 |
70.53E0,
|
|
Packit |
67cb25 |
75.70E0,
|
|
Packit |
67cb25 |
89.57E0,
|
|
Packit |
67cb25 |
91.14E0,
|
|
Packit |
67cb25 |
96.40E0,
|
|
Packit |
67cb25 |
97.19E0,
|
|
Packit |
67cb25 |
114.26E0,
|
|
Packit |
67cb25 |
120.25E0,
|
|
Packit |
67cb25 |
127.08E0,
|
|
Packit |
67cb25 |
133.55E0,
|
|
Packit |
67cb25 |
133.61E0,
|
|
Packit |
67cb25 |
158.67E0,
|
|
Packit |
67cb25 |
172.74E0,
|
|
Packit |
67cb25 |
171.31E0,
|
|
Packit |
67cb25 |
202.14E0,
|
|
Packit |
67cb25 |
220.55E0,
|
|
Packit |
67cb25 |
221.05E0,
|
|
Packit |
67cb25 |
221.39E0,
|
|
Packit |
67cb25 |
250.99E0,
|
|
Packit |
67cb25 |
268.99E0,
|
|
Packit |
67cb25 |
271.80E0,
|
|
Packit |
67cb25 |
271.97E0,
|
|
Packit |
67cb25 |
321.31E0,
|
|
Packit |
67cb25 |
321.69E0,
|
|
Packit |
67cb25 |
330.14E0,
|
|
Packit |
67cb25 |
333.03E0,
|
|
Packit |
67cb25 |
333.47E0,
|
|
Packit |
67cb25 |
340.77E0,
|
|
Packit |
67cb25 |
345.65E0,
|
|
Packit |
67cb25 |
373.11E0,
|
|
Packit |
67cb25 |
373.79E0,
|
|
Packit |
67cb25 |
411.82E0,
|
|
Packit |
67cb25 |
419.51E0,
|
|
Packit |
67cb25 |
421.59E0,
|
|
Packit |
67cb25 |
422.02E0,
|
|
Packit |
67cb25 |
422.47E0,
|
|
Packit |
67cb25 |
422.61E0,
|
|
Packit |
67cb25 |
441.75E0,
|
|
Packit |
67cb25 |
447.41E0,
|
|
Packit |
67cb25 |
448.7E0 ,
|
|
Packit |
67cb25 |
472.89E0,
|
|
Packit |
67cb25 |
476.69E0,
|
|
Packit |
67cb25 |
522.47E0,
|
|
Packit |
67cb25 |
522.62E0,
|
|
Packit |
67cb25 |
524.43E0,
|
|
Packit |
67cb25 |
546.75E0,
|
|
Packit |
67cb25 |
549.53E0,
|
|
Packit |
67cb25 |
575.29E0,
|
|
Packit |
67cb25 |
576.00E0,
|
|
Packit |
67cb25 |
625.55E0,
|
|
Packit |
67cb25 |
20.15E0,
|
|
Packit |
67cb25 |
28.78E0,
|
|
Packit |
67cb25 |
29.57E0,
|
|
Packit |
67cb25 |
37.41E0,
|
|
Packit |
67cb25 |
39.12E0,
|
|
Packit |
67cb25 |
50.24E0,
|
|
Packit |
67cb25 |
61.38E0,
|
|
Packit |
67cb25 |
66.25E0,
|
|
Packit |
67cb25 |
73.42E0,
|
|
Packit |
67cb25 |
95.52E0,
|
|
Packit |
67cb25 |
107.32E0,
|
|
Packit |
67cb25 |
122.04E0,
|
|
Packit |
67cb25 |
134.03E0,
|
|
Packit |
67cb25 |
163.19E0,
|
|
Packit |
67cb25 |
163.48E0,
|
|
Packit |
67cb25 |
175.70E0,
|
|
Packit |
67cb25 |
179.86E0,
|
|
Packit |
67cb25 |
211.27E0,
|
|
Packit |
67cb25 |
217.78E0,
|
|
Packit |
67cb25 |
219.14E0,
|
|
Packit |
67cb25 |
262.52E0,
|
|
Packit |
67cb25 |
268.01E0,
|
|
Packit |
67cb25 |
268.62E0,
|
|
Packit |
67cb25 |
336.25E0,
|
|
Packit |
67cb25 |
337.23E0,
|
|
Packit |
67cb25 |
339.33E0,
|
|
Packit |
67cb25 |
427.38E0,
|
|
Packit |
67cb25 |
428.58E0,
|
|
Packit |
67cb25 |
432.68E0,
|
|
Packit |
67cb25 |
528.99E0,
|
|
Packit |
67cb25 |
531.08E0,
|
|
Packit |
67cb25 |
628.34E0,
|
|
Packit |
67cb25 |
253.24E0,
|
|
Packit |
67cb25 |
273.13E0,
|
|
Packit |
67cb25 |
273.66E0,
|
|
Packit |
67cb25 |
282.10E0,
|
|
Packit |
67cb25 |
346.62E0,
|
|
Packit |
67cb25 |
347.19E0,
|
|
Packit |
67cb25 |
348.78E0,
|
|
Packit |
67cb25 |
351.18E0,
|
|
Packit |
67cb25 |
450.10E0,
|
|
Packit |
67cb25 |
450.35E0,
|
|
Packit |
67cb25 |
451.92E0,
|
|
Packit |
67cb25 |
455.56E0,
|
|
Packit |
67cb25 |
552.22E0,
|
|
Packit |
67cb25 |
553.56E0,
|
|
Packit |
67cb25 |
555.74E0,
|
|
Packit |
67cb25 |
652.59E0,
|
|
Packit |
67cb25 |
656.20E0,
|
|
Packit |
67cb25 |
14.13E0,
|
|
Packit |
67cb25 |
20.41E0,
|
|
Packit |
67cb25 |
31.30E0,
|
|
Packit |
67cb25 |
33.84E0,
|
|
Packit |
67cb25 |
39.70E0,
|
|
Packit |
67cb25 |
48.83E0,
|
|
Packit |
67cb25 |
54.50E0,
|
|
Packit |
67cb25 |
60.41E0,
|
|
Packit |
67cb25 |
72.77E0,
|
|
Packit |
67cb25 |
75.25E0,
|
|
Packit |
67cb25 |
86.84E0,
|
|
Packit |
67cb25 |
94.88E0,
|
|
Packit |
67cb25 |
96.40E0,
|
|
Packit |
67cb25 |
117.37E0,
|
|
Packit |
67cb25 |
139.08E0,
|
|
Packit |
67cb25 |
147.73E0,
|
|
Packit |
67cb25 |
158.63E0,
|
|
Packit |
67cb25 |
161.84E0,
|
|
Packit |
67cb25 |
192.11E0,
|
|
Packit |
67cb25 |
206.76E0,
|
|
Packit |
67cb25 |
209.07E0,
|
|
Packit |
67cb25 |
213.32E0,
|
|
Packit |
67cb25 |
226.44E0,
|
|
Packit |
67cb25 |
237.12E0,
|
|
Packit |
67cb25 |
330.90E0,
|
|
Packit |
67cb25 |
358.72E0,
|
|
Packit |
67cb25 |
370.77E0,
|
|
Packit |
67cb25 |
372.72E0,
|
|
Packit |
67cb25 |
396.24E0,
|
|
Packit |
67cb25 |
416.59E0,
|
|
Packit |
67cb25 |
484.02E0,
|
|
Packit |
67cb25 |
495.47E0,
|
|
Packit |
67cb25 |
514.78E0,
|
|
Packit |
67cb25 |
515.65E0,
|
|
Packit |
67cb25 |
519.47E0,
|
|
Packit |
67cb25 |
544.47E0,
|
|
Packit |
67cb25 |
560.11E0,
|
|
Packit |
67cb25 |
620.77E0,
|
|
Packit |
67cb25 |
18.97E0,
|
|
Packit |
67cb25 |
28.93E0,
|
|
Packit |
67cb25 |
33.91E0,
|
|
Packit |
67cb25 |
40.03E0,
|
|
Packit |
67cb25 |
44.66E0,
|
|
Packit |
67cb25 |
49.87E0,
|
|
Packit |
67cb25 |
55.16E0,
|
|
Packit |
67cb25 |
60.90E0,
|
|
Packit |
67cb25 |
72.08E0,
|
|
Packit |
67cb25 |
85.15E0,
|
|
Packit |
67cb25 |
97.06E0,
|
|
Packit |
67cb25 |
119.63E0,
|
|
Packit |
67cb25 |
133.27E0,
|
|
Packit |
67cb25 |
143.84E0,
|
|
Packit |
67cb25 |
161.91E0,
|
|
Packit |
67cb25 |
180.67E0,
|
|
Packit |
67cb25 |
198.44E0,
|
|
Packit |
67cb25 |
226.86E0,
|
|
Packit |
67cb25 |
229.65E0,
|
|
Packit |
67cb25 |
258.27E0,
|
|
Packit |
67cb25 |
273.77E0,
|
|
Packit |
67cb25 |
339.15E0,
|
|
Packit |
67cb25 |
350.13E0,
|
|
Packit |
67cb25 |
362.75E0,
|
|
Packit |
67cb25 |
371.03E0,
|
|
Packit |
67cb25 |
393.32E0,
|
|
Packit |
67cb25 |
448.53E0,
|
|
Packit |
67cb25 |
473.78E0,
|
|
Packit |
67cb25 |
511.12E0,
|
|
Packit |
67cb25 |
524.70E0,
|
|
Packit |
67cb25 |
548.75E0,
|
|
Packit |
67cb25 |
551.64E0,
|
|
Packit |
67cb25 |
574.02E0,
|
|
Packit |
67cb25 |
623.86E0,
|
|
Packit |
67cb25 |
21.46E0,
|
|
Packit |
67cb25 |
24.33E0,
|
|
Packit |
67cb25 |
33.43E0,
|
|
Packit |
67cb25 |
39.22E0,
|
|
Packit |
67cb25 |
44.18E0,
|
|
Packit |
67cb25 |
55.02E0,
|
|
Packit |
67cb25 |
94.33E0,
|
|
Packit |
67cb25 |
96.44E0,
|
|
Packit |
67cb25 |
118.82E0,
|
|
Packit |
67cb25 |
128.48E0,
|
|
Packit |
67cb25 |
141.94E0,
|
|
Packit |
67cb25 |
156.92E0,
|
|
Packit |
67cb25 |
171.65E0,
|
|
Packit |
67cb25 |
190.00E0,
|
|
Packit |
67cb25 |
223.26E0,
|
|
Packit |
67cb25 |
223.88E0,
|
|
Packit |
67cb25 |
231.50E0,
|
|
Packit |
67cb25 |
265.05E0,
|
|
Packit |
67cb25 |
269.44E0,
|
|
Packit |
67cb25 |
271.78E0,
|
|
Packit |
67cb25 |
273.46E0,
|
|
Packit |
67cb25 |
334.61E0,
|
|
Packit |
67cb25 |
339.79E0,
|
|
Packit |
67cb25 |
349.52E0,
|
|
Packit |
67cb25 |
358.18E0,
|
|
Packit |
67cb25 |
377.98E0,
|
|
Packit |
67cb25 |
394.77E0,
|
|
Packit |
67cb25 |
429.66E0,
|
|
Packit |
67cb25 |
468.22E0,
|
|
Packit |
67cb25 |
487.27E0,
|
|
Packit |
67cb25 |
519.54E0,
|
|
Packit |
67cb25 |
523.03E0,
|
|
Packit |
67cb25 |
612.99E0,
|
|
Packit |
67cb25 |
638.59E0,
|
|
Packit |
67cb25 |
641.36E0,
|
|
Packit |
67cb25 |
622.05E0,
|
|
Packit |
67cb25 |
631.50E0,
|
|
Packit |
67cb25 |
663.97E0,
|
|
Packit |
67cb25 |
646.9E0 ,
|
|
Packit |
67cb25 |
748.29E0,
|
|
Packit |
67cb25 |
749.21E0,
|
|
Packit |
67cb25 |
750.14E0,
|
|
Packit |
67cb25 |
647.04E0,
|
|
Packit |
67cb25 |
646.89E0,
|
|
Packit |
67cb25 |
746.9E0 ,
|
|
Packit |
67cb25 |
748.43E0,
|
|
Packit |
67cb25 |
747.35E0,
|
|
Packit |
67cb25 |
749.27E0,
|
|
Packit |
67cb25 |
647.61E0,
|
|
Packit |
67cb25 |
747.78E0,
|
|
Packit |
67cb25 |
750.51E0,
|
|
Packit |
67cb25 |
851.37E0,
|
|
Packit |
67cb25 |
845.97E0,
|
|
Packit |
67cb25 |
847.54E0,
|
|
Packit |
67cb25 |
849.93E0,
|
|
Packit |
67cb25 |
851.61E0,
|
|
Packit |
67cb25 |
849.75E0,
|
|
Packit |
67cb25 |
850.98E0,
|
|
Packit |
67cb25 |
848.23E0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static void
|
|
Packit |
67cb25 |
hahn1_checksol(const double x[], const double sumsq,
|
|
Packit |
67cb25 |
const double epsrel, const char *sname,
|
|
Packit |
67cb25 |
const char *pname)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
const double sumsq_exact = 1.5324382854E+00;
|
|
Packit |
67cb25 |
const double hahn1_x[hahn1_P] = {
|
|
Packit |
67cb25 |
1.0776351733E+00, -1.2269296921E-01, 4.0863750610E-03,
|
|
Packit |
67cb25 |
-1.4262662514E-06, -5.7609940901E-03, 2.4053735503E-04,
|
|
Packit |
67cb25 |
-1.2314450199E-07 };
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test_rel(sumsq, sumsq_exact, epsrel, "%s/%s sumsq",
|
|
Packit |
67cb25 |
sname, pname);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < hahn1_P; ++i)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
gsl_test_rel(x[i], hahn1_x[i], epsrel, "%s/%s i=%zu",
|
|
Packit |
67cb25 |
sname, pname, i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static int
|
|
Packit |
67cb25 |
hahn1_f (const gsl_vector * x, void *params, gsl_vector * f)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double b[hahn1_P];
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < hahn1_P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
b[i] = gsl_vector_get(x, i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < hahn1_N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x = hahn1_F0[i];
|
|
Packit |
67cb25 |
double y = ((b[0] + x* (b[1] + x * (b[2] + x * b[3])))
|
|
Packit |
67cb25 |
/ (1 + x*(b[4] + x *(b[5] + x*b[6]))));
|
|
Packit |
67cb25 |
gsl_vector_set (f, i, hahn1_F1[i] - y);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return GSL_SUCCESS;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static int
|
|
Packit |
67cb25 |
hahn1_df (const gsl_vector * x, void *params, gsl_matrix * df)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double b[hahn1_P];
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < hahn1_P; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
b[i] = gsl_vector_get(x, i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < hahn1_N; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x = hahn1_F0[i];
|
|
Packit |
67cb25 |
double u = (b[0] + x*(b[1] + x*(b[2] + x * b[3])));
|
|
Packit |
67cb25 |
double v = (1 + x*(b[4] + x*(b[5] + x*b[6])));
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 0, -1/v);
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 1, -x/v);
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 2, -x*x/v);
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 3, -x*x*x/v);
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 4, x*u/(v*v));
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 5, x*x*u/(v*v));
|
|
Packit |
67cb25 |
gsl_matrix_set (df, i, 6, x*x*x*u/(v*v));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return GSL_SUCCESS;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static gsl_multifit_function_fdf hahn1_func =
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
&hahn1_f,
|
|
Packit |
67cb25 |
&hahn1_df,
|
|
Packit |
67cb25 |
NULL,
|
|
Packit |
67cb25 |
hahn1_N,
|
|
Packit |
67cb25 |
hahn1_P,
|
|
Packit |
67cb25 |
NULL,
|
|
Packit |
67cb25 |
0,
|
|
Packit |
67cb25 |
0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static test_fdf_problem hahn1_problem =
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
"nist-hahn1",
|
|
Packit |
67cb25 |
hahn1_x0,
|
|
Packit |
67cb25 |
hahn1_sigma,
|
|
Packit |
67cb25 |
&hahn1_epsrel,
|
|
Packit |
67cb25 |
hahn1_NTRIES,
|
|
Packit |
67cb25 |
&hahn1_checksol,
|
|
Packit |
67cb25 |
&hahn1_func
|
|
Packit |
67cb25 |
};
|