genie::BLI2DNonUnifGrid Class Reference

#include <BLI2D.h>

Inheritance diagram for genie::BLI2DNonUnifGrid:

Collaboration diagram for genie::BLI2DNonUnifGrid:

Public Member Functions
	BLI2DNonUnifGrid ()
	BLI2DNonUnifGrid (int nx, double xmin, double xmax, int ny, double ymin, double ymax)
	BLI2DNonUnifGrid (int nx, int ny, double *x, double *y, double *z)
bool	AddPoint (double x, double y, double z)
double	Evaluate (double x, double y) const
Public Member Functions inherited from genie::BLI2DGrid
	BLI2DGrid ()
virtual	~BLI2DGrid ()
double	XMin (void) const
double	XMax (void) const
double	YMin (void) const
double	YMax (void) const
double	ZMin (void) const
double	ZMax (void) const

Private Member Functions
void	Init (int nx=0, double xmin=0, double xmax=0, int ny=0, double ymin=0, double ymax=0)

Private Attributes
int	fNFillX
int	fNFillY

Additional Inherited Members
Protected Member Functions inherited from genie::BLI2DGrid
int	IdxZ (int ix, int iy) const
Protected Attributes inherited from genie::BLI2DGrid
int	fNX
int	fNY
int	fNZ
double *	fX
double *	fY
double *	fZ
double	fDX
double	fDY
double	fXmin
double	fXmax
double	fYmin
double	fYmax
double	fZmin
double	fZmax

Detailed Description

Definition at line 99 of file BLI2D.h.

Constructor & Destructor Documentation

BLI2DNonUnifGrid::BLI2DNonUnifGrid

(

)

Definition at line 177 of file BLI2D.cxx.

References Init().

{
  this->Init();
}

BLI2DNonUnifGrid::BLI2DNonUnifGrid	(	int	nx,
		double	xmin,
		double	xmax,
		int	ny,
		double	ymin,
		double	ymax
	)

Definition at line 182 of file BLI2D.cxx.

References Init().

{
  this->Init(nx, xmin, xmax, ny, ymin, ymax);
}

BLI2DNonUnifGrid::BLI2DNonUnifGrid	(	int	nx,
		int	ny,
		double *	x,
		double *	y,
		double *	z
	)

Definition at line 188 of file BLI2D.cxx.

References AddPoint(), genie::BLI2DGrid::IdxZ(), and Init().

{
  double xmin = x[0];
  double xmax = x[nx-1];
  double ymin = y[0];
  double ymax = y[ny-1];

  this->Init(nx, xmin, xmax, ny, ymin, ymax);

  for(int ix=0; ix<nx; ix++) {
    for(int iy=0; iy<ny; iy++) {
      this->AddPoint(x[ix], y[iy], z[this->IdxZ(ix,iy)]);
    }
  }
}

Member Function Documentation

bool BLI2DNonUnifGrid::AddPoint ( double  x, double  y, double  z  )

virtual

Implements genie::BLI2DGrid.

Definition at line 205 of file BLI2D.cxx.

References fNFillX, fNFillY, genie::BLI2DGrid::fNX, genie::BLI2DGrid::fNY, genie::BLI2DGrid::fX, genie::BLI2DGrid::fXmax, genie::BLI2DGrid::fXmin, genie::BLI2DGrid::fY, genie::BLI2DGrid::fYmax, genie::BLI2DGrid::fYmin, genie::BLI2DGrid::fZ, genie::BLI2DGrid::fZmax, genie::BLI2DGrid::fZmin, genie::BLI2DGrid::IdxZ(), LOG, pDEBUG, and pWARN.

Referenced by BLI2DNonUnifGrid().

