Struct Mg 

struct Mg<T> {
    levels: Vec<Level<T>>,
    depth: usize,
    aux3: [T; 3],
    aux2: T,
    tol: T,
    n1: usize,
    n2: usize,
}

Fields

levels: Vec<Level<T>>

depth: usize

aux3: [T; 3]

aux3: 1/dx², 1/dy², 1/dz² of the finest grid, used on every level

aux2: T

aux2: 1 / (2 Σ aux3)

tol: T

n1: usize

n2: usize

Implementations

impl<T: Real> Mg<T>

fn run( ifreq: &[f64], mask_in: &[u8], grid: &Grid, params: &LbvParams, progress: impl FnMut(usize, usize), ) -> (Vec<f64>, Vec<u8>)

fn fmg( &mut self, d: usize, progress: &mut impl FnMut(usize, usize), steps: usize, )

MG::FMG

fn v_cycle(&mut self, d: usize, f: Rhs)

MG::V_cycle on level d with right-hand side f

fn relax(&mut self, d: usize, f: Rhs, n: usize)

MG::relax: lexicographic Gauss–Seidel on the interior, n sweeps, or on the coarsest level until the summed absolute update is at most tol.

fn update_res(&mut self, d: usize, f: Rhs)

MG::updateRes: res_d = f − L v on the interior. With f = res_d this overwrites the right-hand side, as in the original.

fn f2c(&mut self, d: usize)

MG::F2C: inject the residual of level d onto level d + 1 (the even-indexed voxels).

fn c2f(&mut self, d: usize, cor: bool)

MG::C2F: prolong level d onto level d − 1, replacing (cor false) or adding the correction scaled by prefac = 4^−d (cor true).