PFuncTest.cpp

Go to the documentation of this file.

// ============================================================================
// Include files
// ============================================================================
// STD & STL
// ============================================================================
#include <cmath>
#include <cstdio>
#include <vector>
// ============================================================================
// GaudiMath
// ============================================================================
#include "GaudiMath/Adapters.h"
#include "GaudiMath/GaudiMath.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 ) );
    }
  }
}