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