diff --git a/tdgl/finite_volume/operators.py b/tdgl/finite_volume/operators.py
index 7272439..91419af 100644
--- a/tdgl/finite_volume/operators.py
+++ b/tdgl/finite_volume/operators.py
@@ -353,7 +353,7 @@ def set_link_exponents(self, link_exponents: np.ndarray) -> None:
             )
             if self.sparse_solver is SparseSolver.CUPY:
                 self.psi_gradient = csr_matrix(self.psi_gradient)
-                self.psi_laplacian = csc_matrix(self.psi_laplacian)
+                self.psi_laplacian = csr_matrix(self.psi_laplacian)
                 self.gradient_weights = cupy.asarray(self.gradient_weights)
                 self.laplacian_weights = cupy.asarray(self.laplacian_weights)
             return