@metago-ai/engine 2.1.2 → 2.1.4

package/ENGINES/yili/lib/hetu-luoshu.js

@@ -0,0 +1,123 @@
+/**
+ * 易理工程学 · 河图洛书引擎（P2）
+ * 幻方构造（奇数阶/双偶数阶/单偶数阶——算学宝鉴纵横图法）+ 九宫配卦配数
+ */
+
+/** 洛书九宫数（传统排布：戴九履一，左三右七，二四为肩，六八为足，五居中央） */
+export const LUOSHU_GRID = [
+  [4, 9, 2],
+  [3, 5, 7],
+  [8, 1, 6]
+];
+
+/** 洛书九宫配后天八卦（坎一/坤二/震三/巽四/中五/乾六/兑七/艮八/离九） */
+export const LUOSHU_GUA = {
+  1: '坎', 2: '坤', 3: '震', 4: '巽', 5: '中', 6: '乾', 7: '兑', 8: '艮', 9: '离'
+};
+
+/** 河图生成数：一六共宗水、二七同道火、三八为朋木、四九为友金、五十同途土 */
+export const HETU_PAIRS = {
+  '水': [1, 6], '火': [2, 7], '木': [3, 8], '金': [4, 9], '土': [5, 10]
+};
+
+/** 奇数阶幻方（Siamese 法，n 为奇数） */
+export function magicSquareOdd(n) {
+  if (n % 2 === 0) throw new Error('magicSquareOdd 需要奇数阶');
+  const sq = Array.from({ length: n }, () => new Array(n).fill(0));
+  let r = 0, c = Math.floor(n / 2);
+  for (let k = 1; k <= n * n; k++) {
+    sq[r][c] = k;
+    const nr = (r - 1 + n) % n, nc = (c + 1) % n;
+    if (sq[nr][nc] !== 0) { r = (r + 1) % n; }
+    else { r = nr; c = nc; }
+  }
+  return sq;
+}
+
+/** 双偶数阶幻方（n=4k，对称交换法） */
+export function magicSquareDoubleEven(n) {
+  if (n % 4 !== 0) throw new Error('magicSquareDoubleEven 需要 4 的倍数阶');
+  const sq = Array.from({ length: n }, (_, i) =>
+    Array.from({ length: n }, (_, j) => i * n + j + 1));
+  const k = n / 4;
+  for (let i = 0; i < n; i++) {
+    for (let j = 0; j < n; j++) {
+      const inA = (i % 4 < 2) === (j % 4 < 2) ? false : (i % 4 >= 2) === (j % 4 >= 2);
+      // 对角区（4×4 子块的对角线）需要交换
+      const ii = Math.floor(i / 4), jj = Math.floor(j / 4);
+      const diag = (i % 4 === j % 4) || (i % 4 + j % 4 === 3);
+      const corner = (ii + jj) % 2 === 0;
+      if (diag === corner) {
+        sq[i][j] = n * n + 1 - sq[i][j];
+      }
+    }
+  }
+  return sq;
+}
+
+/**
+ * 单偶数阶幻方（n=4k+2）。
+ * 当前实现：康威 LUX 变体，已验证 n=6（k=1）行/列/对角全等（算学宝鉴纵横图法同族）。
+ * 交换规则（0/2/3/1 布局 + A↔C 左 k 列除中间行 + A↔C 中间行第 k+1 列）。
+ */
+export function magicSquareSingleEven(n) {
+  if (n % 4 !== 2) throw new Error('magicSquareSingleEven 需要 4k+2 阶');
+  if (n !== 6) throw new Error('单偶数阶幻方当前支持 k=1（六阶），已验证；更高阶请扩展 LUX 表');
+  const m = n / 2;
+  const k = (n - 2) / 4;
+  const mid = Math.floor(m / 2);
+  const A = magicSquareOdd(m);
+  const sq = Array.from({ length: n }, () => new Array(n).fill(0));
+  for (let i = 0; i < m; i++) {
+    for (let j = 0; j < m; j++) {
+      const v = A[i][j];
+      sq[i][j] = v;                    // 左上 A
+      sq[i][j + m] = v + 2 * m * m;    // 右上 B
+      sq[i + m][j] = v + 3 * m * m;    // 左下 C
+      sq[i + m][j + m] = v + m * m;    // 右下 D
+    }
+  }
+  // 交换 A ↔ C：左 k 列（除中间行）
+  for (let i = 0; i < m; i++) {
+    if (i === mid) continue;
+    for (let j = 0; j < k; j++) {
+      const t = sq[i][j]; sq[i][j] = sq[i + m][j]; sq[i + m][j] = t;
+    }
+  }
+  // 交换 A ↔ C：中间行第 k+1 列
+  const t = sq[mid][k]; sq[mid][k] = sq[mid + m][k]; sq[mid + m][k] = t;
+  return sq;
+}
+
+/** 通用幻方构造（算学宝鉴纵横图法：奇/双偶/单偶） */
+export function magicSquare(n) {
+  if (n % 2 === 1) return magicSquareOdd(n);
+  if (n % 4 === 0) return magicSquareDoubleEven(n);
+  return magicSquareSingleEven(n);
+}
+
+/** 幻方性质校验：行/列/对角线和相等 */
+export function verifyMagicSquare(sq) {
+  const n = sq.length;
+  const target = (n * (n * n + 1)) / 2;
+  const sums = [];
+  for (let i = 0; i < n; i++) sums.push(sq[i].reduce((a, b) => a + b, 0));
+  for (let j = 0; j < n; j++) { let s = 0; for (let i = 0; i < n; i++) s += sq[i][j]; sums.push(s); }
+  let d1 = 0, d2 = 0;
+  for (let i = 0; i < n; i++) { d1 += sq[i][i]; d2 += sq[i][n - 1 - i]; }
+  sums.push(d1, d2);
+  return { target, ok: sums.every(s => s === target), sums };
+}
+
+/** 九宫配卦：3×3 洛书 → 八卦名 */
+export function luoshuWithGua() {
+  return LUOSHU_GRID.map(row => row.map(v => ({ number: v, gua: LUOSHU_GUA[v] })));
+}
+
+/** 九宫方位：给定数字 → 方位（坎北/离南/震东/兑西/巽东南/乾西北/坤西南/艮东北/中） */
+export function luoshuDirection(number) {
+  const dir = {
+    1: '北', 2: '西南', 3: '东', 4: '东南', 5: '中', 6: '西北', 7: '西', 8: '东北', 9: '南'
+  };
+  return dir[number] || null;
+}

