pub fn linear_residuals( x: &MatrixF64, y: &VectorF64, c: &VectorF64, r: &mut VectorF64 ) -> Result<(), Value>