Function lap1_forward 

fn lap1_forward(
    f: &[f64],
    mask: &[u8],
    l: usize,
    a: usize,
    n_axis: usize,
    coord: usize,
    _n_total: usize,
) -> f64

Forward one-sided second derivative (matching MATLAB’s fd/ff functions)