Function apply_at 

fn apply_at(
    u: &[f64],
    brain_mask: &[f64],
    d_kernel: &[f64],
    bg_mask: &[f64],
    nx: usize,
    ny: usize,
    nz: usize,
) -> Vec<f64>

Apply A^T = M_bg * D * W