pub fn contract_vxc_0( mat_a: &mut MatrixFull<f64>, mat_b: &MatrixFullSlice<'_, f64>, slice_c: &[f64], scaling_factor: Option<f64>)