gtestBLI2DUnifGrid.cxx File Reference

#include <TMath.h>
#include <TFile.h>
#include <TNtuple.h>
#include "Framework/Messenger/Messenger.h"
#include "Framework/Numerical/BLI2D.h"
#include "Framework/Numerical/RandomGen.h"

Include dependency graph for gtestBLI2DUnifGrid.cxx:

Go to the source code of this file.

Functions
double	func (double x, double y)
int	main (int, char **)

Function Documentation

func()

double func	(	double	x,
		double	y )

Definition at line 81 of file gtestBLI2DUnifGrid.cxx.

{
  //return x*y+10*x-7*y+1;
  return TMath::Sin(x)/x * TMath::Sin(y)/y;
}

Referenced by genie::DISXSec::CacheFreeNucleonXSec(), genie::DMDISXSec::CacheFreeNucleonXSec(), genie::ReinSehgalRESXSecWithCache::CacheResExcitationXSec(), genie::ReinSehgalRESXSecWithCacheFast::CacheResExcitationXSec(), genie::SPPXSecWithCache::CacheResExcitationXSec(), genie::utils::kinematics::COHImportanceSamplingEnvelope(), genie::utils::kinematics::DISImportanceSamplingEnvelope(), genie::AlamSimoAtharVacasSKXSec::Integrate(), genie::CEvNSXSec::Integrate(), genie::COHDNuXSec::Integrate(), genie::COHXSec::Integrate(), genie::COHXSecAR::Integrate(), genie::DFRXSec::Integrate(), genie::DISXSec::Integrate(), genie::DMDISXSec::Integrate(), genie::DMElectronXSec::Integrate(), genie::DMELXSec::Integrate(), genie::HEDISXSec::Integrate(), genie::HELeptonXSec::Integrate(), genie::IMDXSec::Integrate(), genie::MECXSec::Integrate(), genie::NewQELXSec::Integrate(), genie::NuElectronXSec::Integrate(), genie::QELXSec::Integrate(), genie::ReinSehgalRESXSec::Integrate(), genie::ReinSehgalRESXSecFast::Integrate(), genie::RESXSec::Integrate(), genie::SmithMonizQELCCXSec::Integrate(), genie::SPPXSec::Integrate(), main(), genie::utils::kinematics::RESImportanceSamplingEnvelope(), and genie::NievesQELCCPXSec::vcr().

main()

int main	(	int	,
		char **	)

Definition at line 31 of file gtestBLI2DUnifGrid.cxx.

{
  int npoints=10000;
 
  int nx = 100;
  int ny = 100;
 
  double xmin = -5;
  double xmax =  5;
  double ymin = -5;
  double ymax =  5;
 
  double dx = (xmax-xmin)/(nx-1);
  double dy = (ymax-ymin)/(ny-1);
 
  BLI2DUnifGrid biln(nx,xmin,xmax,ny,ymin,ymax);
 
  for(int ix=0; ix<nx; ix++) {
      double x = xmin + ix * dx;
      for(int iy=0; iy<ny; iy++) {
         double y = ymin + iy * dy;      
         double z = func(x,y);
         biln.AddPoint(x,y,z);
      }
  }
 
  RandomGen * rnd = RandomGen::Instance();
 
  TNtuple * nt = new TNtuple("nt","billinear interpolation validation","x:y:ztrue:zeval");
 
  for(int ip=0; ip<npoints; ip++) {
   double rx = rnd->RndGen().Uniform();
   double ry = rnd->RndGen().Uniform();
   double x  = xmin + (xmax-xmin)*rx;
   double y  = ymin + (ymax-ymin)*ry;
   double zt = func(x,y);
   double ze = biln.Evaluate(x,y);
 
   nt->Fill(x,y,zt,ze);
  }
 
  TFile f("./bli2dug.root","recreate");
  nt->Write();
  f.Close();
 
  LOG("test", pINFO)  << "Done!";
  return 0;
}

References genie::BLI2DUnifGrid::AddPoint(), genie::BLI2DUnifGrid::Evaluate(), func(), genie::RandomGen::Instance(), LOG, pINFO, and genie::RandomGen::RndGen().