PFuncTest.cpp

Go to the documentation of this file.

// ============================================================================
// Include files
// ============================================================================
// STD & STL
// ============================================================================
#include <cmath>
#include <cstdio>
#include <vector>
// ============================================================================
// GaudiMath
// ============================================================================
#include "GaudiMath/GaudiMath.h"
#include "GaudiMath/Adapters.h"
// ============================================================================

// ============================================================================
// ============================================================================

double fun ( const std::vector<double>& x)
{
  return 3 * x[0] * x[0] + 5 * x[1] * x[1] + 8;
}

int main()
{
  const GaudiMath::Function& myfun = GaudiMath::adapter( 2 , &fun ) ;
  Genfun::Argument arg(myfun.dimensionality());

  std::string format( "x[0]=%6.5f x[1]=%6.5f myfun=%+.10f");
  format +=" fun=%+.10f fun-myfun=%+.10f \n" ;
  std::vector<double> x(myfun.dimensionality());

  for ( double a = 0; a < 3; ++a) {
      for (double b = 0; b < 3; ++b) {
          arg[0] = x[0] = a;
          arg[1] = x[1] = b;

          printf ( format.c_str() ,
                   x[0], x[1], myfun(arg), fun(x),
                   fun(x) - myfun(arg) );
      }
  }
}