Reworked ODE test with linear interpolator

testing/CMakeLists.txt

@@ -14,8 +14,10 @@ dorie_add_system_test(dorie ode_homogeneous_sand.mini)
 dorie_add_system_test(dorie ode_homogeneous_silt.mini)
 dorie_add_system_test(dorie ode_layered.mini)
 
+dorie_add_system_test(dorie ode_linear_interpolator.mini)
+
 # dorie_add_system_test(dorie ode_evaporation.mini)
-# dorie_add_system_test(dorie ode_linear_interpolator.mini)
+
 # dorie_add_system_test(dorie muphi.mini)
 # dorie_add_system_test_dependency(dorie_muphi dorie-rfg_muphi_pfg)
 # dorie_add_system_test(dorie reference_2d.mini)

testing/ode_linear_interpolator.mini.in

@@ -13,11 +13,15 @@ time.end = 1E7
 time.maxTimestep = 1E7
 time.startTimestep = 1E4
 
-adaptivity.threshold = 1E-9
 adaptivity.useAdaptivity = false, true | expand gridtype
+adaptivity.markingStrategy = threshold
+adaptivity.refinementFraction = 1E-8
+adaptivity.coarseningFraction = 1E-9
+adaptivity.threshold = 5E-8
+
 grid.gridType = rectangular, gmsh | expand gridtype
 grid.cells = 1 100
-grid.initialLevel = 2, 0 | expand gridtype
+grid.initialLevel = 0
 grid.gridFile = "{_asset_path}/meshes/square.msh"
 
 boundary.file = "{_asset_path}/bcs/infiltration_2d.dat"
@@ -26,10 +30,10 @@ parameters.arrayFile = "{_asset_path}/parfields/corners.h5"
 parameters.interpolation = linear
 
 [_ode]
-flux = -3e-6
-head_abstol = 1E-5, 5E-4 | expand gridtype
-flux_abstol = 1E-8, 2E-7 | expand gridtype
+flux = -5.55e-6
+head_abstol = 4E-5, 2E-4 | expand gridtype
+flux_abstol = 3E-9, 7E-8 | expand gridtype
 
 [NewtonParameters]
-AbsoluteLimit = 1E-8
-Reduction = 1E-8
+AbsoluteLimit = 1E-10
+Reduction = 1E-10