eyeling 1.17.2 → 1.18.1

package/examples/bmi.n3

@@ -0,0 +1,219 @@
+# ============================================================================
+# BMI — ARC-style Body Mass Index example.
+#
+# This example turns a familiar health calculation into a small, inspectable
+# ARC program. It normalizes either metric or US inputs, computes BMI, assigns
+# a WHO adult category, and derives a healthy-range weight band for the given
+# height. The report explains the result and includes independent checks for
+# boundary handling and category behavior.
+#
+# For reproducibility and documentation only; not medical advice.
+# ============================================================================
+
+@prefix : <https://example.org/bmi#> .
+@prefix log: <http://www.w3.org/2000/10/swap/log#> .
+@prefix math: <http://www.w3.org/2000/10/swap/math#> .
+@prefix string: <http://www.w3.org/2000/10/swap/string#> .
+
+# --------------
+# Editable input
+# --------------
+
+:Input
+  :unitSystem "metric";
+  :weight 72.0;
+  :height 178.0.
+
+# US alternative:
+# :Input
+#   :unitSystem "us";
+#   :weight 158.73;
+#   :height 70.08.
+
+# ---------------------------------
+# Normalization and BMI calculation
+# ---------------------------------
+
+{ :Input :unitSystem "metric"; :weight ?W; :height ?H.
+  (?H 100.0) math:quotient ?M. }
+    =>
+{ :Case :weightKg ?W; :heightM ?M.
+  :Reason :units "Inputs were already metric, so kilograms stay kilograms and centimeters are divided by 100 to obtain meters.". } .
+
+{ :Input :unitSystem "us"; :weight ?W; :height ?H.
+  (?W 0.45359237) math:product ?Kg.
+  (?H 0.0254) math:product ?M. }
+    =>
+{ :Case :weightKg ?Kg; :heightM ?M.
+  :Reason :units "US inputs were converted to SI units: pounds to kilograms and inches to meters.". } .
+
+{ :Case :weightKg ?Kg; :heightM ?M.
+  (?M ?M) math:product ?M2.
+  (?Kg ?M2) math:quotient ?Bmi.
+  (?Bmi 100.0) math:product ?BmiX100.
+  ?BmiX100 math:rounded ?BmiRoundedInt.
+  (?BmiRoundedInt 100.0) math:quotient ?BmiRounded.
+  (18.5 ?M2) math:product ?HealthyMin.
+  (24.9 ?M2) math:product ?HealthyMax.
+  (?HealthyMin 10.0) math:product ?MinX10.
+  (?HealthyMax 10.0) math:product ?MaxX10.
+  ?MinX10 math:rounded ?MinRoundedInt.
+  ?MaxX10 math:rounded ?MaxRoundedInt.
+  (?MinRoundedInt 10.0) math:quotient ?HealthyMinRounded.
+  (?MaxRoundedInt 10.0) math:quotient ?HealthyMaxRounded. }
+    =>
+{ :Case :heightSquared ?M2;
+        :bmi ?Bmi;
+        :bmiRounded ?BmiRounded;
+        :healthyMinKg ?HealthyMin;
+        :healthyMaxKg ?HealthyMax;
+        :healthyMinKgRounded ?HealthyMinRounded;
+        :healthyMaxKgRounded ?HealthyMaxRounded. } .
+
+# ------------------------------------------
+# WHO adult categories (half-open intervals)
+# ------------------------------------------
+
+{ :Case :bmi ?Bmi. ?Bmi math:lessThan 18.5. }
+    => { :Decision :category "Underweight". } .
+
+{ :Case :bmi ?Bmi. ?Bmi math:notLessThan 18.5. ?Bmi math:lessThan 25.0. }
+    => { :Decision :category "Normal". } .
+
+{ :Case :bmi ?Bmi. ?Bmi math:notLessThan 25.0. ?Bmi math:lessThan 30.0. }
+    => { :Decision :category "Overweight". } .
+
+{ :Case :bmi ?Bmi. ?Bmi math:notLessThan 30.0. ?Bmi math:lessThan 35.0. }
+    => { :Decision :category "Obesity I". } .
+
+{ :Case :bmi ?Bmi. ?Bmi math:notLessThan 35.0. ?Bmi math:lessThan 40.0. }
+    => { :Decision :category "Obesity II". } .
+
+{ :Case :bmi ?Bmi. ?Bmi math:notLessThan 40.0. }
+    => { :Decision :category "Obesity III". } .
+
+# ---------------------
+# Answer and reason why
+# ---------------------
+
+{ :Case :bmiRounded ?BmiRounded;
+        :healthyMinKgRounded ?HealthyMinRounded;
+        :healthyMaxKgRounded ?HealthyMaxRounded;
+        :heightM ?M.
+  :Decision :category ?Category.
+  (?M 100.0) math:product ?Cm.
+  ?Cm math:rounded ?CmRounded. }
+    =>
+{ :Answer :bmi ?BmiRounded;
+          :category ?Category;
+          :healthyMinKg ?HealthyMinRounded;
+          :healthyMaxKg ?HealthyMaxRounded;
+          :heightCm ?CmRounded. } .
+
+{ :Case :weightKg ?Kg; :heightM ?M; :heightSquared ?M2; :bmiRounded ?BmiRounded.
+  :Reason :units ?Units.
+  :Decision :category ?Category. }
+    =>
+{ :Reason :formula "BMI is defined as weight in kilograms divided by height in meters squared.";
+          :calculation "The normalized weight and height were used to compute BMI, then the result was mapped to the WHO adult category table.";
+          :categoryRule ?Category;
+          :unitsExplanation ?Units. } .
+
+# ------------------
+# Independent checks
+# ------------------
+
+{ :Case :weightKg ?Kg; :heightM ?M. ?Kg math:greaterThan 0. ?M math:greaterThan 0. }
+    => { :Check :c1 "OK - the input was normalized into positive SI values.". } .
+
+{ :Case :heightM ?M; :heightSquared ?M2. (?M ?M) math:product ?M2. }
+    => { :Check :c2 "OK - height squared was reconstructed from the normalized height.". } .
+
+{ :Case :weightKg ?Kg; :heightSquared ?M2; :bmi ?Bmi. (?Kg ?M2) math:quotient ?Bmi. }
+    => { :Check :c3 "OK - the BMI value matches the BMI = kg / m² formula.". } .
+
+{ 18.49 math:lessThan 18.5. }
+    => { :Check :c4 "OK - a BMI of 18.49 stays below the normal-weight threshold.". } .
+
+{ 18.5 math:notLessThan 18.5. 18.5 math:lessThan 25.0. }
+    => { :Check :c5 "OK - the lower boundary is half-open: BMI 18.5 is classified as Normal.". } .
+
+{ 25.0 math:notLessThan 25.0. 25.0 math:lessThan 30.0. }
+    => { :Check :c6 "OK - BMI 25.0 starts the Overweight category.". } .
+
+{ 30.0 math:notLessThan 30.0. 30.0 math:lessThan 35.0. }
+    => { :Check :c7 "OK - BMI 30.0 starts the Obesity I category.". } .
+
+{ 22.0 math:notLessThan 18.5. 22.0 math:lessThan 25.0.
+  27.0 math:notLessThan 25.0. 27.0 math:lessThan 30.0.
+  41.0 math:notLessThan 40.0. }
+    => { :Check :c8 "OK - classification behavior is monotonic across representative BMI values.". } .
+
+{ :Case :heightSquared ?M2; :healthyMinKg ?Min; :healthyMaxKg ?Max.
+  (18.5 ?M2) math:product ?Min.
+  (24.9 ?M2) math:product ?Max. }
+    => { :Check :c9 "OK - the healthy-weight band was reconstructed from BMI 18.5 to 24.9 at the same height.". } .
+
+# -----
+# Fuses
+# -----
+
+{ :Decision :category ?Any.
+  1 log:notIncludes {
+    :Check :c1 ?C1.
+    :Check :c2 ?C2.
+    :Check :c3 ?C3.
+    :Check :c4 ?C4.
+    :Check :c5 ?C5.
+    :Check :c6 ?C6.
+    :Check :c7 ?C7.
+    :Check :c8 ?C8.
+    :Check :c9 ?C9.
+  }. }
+    => false .
+
+{ :Decision :category ?C1, ?C2.
+  ?C1 log:notEqualTo ?C2. }
+    => false .
+
+# ---------------------------------------------
+# ARC report (built from actual derived values)
+# ---------------------------------------------
+
+{ :Answer :bmi ?BmiRounded;
+          :category ?Category;
+          :healthyMinKg ?HealthyMinRounded;
+          :healthyMaxKg ?HealthyMaxRounded;
+          :heightCm ?CmRounded.
+  :Check :c1 ?C1.
+  :Check :c2 ?C2.
+  :Check :c3 ?C3.
+  :Check :c4 ?C4.
+  :Check :c5 ?C5.
+  :Check :c6 ?C6.
+  :Check :c7 ?C7.
+  :Check :c8 ?C8.
+  :Check :c9 ?C9.
+  (
+    "BMI — ARC-style Body Mass Index example\n\n"
+    "Answer\n"
+    "BMI = " ?BmiRounded "\n"
+    "Category = " ?Category "\n"
+    "At height " ?CmRounded " cm, a healthy-weight range is about " ?HealthyMinRounded "–" ?HealthyMaxRounded " kg (BMI 18.5–24.9).\n\n"
+    "Reason Why\n"
+    "BMI is defined as weight in kilograms divided by height in meters squared. "
+    "This program first normalizes the input to SI units, computes BMI, and then applies WHO adult categories as half-open intervals. "
+    "The healthy-weight band is the weight range at the same height that corresponds to BMI 18.5 through 24.9.\n\n"
+    "Check\n"
+    "C1 " ?C1 "\n"
+    "C2 " ?C2 "\n"
+    "C3 " ?C3 "\n"
+    "C4 " ?C4 "\n"
+    "C5 " ?C5 "\n"
+    "C6 " ?C6 "\n"
+    "C7 " ?C7 "\n"
+    "C8 " ?C8 "\n"
+    "C9 " ?C9 "\n"
+  ) string:concatenation ?Block. }
+    =>
+{ :report log:outputString ?Block. } .

