fix: adding check to see if a feature is valid.

If something doesn't interpolate it will break further steps where
it is needed.

LoopStructural/interpolators/discrete_interpolator.py

@@ -531,11 +531,15 @@ def _solve(self, solver='cg', **kwargs):
         # check solution is not nan
         # self.support.properties[self.propertyname] = self.c
         if np.all(self.c == np.nan):
+            self.valid = False
             logger.warning("Solver not run, no scalar field")
+            return
         # if solution is all 0, probably didn't work
         if np.all(self.c[self.region] == 0):
+            self.valid = False
             logger.warning("No solution, {} scalar field 0. Add more data.".format(self.propertyname))
-
+            return
+        self.valid = True
     def update(self):
         """
         Check if the solver is up to date, if not rerun interpolation using

LoopStructural/interpolators/geological_interpolator.py

@@ -39,7 +39,8 @@ def __init__(self):
         self.constraints = []
         self.propertyname = 'defaultproperty'
         self.__str = 'Base Geological Interpolator'
-
+        self.valid = False
+
     def __str__(self):
 
         return self.__str

LoopStructural/modelling/features/geological_feature.py

@@ -63,6 +63,9 @@ def __init__(self, name, interpolator, builder=None, data=None, region=None, typ
             self.region = 'everywhere'
         self.model = None
 
+    def is_valid(self):
+        return self.interpolator.valid
+
     def __str__(self):
         return self.name
 

docs/source/_static/loopstructural_design.rst

@@ -0,0 +1 @@
+LoopStructural design