diagonal dummy identity preconditioner

Author: otvam

pypeec/lib_matrix/matrix_factorization.py

@@ -117,18 +117,13 @@ def factor(rhs):
     return factor
 
 
-def _get_fact_diagonal(mat):
+def _get_fact_diagonal():
     """
     Dummy factorization with a diagonal matrix.
     """
 
-    # get the inverse of the diagonal
-    inv = IMPORTLIB.diags(1/mat.diagonal(), format="csc")
-
-    # apply the diagonal inverse
     def factor(rhs):
-        sol = inv*rhs
-        return sol
+        return rhs
 
     return factor
 
@@ -167,7 +162,7 @@ def factor_dummy(rhs):
     elif LIBRARY == "PyAMG":
         factor = _get_fact_pyamg(mat)
     elif LIBRARY == "Diagonal":
-        factor = _get_fact_diagonal(mat)
+        factor = _get_fact_diagonal()
     else:
         raise ValueError("invalid matrix factorization library")