eyeling 1.5.42 → 1.6.0

package/examples/cobalt-kepler-kitchen.n3

@@ -0,0 +1,220 @@
+# ================================================================
+# Heavy-Math N3 Demo: "Cobalt Kepler Kitchen"
+# ------------------------------------------------
+# A math-stress N3 ruleset for orbital mechanics (Kepler-ish) plus
+# a Hohmann transfer, intended to exercise strict W3C/N3 math builtins.
+#
+# What it does:
+#  1) For each orbit (a,e,M):
+#     - periapsis/apopsis
+#     - period (Kepler's 3rd law)
+#     - specific orbital energy
+#     - specific angular momentum
+#     - approximate E (eccentric anomaly) from M using a series
+#     - true anomaly nu via tan-half-angle
+#     - radius r and (x,y) in orbital plane
+#  2) Computes average period across the orbit list
+#  3) Hohmann transfer between two circular radii: dv1, dv2, total dv, time
+#
+# Notes:
+#  - Uses only: math:, list:, log:
+#  - Numbers are written as numeric literals (strict typing friendly).
+# ================================================================
+
+@prefix :     <http://example.org/cobalt-kepler#> .
+@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#> .
+
+############
+# DATA (units: km, s; mu in km^3/s^2)
+############
+
+:Earth :mu 398600.4418 .
+
+:OrbitSet :orbits (
+  [ :name "LEO" ; :a 7000.0  ; :e 0.001  ; :M 1.0 ]
+  [ :name "MEO" ; :a 26560.0 ; :e 0.01   ; :M 0.3 ]
+  [ :name "GEO" ; :a 42164.0 ; :e 0.0001 ; :M 2.0 ]
+) .
+
+:Transfer1 :centralBody :Earth ;
+           :r1 7000.0 ;
+           :r2 42164.0 .
+
+############
+# (1) ORBIT METRICS + STATE AT MEAN ANOMALY
+############
+
+{
+  :Earth :mu ?mu .
+  :OrbitSet :orbits ?os .
+  ?os list:member ?orb .
+  ?orb :a ?a ; :e ?e ; :M ?M .
+
+  # periapsis/apopsis: rp=a(1-e), ra=a(1+e)
+  (1.0 ?e) math:difference ?oneMinusE .
+  (1.0 ?e) math:sum ?onePlusE .
+  (?a ?oneMinusE) math:product ?rp .
+  (?a ?onePlusE)  math:product ?ra .
+
+  # period T = 2*pi*sqrt(a^3/mu)
+  (?a 3.0) math:exponentiation ?a3 .
+  (?a3 ?mu) math:quotient ?a3OverMu .
+  (?a3OverMu 0.5) math:exponentiation ?sqrtTerm .
+  (2.0 3.141592653589793) math:product ?twoPi .
+  (?twoPi ?sqrtTerm) math:product ?T .
+
+  # specific orbital energy: eps = -mu/(2a)
+  (2.0 ?a) math:product ?twoA .
+  (?mu ?twoA) math:quotient ?muOver2a .
+  (0.0 ?muOver2a) math:difference ?eps .
+
+  # specific angular momentum: h = sqrt(mu*a*(1-e^2))
+  (?e 2.0) math:exponentiation ?e2 .
+  (1.0 ?e2) math:difference ?oneMinusE2 .
+  (?mu ?a) math:product ?muA .
+  (?muA ?oneMinusE2) math:product ?muAterm .
+  (?muAterm 0.5) math:exponentiation ?h .
+
+  # Approx eccentric anomaly from mean anomaly (series):
+  # E ≈ M + e*sin(M) + 0.5*e^2*sin(2M)
+  ?M math:sin ?sinM .
+  (?e ?sinM) math:product ?eSinM .
+  (2.0 ?M) math:product ?twoM .
+  ?twoM math:sin ?sin2M .
+  (?e2 ?sin2M) math:product ?e2Sin2M .
+  (0.5 ?e2Sin2M) math:product ?halfE2Sin2M .
+  (?M ?eSinM) math:sum ?Mplus .
+  (?Mplus ?halfE2Sin2M) math:sum ?E .
+
+  # True anomaly via tan-half-angle:
+  # tan(nu/2) = sqrt((1+e)/(1-e)) * tan(E/2)
+  (1.0 ?e) math:sum ?onePlusE2a .
+  (1.0 ?e) math:difference ?oneMinusE2a .
+  (?onePlusE2a ?oneMinusE2a) math:quotient ?ratio .
+  (?ratio 0.5) math:exponentiation ?factor .
+
+  (2.0 ?E) math:quotient ?Ehalf .
+  ?Ehalf math:tan ?tanEhalf .
+  (?factor ?tanEhalf) math:product ?tanNuHalf .
+  ?tanNuHalf math:atan ?nuHalf .
+  (2.0 ?nuHalf) math:product ?nu .
+
+  # r = a(1 - e*cos(E))
+  ?E math:cos ?cosE .
+  (?e ?cosE) math:product ?eCosE .
+  (1.0 ?eCosE) math:difference ?oneMinusECosE .
+  (?a ?oneMinusECosE) math:product ?r .
+
+  # position in orbital plane: x=r*cos(nu), y=r*sin(nu)
+  ?nu math:cos ?cosNu .
+  ?nu math:sin ?sinNu .
+  (?r ?cosNu) math:product ?x .
+  (?r ?sinNu) math:product ?y .
+}
+=>
+{
+  ?orb :periapsis ?rp ;
+       :apoapsis ?ra ;
+       :period ?T ;
+       :specificEnergy ?eps ;
+       :angularMomentum ?h ;
+       :stateAtM [
+         :M ?M ;
+         :E ?E ;
+         :nu ?nu ;
+         :r ?r ;
+         :x ?x ;
+         :y ?y
+       ] .
+} .
+
+############
+# (2) AVERAGE PERIOD ACROSS ORBITS
+############
+
+{
+  :OrbitSet :orbits ?os .
+
+  ( ?T { ?os list:member ?orb . ?orb :period ?T . } ?Ts ) log:collectAllIn _:avgScope .
+  ?Ts list:length ?n .
+  ?Ts math:sum ?sumT .
+  (?sumT ?n) math:quotient ?avgT .
+}
+=>
+{
+  :OrbitSet :avgPeriod ?avgT .
+} .
+
+############
+# (3) HOHMANN TRANSFER (circular r1 -> r2)
+#
+# aT = (r1+r2)/2
+# v1 = sqrt(mu/r1), v2 = sqrt(mu/r2)
+# vT1 = sqrt(mu*(2/r1 - 1/aT))
+# vT2 = sqrt(mu*(2/r2 - 1/aT))
+# dv1 = vT1 - v1
+# dv2 = v2 - vT2
+# dvTotal = |dv1| + |dv2|
+# tTransfer = pi * sqrt(aT^3/mu)
+############
+
+{
+  :Transfer1 :centralBody :Earth ; :r1 ?r1 ; :r2 ?r2 .
+  :Earth :mu ?mu .
+
+  (?r1 ?r2) math:sum ?r1r2 .
+  (?r1r2 2.0) math:quotient ?aT .
+
+  # v1, v2
+  (?mu ?r1) math:quotient ?muOverR1 .
+  (?muOverR1 0.5) math:exponentiation ?v1 .
+  (?mu ?r2) math:quotient ?muOverR2 .
+  (?muOverR2 0.5) math:exponentiation ?v2 .
+
+  # vT at r1: sqrt(mu*(2/r1 - 1/aT))
+  (2.0 ?r1) math:quotient ?twoOverR1 .
+  (1.0 ?aT) math:quotient ?oneOverAT .
+  (?twoOverR1 ?oneOverAT) math:difference ?term1 .
+  (?mu ?term1) math:product ?muTerm1 .
+  (?muTerm1 0.5) math:exponentiation ?vT1 .
+
+  # vT at r2
+  (2.0 ?r2) math:quotient ?twoOverR2 .
+  (?twoOverR2 ?oneOverAT) math:difference ?term2 .
+  (?mu ?term2) math:product ?muTerm2 .
+  (?muTerm2 0.5) math:exponentiation ?vT2 .
+
+  # deltas
+  (?vT1 ?v1) math:difference ?dv1 .
+  (?v2 ?vT2) math:difference ?dv2 .
+  ?dv1 math:absoluteValue ?adv1 .
+  ?dv2 math:absoluteValue ?adv2 .
+  (?adv1 ?adv2) math:sum ?dvTotal .
+
+  # transfer time: pi*sqrt(aT^3/mu)
+  (?aT 3.0) math:exponentiation ?aT3 .
+  (?aT3 ?mu) math:quotient ?aT3OverMu .
+  (?aT3OverMu 0.5) math:exponentiation ?sqrtTT .
+  (3.141592653589793 ?sqrtTT) math:product ?tTransfer .
+}
+=>
+{
+  :Transfer1 :transferSemiMajorAxis ?aT ;
+             :v1 ?v1 ; :v2 ?v2 ;
+             :vTransferAtR1 ?vT1 ;
+             :vTransferAtR2 ?vT2 ;
+             :deltaV1 ?dv1 ;
+             :deltaV2 ?dv2 ;
+             :deltaVTotal ?dvTotal ;
+             :timeOfFlight ?tTransfer .
+} .
+
+############
+# STRICTNESS REGRESSION HOOK (should fail in strict math)
+############
+# {
+#   "7000"^^<http://example.org/not-a-number> math:quotient 2.0 .
+# } => { :bad :datatype :accepted . } .
+

