fix: region funct type incorrect causing

interpolator to not use lsqr solution
Loop3D · Feb 14, 2022 · a5d27f1 · a5d27f1
1 parent ad8e9e7
commit a5d27f1
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/LoopStructural/interpolators/_discrete_interpolator.py b/LoopStructural/interpolators/_discrete_interpolator.py
@@ -64,7 +64,7 @@ def nx(self):
 
     @property
     def region(self):
-        return self.region_function(self.support.nodes)
+        return self.region_function(self.support.nodes).astype(bool)
 
     @property
     def region_map(self):
@@ -598,7 +598,9 @@ def _solve(self, solver="cg", **kwargs):
             self.c[self.region] = self._solve_chol(A, B)
         if solver == "lu":
             logger.info("Solving using scipy LU")
+            print(self.region)
             self.c[self.region] = self._solve_lu(A, B)
+            print(self.c)
         if solver == "pyamg":
             try:
                 logger.info("Solving with pyamg solve")