pub fn get_cube(scf_data: &SCF, n_p: usize) -> MatrixFull<f64>