EDepSim::Cubic Class Reference

#include <EDepSimInterpolator.hh>

Public Member Functions
	Cubic ()
double	interpolate (const double *point, const std::vector< std::vector< std::vector< double >>> &g, const double *delta, const double *offset) const
double	interpolate (double x, double y, double z, const std::vector< std::vector< std::vector< double >>> &g, double hx, double hy, double hz, double xo, double yo, double zo) const

Private Member Functions
double	conv_kernel (double s) const

Detailed Description

Definition at line 44 of file EDepSimInterpolator.hh.

Constructor & Destructor Documentation

EDepSim::Cubic::Cubic

(

)

Definition at line 34 of file EDepSimInterpolator.cc.

{
}

Member Function Documentation

double EDepSim::Cubic::conv_kernel ( double  s ) const

private

Definition at line 84 of file EDepSimInterpolator.cc.

{
    double v = 0;
    double z = std::abs(s);

    if(0 <= z && z < 1)
        v = 1 + (0.5) * z * z * (3 * z - 5);

    else if(1 < z && z < 2)
        v = 2 - z * (4 + 0.5 * z * (z - 5));

    return v;
}

double EDepSim::Cubic::interpolate	(	const double *	point,
		const std::vector< std::vector< std::vector< double >>> &	g,
		const double *	delta,
		const double *	offset
	)		const

Definition at line 38 of file EDepSimInterpolator.cc.

{
    return interpolate(point[0], point[1], point[2], g, delta[0], delta[1], delta[2], offset[0], offset[1], offset[2]);
}

double EDepSim::Cubic::interpolate	(	double	x,
		double	y,
		double	z,
		const std::vector< std::vector< std::vector< double >>> &	g,
		double	hx,
		double	hy,
		double	hz,
		double	xo,
		double	yo,
		double	zo
	)		const

Definition at line 44 of file EDepSimInterpolator.cc.

{
    double v = 0;

    const int xsize = g.size();
    const int ysize = g[0].size();
    const int zsize = g[0][0].size();

    const double xp = (x - xo) / hx;
    const double yp = (y - yo) / hy;
    const double zp = (z - zo) / hz;

    const int x_idx = std::floor(xp);
    const int y_idx = std::floor(yp);
    const int z_idx = std::floor(zp);

    for(auto i : {x_idx - 1, x_idx, x_idx + 1, x_idx + 2})
    {
        for(auto j : {y_idx - 1, y_idx, y_idx + 1, y_idx + 2})
        {
            for(auto k : {z_idx - 1, z_idx, z_idx + 1, z_idx + 2})
            {
                if(i < 0 || j < 0 || k < 0)
                    continue;

                if(i >= xsize || j >= ysize || k >= zsize)
                    continue;

                const double x_c = conv_kernel(xp - i);
                const double y_c = conv_kernel(yp - j);
                const double z_c = conv_kernel(zp - k);
                v += g[i][j][k] * x_c * y_c * z_c;
            }
        }
    }

    return v;
}

---

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

• edep-sim/src/EDepSimInterpolator.hh
• edep-sim/src/EDepSimInterpolator.cc