pub fn numerical_density(
    grid: &Grids,
    mol: &Molecule,
    dm: &mut [MatrixFull<f64>; 2]
) -> [f64; 2]