mirror of
https://github.com/OPM/ResInsight.git
synced 2026-10-25 14:40:26 -05:00
Kristian Bendiksen
b78ef0d0a1
Expose the new RimFaultDistanceResult facility to Python by marking the relevant PDM classes and fields scriptable. RimEclipseView's m_faultCollection, RimFaultInViewCollection, RimFaultInView, RimFaultDistanceResult, and RimFaultDistanceResultCollection are now scriptable, so faults can be enumerated from Python via view.fault_collection().faults(). A new Rimc method add_fault_distance_result on RimFaultInViewCollection takes a name and a list of fault objects, creates the result (auto-numbered FAULTDIST<n> when the name is empty), and returns the new RimFaultDistanceResult. An empty fault list falls back to every fault in the view, matching the static-native FAULTDIST behaviour. Add a Python test exercising both a subset call and a parity check that asserts the all-faults result matches the existing static-native FAULTDIST cell-for-cell.
56 lines
1.6 KiB
Python
56 lines
1.6 KiB
Python
import math
|
|
import os
|
|
import sys
|
|
|
|
sys.path.insert(1, os.path.join(sys.path[0], "../../"))
|
|
|
|
import dataroot
|
|
import rips
|
|
|
|
|
|
GRID_PATH = dataroot.PATH + "/TEST10K_FLT_LGR_NNC/TEST10K_FLT_LGR_NNC.EGRID"
|
|
|
|
|
|
def _get_view(case):
|
|
views = case.views()
|
|
if views:
|
|
return views[0]
|
|
return case.create_view()
|
|
|
|
|
|
def test_add_fault_distance_result_subset(rips_instance, initialize_test):
|
|
case = rips_instance.project.load_case(GRID_PATH)
|
|
view = _get_view(case)
|
|
fault_collection = view.fault_collection()
|
|
faults = fault_collection.faults()
|
|
assert len(faults) > 0
|
|
|
|
result = fault_collection.add_fault_distance_result("PY_FAULTDIST", faults[:1])
|
|
assert result is not None
|
|
assert result.result_name == "PY_FAULTDIST"
|
|
|
|
|
|
def test_add_fault_distance_result_all_matches_static(rips_instance, initialize_test):
|
|
"""Passing every fault into add_fault_distance_result must reproduce the
|
|
static-native FAULTDIST values cell-for-cell."""
|
|
case = rips_instance.project.load_case(GRID_PATH)
|
|
view = _get_view(case)
|
|
fault_collection = view.fault_collection()
|
|
all_faults = fault_collection.faults()
|
|
assert len(all_faults) > 0
|
|
|
|
static_values = case.active_cell_property(
|
|
rips.PropertyType.STATIC_NATIVE, "FAULTDIST", 0
|
|
)
|
|
|
|
fault_collection.add_fault_distance_result("PY_FAULTDIST_ALL", all_faults)
|
|
|
|
generated_values = case.active_cell_property(
|
|
rips.PropertyType.GENERATED, "PY_FAULTDIST_ALL", 0
|
|
)
|
|
|
|
assert len(static_values) == len(generated_values)
|
|
assert len(static_values) > 0
|
|
for s, g in zip(static_values, generated_values):
|
|
assert math.isclose(s, g, abs_tol=1e-6)
|