tmmcore 0.3.1 → 0.4.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "tmmcore",
-  "version": "0.3.1",
+  "version": "0.4.0",
   "description": "Transfer-matrix method for multilayer thin-film optics, with exact analytic derivatives and phase dispersion. JavaScript, C and WebAssembly.",
   "type": "module",
   "main": "./src/index.js",
@@ -23,7 +23,7 @@
     "test": "node tests/equivalence.mjs && node tests/growing_equivalence.mjs && node tests/growing_eval_equivalence.mjs && node tests/absorbing_incident.mjs",
     "compare": "node benchmarks/compare.mjs",
     "examples": "node examples/01-single-layer.mjs && node examples/02-ar-coating.mjs && node examples/03-metal-mirror.mjs && node examples/04-refine.mjs && node examples/05-angle-map.mjs && node examples/06-needle.mjs && node examples/07-jacobian-map.mjs && node examples/08-group-delay.mjs",
-    "build:wasm": "emcc src/tmm_kernel.c -O3 --no-entry -sSTANDALONE_WASM=1 -sALLOW_MEMORY_GROWTH=1 -sEXPORTED_FUNCTIONS=_tmm_one,_tmm_spectrum,_tmm_jacobian,_tmm_needle_scan,_tmm_hessian,_tmm_phase_one,_tmm_phase_spectrum,_tmm_phase_jacobian,_tmm_monitor_curve,_tmm_deposition_spectra,_tmm_growing_eval_create,_tmm_growing_eval_set_top,_tmm_growing_eval_sample,_tmm_growing_eval_free,_malloc,_free -o src/tmm_kernel.wasm",
+    "build:wasm": "emcc src/tmm_kernel.c -O3 --no-entry -sSTANDALONE_WASM=1 -sALLOW_MEMORY_GROWTH=1 -sEXPORTED_FUNCTIONS=_tmm_one,_tmm_spectrum,_tmm_jacobian,_tmm_needle_scan,_tmm_hessian,_tmm_phase_one,_tmm_phase_spectrum,_tmm_phase_jacobian,_tmm_phase_jacobian_spectrum,_tmm_monitor_curve,_tmm_deposition_spectra,_tmm_growing_eval_create,_tmm_growing_eval_set_top,_tmm_growing_eval_sample,_tmm_growing_eval_free,_malloc,_free -o src/tmm_kernel.wasm",
     "build:native": "cc -std=c99 -pedantic -Wall -Wextra -O2 -c src/tmm_kernel.c -o tmm_kernel.o"
   },
   "keywords": [

package/src/phase.js

@@ -430,14 +430,17 @@ export function coefficientPhaseDispersion(coefficientJet) {
  */
 export function coefficientPhaseThicknessDerivatives(coefficientJet, thicknessJets) {
     if (!thicknessJets) return null;
-    const result = { phaseDeg: [], gd: [], gdd: [], tod: [] };
+    const result = { phaseDeg: [], gd: [], gdd: [], tod: [], logMagnitudeSquared: [] };
     for (const thicknessJet of thicknessJets) {
         const logarithmicDerivative = jetDivide(thicknessJet, coefficientJet);
         const derivatives = jetDerivatives(logarithmicDerivative);
+        // The imaginary part of d(ln c)/dd is the phase derivative; the real
+        // part is d(ln |c|)/dd, so twice it is the relative derivative of |c|².
         result.phaseDeg.push(-derivatives[0][1] * 180 / Math.PI);
         result.gd.push(derivatives[1][1]);
         result.gdd.push(derivatives[2][1]);
         result.tod.push(derivatives[3][1]);
+        result.logMagnitudeSquared.push(2 * derivatives[0][0]);
     }
     return result;
 }
@@ -500,8 +503,10 @@ export function tmmPhaseDispersion(lambda_nm, theta_deg, pol, n0Jet, nsJet, laye
  * thicknesses are skipped and receive zero derivative entries.
  *
  * @returns {{r, t}} where each side is the phase quantities plus
- *   `{dPhaseDeg, dGd, dGdd, dTod}`, arrays of length `layers.length`. The
- *   derivative arrays are `null` if the matrix product overflowed.
+ *   `{dPhaseDeg, dGd, dGdd, dTod, dLogMagnitudeSquared}`, arrays of length
+ *   `layers.length`. `dLogMagnitudeSquared` is d(ln |coefficient|²)/dd, the
+ *   relative intensity derivative. The derivative arrays are `null` if the
+ *   matrix product overflowed.
  */
 export function tmmPhaseThicknessJacobian(lambda_nm, theta_deg, pol, n0Jet, nsJet, layers, options = {}) {
     const prepared = prepare(lambda_nm, layers, options);
@@ -524,6 +529,7 @@ export function tmmPhaseThicknessJacobian(lambda_nm, theta_deg, pol, n0Jet, nsJe
             dGd: derivatives ? derivatives.gd : null,
             dGdd: derivatives ? derivatives.gdd : null,
             dTod: derivatives ? derivatives.tod : null,
+            dLogMagnitudeSquared: derivatives ? derivatives.logMagnitudeSquared : null,
         };
     };
     return {