package/ENGINES/yili/lib/index.js

@@ -0,0 +1,23 @@
+/**
+ * 易理工程学 · 引擎汇总入口
+ */
+export * from './almanac.js';
+export * from './wuxing.js';
+export * from './hetu-luoshu.js';
+export * from './yijing.js';
+export * from './liuyao.js';
+export * from './bazi.js';
+export * from './ziwei.js';
+export * from './qimen.js';
+export * from './taiyi.js';
+export * from './zeri.js';
+export * from './shushu-calculus.js';
+export * from './diagnose.js';
+
+export const ENGINE = {
+  name: '易理工程学',
+  version: '0.2.0',
+  families: ['历法', '周易', '五行八卦', '河图洛书', '六爻', '四柱', '紫微', '奇门', '太乙', '择日', '文言微积分', '象数诊断'],
+  almanac: '寿星天文历公版公式（Sxing）',
+  confidence: '一切断言均为启发式标签 + 置信度 + 规则编号，可溯源可复算'
+};

package/ENGINES/yili/lib/liuyao.js

@@ -0,0 +1,226 @@
+/**
+ * 易理工程学 · 六爻引擎（P3）
+ * 装卦（纳甲）/ 安世应 / 定六亲 / 六神 / 月建日辰旺衰 / 应期推算
+ */
+import { ZHI, GAN, GANZHI, changsheng, nayin } from './almanac.js';
+import { relation } from './wuxing.js';
+
+/** 八纯卦纳甲表：内卦三爻地支（初→三）、外卦三爻地支（四→上） */
+export const NAJIA = {
+  '乾': { inner: ['子', '寅', '辰'], outer: ['午', '申', '戌'] },
+  '坎': { inner: ['寅', '辰', '午'], outer: ['申', '戌', '子'] },
+  '艮': { inner: ['辰', '午', '申'], outer: ['戌', '子', '寅'] },
+  '震': { inner: ['子', '寅', '辰'], outer: ['午', '申', '戌'] },
+  '巽': { inner: ['丑', '亥', '酉'], outer: ['未', '巳', '卯'] },
+  '离': { inner: ['卯', '丑', '亥'], outer: ['酉', '未', '巳'] },
+  '坤': { inner: ['未', '巳', '卯'], outer: ['丑', '亥', '酉'] },
+  '兑': { inner: ['巳', '卯', '丑'], outer: ['亥', '酉', '未'] }
+};
+
+/** 纳甲天干：乾内甲外壬、坎内戊外戊、艮内丙外丙、震内庚外庚、巽内辛外辛、离内己外己、坤内乙外癸、兑内丁外丁 */
+export const NAJIA_GAN = {
+  '乾': { inner: '甲', outer: '壬' }, '坎': { inner: '戊', outer: '戊' }, '艮': { inner: '丙', outer: '丙' },
+  '震': { inner: '庚', outer: '庚' }, '巽': { inner: '辛', outer: '辛' }, '离': { inner: '己', outer: '己' },
+  '坤': { inner: '乙', outer: '癸' }, '兑': { inner: '丁', outer: '丁' }
+};
+
+/** 八宫世应位置（六十四卦在世爻 1-6 变：宫首世六，一世世一……游魂世四、归魂世三） */
+export const SHI_YAO = {
+  '乾为天': 6, '天风姤': 1, '天山遁': 2, '天地否': 3, '风地观': 4, '山地剥': 5, '火地晋': 4, '火天大有': 3,
+  '兑为泽': 6, '泽水困': 1, '泽地萃': 2, '泽山咸': 3, '水山蹇': 4, '地山谦': 5, '雷山小过': 4, '雷泽归妹': 3,
+  '离为火': 6, '火山旅': 1, '火风鼎': 2, '火水未济': 3, '山水蒙': 4, '风水涣': 5, '天水讼': 4, '天火同人': 3,
+  '震为雷': 6, '雷地豫': 1, '雷水解': 2, '雷风恒': 3, '地风升': 4, '水风井': 5, '泽风大过': 4, '泽雷随': 3,
+  '巽为风': 6, '风天小畜': 1, '风火家人': 2, '风雷益': 3, '天雷无妄': 4, '火雷噬嗑': 5, '山雷颐': 4, '山风蛊': 3,
+  '坎为水': 6, '水泽节': 1, '水雷屯': 2, '水火既济': 3, '泽火革': 4, '雷火丰': 5, '地火明夷': 4, '地水师': 3,
+  '艮为山': 6, '山火贲': 1, '山天大畜': 2, '山泽损': 3, '火泽睽': 4, '天泽履': 5, '风泽中孚': 4, '风山渐': 3,
+  '坤为地': 6, '地雷复': 1, '地泽临': 2, '地天泰': 3, '雷天大壮': 4, '泽天夬': 5, '水天需': 4, '水地比': 3
+};
+
+/** 宫五行 */
+const GONG_WX = { '乾': '金', '兑': '金', '离': '火', '震': '木', '巽': '木', '坎': '水', '艮': '土', '坤': '土' };
+
+/** 六亲：宫五行 × 爻支五行 → 父母/兄弟/子孙/妻财/官鬼 */
+export function liuqin(gongWuxing, zhiWuxing) {
+  const r = relation(gongWuxing, zhiWuxing);
+  switch (r) {
+    case '比和': return '兄弟';
+    case '生': return '子孙';   // 宫生爻 = 子孙（我所生）
+    case '被生': return '父母'; // 爻生宫 = 父母（生我者）
+    case '克': return '妻财';   // 宫克爻 = 妻财（我克者）
+    case '被克': return '官鬼'; // 爻克宫 = 官鬼（克我者）
+    default: return null;
+  }
+}
+
+/** 六神（按日干起青龙：甲乙起青龙，丙丁起朱雀，戊起勾陈，己起腾蛇，庚辛起白虎，壬癸起玄武） */
+export const LIUSHEN_ORDER = ['青龙', '朱雀', '勾陈', '腾蛇', '白虎', '玄武'];
+export function liuShenStart(dayGan) {
+  const G = { '甲': 0, '乙': 0, '丙': 1, '丁': 1, '戊': 2, '己': 3, '庚': 4, '辛': 4, '壬': 5, '癸': 5 };
+  return G[dayGan] ?? 0;
+}
+export function liuShen(dayGan, yaoIdx) {
+  const start = liuShenStart(dayGan);
+  return LIUSHEN_ORDER[(start + yaoIdx) % 6];
+}
+
+/**
+ * 六爻装卦全流程
+ * @param {string} fullName 卦全名（如 '水雷屯'）
+ * @param {string} dayGanZhi 日干支（如 '甲子'）
+ * @param {string} monthZhi 月建（如 '寅'）
+ * @param {string} dayGan 日干（'甲'）
+ * @returns 六爻明细
+ */
+export function zhuanggua(fullName, dayGanZhi, monthZhi, dayGan) {
+  // 卦宫
+  let gong = null;
+  for (const [g, names] of Object.entries(SHI_GONG_MAP)) {
+    if (names.includes(fullName)) { gong = g; break; }
+  }
+  const gongWx = GONG_WX[gong];
+  const shi = SHI_YAO[fullName];
+  const ying = shi === 1 ? 4 : shi === 2 ? 5 : shi === 3 ? 6 : shi === 4 ? 1 : shi === 5 ? 2 : shi === 6 ? 3 : null;
+  // 上下卦
+  const { upper, lower } = upperLower(fullName);
+  const rows = [];
+  const innerZhis = NAJIA[lower].inner;
+  const outerZhis = NAJIA[upper].outer;
+  const innerGan = NAJIA_GAN[lower].inner;
+  const outerGan = NAJIA_GAN[upper].outer;
+  const zhis = [...innerZhis, ...outerZhis];
+  const gans = [innerGan, innerGan, innerGan, outerGan, outerGan, outerGan];
+  for (let i = 0; i < 6; i++) {
+    const yaoIdx = i; // 0=初
+    const zhi = zhis[i];
+    const gan = gans[i];
+    const zhiWx = ZHI_WUXING_MAP[zhi];
+    rows.push({
+      yao: i + 1,
+      position: ['初', '二', '三', '四', '五', '上'][i],
+      ganZhi: gan + zhi,
+      zhi,
+      wuxing: zhiWx,
+      liuqin: liuqin(gongWx, zhiWx),
+      liuShen: liuShen(dayGan, i),
+      shi: i + 1 === shi,
+      ying: i + 1 === ying,
+      // 旺衰（月建）：月支五行对爻支五行
+      yueJian: monthRelation(monthZhi, zhi),
+      // 日辰生克
+      riChen: dayZhiRelation(dayGanZhi, zhi),
+      // 旬空
+      xunKong: isXunKong(dayGanZhi, zhi)
+    });
+  }
+  return { gong, gongWuxing: gongWx, shi, ying, rows, fullName };
+}
+
+/** 八宫→卦名表（与 yijing.js GUA_FULL 对应） */
+const SHI_GONG_MAP = {
+  '乾': ['乾为天', '天风姤', '天山遁', '天地否', '风地观', '山地剥', '火地晋', '火天大有'],
+  '兑': ['兑为泽', '泽水困', '泽地萃', '泽山咸', '水山蹇', '地山谦', '雷山小过', '雷泽归妹'],
+  '离': ['离为火', '火山旅', '火风鼎', '火水未济', '山水蒙', '风水涣', '天水讼', '天火同人'],
+  '震': ['震为雷', '雷地豫', '雷水解', '雷风恒', '地风升', '水风井', '泽风大过', '泽雷随'],
+  '巽': ['巽为风', '风天小畜', '风火家人', '风雷益', '天雷无妄', '火雷噬嗑', '山雷颐', '山风蛊'],
+  '坎': ['坎为水', '水泽节', '水雷屯', '水火既济', '泽火革', '雷火丰', '地火明夷', '地水师'],
+  '艮': ['艮为山', '山火贲', '山天大畜', '山泽损', '火泽睽', '天泽履', '风泽中孚', '风山渐'],
+  '坤': ['坤为地', '地雷复', '地泽临', '地天泰', '雷天大壮', '泽天夬', '水天需', '水地比']
+};
+
+const ZHI_WUXING_MAP = { '子': '水', '丑': '土', '寅': '木', '卯': '木', '辰': '土', '巳': '火', '午': '火', '未': '土', '申': '金', '酉': '金', '戌': '土', '亥': '水' };
+
+/** 自然字 → 经卦名 */
+const NATURE_TO_GUA = { '天': '乾', '泽': '兑', '火': '离', '雷': '震', '风': '巽', '水': '坎', '山': '艮', '地': '坤' };
+
+/** 卦全名 → 上下经卦名（全名 = 上卦字 + 下卦字 + 卦名；纯卦 'X为X' 特例） */
+function upperLower(fullName) {
+  if (fullName.includes('为')) return { upper: fullName[0], lower: fullName[0] };
+  return { upper: NATURE_TO_GUA[fullName[0]], lower: NATURE_TO_GUA[fullName[1]] };
+}
+
+/** 月建对爻支：月建五行与爻支五行关系（旺相休囚死简化：旺/相/休/囚/死） */
+export function monthRelation(monthZhi, zhi) {
+  const m = ZHI_WUXING_MAP[monthZhi];
+  const z = ZHI_WUXING_MAP[zhi];
+  const r = relation(m, z);
+  if (r === '比和') return '旺';
+  if (r === '生') return '相';    // 月生爻
+  if (r === '被生') return '休';  // 爻生月
+  if (r === '克') return '死';    // 月克爻
+  return '囚';                    // 爻克月
+}
+
+/** 日辰与爻支关系（生/克/比） */
+export function dayZhiRelation(dayGanZhi, zhi) {
+  const dz = dayGanZhi.slice(1);
+  if (dz === zhi) return '比';
+  const d = ZHI_WUXING_MAP[dz];
+  const z = ZHI_WUXING_MAP[zhi];
+  const r = relation(d, z);
+  if (r === '生') return '生';
+  if (r === '克') return '克';
+  if (r === '被生') return '受生';
+  if (r === '被克') return '受克';
+  return '比';
+}
+
+/** 旬空判定：日柱所在旬的空亡支 */
+export function isXunKong(dayGanZhi, zhi) {
+  const idx = GANZHI.indexOf(dayGanZhi);
+  const head = idx - (idx % 10);
+  const kong1 = ZHI[(head + 10) % 12];
+  const kong2 = ZHI[(head + 11) % 12];
+  return zhi === kong1 || zhi === kong2;
+}
+
+/** 用神库：问事 → 六亲 */
+export const YONGSHEN_LIB = {
+  '财': '妻财', '生意': '妻财', '交易': '妻财', '求财': '妻财', '投资': '妻财', '债务': '妻财',
+  '官': '官鬼', '事业': '官鬼', '工作': '官鬼', '升迁': '官鬼', '官司': '官鬼', '名誉': '官鬼', '丈夫': '官鬼',
+  '文书': '父母', '考试': '父母', '学业': '父母', '长辈': '父母', '房宅': '父母', '车': '父母', '印章': '父母',
+  '子孙': '子孙', '子女': '子孙', '医药': '子孙', '下属': '子孙', '晚辈': '子孙', '平安': '子孙',
+  '兄弟': '兄弟', '朋友': '兄弟', '同辈': '兄弟', '合伙': '兄弟', '竞争': '兄弟'
+};
+export function yongshenFor(topic) {
+  for (const [k, v] of Object.entries(YONGSHEN_LIB)) {
+    if (topic.includes(k)) return v;
+  }
+  return null;
+}
+
+/** 应期推算：动爻/用神 + 冲合空破规则 → 候选应期 */
+export function yingqiHints(rows, targetLiqin) {
+  const hints = [];
+  const target = rows.find(r => r.liuqin === targetLiqin);
+  if (!target) return hints;
+  const z = target.zhi;
+  const chong = { '子': '午', '午': '子', '丑': '未', '未': '丑', '寅': '申', '申': '寅', '卯': '酉', '酉': '卯', '辰': '戌', '戌': '辰', '巳': '亥', '亥': '巳' };
+  const he = { '子': '丑', '丑': '子', '寅': '亥', '亥': '寅', '卯': '戌', '戌': '卯', '辰': '酉', '酉': '辰', '巳': '申', '申': '巳', '午': '未', '未': '午' };
+  hints.push({ rule: '用神值日', desc: `${targetLiqin}爻${z}日应`, zhi: z });
+  hints.push({ rule: '冲日应', desc: `冲${z}之${chong[z]}日（逢冲则动）`, zhi: chong[z] });
+  hints.push({ rule: '合日应', desc: `${z}与${he[z]}合日`, zhi: he[z] });
+  if (target.xunKong) {
+    hints.push({ rule: '出空', desc: `${z}旬空，出空之日应`, zhi: z });
+    hints.push({ rule: '冲空', desc: `冲${z}之${chong[z]}日（冲空则实）`, zhi: chong[z] });
+  }
+  return hints;
+}
+
+/** 完整六爻决策分析 */
+export function liuyaoAnalysis({ fullName, dayGanZhi, monthZhi, dayGan, topic, hour = 12 }) {
+  const gua = zhuanggua(fullName, dayGanZhi, monthZhi, dayGan);
+  const yongshen = topic ? yongshenFor(topic) : null;
+  const yongshenRow = yongshen ? gua.rows.find(r => r.liuqin === yongshen) : null;
+  const shiRow = gua.rows.find(r => r.shi);
+  const yingRow = gua.rows.find(r => r.ying);
+  return {
+    ...gua,
+    yongshen,
+    yongshenRow: yongshenRow ? { position: yongshenRow.position, zhi: yongshenRow.zhi, wuxing: yongshenRow.wuxing, yueJian: yongshenRow.yueJian, riChen: yongshenRow.riChen, xunKong: yongshenRow.xunKong } : null,
+    shiVsYing: shiRow && yingRow ? {
+      shiWuxing: shiRow.wuxing, yingWuxing: yingRow.wuxing,
+      relation: relation(shiRow.wuxing, yingRow.wuxing)
+    } : null,
+    yingqi: yongshen ? yingqiHints(gua.rows, yongshen) : []
+  };
+}

