pub fn gen_c(scf_data: &SCF, n_ele: usize) -> MatrixFull<f64>