eyelang 1.7.11 → 1.7.12

package/README.md

@@ -3,7 +3,7 @@
 [![npm version](https://img.shields.io/npm/v/eyelang.svg)](https://www.npmjs.com/package/eyelang)
 [![DOI](https://img.shields.io/badge/DOI-10.5281%2Fzenodo.20761726-blue.svg)](https://doi.org/10.5281/zenodo.20761726)
 
-Eyelang is a small logic programming language for rules, goals, answers, and proofs.
+Eyelang is a small logic programming language for facts, rules, goals, answers, and proofs.
 Its source syntax is Prolog-like Horn-clause syntax with deliberate eyelang choices, including `?x` variables for N3/SPARQL-style readability, explicit `table(path, 2).` declarations for tabled predicates, advisory `mode/3` declarations for host tooling, explicit Herbrand witness terms for executable existential-style consequences, and stratified-negation diagnostics for portable `not/1` usage.
 It grew out of logic-language experiments in the EYE/N3 reasoning tradition, but is packaged here as its own project.
 

package/docs/guide.md

@@ -1,6 +1,6 @@
 # Eyelang Guide
 
-This guide introduces Eyelang, a small Horn-clause language and engine whose source syntax is Prolog-like but deliberately its own compact language for rules, goals, answers, and proofs. Eyelang works over ordinary terms, lists, arithmetic, strings, and finite search. Run it with the `eyelang` CLI, or use `node bin/eyelang.js` when working directly from a source checkout.
+This guide introduces Eyelang, a small Horn-clause language and engine whose source syntax is Prolog-like but deliberately its own compact language for facts, rules, goals, answers, and proofs. Eyelang works over ordinary terms, lists, arithmetic, strings, and finite search. Run it with the `eyelang` CLI, or use `node bin/eyelang.js` when working directly from a source checkout.
 
 Programs write relations directly, for example `ancestor(pat, emma)` or `status(case1, accepted)`. Absolute IRI atoms can be written explicitly with angle brackets, for example `<https://schema.org/name>`, when a program needs web identifiers without prefix declarations. Eyelang output is ordinary Eyelang syntax: by default, the CLI materializes selected answer facts and prints those facts only. Pass `--proof` (or `-p`) when you also want each answer followed by a `why/2` explanation fact that records the proof. Programs may add `materialize/2` declarations such as `materialize(answer, 2).` to focus output on selected predicates.
 

package/examples/stirling-bell-numbers.eye

@@ -1,30 +1,67 @@
 % Stirling numbers of the second kind and Bell numbers.
 %
-% stirling2(N, K, Count) counts partitions of N labelled elements into K nonempty blocks.
-% The recurrence either places the newest labelled element into a new singleton block or into
-% one of the K existing blocks.  bell/2 then sums a whole row of Stirling numbers.
-% Memoization turns the overlapping recurrence into a small dynamic-programming table.
+% The Stirling count S(N,K) is computed with the inclusion-exclusion formula
+%   S(N,K) = (1/K!) * sum(I=0..K, (-1)^(K-I) * C(K,I) * I^N)
+% instead of the overlapping two-branch recurrence.  Bell numbers use
+%   B(0) = 1, B(N) = sum(K=0..N-1, C(N-1,K) * B(K)).
+% The table declarations memoize the smaller helper relations used by both formulas.
 materialize(stirling_bell_answer, 2).
 
+table(factorial, 2).
+table(binomial, 3).
 table(stirling2, 3).
+table(bell, 2).
+
+factorial(0, 1).
+factorial(?n, ?value) :-
+  gt(?n, 0),
+  sub(?n, 1, ?n1),
+  factorial(?n1, ?previous),
+  mul(?n, ?previous, ?value).
+
+binomial(0, 0, 1).
+binomial(?n, 0, 1) :- gt(?n, 0).
+binomial(?n, ?n, 1) :- gt(?n, 0).
+binomial(?n, ?k, ?value) :-
+  gt(?n, 0),
+  gt(?k, 0),
+  lt(?k, ?n),
+  sub(?n, 1, ?n1),
+  sub(?k, 1, ?k1),
+  binomial(?n1, ?k1, ?left),
+  binomial(?n1, ?k, ?right),
+  add(?left, ?right, ?value).
+
+signed_term(?n, ?k, ?i, ?term) :-
+  binomial(?k, ?i, ?c),
+  pow(?i, ?n, ?p),
+  mul(?c, ?p, ?unsigned),
+  sub(?k, ?i, ?d),
+  mod(?d, 2, 0),
+  eq(?term, ?unsigned).
+signed_term(?n, ?k, ?i, ?term) :-
+  binomial(?k, ?i, ?c),
+  pow(?i, ?n, ?p),
+  mul(?c, ?p, ?unsigned),
+  sub(?k, ?i, ?d),
+  mod(?d, 2, 1),
+  neg(?unsigned, ?term).
 
-% Boundary cases define the empty partition and the impossible zero-block columns.
 stirling2(0, 0, 1).
 stirling2(?n, 0, 0) :- gt(?n, 0).
 stirling2(0, ?k, 0) :- gt(?k, 0).
 stirling2(?n, ?k, ?count) :-
   gt(?n, 0),
   gt(?k, 0),
-  sub(?n, 1, ?n1),
-  sub(?k, 1, ?k1),
-  stirling2(?n1, ?k1, ?newblock),
-  stirling2(?n1, ?k, ?existingblocks),
-  mul(?k, ?existingblocks, ?extended),
-  add(?newblock, ?extended, ?count).
+  sumall(?term, (between(0, ?k, ?i), signed_term(?n, ?k, ?i, ?term)), ?sum),
+  factorial(?k, ?factorial),
+  div(?sum, ?factorial, ?count).
 
-% Bell(N) is the row sum S(N,0)+...+S(N,N).
+bell(0, 1).
 bell(?n, ?count) :-
-  sumall(?s, (between(0, ?n, ?k), stirling2(?n, ?k, ?s)), ?count).
+  gt(?n, 0),
+  sub(?n, 1, ?n1),
+  sumall(?term, (between(0, ?n1, ?k), binomial(?n1, ?k, ?choose), bell(?k, ?bell), mul(?choose, ?bell, ?term)), ?count).
 
 stirling_bell_answer(stirling_10_4, ?count) :- stirling2(10, 4, ?count).
 stirling_bell_answer(stirling_12_5, ?count) :- stirling2(12, 5, ?count).

package/package.json

@@ -1,7 +1,7 @@
 {
   "name": "eyelang",
-  "version": "1.7.11",
-  "description": "A small Prolog-like logic programming language for rules, goals, answers, and proofs.",
+  "version": "1.7.12",
+  "description": "A small Prolog-like logic programming language for facts, rules, goals, answers, and proofs.",
   "type": "module",
   "main": "./index.js",
   "types": "./index.d.ts",

package/test/run-regression.mjs

@@ -800,7 +800,7 @@ function whiteBoxCases() {
           ['matrix-chain-order.eye', 'cost', 3, false],
           ['modular-exponentiation.eye', 'pow_mod', 4, true],
           ['pell-equation.eye', 'pell', 3, true],
-          ['stirling-bell-numbers.eye', 'stirling2', 3, true],
+          ['stirling-bell-numbers.eye', 'stirling2', 3, false],
           ['totient-summatory.eye', 'gcd', 3, true],
           ['totient-summatory.eye', 'totient', 2, false],
           ['weighted-interval-scheduling.eye', 'best_from', 2, true],