Struct InterpolatedKernel 

pub struct InterpolatedKernel<T> {
    pub segments_x: Vec<T>,
    pub segments_y: Vec<T>,
    pub domain_x: (T, T),
    pub domain_y: (T, T),
    pub interpolators: DTensor<Interpolate2D<T>, 2>,
    pub n_cells_x: usize,
    pub n_cells_y: usize,
}

2D interpolation kernel for efficient evaluation at arbitrary points

This structure manages a grid of Interpolate2D objects for piecewise polynomial interpolation across the entire kernel domain.

Fields

segments_x: Vec<T>

X-axis segment boundaries (from SVEHints)

segments_y: Vec<T>

Y-axis segment boundaries (from SVEHints)

domain_x: (T, T)

Domain boundaries

domain_y: (T, T)

interpolators: DTensor<Interpolate2D<T>, 2>

Interpolators for each cell ((segments_x.len()-1) × (segments_y.len()-1))

n_cells_x: usize

Number of cells (for efficiency)

n_cells_y: usize

Implementations

impl<T: CustomNumeric + Debug + Clone + 'static> InterpolatedKernel<T>

pub fn from_kernel_and_segments<K: CentrosymmKernel + KernelProperties>( kernel: &K, segments_x: Vec<T>, segments_y: Vec<T>, gauss_per_cell: usize, symmetry: SymmetryType, ) -> Self

Create InterpolatedKernel from kernel and segments

This function creates a grid of Interpolate2D objects, one for each cell defined by the segments. Each cell uses independent Gauss rules and kernel evaluation for optimal interpolation.

Arguments

• kernel - Kernel to interpolate
• segments_x - X-axis segment boundaries
• segments_y - Y-axis segment boundaries
• gauss_per_cell - Number of Gauss points per cell (e.g., 4 for degree 3)
• symmetry - Symmetry type for kernel evaluation

Returns

New InterpolatedKernel instance

pub fn find_cell(&self, x: T, y: T) -> Option<(usize, usize)>

Find the cell containing point (x, y) using binary search

Arguments

• x - x-coordinate
• y - y-coordinate

Returns

Some((i, j)) if point is in domain, None otherwise

pub fn evaluate(&self, x: T, y: T) -> T

Evaluate interpolated kernel at point (x, y)

Arguments

• x - x-coordinate
• y - y-coordinate

Returns

Interpolated kernel value at (x, y)

Panics

Panics if (x, y) is outside the interpolation domain

pub fn domain(&self) -> ((T, T), (T, T))

Get domain boundaries

pub fn n_cells_x(&self) -> usize

Get number of cells in x direction

pub fn n_cells_y(&self) -> usize

Get number of cells in y direction

Trait Implementations

impl<T: Clone> Clone for InterpolatedKernel<T>

fn clone(&self) -> InterpolatedKernel<T>

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<T: Debug> Debug for InterpolatedKernel<T>

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.