package/ENGINES/yili/lib/qimen.js

@@ -0,0 +1,167 @@
+/**
+ * 易理工程学 · 奇门遁甲引擎（P3 · 工程最重）
+ * 时家奇门：节气定局（阴阳遁 18 局）→ 三元 → 三奇六仪排布 → 值符值使 → 九星/八门/八神落宫
+ * 状态空间：1080 局全枚举可生成
+ */
+import { GANZHI, ganzhiDayOfDate, jieqiOfYear } from './almanac.js';
+
+/** 九宫（洛书序：坎一坤二震三巽四中五乾六兑七艮八离九） */
+export const NINE_PALACES = ['坎', '坤', '震', '巽', '中', '乾', '兑', '艮', '离'];
+export const PALACE_NUMBERS = [1, 2, 3, 4, 5, 6, 7, 8, 9];
+
+/** 节气 → 上中下元局数口诀（阳遁 9 局/阴遁 9 局） */
+export const DUN_RULES = {
+  // 阳遁（冬至→夏至前）
+  yang: {
+    '冬至': [1, 7, 4], '小寒': [2, 8, 5], '大寒': [3, 9, 6],
+    '立春': [8, 5, 2], '雨水': [9, 6, 3], '惊蛰': [1, 7, 4],
+    '春分': [3, 9, 6], '清明': [4, 1, 7], '谷雨': [5, 2, 8],
+    '立夏': [4, 1, 7], '小满': [5, 2, 8], '芒种': [6, 3, 9]
+  },
+  // 阴遁（夏至→冬至前）
+  yin: {
+    '夏至': [9, 3, 6], '小暑': [8, 2, 5], '大暑': [7, 1, 4],
+    '立秋': [2, 5, 8], '处暑': [1, 4, 7], '白露': [9, 3, 6],
+    '秋分': [7, 1, 4], '寒露': [6, 9, 3], '霜降': [5, 8, 2],
+    '立冬': [6, 9, 3], '小雪': [5, 8, 2], '大雪': [4, 7, 1]
+  }
+};
+
+/** 三元定元：日柱序 mod 15 = 0 上元、5 下元、10 中元 */
+export function yuanOfDay(dayGanZhi) {
+  const idx = GANZHI.indexOf(dayGanZhi);
+  const r = idx % 15;
+  if (r === 0) return '上元';
+  if (r === 5) return '下元';
+  if (r === 10) return '中元';
+  // 非符头日：归入前一个符头之元
+  const back = (idx - r) % 60; // 最近符头
+  const rb = back % 15;
+  if (rb === 0) return '上元';
+  if (rb === 5) return '下元';
+  return '中元';
+}
+
+/** 定局：日期 → { dunType: 'yang'|'yin', jieqi, yuan, ju } */
+export function dingJu(y, m, d) {
+  const jq = jieqiOfYear(y);
+  const prevJq = jieqiOfYear(y - 1);
+  const ts = Date.UTC(y, m - 1, d);
+  // 最近的冬至/夏至（跨年）：日期在 [冬至, 夏至) → 阳遁；[夏至, 冬至) → 阴遁
+  const dzThis = jq.find(j => j.name === '冬至');
+  const xzThis = jq.find(j => j.name === '夏至');
+  const dzPrev = prevJq.find(j => j.name === '冬至');
+  const candidates = [
+    { name: '冬至', ts: Date.parse(dzPrev.solar.replace(' ', 'T') + '+08:00'), type: 'yang' },
+    { name: '夏至', ts: Date.parse(xzThis.solar.replace(' ', 'T') + '+08:00'), type: 'yin' },
+    { name: '冬至', ts: Date.parse(dzThis.solar.replace(' ', 'T') + '+08:00'), type: 'yang' }
+  ];
+  let latest = candidates[0];
+  for (const c of candidates) if (c.ts <= ts && c.ts >= latest.ts) latest = c;
+  const dunType = latest.type;
+  // 当前节（立春惊蛰清明立夏芒种小暑立秋白露寒露立冬大雪小寒 + 冬至/夏至）
+  const JIE_NAMES = ['立春', '惊蛰', '清明', '立夏', '芒种', '小暑', '立秋', '白露', '寒露', '立冬', '大雪', '小寒'];
+  let jieqiName = '冬至';
+  for (const j of jq) {
+    const jt = Date.parse(j.solar.replace(' ', 'T') + '+08:00');
+    if (jt <= ts && (JIE_NAMES.includes(j.name) || j.name === '冬至' || j.name === '夏至')) jieqiName = j.name;
+  }
+  const yuan = yuanOfDay(ganzhiDayOfDate(y, m, d));
+  const rules = DUN_RULES[dunType];
+  const ju = rules[jieqiName] ? rules[jieqiName][{ '上元': 0, '中元': 1, '下元': 2 }[yuan]] : 1;
+  return { dunType, jieqi: jieqiName, yuan, ju };
+}
+
+/** 三奇六仪排布：六仪戊己庚辛壬癸 + 三奇丁丙乙（乙丙丁）。阳遁顺布、阴遁逆布，从局数宫起 */
+export function sanqiLiuyi(dunType, ju) {
+  const seq = ['戊', '己', '庚', '辛', '壬', '癸', '丁', '丙', '乙'];
+  const start = ju; // 局数宫号（1-9）
+  const out = {};
+  for (let i = 0; i < 9; i++) {
+    let palace = start - 1;
+    if (dunType === 'yang') palace = (palace + i) % 9;
+    else palace = ((palace - i) % 9 + 9) % 9;
+    out[NINE_PALACES[palace]] = seq[i];
+  }
+  return out;
+}
+
+/** 九星（中宫天禽寄坤二宫） */
+export const JIU_XING = { '坎': '天蓬', '坤': '天芮', '震': '天冲', '巽': '天辅', '中': '天禽', '乾': '天心', '兑': '天柱', '艮': '天任', '离': '天英' };
+/** 八门 */
+export const BA_MEN = { '坎': '休门', '坤': '死门', '震': '伤门', '巽': '杜门', '中': '中门', '乾': '开门', '兑': '惊门', '艮': '生门', '离': '景门' };
+/** 八神（阳遁顺排、阴遁逆排，值符起） */
+export const BA_SHEN = ['值符', '腾蛇', '太阴', '六合', '白虎', '玄武', '九地', '九天'];
+
+/** 旬首：时柱 → 六甲旬首 */
+export function xunShou(hourGanZhi) {
+  const idx = GANZHI.indexOf(hourGanZhi);
+  const head = idx - (idx % 10);
+  const map = { 0: '甲子', 10: '甲戌', 20: '甲申', 30: '甲午', 40: '甲辰', 50: '甲寅' };
+  return map[head];
+}
+
+/** 六甲旬首对应六仪：甲子戊、甲戌己、甲申庚、甲午辛、甲辰壬、甲寅癸 */
+const XUNSHOU_YI = { '甲子': '戊', '甲戌': '己', '甲申': '庚', '甲午': '辛', '甲辰': '壬', '甲寅': '癸' };
+
+/**
+ * 时家奇门排盘
+ * @returns 完整盘面
+ */
+export function qimenCalc(y, m, d, hour) {
+  const { dunType, jieqi, yuan, ju } = dingJu(y, m, d);
+  const dayGz = ganzhiDayOfDate(y, m, d);
+  const hourZhi = Math.floor(((hour + 1) % 24) / 2);
+  // 时柱干支（五鼠遁简化：按日干起）
+  const hourGz = GANZHI[((GANZHI.indexOf(dayGz) % 10 % 5) * 2 + hourZhi) % 60]; // 简化：日干定子时干
+  const xun = xunShou(hourGz);
+  const yi = XUNSHOU_YI[xun];
+  const liuyi = sanqiLiuyi(dunType, ju);
+  // 值符宫 = 旬首六仪落宫
+  const fuGong = Object.keys(liuyi).find(p => liuyi[p] === yi) || '坎';
+  // 九星八门落宫：值符星随六仪落宫（简化：值符宫起布九星/八门）
+  const xingOrder = ['天蓬', '天任', '天冲', '天辅', '天英', '天芮', '天柱', '天心', '天禽'];
+  const menOrder = ['休门', '生门', '伤门', '杜门', '景门', '死门', '惊门', '开门'];
+  const board = {};
+  const startIdx = PALACE_NUMBERS.indexOf(ju) >= 0 ? ju - 1 : 0;
+  for (let i = 0; i < 9; i++) {
+    const p = NINE_PALACES[(startIdx + (dunType === 'yang' ? i : -i + 8)) % 9];
+    board[p] = { yi: liuyi[p] };
+  }
+  // 九星落宫（值符星在值符宫，其余按宫序）
+  const xingBoard = {};
+  for (let i = 0; i < 9; i++) {
+    const p = NINE_PALACES[i];
+    xingBoard[p] = i === 4 ? '天禽' : xingOrder[(i + (NINE_PALACES.indexOf(fuGong) >= 0 ? NINE_PALACES.indexOf(fuGong) : 0)) % 9];
+  }
+  return {
+    dunType, jieqi, yuan, ju, dayGz, hourGz, hourZhi, xunShou: xun,
+    liuyi, fuGong,
+    board,
+    xingBoard,
+    menBoard: menOrderBoard(dunType, ju),
+    shenBoard: shenBoard(dunType, fuGong),
+    fullEnumeration: 1080 // 状态空间：18 局 × 60 时
+  };
+}
+
+function menOrderBoard(dunType, ju) {
+  const men = ['休门', '生门', '伤门', '杜门', '景门', '死门', '惊门', '开门'];
+  const out = {};
+  const start = ju - 1;
+  for (let i = 0; i < 9; i++) {
+    const p = NINE_PALACES[i];
+    out[p] = i === 4 ? '中门' : men[(i - start + 8) % 8];
+  }
+  return out;
+}
+
+function shenBoard(dunType, fuGong) {
+  const out = {};
+  const start = NINE_PALACES.indexOf(fuGong);
+  for (let i = 0; i < 8; i++) {
+    const idx = dunType === 'yang' ? (start + i) % 9 : ((start - i) % 9 + 9) % 9;
+    out[NINE_PALACES[idx]] = BA_SHEN[i];
+  }
+  return out;
+}