{

  // check the x,y values' existence before moving anything
  //   if they do, use them
  //   if they don't, add them in the appropriate place
  //
  int xidx = -1;
  int yidx = -1;
  bool changex = true;
  bool changey = true;
  // first, check and see if the x,y values exist
  // NOTE: == should not be used with double values:
  //   instead, a tolerance of ~5 significant digits is checked
  for(int i=0;i<fNFillX;i++)
    {
      if (!(TMath::Abs(x-fX[i])<=.5*.0000001*(TMath::Abs(x)+TMath::Abs(fX[i]))) && x<fX[i])
        {
          // missed one, make sure there aren't too many x and move everything up
          if (fNFillX==fNX)
            {
              LOG("BLI2DNonUnifGrid", pWARN) << "Warning: too many x values, cannot add x = "<<x;
              return false;
            }
          else
            {
              xidx=i;
              changex=true;
              break;
            }
        }
      else
        {
          if (TMath::Abs(x-fX[i])<=.5*.0000001*(TMath::Abs(x)+TMath::Abs(fX[i]))) {
            xidx=i;
            LOG("BLI2DNonUnifGrid", pDEBUG) << "x value found at index "<<i;
            changex=false; break;}
          changex=true;
        }
    }
  if (changex && xidx<0) xidx=fNFillX;
  if (changex) LOG("BLI2DNonUnifGrid", pDEBUG) << "Adding x value at index "<<xidx;

  for(int i=0;i<fNFillY;i++)
    {
      if (!(TMath::Abs(y-fY[i])<=.5*.0000001*(TMath::Abs(y)+TMath::Abs(fY[i]))) && y<fY[i])
        {
          if (fNFillY==fNY)
            {
              LOG("BLI2DNonUnifGrid", pWARN) << "Warning: too many y values, cannot add y = "<<y;
              return false;
            }
          else
            {
              yidx=i;
              changey=true;
              break;
            }
        }
      else
        {
          if (TMath::Abs(y-fY[i])<=.5*.0000001*(TMath::Abs(y)+TMath::Abs(fY[i]))) {
            yidx=i;
            LOG("BLI2DNonUnifGrid", pDEBUG) << "y value found at index "<<i;
            changey=false; break;}
          changey=true;
        }
    }
  if (changey && yidx<0) yidx=fNFillY;
  if (changey) LOG("BLI2DNonUnifGrid", pDEBUG) << "Adding y value at index "<<yidx;

  // make new entries if needed
  if (changex && xidx>=0)
    {
      for (int j=0;j<fNFillX-xidx;j++)
        {
          fX[fNFillX-j]=fX[fNFillX-j-1];
          for (int k=0;k<fNFillY;k++)
            {
              fZ[ this->IdxZ(fNFillX-j,k) ]
                = fZ[ this->IdxZ(fNFillX-j-1,k) ];
            }
        }
      fX[xidx]=x;
      fNFillX++;
      fXmin = TMath::Min(x,fXmin);
      fXmax = TMath::Max(x,fXmax);
    }
  if (changey && yidx>=0)
    {
      for (int j=0;j<fNFillY-yidx;j++)
        {
          fY[fNFillY-j]=fY[fNFillY-j-1];
          for (int k=0;k<fNFillX;k++)
            {
              fZ[ this->IdxZ(k,fNFillY-j) ]
                = fZ[ this->IdxZ(k,fNFillY-j-1) ];
            }
        }
      fY[yidx]=y;
      fNFillY++;
      fYmin = TMath::Min(y,fYmin);
      fYmax = TMath::Max(y,fYmax);
    }

  int iz = this->IdxZ(xidx,yidx);

  fZ[iz] = z;

  fZmin = TMath::Min(z, fZmin);
  fZmax = TMath::Max(z, fZmax);

  LOG("BLI2DNonUnifGrid", pDEBUG)
    << "Added x = " << x << " (ix = " << xidx << ")"
         << " y = " << y << " (iy = " << yidx << ") -> "
         << " z = " << z << " (iz = " << iz   << ")";

  return true;
}

double BLI2DNonUnifGrid::Evaluate ( double  x, double  y  ) const

virtual

Implements genie::BLI2DGrid.

Definition at line 325 of file BLI2D.cxx.

References fNFillX, fNFillY, genie::BLI2DGrid::fX, genie::BLI2DGrid::fXmax, genie::BLI2DGrid::fXmin, genie::BLI2DGrid::fY, genie::BLI2DGrid::fYmax, genie::BLI2DGrid::fYmin, genie::BLI2DGrid::fZ, and genie::BLI2DGrid::IdxZ().

