SolvePeriodicInterp1D< double > Struct Template Reference

Collaboration diagram for SolvePeriodicInterp1D< double >:

Static Public Member Functions
static void	apply (const std::vector< double > &data, std::vector< double > &p)
template<typename CT >
static void	apply (const CT &data, CT &p, int N)

Detailed Description

template<>
struct SolvePeriodicInterp1D< double >

Definition at line 64 of file BsplineOneDimSolvers.h.

Member Function Documentation

apply() [1/2]

static void apply ( const std::vector< double > &  data, std::vector< double > &  p  )

inlinestatic

Definition at line 66 of file BsplineOneDimSolvers.h.

References FuncSolvePeriodicInterp1D().

  {
    FuncSolvePeriodicInterp1D(data, p);
  }

apply() [2/2]

static void apply ( const CT &  data, CT &  p, int  N  )

inlinestatic

Definition at line 72 of file BsplineOneDimSolvers.h.

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

  {
    const double ratio = 0.25;
    CT d(N), gamma(N), mu(N);
    for (int i = 0; i < N; i++)
      d[i] = 1.5 * data[i];
    // First, eliminate leading coefficients
    gamma[0]     = ratio;
    mu[0]        = ratio;
    mu[N - 1]    = ratio;
    gamma[N - 1] = 1.0;
    for (int row = 1; row < (N - 1); row++)
    {
      double diag    = 1.0 - mu[row - 1] * ratio;
      double diagInv = 1.0 / diag;
      gamma[row]     = -ratio * gamma[row - 1] * diagInv;
      mu[row]        = diagInv * ratio;
      d[row]         = diagInv * (d[row] - ratio * d[row - 1]);
      // Last row
      d[N - 1] -= mu[N - 1] * d[row - 1];
      gamma[N - 1] -= mu[N - 1] * gamma[row - 1];
      mu[N - 1] = -mu[N - 1] * mu[row - 1];
    }
    // Last row:  gamma(N-1) hold diagonal element
    mu[N - 1] += ratio;
    gamma[N - 1] -= mu[N - 1] * (mu[N - 2] + gamma[N - 2]);
    d[N - 1] -= mu[N - 1] * d[N - 2];
    p[N] = d[N - 1] / gamma[N - 1];
    // Now go back upward, back substituting
    for (int row = N - 2; row >= 0; row--)
      p[row + 1] = d[row] - mu[row] * p[row + 2] - gamma[row] * p[N];
  }

---

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

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