Function solve 

pub fn solve(
    operator: &DMatrix<Complex>,
    rhs: DMatrix<Complex>,
) -> Result<SolveResidual>

Solve A X = B with pivoted LU; B is a matrix of right-hand sides.