Inheritance diagram for CenterGrid:

Collaboration diagram for CenterGrid:

Public Member Functions
GridType	Type () override
	Returns the type of the grid (i.e. linear, optimal, etc) More...

int	ReverseMap (double x) override
	Returns the index of the nearest point below r. More...

void	Write (IOSectionClass &out) override

void	Read (IOSectionClass &in) override

void	Init (double start, double end, double ratio, int numPoints)

Public Member Functions inherited from Grid
double	operator() (int i) const
	The i'th point in the grid. More...

double *	data ()

Array< double, 1 > &	Points ()

virtual	~Grid ()

Private Attributes
double	a

double	aInv

double	b

double	bInv

double	center

int	HalfPoints

bool	Odd

double	EvenHalf

int	OddOne

Additional Inherited Members
Public Attributes inherited from Grid
double	Start
	First and last grid points. More...

double	End

int	NumPoints
	Number of points in the grid. More...

Protected Attributes inherited from Grid
Array< double, 1 >	grid
	Contains the grid points. More...

Detailed Description

Definition at line 442 of file Grid.h.

Member Function Documentation

◆ Init()

void Init	(	double	start,
		double	end,
		double	ratio,
		int	numPoints
	)

inline

Definition at line 464 of file Grid.h.

References a, qmcplusplus::abs(), aInv, b, bInv, center, Grid::End, EvenHalf, Grid::grid, HalfPoints, qmcplusplus::log(), Grid::NumPoints, Odd, OddOne, Array< T, D, ALLOC >::resize(), sign(), and Grid::Start.

   {
     assert(ratio > 1.0);
     Start      = start;
     End        = end;
     center     = 0.5 * (start + end);
     NumPoints  = numPoints;
     HalfPoints = numPoints / 2;
     Odd        = ((numPoints % 2) == 1);
     b          = log(ratio) / (double)(HalfPoints - 1);
     bInv       = 1.0 / b;
     grid.resize(numPoints);
     if (Odd)
     {
       EvenHalf = 0.0;
       OddOne   = 1;
       a        = 0.5 * (end - start) / expm1(b * HalfPoints);
       aInv     = 1.0 / a;
       for (int i = -HalfPoints; i <= HalfPoints; i++)
       {
         double sign          = (i < 0) ? -1.0 : 1.0;
         grid(i + HalfPoints) = sign * a * expm1(b * std::abs(i)) + center;
       }
     }
     else
     {
       EvenHalf = 0.5;
       OddOne   = 0;
       a        = 0.5 * (end - start) / expm1(b * (-0.5 + HalfPoints));
       aInv     = 1.0 / a;
       for (int i = -HalfPoints; i < HalfPoints; i++)
       {
         double sign          = (i < 0) ? -1.0 : 1.0;
         grid(i + HalfPoints) = sign * a * expm1(b * std::abs(0.5 + i)) + center;
       }
     }
   }

◆ Read()

void Read ( IOSectionClass & in )

inlineoverridevirtual

Implements Grid.

Definition at line 463 of file Grid.h.

463 {}

◆ ReverseMap()

int ReverseMap ( double r )

inlineoverridevirtual

Returns the index of the nearest point below r.

Implements Grid.

Definition at line 456 of file Grid.h.

References qmcplusplus::abs(), aInv, bInv, center, EvenHalf, qmcplusplus::floor(), and HalfPoints.

   {
     x -= center;
     double index = copysign(log1p(std::abs(x) * aInv) * bInv, x);
     return (int)floor(HalfPoints + index - EvenHalf);
   }

◆ Type()

GridType Type ( )

inlineoverridevirtual

Returns the type of the grid (i.e. linear, optimal, etc)

Implements Grid.

Definition at line 454 of file Grid.h.

References CENTER.

454 { return CENTER; }

CENTER

Definition: Grid.h:38

◆ Write()

void Write ( IOSectionClass & out )

inlineoverridevirtual

Implements Grid.

Definition at line 462 of file Grid.h.

462 {}