SolveFirstDerivInterp1D< double > Struct Template Reference

Collaboration diagram for SolveFirstDerivInterp1D< double >:

Static Public Member Functions
template<class CT >
static void	apply (const CT &data, CT &p, int N, double *bcLower, double *bcUpper)

Detailed Description

template<>
struct SolveFirstDerivInterp1D< double >

Definition at line 134 of file BsplineOneDimSolvers.h.

Member Function Documentation

apply()

static void apply ( const CT &  data, CT &  p, int  N, double *  bcLower, double *  bcUpper  )

inlinestatic

Definition at line 137 of file BsplineOneDimSolvers.h.

References qmcplusplus::Units::force::N.

  {
    const double ratio = 0.25;
    CT d(N + 2), mu(N + 2, ratio);
    for (int i = 1; i <= N; i++)
      d[i] = 1.5 * data[i - 1];
    double al = 0.25 * bcLower[0];
    double bl = 0.25 * bcLower[1];
    double cl = 0.25 * bcLower[2];
    double dl = 1.5 * bcLower[3];
    double ar = 0.25 * bcUpper[0];
    double br = 0.25 * bcUpper[1];
    double cr = 0.25 * bcUpper[2];
    double dr = 1.5 * bcUpper[3];
    // First, eliminate leading coefficients
    double alInv = 1.0 / al;
    bl *= alInv;
    cl *= alInv;
    dl *= alInv;
    d[0]  = dl;
    mu[0] = bl;
    mu[1] = ratio - ratio * cl;
    for (int row = 1; row <= N; row++)
    {
      double diag    = 1.0 - mu[row - 1] * ratio;
      double diagInv = 1.0 / diag;
      mu[row] *= diagInv;
      d[row] = diagInv * (d[row] - ratio * d[row - 1]);
    }
    br -= ar * mu[N - 1];
    dr -= ar * d[N - 1];
    cr -= br * mu[N];
    dr -= br * d[N];
    p[N + 1] = dr / cr;
    // Now go back upward, back substituting
    for (int row = N; row >= 1; row--)
      p[row] = d[row] - mu[row] * p[row + 1];
    // And do 0th row
    p[0] = dl - bl * p[1] - cl * p[2];
  }

---

The documentation for this struct was generated from the following file:

• /home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/Numerics/BsplineOneDimSolvers.h