Convert to scipy.sparse.csc_matrix for PyPardiso

#74
loganbvh · Mar 5, 2024 · c09cafe · c09cafe
1 parent ac8b2d9
commit c09cafe
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/tdgl/finite_volume/operators.py b/tdgl/finite_volume/operators.py
@@ -295,9 +295,12 @@ def build_operators(self) -> None:
             self.edge_directions = cupy.array(self.edge_directions)
         if self.sparse_solver is SparseSolver.CUPY:
             assert cupy is not None
+            # https://github.com/loganbvh/py-tdgl/issues/74
+            # https://github.com/haasad/PyPardiso/issues/68
             self.mu_laplacian = csc_matrix(self.mu_laplacian)
             self.mu_laplacian_lu = factorized(self.mu_laplacian)
         elif self.sparse_solver is SparseSolver.PARDISO:
+            self.mu_laplacian = sp.csc_matrix(self.mu_laplacian)
             self.mu_laplacian_lu = None
         else:
             use_umfpack = self.sparse_solver is SparseSolver.UMFPACK

diff --git a/tdgl/version.py b/tdgl/version.py
@@ -1,4 +1,4 @@
-__version_info__ = (0, 8, 1)
+__version_info__ = (0, 8, 2)
 __version__ = ".".join(map(str, __version_info__))