fn evaluate_response( scf_data: &SCF, rimo: &mut Vec<RIFull<f64>>, freq: f64) -> Result<MatrixFull<f64>>