fn compute_laplacian_eigenvalues( nx: usize, ny: usize, nz: usize, vsx: f64, vsy: f64, vsz: f64, ) -> Vec<f64>
Compute discrete Laplacian eigenvalues for inverse Laplacian