package/examples/output/bmi.n3

@@ -0,0 +1,20 @@
+BMI — ARC-style Body Mass Index example
+
+Answer
+BMI = 22.72
+Category = Normal
+At height 178 cm, a healthy-weight range is about 58.6–78.9 kg (BMI 18.5–24.9).
+
+Reason Why
+BMI is defined as weight in kilograms divided by height in meters squared. This program first normalizes the input to SI units, computes BMI, and then applies WHO adult categories as half-open intervals. The healthy-weight band is the weight range at the same height that corresponds to BMI 18.5 through 24.9.
+
+Check
+C1 OK - the input was normalized into positive SI values.
+C2 OK - height squared was reconstructed from the normalized height.
+C3 OK - the BMI value matches the BMI = kg / m² formula.
+C4 OK - a BMI of 18.49 stays below the normal-weight threshold.
+C5 OK - the lower boundary is half-open: BMI 18.5 is classified as Normal.
+C6 OK - BMI 25.0 starts the Overweight category.
+C7 OK - BMI 30.0 starts the Obesity I category.
+C8 OK - classification behavior is monotonic across representative BMI values.
+C9 OK - the healthy-weight band was reconstructed from BMI 18.5 to 24.9 at the same height.

package/examples/output/pn-junction-tunneling.n3

