// @(#)root/minuit2:$Id$
// Authors: M. Winkler, F. James, L. Moneta, A. Zsenei   2003-2005  

/**********************************************************************
 *                                                                    *
 * Copyright (c) 2005 LCG ROOT Math team,  CERN/PH-SFT                *
 *                                                                    *
 **********************************************************************/

#include "Minuit2/FCNBase.h"
#include "Minuit2/FunctionMinimum.h"
#include "Minuit2/MnMigrad.h"
#include "Minuit2/MnUserParameters.h"
#include "Minuit2/MnPrint.h"
namespace ROOT {

   namespace Minuit2 {


class Quad8F : public FCNBase {

public:

  double operator()(const std::vector<double>& par) const {

    double x = par[0];
    double y = par[1];
    double z = par[2];
    double w = par[3];
    double x0 = par[4];
    double y0 = par[5];
    double z0 = par[6];
    double w0 = par[7];


    return ( (1./70.)*(21*x*x + 20*y*y + 19*z*z - 14*x*z - 20*y*z) + w*w + (1./70.)*(21*x0*x0 + 20*y0*y0 + 19*z0*z0 - 14*x0*z0 - 20*y0*z0) + w0*w0 );
  }

  double Up() const {return 1.;}

private:

};

  }  // namespace Minuit2

}  // namespace ROOT

void FaiQuad8F()
{
  ROOT::Minuit2::Quad8F fcn;
  
//   PentiumTimer stopwatch;
//   stopwatch.start();

//   long long int start = stopwatch.lap().ticks();
//   long long int stop = stopwatch.lap().ticks();
//   std::cout<<"stop-start: "<<stop - start<<std::endl;
//   start = stopwatch.lap().ticks();
  {
  //test constructor
   ROOT::Minuit2::MnUserParameters upar;
  upar.Add("x", 10., 0.1);
  upar.Add("y", 10., 0.1);
  upar.Add("z", 10., 0.1);
  upar.Add("w", 10., 0.1);
  upar.Add("x0", 10., 0.1);
  upar.Add("y0", 10., 0.1);
  upar.Add("z0", 10., 0.1);
  upar.Add("w0", 10000., 0.1);

   ROOT::Minuit2::MnMigrad migrad(fcn, upar);
   ROOT::Minuit2::FunctionMinimum min = migrad();
  std::cout<<"minimum: "<<min<<std::endl;
  }
}