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