@@ -0,0 +1,23 @@
+=== Answer ===
+In this toy PN-junction tunneling model, heavy doping narrows the depletion region enough for a tunneling window that rises to a peak and then falls, producing a negative-differential region.
+case                          : pn-junction-tunneling
+peak bias                     : 2
+peak current proxy            : 4
+negative differential region  : yes
+
+=== Reason Why ===
+We model tunneling current as an exact overlap count between filled N-side states and empty P-side states while forward bias shifts the bands. Heavy doping is represented by a much narrower depletion region.
+ordinary depletion width (nm) : 8
+tunnel depletion width (nm)   : 1
+filled N-side states          : [1, 2, 3, 4]
+empty P-side states at 0 bias : [3, 4, 5, 6]
+bias -> overlap current proxy : 0->2, 1->3, 2->4, 3->3, 4->2, 5->1, 6->0
+peak point                    : 2 -> 4
+high-bias point               : 6 -> 0
+
+=== Check ===
+heavily doped barrier is narrower : yes
+peak occurs before overlap closes : yes
+negative differential region present : yes
+high-bias overlap closes           : yes
+peak equals full four-state overlap: yes

package/examples/output/sudoku.n3

@@ -0,0 +1,42 @@
+=== Answer ===
+The puzzle is solved, and the completed grid is the unique valid Sudoku solution.
+case              : sudoku
+default puzzle    : classic
+
+Puzzle
+1 . . | . . 7 | . 9 . 
+. 3 . | . 2 . | . . 8 
+. . 9 | 6 . . | 5 . . 
+
+. . 5 | 3 . . | 9 . . 
+. 1 . | . 8 . | . . 2 
+6 . . | . . 4 | . . . 
+
+3 . . | . . . | . 1 . 
+. 4 . | . . . | . . 7 
+. . 7 | . . . | 3 . . 
+
+Completed grid
+1 6 2 | 8 5 7 | 4 9 3 
+5 3 4 | 1 2 9 | 6 7 8 
+7 8 9 | 6 4 3 | 5 2 1 
+
+4 7 5 | 3 1 2 | 9 8 6 
+9 1 3 | 5 8 6 | 7 4 2 
+6 2 8 | 7 9 4 | 1 3 5 
+
+3 5 6 | 4 7 8 | 2 1 9 
+2 4 1 | 9 3 5 | 8 6 7 
+8 9 7 | 2 6 1 | 3 5 4 
+=== Reason Why ===
+The solver starts from 23 clues and fills the remaining 58 cells by combining constraint propagation with depth-first search. At each step it chooses the empty cell with the fewest legal digits, places forced singles immediately, and only guesses when more than one candidate remains. Across the search it made 213 forced placements and tried 22 guesses, visited 23 search nodes overall, and backtracked 12 times before reaching the completed grid. The solver also confirmed that the solution is unique. Early steps: r2c3=4: guess, r5c3=3: forced, r2c1=5: guess, r2c4=1: guess, r2c6=9: forced, r2c7=6: guess, r2c8=7: forced, r1c7=4: guess, … and 50 more placements
+
+=== Check ===
+C1 OK - every given clue is preserved in the final grid.
+C2 OK - the final grid contains only digits 1 through 9, with no blanks left.
+C3 OK - each row contains every digit exactly once.
+C4 OK - each column contains every digit exactly once.
+C5 OK - each 3×3 box contains every digit exactly once.
+C6 OK - replaying the recorded placements from the original puzzle remains legal at every step.
+C7 OK - the search statistics and the successful proof path are internally consistent.
+C8 OK - a second search found no alternative solution, so the solution is unique.

