2DoubleFuncTest.cpp

Go to the documentation of this file.

// ============================================================================
// Include files 
// ============================================================================
// STD & STL 
// ============================================================================
#include <math.h>
#include <stdio.h>
// ============================================================================
// GaudiMath
// ============================================================================
#include "GaudiMath/GaudiMath.h"
#include "GaudiMath/Adapters.h"
// ============================================================================

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

int main()
{
  const GaudiMath::Function& mypow = GaudiMath::adapter ( pow );
  Genfun::Argument arg(2);
  for ( double x = 0; x < 10; ++x)
  {
    for (double y = 0; y < 10; ++y)
    {
      arg[0] = x;
      arg[1] = y;
      printf ("x=%6.5f y=%6.5f mypow=%+.10f pow=%+.10f pow-mypow=%+.10f \n",
              x, y, mypow(arg), pow(x, y), 
              pow(x, y) - mypow(arg) );
    }
  }
  
  exit(0);
  
}