G4INCL::InverseInterpolationTable Class Reference

Class for interpolating the inverse of a 1-dimensional function. More...

#include <G4INCLInverseInterpolationTable.hh>

Inheritance diagram for G4INCL::InverseInterpolationTable:

Public Member Functions
	InverseInterpolationTable (IFunction1D const &f, const unsigned int nNodes=30)
	InverseInterpolationTable (std::vector< G4double > const &x, std::vector< G4double > const &y)
virtual	~InverseInterpolationTable ()
unsigned int	getNumberOfNodes () const
std::vector< G4double >	getNodeAbscissae () const
std::vector< G4double >	getNodeValues () const
G4double	operator() (const G4double x) const
	Compute the value of the function.
std::string	print () const

---

Detailed Description

Class for interpolating the inverse of a 1-dimensional function.

Definition at line 113 of file G4INCLInverseInterpolationTable.hh.

---

Constructor & Destructor Documentation

G4INCL::InverseInterpolationTable::InverseInterpolationTable	(	IFunction1D const &	f,
		const unsigned int	nNodes = 30
	)

Definition at line 51 of file G4INCLInverseInterpolationTable.cc.

References G4INCL::IFunction1D::getXMaximum(), and G4INCL::IFunction1D::getXMinimum().

                                                                                                      {
// assert(nNodes>2);

    const G4double x0 = f.getXMinimum();
    const G4double x1 = f.getXMaximum();

    // Build the nodes
    G4double last = f(x0);
    InterpolationNode firstNode(last, x0, 0.);
    nodes.push_back(firstNode);
    G4int skippedNodes = 0;
    for(unsigned i = 1; i < nNodes; i++) {
      const G4double xi = x0 + i*(x1-x0)/((G4double)(nNodes-1));
      // Make sure that the x vector is sorted (corresponding to a monotonous
      // function)
      const G4double value = f(xi);
      if(value <= last) {
        ++skippedNodes;
        continue;
      }
      InterpolationNode node(value, xi, 0.);
      nodes.push_back(node);
      last = value;
    }

// assert(nNodes==nodes.size()+skippedNodes);

    // Initialise the "derivative" values
    initDerivatives();
    setFunctionDomain();
  }

G4INCL::InverseInterpolationTable::InverseInterpolationTable	(	std::vector< G4double > const &	x,
		std::vector< G4double > const &	y
	)

Definition at line 83 of file G4INCLInverseInterpolationTable.cc.

                                                                                                                 {
// assert(x.size()==y.size());
    // Assert that the x vector is sorted (corresponding to a monotonous
    // function
// assert(std::adjacent_find(nodes.begin(), nodes.end(), std::greater<InterpolationNode>()) == nodes.end());

    for(unsigned i = 0; i < x.size(); ++i)
      nodes.push_back(InterpolationNode(x.at(i), y.at(i), 0.));

    initDerivatives();
    setFunctionDomain();
  }

virtual G4INCL::InverseInterpolationTable::~InverseInterpolationTable

(

)

[inline, virtual]

Definition at line 117 of file G4INCLInverseInterpolationTable.hh.

{}

---

Member Function Documentation

std::vector<G4double> G4INCL::InverseInterpolationTable::getNodeAbscissae

(

)

const [inline]

Definition at line 121 of file G4INCLInverseInterpolationTable.hh.

References G4INCL::InterpolationNode::getX().

                                                   {
        std::vector<G4double> x(nodes.size());
        std::transform(nodes.begin(), nodes.end(), x.begin(),
            std::mem_fun_ref(&InterpolationNode::getX));
        return x;
      }

std::vector<G4double> G4INCL::InverseInterpolationTable::getNodeValues

(

)

const [inline]

Definition at line 128 of file G4INCLInverseInterpolationTable.hh.

References G4INCL::InterpolationNode::getY().

                                                {
        std::vector<G4double> y(nodes.size());
        std::transform(nodes.begin(), nodes.end(), y.begin(),
            std::mem_fun_ref(&InterpolationNode::getY));
        return y;
      }

unsigned int G4INCL::InverseInterpolationTable::getNumberOfNodes

(

)

const [inline]

Definition at line 119 of file G4INCLInverseInterpolationTable.hh.

{ return nodes.size(); }

G4double G4INCL::InverseInterpolationTable::operator()

(

const G4double

x

)

const [virtual]

Compute the value of the function.

Implements G4INCL::IFunction1D.

Definition at line 117 of file G4INCLInverseInterpolationTable.cc.

                                                                       {
    // Find the relevant interpolation bin
    std::vector<InterpolationNode>::const_iterator iter =
      std::lower_bound(nodes.begin(), nodes.end(), x);

    if(iter==nodes.begin())
      return nodes.front().getY();

    if(iter==nodes.end())
      return nodes.back().getY();

    std::vector<InterpolationNode>::const_iterator previousIter = iter - 1;
    const G4double dx = x - previousIter->getX();
    return previousIter->getY() + previousIter->getYPrime()*dx;
  }

std::string G4INCL::InverseInterpolationTable::print

(

)

const

Definition at line 133 of file G4INCLInverseInterpolationTable.cc.

References CLHEP::detail::n.

Referenced by G4INCL::NuclearDensityFactory::createPCDFTable(), G4INCL::NuclearDensityFactory::createRCDFTable(), and G4INCL::NuclearDensity::NuclearDensity().

                                                   {
    std::string message;
    for(std::vector<InterpolationNode>::const_iterator n=nodes.begin(); n!=nodes.end(); ++n)
      message += n->print();
    return message;
  }

---

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

• G4INCLInverseInterpolationTable.hh
• G4INCLInverseInterpolationTable.cc

---