#F-239 In Progress

Inverse solution for three-layer model: ашч H1 and H2 variations

Author: ak-git

Util/src/main/java/com/ak/rsm/inverse/DynamicInverse.java

@@ -43,7 +43,7 @@ public double applyAsDouble(TetrapolarSystem s, double[] kw) {
     };
   }
 
-  private static double dH(Collection<? extends DerivativeResistivity> r) {
+  static double dH(Collection<? extends DerivativeResistivity> r) {
     return r.stream().mapToDouble(DerivativeResistivity::dh)
         .reduce((left, right) -> Double.compare(left, right) == 0 ? left : Double.NaN).orElse(Double.NaN);
   }

Util/src/test/java/com/ak/rsm/inverse/Inverse3Test.java

@@ -93,7 +93,8 @@ private static void testSingle(Collection<? extends DerivativeMeasurement> dm) {
 
     assertThatNoException().isThrownBy(() -> {
       int scale = 2;
-      double hStep = dm.stream().flatMapToDouble(m -> DoubleStream.of(Math.abs(m.dh()))).min().orElseThrow() / scale;
+      double dh1 = DynamicInverse.dH(dm);
+      double hStep = Math.abs(dh1 / scale);
       LOGGER.info("3-layer step {}", ValuePair.Name.H.of(hStep, 0.0));
 
           ValuePair.Name[] keys = {
@@ -110,7 +111,9 @@ private static void testSingle(Collection<? extends DerivativeMeasurement> dm) {
                 double[] p = Stream.of(h1, dh2).mapToDouble(x -> x.value() / hStep).toArray();
 
                 PointValuePair optimizedK = Simplex.optimizeAll(
-                    k -> DynamicInverse.ofH1Changed(dm, hStep).applyAsDouble(DoubleStream.concat(Arrays.stream(k), Arrays.stream(p)).toArray()),
+                    k -> DynamicInverse.of(dm, hStep, DeltaH.H1.apply(dh1)).applyAsDouble(
+                        DoubleStream.concat(Arrays.stream(k), Arrays.stream(p)).toArray()
+                    ),
                     new Simplex.Bounds(0.0, 1.0), new Simplex.Bounds(-1.0, 0.0)
                 );