docs: def factorize(self,A):
          """
          Factorizes A.
  
          Parameters
          ----------
          A : matrix
             For symmetric systems, should contain only lower diagonal part.
          """
  
          A = coo_matrix(A)
  
          self.mumps.set_centralized_assembled_values(A.data)
          self.mumps.run(job=2)