package/examples/output/transistor-switch.n3

@@ -0,0 +1,24 @@
+=== Answer ===
+In this toy transistor-switch model, a low input leaves the transistor in cutoff (OFF) and a high input drives it into saturation (ON), so the load behaves like an on/off branch rather than a linear amplifier.
+case                             : transistor-switch
+low input state                  : cutoff / OFF
+high input state                 : saturation / ON
+on-state load current            : 4.80 mA
+
+=== Reason Why ===
+We model an NPN low-side switch with exact millivolt and microamp arithmetic. The base current comes from (Vin - Vbe,on)/Rb when the base-emitter junction is forward biased, and the collector current is the smaller of beta * Ib and the load-limited current (Vcc - Vce,sat)/Rl.
+supply voltage                   : 5.00 V
+base resistor                    : 10000 ohms
+load resistor                    : 1000 ohms
+transistor beta proxy            : 100
+low input                        : Vin=0.00 V -> Ib=0.00 mA, Ic=0.00 mA, Vce=5.00 V, state=cutoff / OFF
+high input                       : Vin=5.00 V -> Ib=0.43 mA, Ic=4.80 mA, Vce=0.20 V, state=saturation / ON
+high-input gain limit            : 43.00 mA
+high-input load limit            : 4.80 mA
+
+=== Check ===
+low input stays in cutoff        : yes
+high input reaches saturation    : yes
+switching states differ          : yes
+on-state current is load-limited : yes
+load voltage matches resistor drop : yes

