#define kirby2_N       151

#define kirby2_P       5

#define kirby2_NTRIES  1

/* double kirby2_x0[kirby2_P] = { 2, -0.1, 0.003, -0.001, 0.00001 }; */

static double kirby2_x0[kirby2_P] = { 1.5, -0.15, 0.0025, -0.0015, 0.00002 }; 

static double kirby2_epsrel = 1.0e-5;

static double kirby2_sigma[kirby2_P] = {

  8.7989634338E-02,

  4.1182041386E-03,

  4.1856520458E-05,

  5.8931897355E-05,

  2.0129761919E-07

};

static double kirby2_F1[kirby2_N] = {

       0.0082E0,

       0.0112E0,

       0.0149E0,

       0.0198E0,

       0.0248E0,

       0.0324E0,

       0.0420E0,

       0.0549E0,

       0.0719E0,

       0.0963E0,

       0.1291E0,

       0.1710E0,

       0.2314E0,

       0.3227E0,

       0.4809E0,

       0.7084E0,

       1.0220E0,

       1.4580E0,

       1.9520E0,

       2.5410E0,

       3.2230E0,

       3.9990E0,

       4.8520E0,

       5.7320E0,

       6.7270E0,

       7.8350E0,

       9.0250E0,

      10.2670E0,

      11.5780E0,

      12.9440E0,

      14.3770E0,

      15.8560E0,

      17.3310E0,

      18.8850E0,

      20.5750E0,

      22.3200E0,

      22.3030E0,

      23.4600E0,

      24.0600E0,

      25.2720E0,

      25.8530E0,

      27.1100E0,

      27.6580E0,

      28.9240E0,

      29.5110E0,

      30.7100E0,

      31.3500E0,

      32.5200E0,

      33.2300E0,

      34.3300E0,

      35.0600E0,

      36.1700E0,

      36.8400E0,

      38.0100E0,

      38.6700E0,

      39.8700E0,

      40.0300E0,

      40.5000E0,

      41.3700E0,

      41.6700E0,

      42.3100E0,

      42.7300E0,

      43.4600E0,

      44.1400E0,

      44.5500E0,

      45.2200E0,

      45.9200E0,

      46.3000E0,

      47.0000E0,

      47.6800E0,

      48.0600E0,

      48.7400E0,

      49.4100E0,

      49.7600E0,

      50.4300E0,

      51.1100E0,

      51.5000E0,

      52.1200E0,

      52.7600E0,

      53.1800E0,

      53.7800E0,

      54.4600E0,

      54.8300E0,

      55.4000E0,

      56.4300E0,

      57.0300E0,

      58.0000E0,

      58.6100E0,

      59.5800E0,

      60.1100E0,

      61.1000E0,

      61.6500E0,

      62.5900E0,

      63.1200E0,

      64.0300E0,

      64.6200E0,

      65.4900E0,

      66.0300E0,

      66.8900E0,

      67.4200E0,

      68.2300E0,

      68.7700E0,

      69.5900E0,

      70.1100E0,

      70.8600E0,

      71.4300E0,

      72.1600E0,

      72.7000E0,

      73.4000E0,

      73.9300E0,

      74.6000E0,

      75.1600E0,

      75.8200E0,

      76.3400E0,

      76.9800E0,

      77.4800E0,

      78.0800E0,

      78.6000E0,

      79.1700E0,

      79.6200E0,

      79.8800E0,

      80.1900E0,

      80.6600E0,

      81.2200E0,

      81.6600E0,

      82.1600E0,

      82.5900E0,

      83.1400E0,

      83.5000E0,

      84.0000E0,

      84.4000E0,

      84.8900E0,

      85.2600E0,

      85.7400E0,

      86.0700E0,

      86.5400E0,

      86.8900E0,

      87.3200E0,

      87.6500E0,

      88.1000E0,

      88.4300E0,

      88.8300E0,

      89.1200E0,

      89.5400E0,

      89.8500E0,

      90.2500E0,

      90.5500E0,

      90.9300E0,

      91.2000E0,

      91.5500E0,

      92.2000E0

};

