{
  "passed": true,
  "critical_angle_deg": 40.49266149200838,
  "cases": {
    "30": {
      "run": {
        "passed": true,
        "source_sha256": "ca2c371dc80d4eb72e6a4be4801ca2261abed4e05514899e3e706783291caccf",
        "upstream_commit": "ca145119e5bd0f0804d846ea42e35e4c3daa91de",
        "unchanged_original_script": true,
        "changes": "Wrapper adds read-only final-two-period field samples only",
        "scheme_unavailable": "No meep, meep-mpi or guile executable in the existing environment",
        "physical": {
          "interface": "planar",
          "s_pol": true,
          "ref_medium": 0,
          "n1": 1.54,
          "n2": 1.0,
          "kw_0": 8,
          "kr_w": 60,
          "kr_c": 150,
          "chi_deg": 30.0
        },
        "numerical": {
          "cell": [
            5,
            5
          ],
          "pml": 0.25,
          "vacuum_frequency": 12,
          "resolution": 184.8,
          "pixels_per_dense_wavelength": 10,
          "courant": 0.5,
          "runtime": 10,
          "optical_periods": 120,
          "unchanged_default_numerics": true,
          "grid_shape": [
            924,
            924
          ]
        },
        "sampling": {
          "count": 62,
          "start": 9.83495670995671,
          "end": 10.0,
          "harmonic_fit_relative_l2_residual": 0.0008222789318154763,
          "formula": "Ez(t)=Re(phasor*exp(-i*2*pi*12*t))+dc",
          "intensity": "mean(Ez^2)=abs(phasor)^2/2, fitted over final two periods; not Poynting flux",
          "animation": "actual FDTD snapshots; every third grid point for display only"
        },
        "environment": {
          "python": "3.11.15 | packaged by conda-forge | (main, Jun 11 2026, 03:34:02) [GCC 14.3.0]",
          "meep": "1.33.0",
          "optbeam": "2.1.2",
          "numpy": "2.4.6",
          "scipy": "1.17.1",
          "h5py": "3.16.0",
          "platform": "Linux-6.6.114.1-microsoft-standard-WSL2-x86_64-with-glibc2.43",
          "mpi_processes": 1,
          "threads": {
            "OMP_NUM_THREADS": "1",
            "OPENBLAS_NUM_THREADS": "1",
            "MKL_NUM_THREADS": "1"
          }
        },
        "elapsed_seconds": 22.879783203999978,
        "command": [
          "/mnt/c/Users/bgs43/Documents/Codex/2026-10-23/https-jin-castle-github-io-meep-4/outputs/optical-beams-run/run_original.py",
          "--upstream",
          "/mnt/c/Users/bgs43/Documents/Codex/2026-10-23/https-jin-castle-github-io-meep-4/work/optical-beams-upstream",
          "--output",
          "/mnt/c/Users/bgs43/Documents/Codex/2026-10-23/https-jin-castle-github-io-meep-4/outputs/optical-beams-run/planar30",
          "--chi",
          "30"
        ]
      },
      "original_instantaneous_e2_max_error": 7.105427357601002e-15,
      "epsilon_min": 1.0,
      "epsilon_max": 2.3716,
      "all_arrays_finite": true,
      "reflected_direction": {
        "angle_from_positive_x_deg": 119.97571438340682,
        "spatial_wavevector": [
          -56.020848416014005,
          97.12600640410747
        ],
        "method": "intensity-weighted local phasor phase gradient",
        "mask": {
          "along_ray": [
            0.7,
            1.7
          ],
          "transverse_half_width": 0.12
        }
      },
      "expected_reflected_direction_from_x_deg": 120,
      "transmitted_direction": {
        "angle_from_positive_x_deg": 20.661840106239556,
        "spatial_wavevector": [
          69.11996014290035,
          26.065661629591972
        ],
        "method": "intensity-weighted local phasor phase gradient",
        "mask": {
          "along_ray": [
            0.6,
            1.6
          ],
          "transverse_half_width": 0.12
        }
      },
      "plane_wave_snell_angle_from_normal_deg": 50.35388885302628,
      "measured_angle_from_normal_deg": 50.66184010623955,
      "plane_wave_snell_angle_error_deg": 0.3079512532132753,
      "passed": true
    },
    "45": {
      "run": {
        "passed": true,
        "source_sha256": "ca2c371dc80d4eb72e6a4be4801ca2261abed4e05514899e3e706783291caccf",
        "upstream_commit": "ca145119e5bd0f0804d846ea42e35e4c3daa91de",
        "unchanged_original_script": true,
        "changes": "Wrapper adds read-only final-two-period field samples only",
        "scheme_unavailable": "No meep, meep-mpi or guile executable in the existing environment",
        "physical": {
          "interface": "planar",
          "s_pol": true,
          "ref_medium": 0,
          "n1": 1.54,
          "n2": 1.0,
          "kw_0": 8,
          "kr_w": 60,
          "kr_c": 150,
          "chi_deg": 45.0
        },
        "numerical": {
          "cell": [
            5,
            5
          ],
          "pml": 0.25,
          "vacuum_frequency": 12,
          "resolution": 184.8,
          "pixels_per_dense_wavelength": 10,
          "courant": 0.5,
          "runtime": 10,
          "optical_periods": 120,
          "unchanged_default_numerics": true,
          "grid_shape": [
            924,
            924
          ]
        },
        "sampling": {
          "count": 62,
          "start": 9.83495670995671,
          "end": 10.0,
          "harmonic_fit_relative_l2_residual": 0.00020694190276210927,
          "formula": "Ez(t)=Re(phasor*exp(-i*2*pi*12*t))+dc",
          "intensity": "mean(Ez^2)=abs(phasor)^2/2, fitted over final two periods; not Poynting flux",
          "animation": "actual FDTD snapshots; every third grid point for display only"
        },
        "environment": {
          "python": "3.11.15 | packaged by conda-forge | (main, Jun 11 2026, 03:34:02) [GCC 14.3.0]",
          "meep": "1.33.0",
          "optbeam": "2.1.2",
          "numpy": "2.4.6",
          "scipy": "1.17.1",
          "h5py": "3.16.0",
          "platform": "Linux-6.6.114.1-microsoft-standard-WSL2-x86_64-with-glibc2.43",
          "mpi_processes": 1,
          "threads": {
            "OMP_NUM_THREADS": "1",
            "OPENBLAS_NUM_THREADS": "1",
            "MKL_NUM_THREADS": "1"
          }
        },
        "elapsed_seconds": 21.290379265999945,
        "command": [
          "/mnt/c/Users/bgs43/Documents/Codex/2026-10-23/https-jin-castle-github-io-meep-4/outputs/optical-beams-run/run_original.py",
          "--upstream",
          "/mnt/c/Users/bgs43/Documents/Codex/2026-10-23/https-jin-castle-github-io-meep-4/work/optical-beams-upstream",
          "--output",
          "/mnt/c/Users/bgs43/Documents/Codex/2026-10-23/https-jin-castle-github-io-meep-4/outputs/optical-beams-run/planar45",
          "--chi",
          "45"
        ]
      },
      "original_instantaneous_e2_max_error": 1.4210854715202004e-14,
      "epsilon_min": 1.0,
      "epsilon_max": 2.3716,
      "all_arrays_finite": true,
      "reflected_direction": {
        "angle_from_positive_x_deg": 90.85074682046765,
        "spatial_wavevector": [
          -1.6328006927503544,
          109.95717040116263
        ],
        "method": "intensity-weighted local phasor phase gradient",
        "mask": {
          "along_ray": [
            0.7,
            1.7
          ],
          "transverse_half_width": 0.12
        }
      },
      "expected_reflected_direction_from_x_deg": 90,
      "passed": true
    }
  },
  "rendering": {
    "shared_wavefront_max": 14.864344142343416,
    "shared_cycle_mean_ez2_max": 112.30480515180463,
    "spatial_display_region": [
      -2.2,
      2.2,
      -2.2,
      2.2
    ],
    "field_opacity": "sqrt(abs(Ez)/shared_max), display only",
    "intensity_scaling": "linear, common limit across both cases"
  },
  "limitations": [
    "Default 10 grid points per denser-medium wavelength; no resolution convergence study",
    "Finite-width beam contains an angular spectrum; the central 45-degree ray is above critical but zero total transmission is not asserted",
    "No flux monitors or incident-only subtraction; no measured R/T efficiency claimed",
    "Phase-gradient direction is a local diagnostic, not a flux or beam-centroid measurement",
    "The length unit a is arbitrary; no micrometer mapping is assumed"
  ]
}