eyeling 1.15.6 → 1.15.7

package/examples/fft32-numeric.n3

@@ -0,0 +1,302 @@
+# ===================================================
+# Fast Fourier Transform (FFT)
+# Numeric version with explicit complex pairs (re im)
+# 32nd roots of unity are baked in as constants
+# Ccompute the whole spectrum once per waveform
+# ===================================================
+
+@prefix : <http://example.org/fft32#>.
+@prefix math: <http://www.w3.org/2000/10/swap/math#>.
+@prefix log: <http://www.w3.org/2000/10/swap/log#>.
+
+# ----------------------------------------------
+# 32nd roots of unity: W32^k = exp(-2*pi*i*k/32)
+# represented as (re im)
+# ----------------------------------------------
+0 :w32 (1.0 0.0).
+1 :w32 (0.98078528040323043 -0.19509032201612825).
+2 :w32 (0.92387953251128674 -0.38268343236508978).
+3 :w32 (0.83146961230254524 -0.55557023301960218).
+4 :w32 (0.70710678118654757 -0.70710678118654746).
+5 :w32 (0.55557023301960229 -0.83146961230254524).
+6 :w32 (0.38268343236508984 -0.92387953251128674).
+7 :w32 (0.19509032201612833 -0.98078528040323043).
+8 :w32 (0.0 -1.0).
+9 :w32 (-0.19509032201612819 -0.98078528040323043).
+10 :w32 (-0.38268343236508973 -0.92387953251128674).
+11 :w32 (-0.55557023301960196 -0.83146961230254546).
+12 :w32 (-0.70710678118654746 -0.70710678118654757).
+13 :w32 (-0.83146961230254535 -0.55557023301960218).
+14 :w32 (-0.92387953251128674 -0.38268343236508989).
+15 :w32 (-0.98078528040323043 -0.19509032201612861).
+16 :w32 (-1.0 0.0).
+17 :w32 (-0.98078528040323043 0.19509032201612836).
+18 :w32 (-0.92387953251128685 0.38268343236508967).
+19 :w32 (-0.83146961230254546 0.55557023301960196).
+20 :w32 (-0.70710678118654768 0.70710678118654746).
+21 :w32 (-0.55557023301960218 0.83146961230254524).
+22 :w32 (-0.38268343236509034 0.92387953251128652).
+23 :w32 (-0.19509032201612866 0.98078528040323032).
+24 :w32 (0.0 1.0).
+25 :w32 (0.1950903220161283 0.98078528040323043).
+26 :w32 (0.38268343236509 0.92387953251128663).
+27 :w32 (0.55557023301960184 0.83146961230254546).
+28 :w32 (0.70710678118654735 0.70710678118654768).
+29 :w32 (0.83146961230254524 0.55557023301960218).
+30 :w32 (0.92387953251128652 0.38268343236509039).
+31 :w32 (0.98078528040323032 0.19509032201612872).
+
+# --------------------------------------------
+# Complex arithmetic on explicit pairs (re im)
+# --------------------------------------------
+
+# (a+bi) + (c+di) = (a+c) + (b+d)i
+{ ((?ar ?ai) (?br ?bi)) :cAdd (?cr ?ci). } <= { (?ar ?br) math:sum ?cr. (?ai ?bi) math:sum ?ci. }.
+
+# (a+bi) - (c+di) = (a-c) + (b-d)i
+{ ((?ar ?ai) (?br ?bi)) :cSub (?cr ?ci). } <= { (?ar ?br) math:difference ?cr. (?ai ?bi) math:difference ?ci. }.
+
+# (a+bi) * (c+di) = (ac-bd) + (ad+bc)i
+{ ((?ar ?ai) (?br ?bi)) :cMul (?cr ?ci). } <= { (?ar ?br) math:product ?ac. (?ai ?bi) math:product ?bd. (?ac ?bd) math:difference ?cr. (?ar ?bi) math:product ?ad. (?ai ?br) math:product ?bc. (?ad ?bc) math:sum ?ci. }.
+
+# -------------------------
+# Numeric FFT vector rules
+# -------------------------
+
+# Base case: one real sample x becomes the complex number (x 0)
+{ (?x) :fft1 (?x 0.0). } <= {}.
+
+# ---------
+# Length 2
+# ---------
+
+{
+  (?x0 ?x1) :fft2 (?y0 ?y1).
+} <= {
+  (?x0) :fft1 ?e0.
+  (?x1) :fft1 ?o0.
+
+  (?e0 ?o0) :cAdd ?y0.
+  (?e0 ?o0) :cSub ?y1.
+}.
+
+# ---------
+# Length 4
+# ---------
+
+{
+  (?x0 ?x1 ?x2 ?x3) :fft4 (?y0 ?y1 ?y2 ?y3).
+} <= {
+  (?x0 ?x2) :fft2 (?e0 ?e1).
+  (?x1 ?x3) :fft2 (?o0 ?o1).
+
+  (?e0 ?o0) :cAdd ?y0.
+  (?e0 ?o0) :cSub ?y2.
+
+  8 :w32 ?w1.
+  (?w1 ?o1) :cMul ?t1.
+  (?e1 ?t1) :cAdd ?y1.
+  (?e1 ?t1) :cSub ?y3.
+}.
+
+# ---------
+# Length 8
+# ---------
+
+{
+  (?x0 ?x1 ?x2 ?x3 ?x4 ?x5 ?x6 ?x7) :fft8 (?y0 ?y1 ?y2 ?y3 ?y4 ?y5 ?y6 ?y7).
+} <= {
+  (?x0 ?x2 ?x4 ?x6) :fft4 (?e0 ?e1 ?e2 ?e3).
+  (?x1 ?x3 ?x5 ?x7) :fft4 (?o0 ?o1 ?o2 ?o3).
+
+  (?e0 ?o0) :cAdd ?y0.
+  (?e0 ?o0) :cSub ?y4.
+
+  4 :w32 ?w1.
+  (?w1 ?o1) :cMul ?t1.
+  (?e1 ?t1) :cAdd ?y1.
+  (?e1 ?t1) :cSub ?y5.
+
+  8 :w32 ?w2.
+  (?w2 ?o2) :cMul ?t2.
+  (?e2 ?t2) :cAdd ?y2.
+  (?e2 ?t2) :cSub ?y6.
+
+  12 :w32 ?w3.
+  (?w3 ?o3) :cMul ?t3.
+  (?e3 ?t3) :cAdd ?y3.
+  (?e3 ?t3) :cSub ?y7.
+}.
+
+# ---------
+# Length 16
+# ---------
+
+{
+  (?x0 ?x1 ?x2 ?x3 ?x4 ?x5 ?x6 ?x7 ?x8 ?x9 ?x10 ?x11 ?x12 ?x13 ?x14 ?x15) :fft16 (?y0 ?y1 ?y2 ?y3 ?y4 ?y5 ?y6 ?y7 ?y8 ?y9 ?y10 ?y11 ?y12 ?y13 ?y14 ?y15).
+} <= {
+  (?x0 ?x2 ?x4 ?x6 ?x8 ?x10 ?x12 ?x14) :fft8 (?e0 ?e1 ?e2 ?e3 ?e4 ?e5 ?e6 ?e7).
+  (?x1 ?x3 ?x5 ?x7 ?x9 ?x11 ?x13 ?x15) :fft8 (?o0 ?o1 ?o2 ?o3 ?o4 ?o5 ?o6 ?o7).
+
+  (?e0 ?o0) :cAdd ?y0.
+  (?e0 ?o0) :cSub ?y8.
+
+  2 :w32 ?w1.
+  (?w1 ?o1) :cMul ?t1.
+  (?e1 ?t1) :cAdd ?y1.
+  (?e1 ?t1) :cSub ?y9.
+
+  4 :w32 ?w2.
+  (?w2 ?o2) :cMul ?t2.
+  (?e2 ?t2) :cAdd ?y2.
+  (?e2 ?t2) :cSub ?y10.
+
+  6 :w32 ?w3.
+  (?w3 ?o3) :cMul ?t3.
+  (?e3 ?t3) :cAdd ?y3.
+  (?e3 ?t3) :cSub ?y11.
+
+  8 :w32 ?w4.
+  (?w4 ?o4) :cMul ?t4.
+  (?e4 ?t4) :cAdd ?y4.
+  (?e4 ?t4) :cSub ?y12.
+
+  10 :w32 ?w5.
+  (?w5 ?o5) :cMul ?t5.
+  (?e5 ?t5) :cAdd ?y5.
+  (?e5 ?t5) :cSub ?y13.
+
+  12 :w32 ?w6.
+  (?w6 ?o6) :cMul ?t6.
+  (?e6 ?t6) :cAdd ?y6.
+  (?e6 ?t6) :cSub ?y14.
+
+  14 :w32 ?w7.
+  (?w7 ?o7) :cMul ?t7.
+  (?e7 ?t7) :cAdd ?y7.
+  (?e7 ?t7) :cSub ?y15.
+}.
+
+# ---------
+# Length 32
+# ---------
+
+{
+  (?x0 ?x1 ?x2 ?x3 ?x4 ?x5 ?x6 ?x7 ?x8 ?x9 ?x10 ?x11 ?x12 ?x13 ?x14 ?x15 ?x16 ?x17 ?x18 ?x19 ?x20 ?x21 ?x22 ?x23 ?x24 ?x25 ?x26 ?x27 ?x28 ?x29 ?x30 ?x31) :fft32 (?y0 ?y1 ?y2 ?y3 ?y4 ?y5 ?y6 ?y7 ?y8 ?y9 ?y10 ?y11 ?y12 ?y13 ?y14 ?y15 ?y16 ?y17 ?y18 ?y19 ?y20 ?y21 ?y22 ?y23 ?y24 ?y25 ?y26 ?y27 ?y28 ?y29 ?y30 ?y31).
+} <= {
+  (?x0 ?x2 ?x4 ?x6 ?x8 ?x10 ?x12 ?x14 ?x16 ?x18 ?x20 ?x22 ?x24 ?x26 ?x28 ?x30) :fft16 (?e0 ?e1 ?e2 ?e3 ?e4 ?e5 ?e6 ?e7 ?e8 ?e9 ?e10 ?e11 ?e12 ?e13 ?e14 ?e15).
+  (?x1 ?x3 ?x5 ?x7 ?x9 ?x11 ?x13 ?x15 ?x17 ?x19 ?x21 ?x23 ?x25 ?x27 ?x29 ?x31) :fft16 (?o0 ?o1 ?o2 ?o3 ?o4 ?o5 ?o6 ?o7 ?o8 ?o9 ?o10 ?o11 ?o12 ?o13 ?o14 ?o15).
+
+  (?e0 ?o0) :cAdd ?y0.
+  (?e0 ?o0) :cSub ?y16.
+
+  1 :w32 ?w1.
+  (?w1 ?o1) :cMul ?t1.
+  (?e1 ?t1) :cAdd ?y1.
+  (?e1 ?t1) :cSub ?y17.
+
+  2 :w32 ?w2.
+  (?w2 ?o2) :cMul ?t2.
+  (?e2 ?t2) :cAdd ?y2.
+  (?e2 ?t2) :cSub ?y18.
+
+  3 :w32 ?w3.
+  (?w3 ?o3) :cMul ?t3.
+  (?e3 ?t3) :cAdd ?y3.
+  (?e3 ?t3) :cSub ?y19.
+
+  4 :w32 ?w4.
+  (?w4 ?o4) :cMul ?t4.
+  (?e4 ?t4) :cAdd ?y4.
+  (?e4 ?t4) :cSub ?y20.
+
+  5 :w32 ?w5.
+  (?w5 ?o5) :cMul ?t5.
+  (?e5 ?t5) :cAdd ?y5.
+  (?e5 ?t5) :cSub ?y21.
+
+  6 :w32 ?w6.
+  (?w6 ?o6) :cMul ?t6.
+  (?e6 ?t6) :cAdd ?y6.
+  (?e6 ?t6) :cSub ?y22.
+
+  7 :w32 ?w7.
+  (?w7 ?o7) :cMul ?t7.
+  (?e7 ?t7) :cAdd ?y7.
+  (?e7 ?t7) :cSub ?y23.
+
+  8 :w32 ?w8.
+  (?w8 ?o8) :cMul ?t8.
+  (?e8 ?t8) :cAdd ?y8.
+  (?e8 ?t8) :cSub ?y24.
+
+  9 :w32 ?w9.
+  (?w9 ?o9) :cMul ?t9.
+  (?e9 ?t9) :cAdd ?y9.
+  (?e9 ?t9) :cSub ?y25.
+
+  10 :w32 ?w10.
+  (?w10 ?o10) :cMul ?t10.
+  (?e10 ?t10) :cAdd ?y10.
+  (?e10 ?t10) :cSub ?y26.
+
+  11 :w32 ?w11.
+  (?w11 ?o11) :cMul ?t11.
+  (?e11 ?t11) :cAdd ?y11.
+  (?e11 ?t11) :cSub ?y27.
+
+  12 :w32 ?w12.
+  (?w12 ?o12) :cMul ?t12.
+  (?e12 ?t12) :cAdd ?y12.
+  (?e12 ?t12) :cSub ?y28.
+
+  13 :w32 ?w13.
+  (?w13 ?o13) :cMul ?t13.
+  (?e13 ?t13) :cAdd ?y13.
+  (?e13 ?t13) :cSub ?y29.
+
+  14 :w32 ?w14.
+  (?w14 ?o14) :cMul ?t14.
+  (?e14 ?t14) :cAdd ?y14.
+  (?e14 ?t14) :cSub ?y30.
+
+  15 :w32 ?w15.
+  (?w15 ?o15) :cMul ?t15.
+  (?e15 ?t15) :cAdd ?y15.
+  (?e15 ?t15) :cSub ?y31.
+}.
+
+# ---------
+# waveforms
+# ---------
+
+# alternating
+(1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1) a :WaveForm.
+
+# constant ones
+(1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1) a :WaveForm.
+
+# cosine bin3
+(1 0.83146961230254501 0.38268343236509 -0.195090322016128 -0.70710678118654702 -0.98078528040322999 -0.92387953251128696 -0.55557023301960196 0 0.55557023301960196 0.92387953251128696 0.98078528040322999 0.70710678118654802 0.195090322016129 -0.38268343236509 -0.83146961230254501 -1 -0.83146961230254601 -0.382683432365091 0.195090322016127 0.70710678118654702 0.98078528040322999 0.92387953251128696 0.55557023301960295 0.000000000000001 -0.55557023301960196 -0.92387953251128596 -0.98078528040323099 -0.70710678118654702 -0.19509032201613 0.38268343236509 0.83146961230254401) a :WaveForm.
+
+# impulse
+(1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0) a :WaveForm.
+
+# ramp 0-32
+(0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31) a :WaveForm.
+
+# sine bin5
+(0 0.83146961230254501 0.92387953251128696 0.195090322016129 -0.70710678118654702 -0.98078528040322999 -0.38268343236509 0.55557023301960196 1 0.55557023301960196 -0.38268343236509 -0.98078528040323099 -0.70710678118654802 0.195090322016127 0.92387953251128596 0.83146961230254601 0.000000000000001 -0.83146961230254501 -0.92387953251128696 -0.195090322016128 0.70710678118654802 0.98078528040322999 0.382683432365089 -0.55557023301960295 -1 -0.55557023301960395 0.38268343236508801 0.98078528040322999 0.70710678118654902 -0.195090322016127 -0.92387953251128596 -0.83146961230254601) a :WaveForm.
+
+# -----
+# Query
+# -----
+
+{
+  ?in a :WaveForm.
+  ?in :fft32 ?out.
+}
+log:query
+{
+  ?in :fft ?out.
+}.