{
  double evalx=TMath::Min(x,fXmax);
  evalx=TMath::Max(evalx,fXmin);
  double evaly=TMath::Min(y,fYmax);
  evaly=TMath::Max(evaly,fYmin);

  int ix_lo  = -2;
  int iy_lo  = -2;
  for (int i=0;i<fNFillX;i++)
    {
      if (evalx<=fX[fNFillX-1-i]) ix_lo=fNFillX-2-i;
      else break;
    }
  for (int i=0;i<fNFillY;i++)
    {
      if (evaly<=fY[fNFillY-1-i]) iy_lo=fNFillY-2-i;
      else break;
    }
  int ix_hi  = ix_lo + 1;
  int iy_hi  = iy_lo + 1;

  // in case x = xmin
  if (ix_lo==-1) {ix_lo++; ix_hi++;}
  if (iy_lo==-1) {iy_lo++; iy_hi++;}
  // in case x = xmax
  if (ix_lo==-2) {ix_lo=fNFillX-2; ix_hi=fNFillX-1;}
  if (iy_lo==-2) {iy_lo=fNFillY-2; iy_hi=fNFillY-1;}
  // if an error occurs
  if (ix_lo<0       || iy_lo<0      ) return 0.;
  if (ix_hi>fNFillX || iy_hi>fNFillY) return 0.;

  double x1  = fX[ix_lo];
  double x2  = fX[ix_hi];
  double y1  = fY[iy_lo];
  double y2  = fY[iy_hi];

  double z11 = fZ[ this->IdxZ(ix_lo,iy_lo) ];
  double z21 = fZ[ this->IdxZ(ix_hi,iy_lo) ];
  double z12 = fZ[ this->IdxZ(ix_lo,iy_hi) ];
  double z22 = fZ[ this->IdxZ(ix_hi,iy_hi) ];

  double z1  = z11 * (x2-evalx)/(x2-x1) + z21 * (evalx-x1)/(x2-x1);
  double z2  = z12 * (x2-evalx)/(x2-x1) + z22 * (evalx-x1)/(x2-x1);
  double z   = z1  * (y2-evaly)/(y2-y1) + z2  * (evaly-y1)/(y2-y1);

  /*
  LOG("BLI2DNonUnifGrid", pINFO) << "x = " << evalx << " -> nearby nodes: " << x1 << ", " << x2;
  LOG("BLI2DNonUnifGrid", pINFO) << "y = " << evaly << " -> nearby nodes: " << y1 << ", " << y2;
  LOG("BLI2DNonUnifGrid", pINFO) << "xmin = " << fXmin << ", xmax = " << fXmax;
  LOG("BLI2DNonUnifGrid", pINFO) << "ymin = " << fYmin << ", ymax = " << fYmax;
  LOG("BLI2DNonUnifGrid", pINFO) << "z11 := z(" << this->IdxZ(ix_lo,iy_lo) << ")  = " << z11;
  LOG("BLI2DNonUnifGrid", pINFO) << "z21 := z(" << this->IdxZ(ix_hi,iy_lo) << ")  = " << z21;
  LOG("BLI2DNonUnifGrid", pINFO) << "z12 := z(" << this->IdxZ(ix_lo,iy_hi) << ")  = " << z12;
  LOG("BLI2DNonUnifGrid", pINFO) << "z22 := z(" << this->IdxZ(ix_hi,iy_hi) << ")  = " << z22;
  LOG("BLI2DNonUnifGrid", pINFO) << "z1 = " << z1 << ", z2 = " << z2;
  LOG("BLI2DNonUnifGrid", pINFO) << "interpolated z(x,y) = " << z;
  */

  return z;
}

void BLI2DNonUnifGrid::Init ( int  nx = 0, double  xmin = 0, double  xmax = 0, int  ny = 0, double  ymin = 0, double  ymax = 0  )

privatevirtual

Implements genie::BLI2DGrid.

Definition at line 387 of file BLI2D.cxx.

References fNFillX, fNFillY, genie::BLI2DGrid::fNX, genie::BLI2DGrid::fNY, genie::BLI2DGrid::fNZ, genie::BLI2DGrid::fX, genie::BLI2DGrid::fXmax, genie::BLI2DGrid::fXmin, genie::BLI2DGrid::fY, genie::BLI2DGrid::fYmax, genie::BLI2DGrid::fYmin, genie::BLI2DGrid::fZ, genie::BLI2DGrid::fZmax, and genie::BLI2DGrid::fZmin.

Referenced by BLI2DNonUnifGrid().

{
  fNX    = 0;
  fNY    = 0;
  fNZ    = 0;
  fNFillX= 0;
  fNFillY= 0;
  fXmin  = 0.;
  fXmax  = 0.;
  fYmin  = 0.;
  fYmax  = 0.;
  fZmin  = std::numeric_limits<double>::max();
  fZmax  = std::numeric_limits<double>::min();
  fX     = 0;
  fY     = 0;
  fZ     = 0;

  if(nx>1 && ny>1) {
    fNX = nx;
    fNY = ny;
    fNZ = nx * ny;

    fXmin = xmin;
    fXmax = xmax;
    fYmin = ymin;
    fYmax = ymax;
    fZmin = std::numeric_limits<double>::max();
    fZmax = std::numeric_limits<double>::min();

    fX = new double[fNX];
    fY = new double[fNY];
    fZ = new double[fNZ];

    for(int i=0; i<fNX; i++) { fX[i] = 0.; }
    for(int i=0; i<fNY; i++) { fY[i] = 0.; }
    for(int i=0; i<fNZ; i++) { fZ[i] = 0.; }
  }
}

Member Data Documentation

int genie::BLI2DNonUnifGrid::fNFillX

private

Definition at line 116 of file BLI2D.h.

Referenced by AddPoint(), Evaluate(), and Init().

int genie::BLI2DNonUnifGrid::fNFillY

private

Definition at line 117 of file BLI2D.h.

Referenced by AddPoint(), Evaluate(), and Init().

---

The documentation for this class was generated from the following files:

• BLI2D.h
• BLI2D.cxx