package/examples/jade-eigen-loom.n3

@@ -0,0 +1,322 @@
+# ================================================================
+# Heavy-Math N3 Demo: "Jade Eigen Loom"
+# ------------------------------------------------
+# A math-stress N3 ruleset for 2D PCA (principal component analysis)
+# using only standard N3 builtins (math:, list:, log:).
+#
+# What it does:
+#  1) Mean-center a cloud of 2D points and compute covariance matrix:
+#        [ covXX  covXY ]
+#        [ covXY  covYY ]
+#  2) Closed-form eigenvalues λ1, λ2 (principal variances)
+#  3) Principal axis angle θ (piecewise, without atan2)
+#  4) Per-point PCA scores (u,v) and Mahalanobis distance squared:
+#        md2 = u^2/λ1 + v^2/λ2
+#     Flags outliers if md2 > 5.991 (≈ 95% ellipse, df=2)
+#  5) Ellipse radii for 1σ and 95% (scaled by sqrt(5.991)), plus area
+#
+# Notes:
+#  - Uses only: math:, list:, log:
+#  - Designed to be strict-datatype friendly: all numbers are numeric literals.
+# ================================================================
+
+@prefix :     <http://example.org/jade-eigen-loom#> .
+@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#> .
+
+############
+# DATA
+############
+
+:PCA1 :points (
+  [ :id 1 ; :x 2.0  ; :y 1.0  ]
+  [ :id 2 ; :x 3.0  ; :y 2.0  ]
+  [ :id 3 ; :x 4.0  ; :y 3.2  ]
+  [ :id 4 ; :x 5.0  ; :y 5.1  ]
+  [ :id 5 ; :x 6.0  ; :y 7.9  ]
+  [ :id 6 ; :x 7.0  ; :y 13.0 ]
+  [ :id 7 ; :x 20.0 ; :y -3.0 ]   # outlier-ish
+) .
+
+############
+# (1) MEANS + COVARIANCE
+# covXX = Σ(dx^2)/n, covYY = Σ(dy^2)/n, covXY = Σ(dx*dy)/n  (population form)
+############
+
+{
+  :PCA1 :points ?pts .
+  ?pts list:length ?n .
+
+  ( ?x { ?pts list:member ?p . ?p :x ?x . } ?xs ) log:collectAllIn _:m1 .
+  ( ?y { ?pts list:member ?p . ?p :y ?y . } ?ys ) log:collectAllIn _:m1 .
+
+  ?xs math:sum ?sumX .
+  ?ys math:sum ?sumY .
+  (?sumX ?n) math:quotient ?meanX .
+  (?sumY ?n) math:quotient ?meanY .
+
+  ( ?dx2 {
+      ?pts list:member ?p . ?p :x ?x .
+      (?x ?meanX) math:difference ?dx .
+      (?dx 2.0) math:exponentiation ?dx2 .
+    } ?dx2s ) log:collectAllIn _:m1 .
+  ?dx2s math:sum ?ssXX .
+  (?ssXX ?n) math:quotient ?covXX .
+
+  ( ?dy2 {
+      ?pts list:member ?p . ?p :y ?y .
+      (?y ?meanY) math:difference ?dy .
+      (?dy 2.0) math:exponentiation ?dy2 .
+    } ?dy2s ) log:collectAllIn _:m1 .
+  ?dy2s math:sum ?ssYY .
+  (?ssYY ?n) math:quotient ?covYY .
+
+  ( ?dxdy {
+      ?pts list:member ?p .
+      ?p :x ?x ; :y ?y .
+      (?x ?meanX) math:difference ?dx .
+      (?y ?meanY) math:difference ?dy .
+      (?dx ?dy) math:product ?dxdy .
+    } ?dxdys ) log:collectAllIn _:m1 .
+  ?dxdys math:sum ?ssXY .
+  (?ssXY ?n) math:quotient ?covXY .
+}
+=>
+{
+  :PCA1 :n ?n ;
+        :meanX ?meanX ;
+        :meanY ?meanY ;
+        :covXX ?covXX ;
+        :covYY ?covYY ;
+        :covXY ?covXY .
+} .
+
+############
+# (2) EIGENVALUES of 2x2 symmetric covariance matrix
+# tr = a + d
+# det = a*d - b^2
+# disc = tr^2 - 4*det
+# λ1 = (tr + sqrt(disc))/2
+# λ2 = (tr - sqrt(disc))/2
+############
+
+{
+  :PCA1 :covXX ?a ; :covYY ?d ; :covXY ?b .
+
+  (?a ?d) math:sum ?tr .
+  (?a ?d) math:product ?ad .
+  (?b 2.0) math:exponentiation ?b2 .
+  (?ad ?b2) math:difference ?det .
+
+  (?tr 2.0) math:exponentiation ?tr2 .
+  (4.0 ?det) math:product ?fourDet .
+  (?tr2 ?fourDet) math:difference ?disc .
+  (?disc 0.5) math:exponentiation ?sqrtDisc .
+
+  (?tr ?sqrtDisc) math:sum ?trPlus .
+  (?tr ?sqrtDisc) math:difference ?trMinus .
+  (?trPlus 2.0) math:quotient ?lambda1 .
+  (?trMinus 2.0) math:quotient ?lambda2 .
+
+  (?lambda1 ?tr) math:quotient ?explained1 .  # λ1 / (λ1+λ2) since tr=λ1+λ2
+}
+=>
+{
+  :PCA1 :lambda1 ?lambda1 ;
+        :lambda2 ?lambda2 ;
+        :explainedVar1 ?explained1 .
+} .
+
+############
+# (3) PRINCIPAL AXIS ANGLE θ
+# θ = 0.5 * atan2(2b, a-d)
+# We emulate atan2 with piecewise rules using atan(ratio) and π adjustments.
+############
+
+# Common helpers (only when diff != 0): phi = atan((2b)/(a-d))
+{
+  :PCA1 :covXX ?a ; :covYY ?d ; :covXY ?b .
+  (?a ?d) math:difference ?diff .
+  ?diff math:notEqualTo 0.0 .
+  (2.0 ?b) math:product ?twoB .
+  (?twoB ?diff) math:quotient ?ratio .
+  ?ratio math:atan ?phi .
+}
+=>
+{
+  :PCA1 :_phi ?phi ;
+        :_diff ?diff ;
+        :_twoB ?twoB ;
+        :_b ?b .
+} .
+
+# Case 1: diff > 0  => theta = 0.5*phi
+{
+  :PCA1 :_phi ?phi ; :_diff ?diff .
+  ?diff math:greaterThan 0.0 .
+  (0.5 ?phi) math:product ?theta .
+  ?theta math:degrees ?thetaDeg .
+}
+=>
+{
+  :PCA1 :thetaRad ?theta ;
+        :thetaDeg ?thetaDeg .
+} .
+
+# Case 2: diff < 0 and b >= 0  => theta = 0.5*(phi + pi)
+{
+  :PCA1 :_phi ?phi ; :_diff ?diff ; :_b ?b .
+  ?diff math:lessThan 0.0 .
+  ?b math:notLessThan 0.0 .
+  (?phi 3.141592653589793) math:sum ?phiP .
+  (0.5 ?phiP) math:product ?theta .
+  ?theta math:degrees ?thetaDeg .
+}
+=>
+{
+  :PCA1 :thetaRad ?theta ;
+        :thetaDeg ?thetaDeg .
+} .
+
+# Case 3: diff < 0 and b < 0  => theta = 0.5*(phi - pi)
+{
+  :PCA1 :_phi ?phi ; :_diff ?diff ; :_b ?b .
+  ?diff math:lessThan 0.0 .
+  ?b math:lessThan 0.0 .
+  (?phi 3.141592653589793) math:difference ?phiM .
+  (0.5 ?phiM) math:product ?theta .
+  ?theta math:degrees ?thetaDeg .
+}
+=>
+{
+  :PCA1 :thetaRad ?theta ;
+        :thetaDeg ?thetaDeg .
+} .
+
+# Case 4: diff == 0 and b > 0 => theta = +pi/4
+{
+  :PCA1 :covXX ?a ; :covYY ?d ; :covXY ?b .
+  (?a ?d) math:difference ?diff .
+  ?diff math:equalTo 0.0 .
+  ?b math:greaterThan 0.0 .
+  0.7853981633974483 math:degrees ?thetaDeg .
+}
+=>
+{
+  :PCA1 :thetaRad 0.7853981633974483 ;
+        :thetaDeg ?thetaDeg .
+} .
+
+# Case 5: diff == 0 and b < 0 => theta = -pi/4
+{
+  :PCA1 :covXX ?a ; :covYY ?d ; :covXY ?b .
+  (?a ?d) math:difference ?diff .
+  ?diff math:equalTo 0.0 .
+  ?b math:lessThan 0.0 .
+  (0.0 0.7853981633974483) math:difference ?theta .
+  ?theta math:degrees ?thetaDeg .
+}
+=>
+{
+  :PCA1 :thetaRad ?theta ;
+        :thetaDeg ?thetaDeg .
+} .
+
+# Case 6: diff == 0 and b == 0 => theta = 0
+{
+  :PCA1 :covXX ?a ; :covYY ?d ; :covXY ?b .
+  (?a ?d) math:difference ?diff .
+  ?diff math:equalTo 0.0 .
+  ?b math:equalTo 0.0 .
+  0.0 math:degrees ?thetaDeg .
+}
+=>
+{
+  :PCA1 :thetaRad 0.0 ;
+        :thetaDeg ?thetaDeg .
+} .
+
+############
+# (4) PER-POINT PCA SCORES + MAHALANOBIS DISTANCE AND OUTLIERS
+# u = dx*cosθ + dy*sinθ
+# v = -dx*sinθ + dy*cosθ
+# md2 = u^2/λ1 + v^2/λ2
+# outlier if md2 > 5.991   (≈ 95% ellipse, df=2)
+############
+
+{
+  :PCA1 :points ?pts ;
+        :meanX ?mx ; :meanY ?my ;
+        :thetaRad ?theta ;
+        :lambda1 ?l1 ; :lambda2 ?l2 .
+
+  ?theta math:cos ?c .
+  ?theta math:sin ?s .
+
+  ?pts list:member ?p .
+  ?p :x ?x ; :y ?y .
+
+  (?x ?mx) math:difference ?dx .
+  (?y ?my) math:difference ?dy .
+
+  (?dx ?c) math:product ?dxC .
+  (?dy ?s) math:product ?dyS .
+  (?dxC ?dyS) math:sum ?u .
+
+  (?dx ?s) math:product ?dxS .
+  (0.0 ?dxS) math:difference ?negDxS .
+  (?dy ?c) math:product ?dyC .
+  (?negDxS ?dyC) math:sum ?v .
+
+  (?u 2.0) math:exponentiation ?u2 .
+  (?v 2.0) math:exponentiation ?v2 .
+  (?u2 ?l1) math:quotient ?u2Over .
+  (?v2 ?l2) math:quotient ?v2Over .
+  (?u2Over ?v2Over) math:sum ?md2 .
+}
+=>
+{
+  :PCA1 :score [ :point ?p ; :u ?u ; :v ?v ; :md2 ?md2 ] .
+} .
+
+{
+  :PCA1 :score [ :point ?p ; :u ?u ; :v ?v ; :md2 ?md2 ] .
+  ?md2 math:greaterThan 5.991 .
+}
+=>
+{
+  :PCA1 :outlier [ :point ?p ; :u ?u ; :v ?v ; :md2 ?md2 ; :threshold 5.991 ] .
+} .
+
+############
+# (5) ELLIPSE RADII (1σ and ~95%) + AREA
+############
+
+{
+  :PCA1 :lambda1 ?l1 ; :lambda2 ?l2 .
+
+  (?l1 0.5) math:exponentiation ?sigma1 .
+  (?l2 0.5) math:exponentiation ?sigma2 .
+
+  (5.991 0.5) math:exponentiation ?k95 .
+  (?k95 ?sigma1) math:product ?a95 .
+  (?k95 ?sigma2) math:product ?b95 .
+
+  (3.141592653589793 ?a95) math:product ?piA .
+  (?piA ?b95) math:product ?area95 .
+}
+=>
+{
+  :PCA1 :sigma1 ?sigma1 ;
+        :sigma2 ?sigma2 ;
+        :ellipse95 [ :k ?k95 ; :a ?a95 ; :b ?b95 ; :area ?area95 ] .
+} .
+
+############
+# STRICTNESS REGRESSION HOOK (should fail in strict math)
+############
+# {
+#   "1"^^<http://example.org/not-a-number> math:sum 2.0 .
+# } => { :bad :datatype :accepted . } .
+

package/examples/not-includes.n3

@@ -0,0 +1,32 @@
+# =======================
+# log:notIncludes example
+# =======================
+
+@prefix : <http://example.org/> .
+@prefix log: <http://www.w3.org/2000/10/swap/log#> .
+@prefix xsd: <http://www.w3.org/2001/XMLSchema#> .
+@prefix math: <http://www.w3.org/2000/10/swap/math#> .
+
+{
+  1 math:equalTo 1.
+}
+=>
+{
+  :x :y 1 .
+  
+  {
+    :x :y 1
+  }
+  => 
+  { 
+    :a :b :c
+  }.
+}.
+
+{
+  ?SCOPE log:notIncludes { :a :b :c }.
+}
+=>
+{
+  :test :is false 
+}.