package/examples/output/fft32-numeric.n3

@@ -0,0 +1,9 @@
+@prefix : <http://example.org/fft32#> .
+@prefix xsd: <http://www.w3.org/2001/XMLSchema#> .
+
+(0 0.83146961230254501 0.92387953251128696 0.195090322016129 -0.70710678118654702 -0.98078528040322999 -0.38268343236509 0.55557023301960196 1 0.55557023301960196 -0.38268343236509 -0.98078528040323099 -0.70710678118654802 0.195090322016127 0.92387953251128596 0.83146961230254601 0.000000000000001 -0.83146961230254501 -0.92387953251128696 -0.195090322016128 0.70710678118654802 0.98078528040322999 0.382683432365089 -0.55557023301960295 -1 -0.55557023301960395 0.38268343236508801 0.98078528040322999 0.70710678118654902 -0.195090322016127 -0.92387953251128596 -0.83146961230254601) :fft (("-2.9968028886505635e-15"^^xsd:decimal "0"^^xsd:decimal) ("2.7410764551198392e-15"^^xsd:decimal "-2.0320206323861304e-15"^^xsd:decimal) ("6.653803231410531e-15"^^xsd:decimal "-2.935201360539696e-15"^^xsd:decimal) ("-3.263161509486249e-15"^^xsd:decimal "5.895742328325398e-15"^^xsd:decimal) ("-3.118026263548378e-15"^^xsd:decimal "1.7057423285703395e-15"^^xsd:decimal) ("-8.327471962526034e-15"^^xsd:decimal "-16"^^xsd:decimal) ("-1.8187743254395905e-15"^^xsd:decimal "-8.649275785265081e-16"^^xsd:decimal) ("-7.698815671630478e-16"^^xsd:decimal "-5.907250600715114e-15"^^xsd:decimal) ("5.996003610813205e-15"^^xsd:decimal "9.992007221626409e-16"^^xsd:decimal) ("-1.046010004595602e-15"^^xsd:decimal "7.32799883525688e-17"^^xsd:decimal) ("-4.2047531300660247e-16"^^xsd:decimal "5.481556342888897e-16"^^xsd:decimal) ("3.4408920985006257e-15"^^xsd:decimal "0"^^xsd:decimal) ("1.1212233748978148e-15"^^xsd:decimal "-2.9265911575494193e-16"^^xsd:decimal) ("1.6825725784634528e-15"^^xsd:decimal "-8.40438346133828e-16"^^xsd:decimal) ("-4.145535929643375e-16"^^xsd:decimal "-1.522118147724299e-15"^^xsd:decimal) ("1.3184849660375944e-15"^^xsd:decimal "-1.5590910214806506e-15"^^xsd:decimal) ("2.998401444325282e-15"^^xsd:decimal "0"^^xsd:decimal) ("9.36240855879889e-16"^^xsd:decimal "1.5259059172126727e-15"^^xsd:decimal) ("-4.145535929643375e-16"^^xsd:decimal "1.5221181477242987e-15"^^xsd:decimal) ("2.151339929186374e-15"^^xsd:decimal "1.2096850292756035e-15"^^xsd:decimal) ("1.1212233748978148e-15"^^xsd:decimal "2.926591157549422e-16"^^xsd:decimal) ("5.543122344752187e-16"^^xsd:decimal "0"^^xsd:decimal) ("-4.2047531300660207e-16"^^xsd:decimal "-5.481556342888898e-16"^^xsd:decimal) ("-5.348202489988052e-16"^^xsd:decimal "8.574383478270009e-16"^^xsd:decimal) ("5.996003610813205e-15"^^xsd:decimal "-9.992007221626409e-16"^^xsd:decimal) ("-1.3022367882033748e-15"^^xsd:decimal "5.761905245021641e-15"^^xsd:decimal) ("-1.8187743254395905e-15"^^xsd:decimal "8.649275785265076e-16"^^xsd:decimal) ("-1.3434497875801754e-14"^^xsd:decimal "16"^^xsd:decimal) ("-3.118026263548378e-15"^^xsd:decimal "-1.7057423285703398e-15"^^xsd:decimal) ("-3.2384833686133903e-15"^^xsd:decimal "-6.2649890114671735e-15"^^xsd:decimal) ("6.653803231410531e-15"^^xsd:decimal "2.935201360539697e-15"^^xsd:decimal) ("3.0916442077252603e-15"^^xsd:decimal "1.2798327561680123e-15"^^xsd:decimal)) .
+(0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31) :fft ((496 "0"^^xsd:decimal) ("-16.000000000000004"^^xsd:decimal "162.45072620174173"^^xsd:decimal) ("-16"^^xsd:decimal "80.43743187401355"^^xsd:decimal) ("-16.000000000000007"^^xsd:decimal "52.74493134301312"^^xsd:decimal) ("-16"^^xsd:decimal "38.62741699796952"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "29.93389458863023"^^xsd:decimal) ("-16.000000000000004"^^xsd:decimal "23.945692202647823"^^xsd:decimal) ("-16.000000000000004"^^xsd:decimal "19.49605640940763"^^xsd:decimal) ("-16"^^xsd:decimal "16"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "13.13086065325857"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "10.690858206708782"^^xsd:decimal) ("-16.000000000000007"^^xsd:decimal "8.552178175212664"^^xsd:decimal) ("-16"^^xsd:decimal "6.627416997969519"^^xsd:decimal) ("-15.999999999999995"^^xsd:decimal "4.853546937717482"^^xsd:decimal) ("-16.000000000000004"^^xsd:decimal "3.1825978780745316"^^xsd:decimal) ("-16.00000000000003"^^xsd:decimal "1.5758624537146346"^^xsd:decimal) (-16 "0"^^xsd:decimal) ("-15.999999999999996"^^xsd:decimal "-1.5758624537146204"^^xsd:decimal) ("-16"^^xsd:decimal "-3.1825978780745245"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "-4.8535469377174785"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "-6.627416997969519"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "-8.552178175212667"^^xsd:decimal) ("-15.999999999999996"^^xsd:decimal "-10.690858206708786"^^xsd:decimal) ("-16.000000000000004"^^xsd:decimal "-13.130860653258566"^^xsd:decimal) ("-16"^^xsd:decimal "-16"^^xsd:decimal) ("-16.000000000000004"^^xsd:decimal "-19.49605640940762"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "-23.94569220264782"^^xsd:decimal) ("-15.999999999999986"^^xsd:decimal "-29.933894588630235"^^xsd:decimal) ("-15.999999999999998"^^xsd:decimal "-38.62741699796952"^^xsd:decimal) ("-16"^^xsd:decimal "-52.74493134301312"^^xsd:decimal) ("-15.999999999999993"^^xsd:decimal "-80.43743187401357"^^xsd:decimal) ("-15.999999999999957"^^xsd:decimal "-162.45072620174176"^^xsd:decimal)) .
+(1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1) :fft ((0 "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) (32 "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal)) .
+(1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0) :fft ((1 "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) (1 "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal) ("1"^^xsd:decimal "0"^^xsd:decimal)) .
+(1 0.83146961230254501 0.38268343236509 -0.195090322016128 -0.70710678118654702 -0.98078528040322999 -0.92387953251128696 -0.55557023301960196 0 0.55557023301960196 0.92387953251128696 0.98078528040322999 0.70710678118654802 0.195090322016129 -0.38268343236509 -0.83146961230254501 -1 -0.83146961230254601 -0.382683432365091 0.195090322016127 0.70710678118654702 0.98078528040322999 0.92387953251128696 0.55557023301960295 0.000000000000001 -0.55557023301960196 -0.92387953251128596 -0.98078528040323099 -0.70710678118654702 -0.19509032201613 0.38268343236509 0.83146961230254401) :fft (("-1.9976021664879226e-15"^^xsd:decimal "0"^^xsd:decimal) ("1.3737538074197586e-15"^^xsd:decimal "-1.946465856829794e-15"^^xsd:decimal) ("-5.564884063023956e-15"^^xsd:decimal "1.1067652769226016e-15"^^xsd:decimal) ("15.999999999999998"^^xsd:decimal "-7.66133814775094e-15"^^xsd:decimal) ("1.4138824906527567e-15"^^xsd:decimal "1.4130832128153977e-15"^^xsd:decimal) ("-2.705341196503391e-16"^^xsd:decimal "-1.5008655956026007e-15"^^xsd:decimal) ("1.890916229553308e-16"^^xsd:decimal "-1.2688104344492827e-16"^^xsd:decimal) ("-1.431798885729609e-16"^^xsd:decimal "6.382299624403909e-17"^^xsd:decimal) ("1.999200722162641e-15"^^xsd:decimal "0"^^xsd:decimal) ("-1.974143070189294e-15"^^xsd:decimal "1.2758998416373139e-16"^^xsd:decimal) ("6.370748026754646e-16"^^xsd:decimal "-9.546460247504417e-16"^^xsd:decimal) ("-2.8857717599356964e-15"^^xsd:decimal "-1.4169280937773732e-15"^^xsd:decimal) ("-1.4122839349780385e-15"^^xsd:decimal "1.4130832128153975e-15"^^xsd:decimal) ("8.881784197001252e-16"^^xsd:decimal "-1.55351295663786e-15"^^xsd:decimal) ("7.38717637393161e-16"^^xsd:decimal "-3.717802593033475e-15"^^xsd:decimal) ("3.3987495982314057e-15"^^xsd:decimal "-3.4349527686803807e-15"^^xsd:decimal) ("5.996003610813205e-15"^^xsd:decimal "0"^^xsd:decimal) ("3.919183002821392e-15"^^xsd:decimal "3.338931990141291e-15"^^xsd:decimal) ("7.387176373931606e-16"^^xsd:decimal "3.717802593033475e-15"^^xsd:decimal) ("1.7763568394002505e-15"^^xsd:decimal "2.9968028886505635e-15"^^xsd:decimal) ("-1.4122839349780383e-15"^^xsd:decimal "-1.4130832128153977e-15"^^xsd:decimal) ("-1.5058227197499114e-15"^^xsd:decimal "1.946553361127382e-15"^^xsd:decimal) ("6.370748026754641e-16"^^xsd:decimal "9.54646024750442e-16"^^xsd:decimal) ("-1.6322401255607569e-15"^^xsd:decimal "-7.241239318976423e-16"^^xsd:decimal) ("1.999200722162641e-15"^^xsd:decimal "0"^^xsd:decimal) ("2.339199387486444e-16"^^xsd:decimal "-6.286805864256665e-16"^^xsd:decimal) ("1.8909162295533075e-16"^^xsd:decimal "1.2688104344492837e-16"^^xsd:decimal) ("-6.669419188648045e-16"^^xsd:decimal "3.0510651347749833e-16"^^xsd:decimal) ("1.413882490652757e-15"^^xsd:decimal "-1.4130832128153975e-15"^^xsd:decimal) ("15.999999999999996"^^xsd:decimal "8.216449660063518e-15"^^xsd:decimal) ("-5.564884063023956e-15"^^xsd:decimal "-1.1067652769226028e-15"^^xsd:decimal) ("1.929384094702813e-15"^^xsd:decimal "1.8716105437342336e-15"^^xsd:decimal)) .
+(1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1) :fft ((32 "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) (0 "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal) ("0"^^xsd:decimal "0"^^xsd:decimal)) .

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "eyeling",
-  "version": "1.15.6",
+  "version": "1.15.7",
   "description": "A minimal Notation3 (N3) reasoner in JavaScript.",
   "main": "./index.js",
   "keywords": [
@@ -42,11 +42,10 @@
     "test:api": "node test/api.test.js",
     "test:n3gen": "node test/n3gen.test.js",
     "test:examples": "node test/examples.test.js",
-    "test:check": "node test/check.test.js",
     "test:manifest": "node test/manifest.test.js",
     "test:playground": "node test/playground.test.js",
     "test:package": "node test/package.test.js",
-    "test:all": "npm run test:api && npm run test:n3gen && npm run test:examples && npm run test:check && npm run test:manifest && npm run test:playground",
+    "test:all": "npm run test:api && npm run test:n3gen && npm run test:examples && npm run && npm run test:manifest && npm run test:playground",
     "pretest": "npm run build && npm run test:packlist",
     "test": "npm run test:all",
     "posttest": "npm run test:package",

package/examples/check/input/ackermann.c

@@ -1,198 +0,0 @@
-#include <stdio.h>
-#include <stdlib.h>
-#include <stdint.h>
-#include <string.h>
-
-#define BIG_BASE 1000000000u
-
-typedef struct {
-    uint32_t *d;
-    size_t len;
-    size_t cap;
-} BigInt;
-
-static void die(const char *msg) {
-    fputs(msg, stderr);
-    fputc('\n', stderr);
-    exit(EXIT_FAILURE);
-}
-
-static void bi_init(BigInt *a, unsigned int value) {
-    a->cap = 4;
-    a->len = 1;
-    a->d = (uint32_t *)calloc(a->cap, sizeof(uint32_t));
-    if (!a->d) die("out of memory");
-    a->d[0] = value;
-}
-
-static void bi_reserve(BigInt *a, size_t need) {
-    if (need <= a->cap) return;
-    size_t cap = a->cap;
-    while (cap < need) cap *= 2;
-    uint32_t *p = (uint32_t *)realloc(a->d, cap * sizeof(uint32_t));
-    if (!p) die("out of memory");
-    memset(p + a->cap, 0, (cap - a->cap) * sizeof(uint32_t));
-    a->d = p;
-    a->cap = cap;
-}
-
-static void bi_trim(BigInt *a) {
-    while (a->len > 1 && a->d[a->len - 1] == 0) a->len--;
-}
-
-static void bi_mul_small(BigInt *a, unsigned int m) {
-    uint64_t carry = 0;
-    for (size_t i = 0; i < a->len; i++) {
-        uint64_t cur = (uint64_t)a->d[i] * m + carry;
-        a->d[i] = (uint32_t)(cur % BIG_BASE);
-        carry = cur / BIG_BASE;
-    }
-    while (carry) {
-        bi_reserve(a, a->len + 1);
-        a->d[a->len++] = (uint32_t)(carry % BIG_BASE);
-        carry /= BIG_BASE;
-    }
-}
-
-static void bi_add_small(BigInt *a, unsigned int add) {
-    uint64_t carry = add;
-    size_t i = 0;
-    while (carry) {
-        if (i == a->len) {
-            bi_reserve(a, a->len + 1);
-            a->d[a->len++] = 0;
-        }
-        uint64_t cur = (uint64_t)a->d[i] + carry;
-        a->d[i] = (uint32_t)(cur % BIG_BASE);
-        carry = cur / BIG_BASE;
-        i++;
-    }
-}
-
-static void bi_sub_small(BigInt *a, unsigned int sub) {
-    uint64_t borrow = sub;
-    size_t i = 0;
-    while (borrow) {
-        if (i >= a->len) die("underflow");
-        uint64_t cur = a->d[i];
-        uint64_t part = borrow % BIG_BASE;
-        borrow /= BIG_BASE;
-        if (cur < part) {
-            a->d[i] = (uint32_t)(cur + BIG_BASE - part);
-            borrow += 1;
-        } else {
-            a->d[i] = (uint32_t)(cur - part);
-        }
-        i++;
-    }
-    bi_trim(a);
-}
-
-static unsigned int bi_to_uint(const BigInt *a) {
-    if (a->len == 0) return 0;
-    if (a->len > 2) die("value too large for unsigned int");
-    uint64_t v = a->d[0];
-    if (a->len == 2) v += (uint64_t)a->d[1] * BIG_BASE;
-    if (v > 0xffffffffu) die("value too large for unsigned int");
-    return (unsigned int)v;
-}
-
-static char *bi_to_string(const BigInt *a) {
-    size_t bytes = a->len * 10 + 1;
-    char *s = (char *)malloc(bytes);
-    if (!s) die("out of memory");
-    int n = snprintf(s, bytes, "%u", a->d[a->len - 1]);
-    if (n < 0) die("snprintf failed");
-    size_t pos = (size_t)n;
-    for (size_t i = a->len - 1; i-- > 0;) {
-        n = snprintf(s + pos, bytes - pos, "%09u", a->d[i]);
-        if (n < 0) die("snprintf failed");
-        pos += (size_t)n;
-    }
-    return s;
-}
-
-static void bi_free(BigInt *a) {
-    free(a->d);
-    a->d = NULL;
-    a->len = 0;
-    a->cap = 0;
-}
-
-static BigInt bi_pow2(unsigned int exp) {
-    BigInt r;
-    bi_init(&r, 1);
-    for (unsigned int i = 0; i < exp; i++) bi_mul_small(&r, 2);
-    return r;
-}
-
-/* Independent evaluator for the fixed-z hyperoperation chain used by
-   examples/ackermann.n3:
-     A(x,y) = H(x, y+3, 2) - 3
-   where H(0,y,2)=y+1, H(1,y,2)=y+2, H(2,y,2)=2y, H(3,y,2)=2^y,
-   and H(x,0,2)=1 for x>3, H(x,y,2)=H(x-1, H(x,y-1,2), 2).
-*/
-static BigInt hyper2(unsigned int x, unsigned int y) {
-    if (x == 0) {
-        BigInt r;
-        bi_init(&r, y + 1);
-        return r;
-    }
-    if (x == 1) {
-        BigInt r;
-        bi_init(&r, y + 2);
-        return r;
-    }
-    if (x == 2) {
-        BigInt r;
-        bi_init(&r, 2u * y);
-        return r;
-    }
-    if (x == 3) {
-        return bi_pow2(y);
-    }
-    if (y == 0) {
-        BigInt r;
-        bi_init(&r, 1);
-        return r;
-    }
-
-    BigInt inner = hyper2(x, y - 1);
-    unsigned int inner_u = bi_to_uint(&inner);
-    bi_free(&inner);
-    return hyper2(x - 1, inner_u);
-}
-
-static BigInt ackermann2(unsigned int x, unsigned int y) {
-    BigInt r = hyper2(x, y + 3);
-    bi_sub_small(&r, 3);
-    return r;
-}
-
-static void print_line(unsigned int x, unsigned int y) {
-    BigInt v = ackermann2(x, y);
-    char *s = bi_to_string(&v);
-    printf("    (%u %u) :ackermann %s .\n", x, y, s);
-    free(s);
-    bi_free(&v);
-}
-
-int main(void) {
-    puts("@prefix : <https://eyereasoner.github.io/ns#> .");
-    puts("");
-    puts(":test :is {");
-    print_line(0, 0);
-    print_line(0, 6);
-    print_line(1, 2);
-    print_line(1, 7);
-    print_line(2, 2);
-    print_line(2, 9);
-    print_line(3, 4);
-    print_line(3, 1000);
-    print_line(4, 0);
-    print_line(4, 1);
-    print_line(4, 2);
-    print_line(5, 0);
-    puts("} .");
-    return 0;
-}

package/examples/check/input/deep-taxonomy-10.c

@@ -1,97 +0,0 @@
-#include <stdio.h>
-#include <stdlib.h>
-#include <string.h>
-
-enum { DEPTH = 10 };
-
-static int emit(const char *s) {
-    return fputs(s, stdout) != EOF;
-}
-
-static int emit_type(const char *kind, int index) {
-    return printf(":ind a :%s%d .\n", kind, index) >= 0;
-}
-
-int main(void) {
-    unsigned char hasN[DEPTH + 1];
-    unsigned char hasI[DEPTH + 1];
-    unsigned char hasJ[DEPTH + 1];
-    unsigned char hasA2 = 0;
-    unsigned char hasTest = 0;
-    int changed = 1;
-
-    memset(hasN, 0, sizeof(hasN));
-    memset(hasI, 0, sizeof(hasI));
-    memset(hasJ, 0, sizeof(hasJ));
-
-    /* fact: :ind a :N0. */
-    hasN[0] = 1;
-
-    /*
-     * Forward-chain the rules to a fixpoint.
-     *
-     * {?X a :N0}      => {?X a :N1,      :I1,      :J1}.
-     * {?X a :N1}      => {?X a :N2,      :I2,      :J2}.
-     * ...
-     * {?X a :N9}  => {?X a :N10, :I10, :J10}.
-     * {?X a :N10} => {?X a :A2}.
-     * {:ind a :A2}    => {:test :is true}.
-     */
-    while (changed) {
-        changed = 0;
-
-        for (int i = 0; i < DEPTH; i++) {
-            if (!hasN[i]) {
-                continue;
-            }
-
-            if (!hasN[i + 1]) {
-                hasN[i + 1] = 1;
-                changed = 1;
-            }
-            if (!hasI[i + 1]) {
-                hasI[i + 1] = 1;
-                changed = 1;
-            }
-            if (!hasJ[i + 1]) {
-                hasJ[i + 1] = 1;
-                changed = 1;
-            }
-        }
-
-        if (hasN[DEPTH] && !hasA2) {
-            hasA2 = 1;
-            changed = 1;
-        }
-
-        if (hasA2 && !hasTest) {
-            hasTest = 1;
-            changed = 1;
-        }
-    }
-
-    if (!emit("@prefix : <http://eulersharp.sourceforge.net/2009/12dtb/test#> .\n\n")) {
-        return EXIT_FAILURE;
-    }
-
-    for (int i = 1; i <= DEPTH; i++) {
-        if (!hasN[i] || !hasI[i] || !hasJ[i]) {
-            fprintf(stderr, "incomplete closure at depth %d\n", i);
-            return EXIT_FAILURE;
-        }
-        if (!emit_type("N", i) || !emit_type("I", i) || !emit_type("J", i)) {
-            return EXIT_FAILURE;
-        }
-    }
-
-    if (!hasA2 || !hasTest) {
-        fputs("expected conclusions were not derived\n", stderr);
-        return EXIT_FAILURE;
-    }
-
-    if (!emit(":ind a :A2 .\n:test :is true .\n")) {
-        return EXIT_FAILURE;
-    }
-
-    return EXIT_SUCCESS;
-}