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