|
Packit |
67cb25 |
#define bard_N 15
|
|
Packit |
67cb25 |
#define bard_P 3
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#define bard_NTRIES 3
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double bard_x0[bard_P] = { 1.0, 1.0, 1.0 };
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double bard_epsrel = 1.0e-8;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static double bard_Y[bard_N] = {
|
|
Packit |
67cb25 |
0.14, 0.18, 0.22, 0.25, 0.29, 0.32, 0.35, 0.39, 0.37,
|
|
Packit |
67cb25 |
0.58, 0.73, 0.96, 1.34, 2.10, 4.39
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static void
|
|
Packit |
67cb25 |
bard_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_exact1 = 8.214877306578963e-03;
|
|
Packit |
67cb25 |
const double bard_x1[bard_P] = { 8.241055975623580e-02,
|
|
Packit |
67cb25 |
1.133036092245175,
|
|
Packit |
67cb25 |
2.343695178435405 };
|
|
Packit |
67cb25 |
const double sumsq_exact2 = 17.42869333333333;
|
|
Packit |
67cb25 |
const double bard_x2[bard_P] = { 8.406666666666666e-01,
|
|
Packit |
67cb25 |
-99999.9, /* -inf */
|
|
Packit |
67cb25 |
-99999.9 }; /* -inf */
|
|
Packit |
67cb25 |
const double *bard_x;
|
|
Packit |
67cb25 |
double sumsq_exact;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
if (fabs(x[1]) < 10.0 && fabs(x[2]) < 10.0)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
bard_x = bard_x1;
|
|
Packit |
67cb25 |
sumsq_exact = sumsq_exact1;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
else
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
bard_x = bard_x2;
|
|
Packit |
67cb25 |
sumsq_exact = sumsq_exact2;
|
|
Packit |
67cb25 |
}
|
|
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 < bard_P; ++i)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (bard_x[i] < -90000.0)
|
|
Packit |
67cb25 |
continue;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_test_rel(x[i], bard_x[i], epsrel, "%s/%s i=%zu",
|
|
Packit |
67cb25 |
sname, pname, i);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static int
|
|
Packit |
67cb25 |
bard_f (const gsl_vector * x, void *params, gsl_vector * f)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x1 = gsl_vector_get(x, 0);
|
|
Packit |
67cb25 |
double x2 = gsl_vector_get(x, 1);
|
|
Packit |
67cb25 |
double x3 = gsl_vector_get(x, 2);
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < bard_N; ++i)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double ui = i + 1.0;
|
|
Packit |
67cb25 |
double vi = 16.0 - i - 1.0;
|
|
Packit |
67cb25 |
double wi = GSL_MIN(ui, vi);
|
|
Packit |
67cb25 |
double yi = bard_Y[i];
|
|
Packit |
67cb25 |
double fi = yi - (x1 + (ui / (x2*vi + x3*wi)));
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_vector_set(f, i, fi);
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return GSL_SUCCESS;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static int
|
|
Packit |
67cb25 |
bard_df (const gsl_vector * x, void *params, gsl_matrix * J)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x2 = gsl_vector_get(x, 1);
|
|
Packit |
67cb25 |
double x3 = gsl_vector_get(x, 2);
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < bard_N; ++i)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double ui = i + 1.0;
|
|
Packit |
67cb25 |
double vi = 16.0 - i - 1.0;
|
|
Packit |
67cb25 |
double wi = GSL_MIN(ui, vi);
|
|
Packit |
67cb25 |
double term = x2 * vi + x3 * wi;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
gsl_matrix_set(J, i, 0, -1.0);
|
|
Packit |
67cb25 |
gsl_matrix_set(J, i, 1, ui * vi / (term * term));
|
|
Packit |
67cb25 |
gsl_matrix_set(J, i, 2, ui * wi / (term * term));
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return GSL_SUCCESS;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static gsl_multifit_function_fdf bard_func =
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
&bard_f,
|
|
Packit |
67cb25 |
&bard_df,
|
|
Packit |
67cb25 |
NULL,
|
|
Packit |
67cb25 |
bard_N,
|
|
Packit |
67cb25 |
bard_P,
|
|
Packit |
67cb25 |
NULL,
|
|
Packit |
67cb25 |
0,
|
|
Packit |
67cb25 |
0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static test_fdf_problem bard_problem =
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
"bard",
|
|
Packit |
67cb25 |
bard_x0,
|
|
Packit |
67cb25 |
NULL,
|
|
Packit |
67cb25 |
&bard_epsrel,
|
|
Packit |
67cb25 |
bard_NTRIES,
|
|
Packit |
67cb25 |
&bard_checksol,
|
|
Packit |
67cb25 |
&bard_func
|
|
Packit |
67cb25 |
};
|