GSLAdaptersTest.cpp

Go to the documentation of this file.

// ============================================================================
// Include file
// ============================================================================
// STD & STL
// ============================================================================
#include <cstdio>
#include <iostream>
// ============================================================================
// GaudiGSL/GaudiMath
// ============================================================================
#include "GaudiMath/GaudiMath.h"
// ============================================================================
// GSL
// ============================================================================
#include "gsl/gsl_sf_bessel.h"
// ============================================================================

int main()
{

  const GaudiMath::Function& I0 = GaudiMath::adapter( gsl_sf_bessel_I0 );
  const GaudiMath::Function& I1 = GaudiMath::adapter( gsl_sf_bessel_I1_e );

  for ( int i = 0; i < 20; ++i ) {
    const double x  = double( i ) * 1 / 20.0;
    const double x0 = I0( x );
    const double x1 = I1( x );
    printf( " x = %+8.4f I0=%15.8f I1=%15.8f \n", x, x0, x1 );
  };
}