[BUG] Handle potential ValueError in scalar field computation and enforce proper cache initialization in tests and fixtures

Author: Leguark

gempy_engine/modules/evaluator/generic_evaluator.py

@@ -72,7 +72,10 @@ def _eval_on(
     eval_kernel = yield_evaluation_kernel(
         solver_input, options.kernel_options, slice_array=slice_array
     )
-    scalar_field = (eval_kernel.T @ weights).reshape(-1)
+    try:
+        scalar_field = (eval_kernel.T @ weights).reshape(-1)
+    except ValueError:
+        pass
 
     gx_field: Optional[np.ndarray] = None
     gy_field: Optional[np.ndarray] = None

tests/fixtures/simple_geometries.py

@@ -49,7 +49,7 @@ def unconformity() -> Tuple[InterpolationInput, InterpolationOptions, InputDataD
     options = InterpolationOptions.from_args(range_, c_o, uni_degree=1, i_res=i_r, gi_res=gi_r,
                                    number_dimensions=3,
                                    kernel_function=AvailableKernelFunctions.cubic)
-
+    options.cache_mode = InterpolationOptions.CacheMode.NO_CACHE
     resolution = [2, 2, 2]
     extent = [0, 1000, 0, 1000, 0, 1000]
 

tests/fixtures/simple_models.py

@@ -169,6 +169,7 @@ def simple_model_interpolation_input_factory():
     ori_i = Orientations(dip_positions, dip_gradients, nugget_effect_grad)
     interpolation_options = InterpolationOptions.from_args(range_, co, 0, number_dimensions=3,
                                                  kernel_function=AvailableKernelFunctions.cubic)
+    interpolation_options.cache_mode = InterpolationOptions.CacheMode.NO_CACHE
     ids = np.array([1, 2])
     interpolation_input = InterpolationInput(spi, ori_i, grid_0_centers, ids)
     tensor_struct = TensorsStructure(np.array([7]))
@@ -279,6 +280,7 @@ def simple_model_3_layers_high_res(simple_grid_3d_more_points_grid) -> Tuple[Int
 
     interpolation_options = InterpolationOptions.from_args(range_, co, 0,
                                                  number_dimensions=3, kernel_function=AvailableKernelFunctions.cubic)
+    interpolation_options.cache_mode = InterpolationOptions.CacheMode.NO_CACHE
 
     ids = np.array([1, 2, 3, 4])
 

tests/test_common/test_integrations/test_interpolate_model.py

@@ -12,6 +12,9 @@
 
 
 def test_interpolate_model(simple_model_interpolation_input, n_oct_levels=3):
+
+    from gempy_engine.modules.weights_cache.weights_cache_interface import WeightCache
+    WeightCache.initialize_cache_dir()
     """Kernel function Cubic"""
     interpolation_input, options, structure = simple_model_interpolation_input
     print(interpolation_input)
@@ -29,6 +32,9 @@ def test_interpolate_model(simple_model_interpolation_input, n_oct_levels=3):
 
 @pytest.mark.skipif(TEST_SPEED.value <= 1, reason="Global test speed below this test value.")
 def test_interpolate_model_no_octtree(simple_model_3_layers_high_res, n_oct_levels=2):
+
+    from gempy_engine.modules.weights_cache.weights_cache_interface import WeightCache
+    WeightCache.initialize_cache_dir()
     interpolation_input, options, structure = simple_model_3_layers_high_res
     print(interpolation_input)
 
@@ -44,6 +50,9 @@ def test_interpolate_model_no_octtree(simple_model_3_layers_high_res, n_oct_leve
 
 
 def test_interpolate_model_several_surfaces(simple_model_3_layers, n_oct_levels=3):
+
+    from gempy_engine.modules.weights_cache.weights_cache_interface import WeightCache
+    WeightCache.initialize_cache_dir()
     interpolation_input, options, structure = simple_model_3_layers
     print(interpolation_input)