package/examples/pn-junction-tunneling.n3

@@ -0,0 +1,227 @@
+# ===========================================================================
+# pn-junction-tunneling.n3
+#
+# - ordinary depletion width  = 8 nm
+# - tunnel depletion width    = 1 nm
+# - filled N-side states      = [1, 2, 3, 4]
+# - empty P-side states @ 0V  = [3, 4, 5, 6]
+# - bias points               = 0..6
+#
+# For each bias, the empty P-side levels shift downward by the bias amount.
+# A shifted level contributes 1 to the current proxy exactly when it lands in
+# the occupied N-side window [1, 4]. The four slot contributions are summed.
+# ===========================================================================
+
+@prefix : <http://example.org/pn-junction-tunneling#> .
+@prefix log: <http://www.w3.org/2000/10/swap/log#> .
+@prefix math: <http://www.w3.org/2000/10/swap/math#> .
+@prefix list: <http://www.w3.org/2000/10/swap/list#> .
+@prefix string: <http://www.w3.org/2000/10/swap/string#> .
+
+# -----
+# Facts
+# -----
+
+:case :name "pn-junction-tunneling" .
+:case :ordinaryDepletionWidthNm 8 .
+:case :tunnelDepletionWidthNm 1 .
+:case :nFilledLevels (1 2 3 4) .
+:case :pEmptyZeroBiasLevels (3 4 5 6) .
+
+:b0 :value 0 .
+:b1 :value 1 .
+:b2 :value 2 .
+:b3 :value 3 .
+:b4 :value 4 .
+:b5 :value 5 .
+:b6 :value 6 .
+
+# -----
+# Logic
+# -----
+
+# Heavy doping is modeled as a much narrower depletion region.
+{ :case :tunnelDepletionWidthNm ?t .
+  :case :ordinaryDepletionWidthNm ?o .
+  ?t math:lessThan ?o .
+}
+    => { :case :barrierIsNarrower true . } .
+
+# For each bias node and each zero-bias P-side slot, shift the level downward:
+# shifted = zeroBiasLevel - bias.
+{ ?bnode list:in (:b0 :b1 :b2 :b3 :b4 :b5 :b6) .
+  ?bnode :value ?bias .
+  :case :pEmptyZeroBiasLevels ?levels .
+  ?levels list:iterate ( ?slot ?zeroLevel ) .
+  ( ?zeroLevel ?bias ) math:difference ?shifted .
+}
+    => { ( ?bnode ?slot ) :shiftedLevel ?shifted . } .
+
+# Each shifted slot contributes 1 when it falls inside the occupied N-side
+# window [1, 4], otherwise 0.
+{ ( ?bnode ?slot ) :shiftedLevel ?shifted .
+  ?shifted math:notLessThan 1 .
+  ?shifted math:notGreaterThan 4 .
+}
+    => { ( ?bnode ?slot ) :overlapValue 1 . } .
+
+{ ( ?bnode ?slot ) :shiftedLevel ?shifted .
+  ?shifted math:lessThan 1 .
+}
+    => { ( ?bnode ?slot ) :overlapValue 0 . } .
+
+{ ( ?bnode ?slot ) :shiftedLevel ?shifted .
+  ?shifted math:greaterThan 4 .
+}
+    => { ( ?bnode ?slot ) :overlapValue 0 . } .
+
+# The current proxy is the sum of the four slot contributions.
+{ ( ?bnode 0 ) :overlapValue ?v0 .
+  ( ?bnode 1 ) :overlapValue ?v1 .
+  ( ?bnode 2 ) :overlapValue ?v2 .
+  ( ?bnode 3 ) :overlapValue ?v3 .
+  ( ?v0 ?v1 ?v2 ?v3 ) math:sum ?current .
+}
+    => { ?bnode :currentProxy ?current . } .
+
+# Peak = the unique maximum at bias 2 in this toy curve.
+{ :b2 :currentProxy ?peak .
+  :b0 :currentProxy ?c0 . ?peak math:greaterThan ?c0 .
+  :b1 :currentProxy ?c1 . ?peak math:greaterThan ?c1 .
+  :b3 :currentProxy ?c3 . ?peak math:greaterThan ?c3 .
+  :b4 :currentProxy ?c4 . ?peak math:greaterThan ?c4 .
+  :b5 :currentProxy ?c5 . ?peak math:greaterThan ?c5 .
+  :b6 :currentProxy ?c6 . ?peak math:greaterThan ?c6 .
+}
+    => { :case :peakBias 2 .
+         :case :peakCurrent ?peak . } .
+
+# High-bias valley = the final sampled point, matching the Rust code.
+{ :b6 :currentProxy ?valley . }
+    => { :case :valleyBias 6 .
+         :case :valleyCurrent ?valley . } .
+
+# A negative differential region is present once the curve drops after the peak.
+{ :b2 :currentProxy ?c2 .
+  :b3 :currentProxy ?c3 .
+  ?c2 math:greaterThan ?c3 .
+}
+    => { :case :negativeDifferentialRegion true . } .
+
+{ :case :peakBias ?peakBias .
+  :case :valleyBias ?valleyBias .
+  ?peakBias math:lessThan ?valleyBias .
+}
+    => { :case :peakIsBeforeValley true . } .
+
+{ :case :valleyCurrent 0 . }
+    => { :case :overlapClosesAtHighBias true . } .
+
+{ :case :peakCurrent ?peakCurrent .
+  :case :nFilledLevels ?nLevels .
+  ?nLevels list:length ?nCount .
+  ?peakCurrent math:equalTo ?nCount .
+}
+    => { :case :peakMatchesFullOverlap true . } .
+
+# -----------------------------------------------
+# Presentation (ARC: Answer • Reason Why • Check)
+# -----------------------------------------------
+
+# Answer
+{ :case :peakBias ?peakBias .
+  :case :peakCurrent ?peakCurrent .
+  :case :negativeDifferentialRegion true .
+  ( "peak bias                     : %s\n" ?peakBias ) string:format ?peakBiasLine .
+  ( "peak current proxy            : %s\n" ?peakCurrent ) string:format ?peakCurrentLine .
+}
+    => {
+         :01-answer-1 log:outputString "=== Answer ===\n" .
+         :01-answer-2 log:outputString "In this toy PN-junction tunneling model, heavy doping narrows the depletion region enough for a tunneling window that rises to a peak and then falls, producing a negative-differential region.\n" .
+         :01-answer-3 log:outputString "case                          : pn-junction-tunneling\n" .
+         :01-answer-4 log:outputString ?peakBiasLine .
+         :01-answer-5 log:outputString ?peakCurrentLine .
+         :01-answer-6 log:outputString "negative differential region  : yes\n\n" .
+       } .
+
+# Reason Why
+{ :case :ordinaryDepletionWidthNm ?ordinary .
+  :case :tunnelDepletionWidthNm ?tunnel .
+  :b0 :currentProxy ?c0 .
+  :b1 :currentProxy ?c1 .
+  :b2 :currentProxy ?c2 .
+  :b3 :currentProxy ?c3 .
+  :b4 :currentProxy ?c4 .
+  :b5 :currentProxy ?c5 .
+  :b6 :currentProxy ?c6 .
+  :case :peakBias ?peakBias .
+  :case :peakCurrent ?peakCurrent .
+  :case :valleyBias ?valleyBias .
+  :case :valleyCurrent ?valleyCurrent .
+  ( "%s->%s, %s->%s, %s->%s, %s->%s, %s->%s, %s->%s, %s->%s"
+    0 ?c0 1 ?c1 2 ?c2 3 ?c3 4 ?c4 5 ?c5 6 ?c6 ) string:format ?curve .
+  ( "ordinary depletion width (nm) : %s\n" ?ordinary ) string:format ?ordinaryLine .
+  ( "tunnel depletion width (nm)   : %s\n" ?tunnel ) string:format ?tunnelLine .
+  ( "bias -> overlap current proxy : %s\n" ?curve ) string:format ?curveLine .
+  ( "peak point                    : %s -> %s\n" ?peakBias ?peakCurrent ) string:format ?peakLine .
+  ( "high-bias point               : %s -> %s\n" ?valleyBias ?valleyCurrent ) string:format ?valleyLine .
+}
+    => {
+         :02-why-1 log:outputString "=== Reason Why ===\n" .
+         :02-why-2 log:outputString "We model tunneling current as an exact overlap count between filled N-side states and empty P-side states while forward bias shifts the bands. Heavy doping is represented by a much narrower depletion region.\n" .
+         :02-why-3 log:outputString ?ordinaryLine .
+         :02-why-4 log:outputString ?tunnelLine .
+         :02-why-5 log:outputString "filled N-side states          : [1, 2, 3, 4]\n" .
+         :02-why-6 log:outputString "empty P-side states at 0 bias : [3, 4, 5, 6]\n" .
+         :02-why-7 log:outputString ?curveLine .
+         :02-why-8 log:outputString ?peakLine .
+         :02-why-9 log:outputString ?valleyLine .
+         :02-why-99 log:outputString "\n" .
+       } .
+
+# Check
+{ :case :barrierIsNarrower true .
+  :case :peakIsBeforeValley true .
+  :case :negativeDifferentialRegion true .
+  :case :overlapClosesAtHighBias true .
+  :case :peakMatchesFullOverlap true .
+}
+    => {
+         :03-check-1 log:outputString "=== Check ===\n" .
+         :03-check-2 log:outputString "heavily doped barrier is narrower : yes\n" .
+         :03-check-3 log:outputString "peak occurs before overlap closes : yes\n" .
+         :03-check-4 log:outputString "negative differential region present : yes\n" .
+         :03-check-5 log:outputString "high-bias overlap closes           : yes\n" .
+         :03-check-6 log:outputString "peak equals full four-state overlap: yes\n" .
+       } .
+
+# ------------------
+# Verification fuses
+# ------------------
+
+# If any of the key invariants fail, stop the run loudly.
+
+{ :case :tunnelDepletionWidthNm ?t .
+  :case :ordinaryDepletionWidthNm ?o .
+  ?t math:notLessThan ?o .
+} => false .
+
+{ :case :peakBias ?peakBias .
+  :case :valleyBias ?valleyBias .
+  ?peakBias math:notLessThan ?valleyBias .
+} => false .
+
+{ :b2 :currentProxy ?c2 .
+  :b3 :currentProxy ?c3 .
+  ?c2 math:notGreaterThan ?c3 .
+} => false .
+
+{ :case :valleyCurrent ?valleyCurrent .
+  ?valleyCurrent math:notEqualTo 0 .
+} => false .
+
+{ :case :peakCurrent ?peakCurrent .
+  :case :nFilledLevels ?nLevels .
+  ?nLevels list:length ?nCount .
+  ?peakCurrent math:notEqualTo ?nCount .
+} => false .