eyeling 1.11.22 → 1.11.23

package/examples/matiyasevich-pi-fib.n3

@@ -0,0 +1,140 @@
+# ==========================================================
+# Matiyasevich (1985): π from Fibonacci numbers
+#
+# Illustrates the identity:
+#   π = 4 * Σ_{n≥1} atan( 1 / F_{2n+1} )
+# where F_k is the k-th Fibonacci number.
+#
+# This N3 program computes a finite partial sum using the list
+# of n values in :run :ns, and compares it with π = 4*atan(1).
+#
+# Notes:
+# - To keep reasoners terminating, Fibonacci computation is gated:
+#   only indices reachable from the required {2n+1} set are derived.
+# ==========================================================
+
+@prefix :     <http://example.org/matiyasevich#>.
+@prefix math: <http://www.w3.org/2000/10/swap/math#>.
+@prefix list: <http://www.w3.org/2000/10/swap/list#>.
+@prefix log:  <http://www.w3.org/2000/10/swap/log#>.
+
+:run :ns (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15).
+
+# --- Determine which Fibonacci indices are needed: k = 2n+1 ----------
+{
+  :run :ns ?Ns0.
+  ?Ns0 list:iterate (?i0 ?n0).
+
+  (2 ?n0) math:product ?twoN0.
+  (?twoN0 1) math:sum ?k0.
+}
+=>
+{
+  :need :fibIndex ?k0.
+}.
+
+# --- Close the needed set downward: if we need n, we also need n-1 and n-2 ----
+{
+  :need :fibIndex ?nN.
+  ?nN math:greaterThan 1.
+
+  (?nN 1) math:difference ?n1N.
+  (?nN 2) math:difference ?n2N.
+}
+=>
+{
+  :need :fibIndex ?n1N.
+  :need :fibIndex ?n2N.
+}.
+
+# --- Fibonacci facts, computed only for indices requested via :need :fibIndex ----
+{
+  :need :fibIndex 0.
+}
+=>
+{
+  0 :fib 0.
+}.
+
+{
+  :need :fibIndex 1.
+}
+=>
+{
+  1 :fib 1.
+}.
+
+{
+  :need :fibIndex ?nF.
+  ?nF math:greaterThan 1.
+
+  (?nF 1) math:difference ?n1F.
+  (?nF 2) math:difference ?n2F.
+
+  ?n1F :fib ?f1F.
+  ?n2F :fib ?f2F.
+  (?f1F ?f2F) math:sum ?fnF.
+}
+=>
+{
+  ?nF :fib ?fnF.
+}.
+
+# --- Final: collect terms, sum angles, approximate pi, and compute error ----
+{
+  :run :ns ?Ns.
+
+  # Rows: (n k F_k theta_n) with theta_n = atan(1/F_k)
+  ( (?n ?k ?F ?theta)
+    {
+      :run :ns ?Ns.
+      ?Ns list:iterate (?i ?n).
+
+      (2 ?n) math:product ?twoN.
+      (?twoN 1) math:sum ?k.
+
+      ?k :fib ?F.
+      (1 ?F) math:quotient ?inv.
+      ?inv math:atan ?theta.
+    }
+    ?rows
+  ) log:collectAllIn _:rowsScope.
+
+  # Angles only, to sum
+  ( ?theta2
+    {
+      :run :ns ?Ns.
+      ?Ns list:iterate (?i2 ?n2).
+
+      (2 ?n2) math:product ?twoN2.
+      (?twoN2 1) math:sum ?k2.
+
+      ?k2 :fib ?F2.
+      (1 ?F2) math:quotient ?inv2.
+      ?inv2 math:atan ?theta2.
+    }
+    ?thetas
+  ) log:collectAllIn _:thetaScope.
+
+  # A small “grounding” guard: list:length requires the list input be present.
+  ?thetas list:length ?len.
+
+  ?thetas math:sum ?sumAngles.
+  (4 ?sumAngles) math:product ?piFromFibs.
+
+  # Reference value: atan(1) = pi/4  =>  pi = 4*atan(1)
+  1 math:atan ?atan1.
+  (4 ?atan1) math:product ?piRef.
+
+  (?piFromFibs ?piRef) math:difference ?delta.
+  ?delta math:absoluteValue ?absError.
+}
+=>
+{
+  :result
+    :rows       ?rows;
+    :piFromFibs ?piFromFibs;
+    :piRef      ?piRef;
+    :absError   ?absError.
+}.
+

package/examples/output/matiyasevich-pi-fib.n3