package/ENGINES/yili/lib/shushu-calculus.js

@@ -0,0 +1,85 @@
+/**
+ * 易理工程学 · 文言微积分引擎（P4）
+ * 察机（势函数驻点）/ 应期（解时方程）/ 势函数 U(象,时)
+ */
+import { jieqiOfYear } from './almanac.js';
+
+/**
+ * 势函数 U(象, 时)：状态 × 时间 → ℝ
+ * 模型：U = Σ 因子贡献 × 时间调制（旺相休囚死系数）+ 周期项
+ * @param {object} state 状态：{ factors: [{name, wuxing, weight}] }
+ * @param {number} t 时间（天，从基准起）
+ * @param {object} seasonCoef 季节调制 {wuxing: coef}
+ */
+export function potential(state, t, seasonCoef) {
+  let U = 0;
+  for (const f of state.factors || []) {
+    const coef = seasonCoef ? (seasonCoef[f.wuxing] ?? 1) : 1;
+    U += (f.weight || 1) * coef;
+  }
+  // 周期项：年周期（365.2422 天）
+  U += (state.annualAmplitude || 0) * Math.cos((2 * Math.PI * t) / 365.2422 + (state.phase || 0));
+  return U;
+}
+
+/**
+ * 察机：求势函数驻点（dU/dt = 0）——数值法（黄金分割/扫描）
+ * @returns {Array<{t, U, kind: '峰'|'谷'|'鞍'}>}
+ */
+export function chaJi(state, seasonCoef, t0 = 0, t1 = 365.2422, step = 1) {
+  const pts = [];
+  for (let t = t0; t <= t1; t += step) {
+    pts.push({ t, U: potential(state, t, seasonCoef) });
+  }
+  const extremes = [];
+  for (let i = 1; i < pts.length - 1; i++) {
+    const a = pts[i - 1].U, b = pts[i].U, c = pts[i + 1].U;
+    if (b > a && b > c) extremes.push({ t: pts[i].t, U: b, kind: '峰' });
+    if (b < a && b < c) extremes.push({ t: pts[i].t, U: c, kind: '谷' });
+  }
+  return extremes;
+}
+
+/** 应期：解 U(象, t) = 阈值 的时（数值扫描） */
+export function yingqi(state, seasonCoef, threshold, t0 = 0, t1 = 365.2422, step = 1) {
+  const hits = [];
+  let prev = potential(state, t0, seasonCoef);
+  for (let t = t0 + step; t <= t1; t += step) {
+    const cur = potential(state, t, seasonCoef);
+    if ((prev - threshold) * (cur - threshold) <= 0) {
+      hits.push({ t, U: cur, approx: true });
+    }
+    prev = cur;
+  }
+  return hits;
+}
+
+/** 禾势（累积势能）：∫U dt（数值积分，梯形法） */
+export function heShi(state, seasonCoef, t0 = 0, t1 = 365.2422, step = 1) {
+  let sum = 0;
+  let prevT = t0, prevU = potential(state, t0, seasonCoef);
+  for (let t = t0 + step; t <= t1; t += step) {
+    const U = potential(state, t, seasonCoef);
+    sum += ((prevU + U) / 2) * (t - prevT);
+    prevT = t; prevU = U;
+  }
+  return sum;
+}
+
+/** 文言微积分报告 */
+export function shushuCalculus(state, seasonCoef, opts = {}) {
+  const extremes = chaJi(state, seasonCoef, opts.t0, opts.t1, opts.step);
+  const jq = [];
+  for (const e of extremes) {
+    const day = Math.floor(e.t);
+    jq.push({ 机: e.kind === '峰' ? '势之极（峰）' : '势之谷', 彳U: 0, 象时: `t=${day.toFixed(1)}天`, U: +e.U.toFixed(3) });
+  }
+  return {
+    象势函数: 'U = Σ 权重×季节系数 + A·cos(2πt/T + φ)',
+    察机: jq,
+    应期: opts.threshold !== undefined ? yingqi(state, seasonCoef, opts.threshold, opts.t0, opts.t1, opts.step) : [],
+    禾势: +heShi(state, seasonCoef, opts.t0, opts.t1, opts.step).toFixed(2),
+    规则: 'R-CAL-01 势函数模型（启发式，非因果断言）',
+    置信度: '结构 100%，结论启发式'
+  };
+}