pub fn solve_implicit_system_multi(
mat: &[Vec<f64>],
rhs_cols: &[Vec<f64>],
) -> OptimizeResult<Vec<Vec<f64>>>Expand description
Solve the system mat * X = rhs_matrix where rhs_matrix has k columns.
Returns the solution matrix (n × k) in row-major form.