@@ -0,0 +1,71 @@
+@prefix : <http://example.org/matiyasevich#> .
+@prefix xsd: <http://www.w3.org/2001/XMLSchema#> .
+
+:need :fibIndex 3 .
+:need :fibIndex 5 .
+:need :fibIndex 7 .
+:need :fibIndex 9 .
+:need :fibIndex 11 .
+:need :fibIndex 13 .
+:need :fibIndex 15 .
+:need :fibIndex 17 .
+:need :fibIndex 19 .
+:need :fibIndex 21 .
+:need :fibIndex 23 .
+:need :fibIndex 25 .
+:need :fibIndex 27 .
+:need :fibIndex 29 .
+:need :fibIndex 31 .
+:need :fibIndex 2 .
+:need :fibIndex 1 .
+:need :fibIndex 4 .
+:need :fibIndex 6 .
+:need :fibIndex 8 .
+:need :fibIndex 10 .
+:need :fibIndex 12 .
+:need :fibIndex 14 .
+:need :fibIndex 16 .
+:need :fibIndex 18 .
+:need :fibIndex 20 .
+:need :fibIndex 22 .
+:need :fibIndex 24 .
+:need :fibIndex 26 .
+:need :fibIndex 28 .
+:need :fibIndex 30 .
+1 :fib 1 .
+:need :fibIndex 0 .
+0 :fib 0 .
+2 :fib 1 .
+3 :fib 2 .
+4 :fib 3 .
+5 :fib 5 .
+6 :fib 8 .
+7 :fib 13 .
+8 :fib 21 .
+9 :fib 34 .
+10 :fib 55 .
+11 :fib 89 .
+12 :fib 144 .
+13 :fib 233 .
+14 :fib 377 .
+15 :fib 610 .
+16 :fib 987 .
+17 :fib 1597 .
+18 :fib 2584 .
+19 :fib 4181 .
+20 :fib 6765 .
+21 :fib 10946 .
+22 :fib 17711 .
+23 :fib 28657 .
+24 :fib 46368 .
+25 :fib 75025 .
+26 :fib 121393 .
+27 :fib 196418 .
+28 :fib 317811 .
+29 :fib 514229 .
+30 :fib 832040 .
+31 :fib 1346269 .
+:result :rows ((1 3 2 "0.4636476090008061"^^xsd:decimal) (2 5 5 "0.19739555984988078"^^xsd:decimal) (3 7 13 "0.07677189126977804"^^xsd:decimal) (4 9 34 "0.029403288204005115"^^xsd:decimal) (5 11 89 "0.011235482257962729"^^xsd:decimal) (6 13 233 "0.004291819142011206"^^xsd:decimal) (7 15 610 "0.0016393427937457503"^^xsd:decimal) (8 17 1597 "0.0006261739945535625"^^xsd:decimal) (9 19 4181 "0.00023917722576690222"^^xsd:decimal) (10 21 10946 "0.00009135757328868366"^^xsd:decimal) (11 23 28657 "0.000034895487999235846"^^xsd:decimal) (12 25 75025 "0.000013328890369087375"^^xsd:decimal) (13 27 196418 "0.0000050911830890822634"^^xsd:decimal) (14 29 514229 "0.000001944658897103702"^^xsd:decimal) (15 31 1346269 "7.427936021700092e-7"^^xsd:decimal)) .
+:result :piFromFibs "3.1415908173030225"^^xsd:decimal .
+:result :piRef "3.141592653589793"^^xsd:decimal .
+:result :absError "0.000001836286770640072"^^xsd:decimal .

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "eyeling",
-  "version": "1.11.22",
+  "version": "1.11.23",
   "description": "A minimal Notation3 (N3) reasoner in JavaScript.",
   "main": "./index.js",
   "keywords": [

package/test/manifest.test.js

@@ -109,15 +109,7 @@ function rmrf(p) {
   }
   ok(`npm ci ${C.dim}(${r.ms} ms)${C.n}`);
 
-  // 3) Build + pack local Eyeling
-  r = run('npm', ['run', 'build'], { cwd: ROOT });
-  if (r.status !== 0) {
-    fail(`npm run build failed (exit ${r.status})`);
-    rmrf(workDir);
-    return;
-  }
-  ok(`built eyeling ${C.dim}(${r.ms} ms)${C.n}`);
-
+  // 3) Pack local Eyeling
   const pack = runCapture('npm', ['pack', '--silent'], {
     cwd: ROOT,
     stdio: ['ignore', 'pipe', 'pipe'],