static double kirby2_F0[kirby2_N] = {

    9.65E0,

   10.74E0,

   11.81E0,

   12.88E0,

   14.06E0,

   15.28E0,

   16.63E0,

   18.19E0,

   19.88E0,

   21.84E0,

   24.00E0,

   26.25E0,

   28.86E0,

   31.85E0,

   35.79E0,

   40.18E0,

   44.74E0,

   49.53E0,

   53.94E0,

   58.29E0,

   62.63E0,

   67.03E0,

   71.25E0,

   75.22E0,

   79.33E0,

   83.56E0,

   87.75E0,

   91.93E0,

   96.10E0,

  100.28E0,

  104.46E0,

  108.66E0,

  112.71E0,

  116.88E0,

  121.33E0,

  125.79E0,

  125.79E0,

  128.74E0,

  130.27E0,

  133.33E0,

  134.79E0,

  137.93E0,

  139.33E0,

  142.46E0,

  143.90E0,

  146.91E0,

  148.51E0,

  151.41E0,

  153.17E0,

  155.97E0,

  157.76E0,

  160.56E0,

  162.30E0,

  165.21E0,

  166.90E0,

  169.92E0,

  170.32E0,

  171.54E0,

  173.79E0,

  174.57E0,

  176.25E0,

  177.34E0,

  179.19E0,

  181.02E0,

  182.08E0,

  183.88E0,

  185.75E0,

  186.80E0,

  188.63E0,

  190.45E0,

  191.48E0,

  193.35E0,

  195.22E0,

  196.23E0,

  198.05E0,

  199.97E0,

  201.06E0,

  202.83E0,

  204.69E0,

  205.86E0,

  207.58E0,

  209.50E0,

  210.65E0,

  212.33E0,

  215.43E0,

  217.16E0,

  220.21E0,

  221.98E0,

  225.06E0,

  226.79E0,

  229.92E0,

  231.69E0,

  234.77E0,

  236.60E0,

  239.63E0,

  241.50E0,

  244.48E0,

  246.40E0,

  249.35E0,

  251.32E0,

  254.22E0,

  256.24E0,

  259.11E0,

  261.18E0,

  264.02E0,

  266.13E0,

  268.94E0,

  271.09E0,

  273.87E0,

  276.08E0,

  278.83E0,

  281.08E0,

  283.81E0,

  286.11E0,

  288.81E0,

  291.08E0,

  293.75E0,

  295.99E0,

  298.64E0,

  300.84E0,

  302.02E0,

  303.48E0,

  305.65E0,

  308.27E0,

  310.41E0,

  313.01E0,

  315.12E0,

  317.71E0,

  319.79E0,

  322.36E0,

  324.42E0,

  326.98E0,

  329.01E0,

  331.56E0,

  333.56E0,

  336.10E0,

  338.08E0,

  340.60E0,

  342.57E0,

  345.08E0,

  347.02E0,

  349.52E0,

  351.44E0,

  353.93E0,

  355.83E0,

  358.32E0,

  360.20E0,

  362.67E0,

  364.53E0,

  367.00E0,

  371.30E0

};

static void

kirby2_checksol(const double x[], const double sumsq,

                const double epsrel, const char *sname,

                const char *pname)

{

  size_t i;

  const double sumsq_exact = 3.9050739624E+00;

  const double kirby2_x[kirby2_P] = {

    1.6745063063E+00, -1.3927397867E-01, 2.5961181191E-03,

   -1.7241811870E-03, 2.1664802578E-05 };

  gsl_test_rel(sumsq, sumsq_exact, epsrel, "%s/%s sumsq",

               sname, pname);

  for (i = 0; i < kirby2_P; ++i)

    {

      gsl_test_rel(x[i], kirby2_x[i], epsrel, "%s/%s i=%zu",

                   sname, pname, i);

    }

}

static int

kirby2_f (const gsl_vector * x, void *params, gsl_vector * f)

{

  double b[kirby2_P];

  size_t i;

  for (i = 0; i < kirby2_P; i++)

    {

      b[i] = gsl_vector_get(x, i);

    }

  for (i = 0; i < kirby2_N; i++)

    {

      double x = kirby2_F0[i];

      double y = ((b[0] + x* (b[1]  + x * b[2]))

                  / (1 + x*(b[3]  + x *b[4])));

      gsl_vector_set (f, i, kirby2_F1[i] - y);

    }

  return GSL_SUCCESS;

}

static int

kirby2_df (const gsl_vector * x, void *params, gsl_matrix * df)

{

  double b[kirby2_P];

  size_t i;

  for (i = 0; i < kirby2_P; i++)

    {

      b[i] = gsl_vector_get(x, i);

    }

  for (i = 0; i < kirby2_N; i++)

    {

      double x = kirby2_F0[i];

      double u = (b[0] + x*(b[1] + x*b[2]));

      double v = (1 + x*(b[3] + x*b[4]));

      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*u/(v*v));

      gsl_matrix_set (df, i, 4, x*x*u/(v*v));

    }

  return GSL_SUCCESS;

}

static gsl_multifit_function_fdf kirby2_func =

{

  &kirby2_f,

  &kirby2_df,

  NULL,

  kirby2_N,

  kirby2_P,

  NULL,

  0,

  0

};

static test_fdf_problem kirby2_problem =

{

  "nist-kirby2",

  kirby2_x0,

  kirby2_sigma,

  &kirby2_epsrel,

  kirby2_NTRIES,

  &kirby2_checksol,

  &kirby2_func

};