Mercurial > repos > pointsource_pde / changeset

--- a/experiments/laser_and_mirrors_aux.py	Wed Oct 07 08:02:41 2026 -0500
+++ b/experiments/laser_and_mirrors_aux.py	Wed Oct 07 08:25:46 2026 -0500
@@ -245,27 +245,27 @@
         true_u_n_list_array=[u_n_list[i].x.array for i in plotter.plot_frames],
     )
 
-    delta = 0.001 * k
-    dk_est = (
-        dat.apply((μ_true, (k + delta, (c1, c2)))) - dat.apply((μ_true, auxtrue))
-    ) / delta
-    print("k-diff-test:", dk, ", ", dk_est)
-    delta = 0.001 * c1
-    dc1_est = (
-        dat.apply((μ_true, (k, (c1 + delta, c2)))) - dat.apply((μ_true, auxtrue))
-    ) / delta
-    print("c1-diff-test:", dc1, ", ", dc1_est)
-    delta = 0.001
-    μ_perturb = DiscreteMeasure_2_f64(
-        [([0.1, 0.3], 0.08 + delta), ([0.4, 0.25], 0.05), ([0.25, 0.13], 0.06)]
-    )
-    dcμ_est = (dat.apply((μ_perturb, auxtrue)) - dat.apply((μ_true, auxtrue))) / delta
-    print(
-        "μ-diff-test:",
-        eval_Function_f64(ω_true._cpp_object, [0.1, 0.3, 0.0]),
-        ", ",
-        dcμ_est,
-    )
+    # delta = 0.001 * k
+    # dk_est = (
+    #     dat.apply((μ_true, (k + delta, (c1, c2)))) - dat.apply((μ_true, auxtrue))
+    # ) / delta
+    # print("k-diff-test:", dk, ", ", dk_est)
+    # delta = 0.001 * c1
+    # dc1_est = (
+    #     dat.apply((μ_true, (k, (c1 + delta, c2)))) - dat.apply((μ_true, auxtrue))
+    # ) / delta
+    # print("c1-diff-test:", dc1, ", ", dc1_est)
+    # delta = 0.001
+    # μ_perturb = DiscreteMeasure_2_f64(
+    #     [([0.1, 0.3], 0.08 + delta), ([0.4, 0.25], 0.05), ([0.25, 0.13], 0.06)]
+    # )
+    # dcμ_est = (dat.apply((μ_perturb, auxtrue)) - dat.apply((μ_true, auxtrue))) / delta
+    # print(
+    #     "μ-diff-test:",
+    #     eval_Function_f64(ω_true._cpp_object, [0.1, 0.3, 0.0]),
+    #     ", ",
+    #     dcμ_est,
+    # )
 
     return dat, auxtrue, aux_noisy, μ_bound, μ_true, plot_factory, pde