InterpTest.cpp

Go to the documentation of this file.

// $Id: Interp.cpp,v 1.3 2006/12/08 12:15:24 hmd Exp $
// ============================================================================
// Include files
// ============================================================================
// STD& STL
// ============================================================================
#include <cmath>
#include <iostream>
#include <cstdio>
#include "GaudiKernel/SystemOfUnits.h"
#include "GaudiMath/GaudiMath.h"

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

int main()
{

  const GaudiMath::Function& mysin = GaudiMath::SimpleFun( sin );
  const GaudiMath::Function& mycos = GaudiMath::SimpleFun( cos );

  std::vector<double> vctx ;
  std::vector<double> vcty ;

  for ( int i = 0 ; i < 20 ; ++i )
  {
    double x = i * 3.141592/20.0 ;
    double y = mysin ( x ) ;
    vctx.push_back   ( x ) ;
    vcty.push_back   ( y ) ;
  }

  using GaudiMath::Interpolation::Cspline_Periodic;

  const GaudiMath::Constant minusOne (-1.0, 1);
  const GaudiMath::Constant plusOne (1.0, 1);

  const GaudiMath::Spline i1(vctx, vcty, Cspline_Periodic);

  const GaudiMath::SplineDeriv i2(vctx, vcty, Cspline_Periodic);

  const GaudiMath::Function& i3 =
      GaudiMath::SplineDeriv2(vctx, vcty, Cspline_Periodic) * minusOne;

  const GaudiMath::Function& i4 = plusOne +
      GaudiMath::SplineInteg(vctx, vcty, Cspline_Periodic, 0.0) * minusOne;

  for ( int i = -10 ; i < 30 ; ++i ) {
    double x = i * 3.141592/30.0 ;
    printf( "x=%+7.1f; Sin-i1=%+14.9f; Cos-i2=%+14.9f; Sin-i3=%+14.9f; Cos-i4=%+14.9f;\n",
            x / Gaudi::Units::degree,
            mysin ( x ) - i1( x ),
            mycos ( x ) - i2( x ),
            mysin ( x ) - i3( x ),
            mycos ( x ) - i4( x ));
  }

  exit(0);

}