def solve(self, rhs): """Solve with the system matrix Args: rhs: Right hand side in system Returns: solution to M*x = rhs """ self.update() return la.lu_solve(self.lupiv, rhs)