docs: def solve(self,b):
          """
          Solves system Ax=b.
          
          Parameters
          ----------
          b : ndarray
          
          Returns
          -------
          x : ndarray
          """
  
          x = b.copy()
          self.mumps.set_rhs(x)
          self.mumps.run(job=3)
  
          return x