← Archive
lm-002939 · 2026-08

跨領域壓力測試與可計算性驗證

下載 MD 檔 ⬇
📎 附件 · Companion files — 隨文交付的程式 / 證明 / 資料,可獨立下載重驗

跨領域壓力測試與可計算性驗證

——從純數學、有限離散域、物理公式到複數相位、Fourier、Excel 與 Python Backend

English Title: Cross-Domain Stress Testing and Computability Validation of Operation Translation: From Pure Mathematics and Finite Decision Domains to Physics, Complex Phase, Fourier Analysis, Excel, and Python Backends

作者: Neo.K
機構: EveMissLab / 一言諾科技有限公司
系列: Operation Translation Series A — Paper 07
版本: v0.1
日期: 2026-08-10


摘要

Series A Paper 01–06 已建立一套由判定域、運算轉譯、加法化、Correction Field、離散 valuation coordinate、局部 chart、tangent linearization、representation 與非交換 BCH correction 所組成的理論框架。本文不再新增母理論,而是回答一個更直接的問題:

上述框架在跨領域、跨數值尺度與跨 backend 的實際運算中是否仍然成立?\boxed{ \text{上述框架在跨領域、跨數值尺度與跨 backend 的實際運算中是否仍然成立?} }

本文整理並重現本系列先前完成的七組壓力測試:

  1. 高動態範圍實數乘除與 residual-log coordinate;
  2. 有限 Collatz-like partial map 與 exact discrete recovery;
  3. 維度感知物理公式;
  4. 複數相位、干涉與 Fourier / DFT;
  5. 純數學巨大乘積、組合數、Vandermonde、Euler product 與多項式乘積;
  6. commuting / noncommuting matrix 對照與 BCH correction hierarchy;
  7. Microsoft Excel 原生函數 backend 驗證。

本文特別建立四層證據制度:

E0: theorem / identity\boxed{ E_0:\text{ theorem / identity} } E1: high-precision computational validation\boxed{ E_1:\text{ high-precision computational validation} } E2: ordinary-backend reproducibility\boxed{ E_2:\text{ ordinary-backend reproducibility} } E3: negative control / structural failure\boxed{ E_3:\text{ negative control / structural failure} }

其目的在於避免把「數學上已證明」「高精度數值吻合」「普通軟體跑得動」與「某次實驗沒有出錯」混成同一種證據。

在 80 decimal-digit Python Decimal 實驗中,表示

x=s2kerx=s\,2^ke^r

成功處理包括

101000,10700,10500,1040010^{1000},\qquad10^{-700},\qquad10^{-500},\qquad10^{400}

等跨越數千數量級的輸入;500 組隨機乘除壓力測試中,所使用 precision 下重建相對誤差為零。零值乘法可由特殊 zero state 處理,而除以零被 Domain Validator 正確拒絕。

有限 partial-map 測試取

N=100000.N=100000.

共檢查 66,667 個合法狀態,其中 50,000 個 even branch 與 16,667 個 odd branch;33,333 個狀態因超出 partial-map 判定域而被標記為 undefined。所有合法狀態的 exact integer recovery failure 為:

0.\boxed{0}.

最大 decode-target 絕對差約為:

1.93×1075,1.93\times10^{-75},

遠低於該有限判定域的 uniform log-error margin:

4.9999875×106.4.9999875\times10^{-6}.

純數學壓力測試涵蓋一個 601 位十進位數的中央二項式係數:

(20001000),\binom{2000}{1000},

284 位數的 factorial ratio:

600!(200!)3,\frac{600!}{(200!)^3},

190 個 multiplicative factors 所構成的 20×2020\times20 Vandermonde determinant、1229 個質數所構成的有限 Euler product,以及 degree-80 root-product polynomial。這些測試均顯示 multiplicative structure 在 additive/log coordinates 中可穩定重建;其中離散整數案例透過 nearest-state / nearest-integer decoding 成功恢復 exact integer。

物理測試以 120 decimal digits 執行,包含:

V=4π3r3,V=\frac{4\pi}{3}r^3, T=2πLg,T=2\pi\sqrt{\frac Lg}, (2π)3=h3,(2\pi\hbar)^3=h^3, σ=2π5kB415h3c2,\sigma=\frac{2\pi^5k_B^4}{15h^3c^2},

blackbody flux、relativistic energy 與極端 Lorentz-factor cancellation。轉譯路徑與直接高精度參考值的相對差大多位於:

101191012010^{-119}\sim10^{-120}

量級。更重要的是,dimension-aware validator 成功拒絕「meter + kelvin」等數學形式可寫、物理維度不合法的加法,顯示判定域不應只檢查數值 domain,也可以包含型別/量綱規則。

複數與 Fourier 測試則展示另一種情況:multiplication 可以被 magnitude-log 與 phase addition 自然處理,但 addition / destructive interference 必須進入 correction / cancellation path。對

1+ei(π1040),1+e^{i(\pi-10^{-40})},

高精度模型得到幅值:

104010^{-40}

且相對誤差約:

2.13×10120.2.13\times10^{-120}.

128 點純 tone DFT 的 target bin 精確得到幅值 128,而其他 bins 僅留下約 1011810^{-118} 的 numerical residues。Gaussian Fourier integral 的 coordinate path 與同一 numerical quadrature 的 relative difference 約:

1.21×10119,1.21\times10^{-119},

而與 analytic full-line answer 的約 3.30×10233.30\times10^{-23} 差距來自有限積分窗口與 quadrature approximation,不是 operation translation 本身。

非交換矩陣測試則構成 negative control。Commuting diagonal matrices 有:

[A,B]=0,[A,B]=0,

且:

log(eAeB)(A+B)F1.46×1016,\|\log(e^Ae^B)-(A+B)\|_F \approx1.46\times10^{-16},

只留下 ordinary floating backend 尺度的誤差。對 noncommuting matrices,naive additive approximation error 為約:

0.042158.0.042158.

加入:

12[A,B]\frac12[A,B]

後降至約:

0.0035874,0.0035874,

再加入三階 BCH commutators 後降至:

0.0004214.0.0004214.

這提供了 Series A 最重要的結構失敗對照:不是所有乘法都能 faithful exact 地變成 ordinary addition;當 source operation 非交換時,correction hierarchy 是結構必要項。

最後,本文使用 Microsoft Excel 原生公式建立 ordinary-backend workbook。Excel 的數值精度規格為 15 digits;因此該測試刻意不要求 symbolic exactness,而只檢查這套方法是否能在一般辦公軟體中被實際執行。結果顯示 multiplication、division、powers、 π3\pi^3π5\pi^5 、Stefan–Boltzmann、有限 partial map、DFT 等皆可用標準 LNEXPPOWERPISQRTCOSSINSUM 等公式重現;其殘餘誤差落在 ordinary double-precision 所預期的尺度。

本文的最終結論不是「operation translation 比傳統運算永遠更快」,也不是「所有運算都能被加法化」,而是:

Series A 所提出的分類,可以被跨數學域與跨 backend 實際區分與驗證。\boxed{ \text{Series A 所提出的分類,可以被跨數學域與跨 backend 實際區分與驗證。} }

其中 exact、corrected、local/approximate、decision-exact 與 structurally non-additivizable 等類別,在實驗中均出現了可重現案例。

關鍵詞: 壓力測試、可計算性、Operation Translation、Logarithmic Coordinate、Exact Recovery、Fourier、BCH、Excel、Python Decimal、Cross-Domain Validation


1. Paper 07 的任務

前六篇回答的是:

應該如何定義與分類 operation translation?

本篇回答:

這套分類碰到真正數值、真正公式、真正 backend 時,會不會立刻失效?

因此 Paper 07 不再追求新的抽象結構。

它是一篇:

validation paper.\boxed{ \text{validation paper}. }

2. 驗證不是證明的替代品

必須先固定一個研究原則:

computational evidencemathematical proof.\boxed{ \text{computational evidence} \neq \text{mathematical proof}. }

例如:

ln(xy)=lnx+lny\ln(xy)=\ln x+\ln y

不需要靠 500 次 random test 證明。

它本身已是 theorem-level identity。

Random test 的用途是確認:

  • implementation 是否正確;
  • chosen chart 是否處理 edge cases;
  • numerical backend 是否符合預期;
  • Runtime 判定域是否正確。

3. 四層 Evidence Model

E0E_0:Theorem / Identity

例如:

ln(xy)=lnx+lny.\ln(xy)=\ln x+\ln y.

或:

Φ(q1q2)=Φ(q1)+Φ(q2).\Phi(q_1q_2)=\Phi(q_1)+\Phi(q_2).

這是數學結構證據。


E1E_1:High-Precision Computational Validation

使用 arbitrary / adjustable precision backend,

檢查:

T^1(ν^+C^)\widehat T^{-1} \left( \widehat\nu+\widehat C \right)

與高精度 reference 是否一致。


E2E_2:Ordinary Backend Reproducibility

例如:

  • binary64;
  • Microsoft Excel;
  • standard spreadsheet functions。

它檢查:

不靠特殊 symbolic engine,方法是否仍可實際執行?


E3E_3:Negative Control / Structural Failure

刻意選一個理論預測「不能 exact additive」的案例。

例如:

[A,B]0.[A,B]\neq0.

若 naive:

X+YX+Y

真的失敗,

而 BCH correction 改善,

則支持我們的分類機制確實有區別能力。


4. Backend 與母理論分離

本系列反覆強調:

Mathematical TransformNumerical Backend.\boxed{ \text{Mathematical Transform} \neq \text{Numerical Backend}. }

Python decimal 支援可調 precision 的 decimal arithmetic。

Microsoft Excel 則是典型 ordinary finite-precision spreadsheet backend。

所以同一:

T,ν,CT,\nu,C

可以在不同 backend 實現。

Paper 07 的目的之一就是確認這一點。


5. Stress Test A:高動態範圍實數乘除

實驗 representation:

x=s2ker,x = s\,2^ke^r,

其中:

s{1,+1},s\in\{-1,+1\}, kZ,k\in\mathbb Z, ln22r<ln22.-\frac{\ln2}{2} \leq r < \frac{\ln2}{2}.

這是 Paper 02 的 engineering chart。


6. 測試尺度

實驗包含:

x=101000,y=10700.x=10^{1000}, \qquad y=10^{-700}.

其乘積:

xy=10300,xy=10^{300},

而商:

x/y=101700.x/y=10^{1700}.

也包含:

x=10500,y=10400.x=10^{-500}, \qquad y=-10^{400}.

此時:

xy=10100,xy=-10^{-100}, x/y=10900.x/y=-10^{-900}.

7. 實數高精度結果

在 80 decimal-digit precision 下,

上述案例的 reconstructed multiplication / division relative error 在輸出 precision 上為:

0.0.

另外進行:

500500

組 random stress tests,

得到:

maxε×=0\boxed{ \max \varepsilon_{\times}=0 }

以及:

maxε÷=0\boxed{ \max \varepsilon_{\div}=0 }

在該 Decimal context 下成立。

注意這裡的「0」表示:

在所設定 decimal precision 與所比較的 Decimal reference 下沒有可見差異。

它不應被誤寫成:

有限精度機器已經變成數學實數。


8. Zero State 與 Domain Validator

實驗:

0×123.450\times123.45

得到:

0.0.

而:

123.45/0123.45/0

被拒絕。

這驗證 Paper 02 的結論:

0 是特殊 absorbing state,不屬於 nonzero logarithmic chart。\boxed{ 0 \text{ 是特殊 absorbing state,不屬於 nonzero logarithmic chart。} }

以及 Paper 01:

operation legality 必須先於 transform execution.\boxed{ \text{operation legality} \text{ 必須先於 transform execution}. }

9. Stress Test B:有限 Collatz-like Partial Map

本文不研究原始 Collatz conjecture。

只借用其兩種 operation branch:

nn/2n\mapsto n/2

與:

n3n+1.n\mapsto3n+1.

限制:

SN={1,,N}.S_N=\{1,\dots,N\}.

定義 partial map:

CN(n)={n/2,n0(mod2),3n+1,n1(mod2), 3n+1N,undefined,otherwise.C_N(n) = \begin{cases} n/2,& n\equiv0\pmod2,\\ 3n+1,& n\equiv1\pmod2,\ 3n+1\leq N,\\ \text{undefined},&\text{otherwise}. \end{cases}

10. 有限 Partial-Map 統計

取:

N=100000.N=100000.

總合法 states:

66667.\boxed{66667}.

其中:

5000050000

個 even branch,

1666716667

個 odd branch。

另外:

3333333333

個 state 為 boundary / undefined。

這些 undefined 狀態不是 failure。

它們正是 partial operation domain 的一部分。


11. Exact Recovery 結果

使用 80 decimal-digit precision,

全部合法 states:

6666766667

個中,

exact integer recovery failure:

0.\boxed{0}.

最大:

m^m|\widehat m-m|

約:

1.9345×1075.1.9345\times10^{-75}.

最大 log update error 約:

1.6868×1080.1.6868\times10^{-80}.

12. Decision Margin

由 Paper 04:

若:

ε<ln(1+12N),|\varepsilon| < \ln\left( 1+\frac1{2N} \right),

則:

round(elnm+ε)=m\operatorname{round} (e^{\ln m+\varepsilon}) = m

對所有:

mNm\leq N

成立。

在:

N=100000N=100000

時:

ln(1+1200000)4.9999875×106.\boxed{ \ln\left( 1+\frac1{200000} \right) \approx 4.9999875\times10^{-6}. }

實際誤差約:

1080.10^{-80}.

因此距離 decision boundary 尚有約:

7474

個十進位數量級以上的安全差距。


13. 這個測試驗證的是 Decision Exactness

它不是說:

L^=L\widehat L = L

成為 exact real equality。

它說:

L^ 沒有跨越 decision cell,\boxed{ \widehat L \text{ 沒有跨越 decision cell}, }

所以:

approximate coordinateexact discrete state.\boxed{ \text{approximate coordinate} \to \text{exact discrete state}. }

這正是 Paper 04 所定義的 Type III exactness。


14. Stress Test C:純數學巨大乘積

純數學測試刻意選擇:

  • 大階乘比;
  • 多因子乘積;
  • prime products;
  • 高次 polynomial product;

避免只測小型 toy arithmetic。


15. 中央二項式係數

測試:

(20001000).\binom{2000}{1000}.

它具有:

601\boxed{601}

位十進位數。

使用:

L=logΓ(2001)2logΓ(1001).L = \log\Gamma(2001)-2\log\Gamma(1001).

再經:

eLe^L

與 nearest integer decoder,

得到:

exact integer recovery = True.\boxed{ \text{exact integer recovery = True}. }

16. Multinomial / Factorial Ratio

測試:

600!(200!)3.\frac{600!}{(200!)^3}.

其有:

284284

位十進位數。

轉成:

L=logΓ(601)3logΓ(201).L = \log\Gamma(601)-3\log\Gamma(201).

nearest integer recovery:

True.\boxed{ \text{True}. }

17. Vandermonde Determinant

取:

xi=i,i=1,,20.x_i=i, \qquad i=1,\dots,20.

Vandermonde determinant:

detV=1i<j20(ji).\det V = \prod_{1\leq i<j\leq20} (j-i).

共有:

(202)=190\binom{20}{2} = 190

個 multiplicative factors。

log-coordinate:

logdetV=i<jlog(ji).\boxed{ \log|\det V| = \sum_{i<j}\log(j-i). }

determinant 有:

138138

位數,

並成功 exact integer recovery。


18. Euler Product

測試:

p10000(1p2)1.\prod_{p\leq10000} (1-p^{-2})^{-1}.

其中 prime 數量:

1229.\boxed{1229}.

直接有限乘積:

1.64491792074628641793395922763.1.64491792074628641793395922763.

log-coordinate reconstruction 得到相同顯示值。

兩條有限乘積路徑的 relative difference 約:

2.28×10900.2.28\times10^{-900}.

19. Truncation Error 與 Translation Error 分離

該 finite Euler product 與:

ζ(2)=π26\zeta(2) = \frac{\pi^2}{6}

仍有約:

9.82×1069.82\times10^{-6}

relative difference。

原因是只取:

p10000.p\leq10000.

因此:

Euler-product truncation erroroperation-translation error.\boxed{ \text{Euler-product truncation error} \neq \text{operation-translation error}. }

這是 Paper 07 的核心證據分類之一。


20. Degree-80 Root-Product Polynomial

測試:

P(x)=r=180(xr)P(x) = \prod_{r=1}^{80}(x-r)

在:

x=80.25x=80.25

處。

直接值約:

7.37088593850695×10116.7.37088593850695\times10^{116}.

representation:

sgnP+rlnxr.\operatorname{sgn}P + \sum_r\ln|x-r|.

relative difference 約:

4.50×10899.4.50\times10^{-899}.

21. Stress Test D:維度感知物理公式

物理測試不是為了證明物理常數。

而是測:

多個 π\pi 、冪次、根號、尺度與物理 dimension metadata 混合後,operation translation 是否仍可執行?


22. Physical Quantity Representation

使用:

Q(s,d,L),Q \longleftrightarrow (s,\mathbf d,L),

其中:

d\mathbf d

是 SI base dimension exponent vector,

而:

L=lnq.L=\ln|q|.

乘法:

d12=d1+d2,\mathbf d_{12} = \mathbf d_1+\mathbf d_2, L12=L1+L2.L_{12} = L_1+L_2.

除法:

d12=d1d2.\mathbf d_{12} = \mathbf d_1-\mathbf d_2.

23. 測試公式

包括:

V=4π3r3,V=\frac{4\pi}{3}r^3, T=2πL/g,T=2\pi\sqrt{L/g}, (2π)3=h3,(2\pi\hbar)^3=h^3, σ=2π5kB415h3c2,\sigma = \frac{2\pi^5k_B^4}{15h^3c^2}, F=σT4,F=\sigma T^4,

以及:

E=(pc)2+(mc2)2.E = \sqrt{(pc)^2+(mc^2)^2}.

24. 120-Digit 結果

在 120 decimal-digit precision 下,

代表性的 relative differences:

sphere:

1.15×10120,1.15\times10^{-120},

phase cell:

1.60×10119,1.60\times10^{-119},

Stefan–Boltzmann:

1.43×10119,1.43\times10^{-119},

blackbody flux:

1.65×10119,1.65\times10^{-119},

relativistic energy:

3.85×10120.3.85\times10^{-120}.

這些數值不是物理精度聲明,

而是 transformed path 與同 backend 高精度 direct reference 的一致性測試。


25. Dimension Illegality Test

測試:

1 meter+1 kelvin.1\ \mathrm{meter} + 1\ \mathrm{kelvin}.

系統拒絕:

Dimensional-illegality test = True.\boxed{ \text{Dimensional-illegality test = True}. }

這說明判定域可以包括:

semantic / type legality,\boxed{ \text{semantic / type legality}, }

而不限於:

x>0,y0x>0,\quad y\neq0

等純數值限制。


26. Extreme Cancellation

測:

β=11040.\beta = 1-10^{-40}.

在 high-precision backend 中:

β1.\beta\neq1.

但 ordinary binary64:

1.010401.0-10^{-40}

會 round 回:

1.0.1.0.

因此:

γ=(1β2)1/2\gamma = (1-\beta^2)^{-1/2}

若完全依賴 binary64 原始 subtraction,

會直接進錯誤 domain。

high-precision log-diff path 則得到:

γ7.0710678118654752440×1019.\gamma \approx 7.0710678118654752440\times10^{19}.

這是一個非常清楚的 backend-vs-theory 分界案例。


27. Stress Test E:複數相位與干涉

非零複數表示:

z=eLeiϕ.z = e^Le^{i\phi}.

multiplication:

LL1+L2,L\to L_1+L_2, ϕϕ1+ϕ2.\phi\to\phi_1+\phi_2.

所以 complex multiplication 仍然高度 additive。

但 addition 需要 phasor sum。


28. 近完全相消

測試:

1+ei(π1040).1+e^{i(\pi-10^{-40})}.

其幅值約:

1040.10^{-40}.

120-digit computation 得:

z1040\boxed{ |z|\approx10^{-40} }

relative difference 約:

2.13×10120.2.13\times10^{-120}.

這顯示 corrected addition 能處理非常強的 destructive interference,

前提是 backend precision 足夠。


29. Symbolic Phase Test 的正確定位

先前 V2 prototype 使用:

ϕ=qπ+δ\phi=q\pi+\delta

保存 rational- π\pi tag,

所以:

1+eiπ1+e^{i\pi}

可被 symbolic rule 判為 exact zero。

但 Paper 07 不把這視為母理論必要條件。

它只是:

one optional exact backend layer.\boxed{ \text{one optional exact backend layer}. }

是否支援 cyclotomic / root-of-unity exact algebra,

屬於 symbolic-engine scope。


30. 128-Point DFT

純 tone:

xn=ei2πk0n/Nx_n = e^{i2\pi k_0n/N}

測試:

N=128,k0=17.N=128, \qquad k_0=17.

target bin:

X17=128.\boxed{ |X_{17}|=128. }

代表性 off-target bins 僅留下:

1011810^{-118}

量級 numerical residual。

因此 DFT 本身不需要被重新發明;

測試只說明 complex phase + summation correction path 可以重建 Fourier cancellation structure。


31. Gaussian Fourier Integral

log-coordinate quadrature:

1.14761804968609211829948358551.1476180496860921182994835855\ldots

direct same-quadrature:

1.14761804968609211829948358551.1476180496860921182994835855\ldots

coordinate relative difference:

1.21×10119.1.21\times10^{-119}.

analytic full-line answer 與有限 numerical quadrature 差約:

3.30×1023.3.30\times10^{-23}.

所以:

coordinate errorquadrature/window error.\boxed{ \text{coordinate error} \ll \text{quadrature/window error}. }

32. Stress Test F:Noncommutative Negative Control

這是整篇最重要的 negative control。

如果所有案例都成功,

我們無法知道分類是不是只是:

只挑適合案例。

所以特意加入 noncommuting matrices。


33. Commuting Matrix Case

取 commuting diagonal matrices。

得到:

[A,B]F=0.\|[A,B]\|_F=0.

並有:

log(eAeB)(A+B)F1.46×1016.\|\log(e^Ae^B)-(A+B)\|_F \approx 1.46\times10^{-16}.

這和 ordinary floating backend noise 同尺度。


34. Noncommuting Matrix Case

取:

A=(00.2000),A= \begin{pmatrix} 0&0.20\\ 0&0 \end{pmatrix}, B=(000.300).B= \begin{pmatrix} 0&0\\ 0.30&0 \end{pmatrix}.

得到:

[A,B]F0.0848528.\|[A,B]\|_F \approx 0.0848528.

35. Naive Addition Failure

只取:

Z1=A+BZ_1=A+B

與:

Z=log(eAeB)Z=\log(e^Ae^B)

比較:

ZZ1F0.0421580.\boxed{ \|Z-Z_1\|_F \approx 0.0421580. }

這不是 floating-point 細微殘差。

它是明顯 structural mismatch。


36. 加入 BCH 二階

取:

Z2=A+B+12[A,B].Z_2 = A+B+\frac12[A,B].

誤差:

ZZ2F0.0035874.\boxed{ \|Z-Z_2\|_F \approx 0.0035874. }

相對 naive model 改善約:

11.7511.75

倍。


37. 加入 BCH 三階

加入:

112[A,[A,B]]+112[B,[B,A]].\frac1{12}[A,[A,B]] + \frac1{12}[B,[B,A]].

誤差降為:

0.0004214.\boxed{ 0.0004214. }

又改善約:

8.518.51

倍。

這正符合 Paper 06:

additive core+commutator correction hierarchy.\boxed{ \text{additive core} + \text{commutator correction hierarchy}. }

38. 這個實驗不證明 BCH 全域收斂

它只驗證:

對此局部矩陣案例,加入預測的 commutator terms 會逐步恢復 non-Abelian structure。

是否收斂仍需 BCH convergence theory 與 domain metadata。

因此:

E1E_1

不能冒充:

E0.E_0.

39. Stress Test G:Microsoft Excel

為了測試:

是否需要特殊 research backend 才能執行這些 transformed formulas?

建立:

excel_log_coordinate_function_stress_test.xlsx

使用原生函數:

  • LN
  • EXP
  • POWER
  • PI
  • SQRT
  • SIGN
  • IF
  • MOD
  • COS
  • SIN
  • SUM
  • SUMPRODUCT

等。


40. Excel 是刻意選擇的低門檻 Backend

Microsoft Excel 官方規格中,

number precision 為:

15 digits.\boxed{15\text{ digits}.}

所以 Excel 測試的目的不是:

證明 symbolic exactness。

而是:

在一般 office spreadsheet 的有限精度下,transformed formulas 是否仍能正確工作到 backend 所允許的程度?


41. Excel Real-Log Results

Real-log multiplication / division transformed formulas 的代表性 difference 落在:

1016101410^{-16} \sim 10^{-14}

量級。

這符合 ordinary finite-precision environment 的預期。


42. Excel π\pi 與物理公式

使用:

π3,π5,\pi^3, \qquad \pi^5,

roots,

Stefan–Boltzmann 型公式,

transformed path 均能以 ordinary formulas 計算。

因此:

πn\boxed{ \pi^n }

本身不是 operation translation 的核心障礙。

在 log-space:

ln(πn)=nlnπ.\ln(\pi^n)=n\ln\pi.

43. Excel Finite Partial Map

有限:

n/2n/2

與:

3n+13n+1

branch 測試可用 ordinary worksheet formulas 執行,

並可由 rounding 回復正確 finite state。

這再次展示:

decision exactness\boxed{ \text{decision exactness} }

不要求 symbolic CAS。


44. Excel DFT

對:

N=32,N=32,

target bin:

k0=5k_0=5

得到:

X5=32.\boxed{ |X_{5}|=32. }

off-target bins 留下:

1015101410^{-15}\sim10^{-14}

量級 residual。

這正是:

E2E_2

ordinary-backend evidence。


45. Excel Exact-Zero Residual

例如:

1+eiπ1+e^{i\pi}

在普通 finite floating representation 中可能留下:

1.224646799×1016.1.224646799\ldots\times10^{-16}.

Paper 07 對此不視為 framework failure。

因為:

symbolic exact zero\boxed{ \text{symbolic exact zero} }

是 CAS / exact phase backend 的責任,

而 ordinary Excel 的工作是 finite-precision numerical realization。


46. Validation Matrix

測試 Source operation Transform type 分類 證據
Positive real multiply/divide ×,÷\times,\div log / residual log Exact E0+E1E_0+E_1
Finite Collatz-like partial affine branches log + correction Corrected + Decision Exact E0+E1E_0+E_1
Binomial / factorial ratio products / ratios log-gamma additive Exact + Decision Exact E1E_1
Vandermonde large product sum of logs Exact + Decision Exact E1E_1
Euler product large product sum of logs Exact finite product E1E_1
Factored polynomial product sign + log magnitude Exact E1E_1
Physics multiplicative formulas products / powers dimension + log Exact on valid dimension domain E1E_1
Positive sum ++ log-sum-exp Corrected E0+E1+E2E_0+E_1+E_2
Positive difference - log-diff-exp Corrected, singular at zero E0+E1E_0+E_1
Complex multiplication ×\times log magnitude + phase Exact modulo phase topology E0+E1E_0+E_1
DFT / phasor sum ++ complex correction Corrected E1+E2E_1+E_2
Commuting matrices multiplication local matrix log Exact local additive E1E_1
Noncommuting matrices multiplication log + BCH Corrected E3E_3
Excel workbook mixed native formulas Backend reproducible E2E_2

47. 通過了什麼?

這一輪測試支持以下命題。

第一:

exact multiplicative additivization\boxed{ \text{exact multiplicative additivization} }

確實能在多種巨大乘積結構上工作。

第二:

correction framework\boxed{ \text{correction framework} }

不是只對 ax+bax+b 有用,

也涵蓋 log-sum-exp、log-diff-exp 與 BCH。

第三:

finite decision recovery\boxed{ \text{finite decision recovery} }

可以把 approximate coordinate 重新投影為 exact state。

第四:

domain legality\boxed{ \text{domain legality} }

可以包含:

  • 除零;
  • positivity;
  • partial map boundary;
  • physical dimensions;
  • logarithm branch;
  • noncommutative structure。

48. 沒有證明什麼?

本文沒有證明:

Operation Translation 永遠更快.\boxed{ \text{Operation Translation 永遠更快}. }

沒有證明:

所有運算都能變成 ordinary addition.\boxed{ \text{所有運算都能變成 ordinary addition}. }

也沒有證明:

log representation 在所有 hardware 上優於 ordinary arithmetic.\boxed{ \text{log representation 在所有 hardware 上優於 ordinary arithmetic}. }

這些都需要指定 cost model。


49. Factorization Cost 仍然存在

Paper 04 valuation model:

Q>0×pZ\mathbb Q_{>0}^{\times} \cong \bigoplus_p\mathbb Z

數學上完全 exact。

但若輸入是一個 arbitrary large integer,

取得:

Φ(n)\Phi(n)

可能需要 factorization。

因此:

operation after encoding cheap\boxed{ \text{operation after encoding cheap} }

不代表:

encoding itself cheap.\boxed{ \text{encoding itself cheap}. }

Paper 07 保留這個重要負面結果。


50. Addition Cost 也只是重新分配

Logarithmic coordinate 讓:

×\times

變:

+,+,

但 ordinary:

++

變成:

LSE.\operatorname{LSE}.

所以:

complexity is redistributed across the operation profile.\boxed{ \text{complexity is redistributed across the operation profile}. }

不能只比較單一 primitive。


51. Backend Failure 不等於 Theory Failure

例如 binary64:

11040=11-10^{-40}=1

的情況。

這是有限 representation 無法分辨兩數。

不是:

110401-10^{-40}

在數學上真的等於:

1.1.

同理 Excel 的:

1+eiπ1+e^{i\pi}

tiny residual,

也不能推翻:

eiπ+1=0.e^{i\pi}+1=0.

52. Theory Failure 也不能推給 Backend

反過來,

noncommuting matrix:

ABBAAB\neq BA

時:

log(eAeY)X+Y\log(e^Ae^Y) \neq X+Y

不是提高到 1000 digits 就會消失。

它是:

structural failure.\boxed{ \text{structural failure}. }

這兩種 failure mode 必須嚴格分開。


53. Error Taxonomy

本文因此整理出:

Type A:Representation Error

有限 bit / decimal digits。

Type B:Approximation Error

Taylor、quadrature、finite truncation。

Type C:Decision Error

approximate result 跨過 decision boundary。

Type D:Domain Error

在不合法輸入上執行 operation。

Type E:Structural Error

錯誤假設 source 可映到不相容 target structure。

Paper 07 的所有壓力測試都可以落入這五類之一。


54. Computability 的正確結論

本篇只主張:

所提出的核心 transforms 與 correction classes 可以在普通現代計算機 backend 上被實際實現。

Python Decimal 提供可調 precision arithmetic。

SciPy 已提供 logsumexp 等成熟 primitive。

Excel 也能用標準 worksheet functions 實現基本 transformed arithmetic。

因此:

這套框架不是只存在於形式符號層。\boxed{ \text{這套框架不是只存在於形式符號層。} }

55. 但「可計算」不等於「值得替換現有硬體」

是否值得建立:

  • dedicated LNS hardware;
  • GPU kernels;
  • SIMD;
  • cache table;
  • exact valuation engine;

屬於 Engineering Whitepaper 的問題。

Paper 07 只證明:

computational realization is feasible.\boxed{ \text{computational realization is feasible}. }

56. Runtime Validation Requirements

由本篇結果,

Runtime 至少需要六類 validator:

  1. Domain Validator
  2. Chart / Transform Validator
  3. Algebra / Commutativity Classifier
  4. Correction Validity Checker
  5. Precision / Error Validator
  6. Decision Recovery Validator

如果缺一層,

某些壓力測試就會失敗。


57. Reproducibility Artifacts

本系列已建立:

  • topological_log_coordinate_mvp.py
  • finite_collatz_additive_coordinates.py
  • dimension_aware_log_physics_stress.py
  • complex_phase_log_coordinate_stress_v2.py
  • pure_math_additive_coordinate_stress.py
  • excel_log_coordinate_function_stress_test.xlsx

以及各自的結果與 README。

這些 artifact 將在 Series A 最終工程 bundle 中統一歸檔。


58. Series A 的實證狀態

至 Paper 07,

各主要類別都已有代表案例:

Exact

xy,x/y,xα,xy,\quad x/y,\quad x^\alpha,

valuation multiplication,

Vandermonde product,

finite Euler product。

Corrected

ax+b,ax+b,

log-sum-exp,

log-diff-exp,

BCH。

Local

complex logarithm,

matrix / Lie logarithm。

Decision Exact

finite integer recovery。

Structurally Impossible as Exact Abelian Model

genuine noncommutative multiplication。

所以分類不是空集合分類。

每一類都實際出現。


59. 本篇主要結論

結論一

Operation Translation 的 exact / corrected / local / impossible classification 可以被真正數值案例區分。

結論二

高精度 backend 證明方法具有 numerical realizability,但不是 theorem 的替代品。

結論三

ordinary Excel backend 仍能執行基本 transformed formulas,顯示方法不依賴特殊硬體才能成立。

結論四

finite decision domains 使 approximate coordinate → exact state 成為可嚴格保證的機制。

結論五

noncommuting matrix negative control 驗證 structural obstruction 分類具有辨識力。

結論六

truncation、quadrature、floating error 與 structural correction 必須分開記錄。


60. 限制

第一,本文不是 performance benchmark paper,因此未比較 CPU cycles、GPU throughput、energy efficiency 或 silicon area。

第二,Python Decimal 的高精度結果不能直接外推到 binary64 hardware performance。

第三,Excel 測試是 ordinary-backend feasibility test,不是高精度 benchmark。

第四,物理案例只驗證 formula-level transformation 與 dimension legality,不宣稱建立新的物理理論。

第五,complex symbolic phase prototype 的 exact-zero 功能屬 optional symbolic backend。

第六,BCH 測試只構成局部 numerical case,不提供 universal convergence theorem。


61. 結論

Series A 從一個非常簡單的問題開始:

log(xy)=logx+logy.\log(xy)=\log x+\log y.

經過七篇論文,

問題已變成:

在什麼判定域、什麼 target structure、 什麼 correction class 與什麼 reconstruction requirement 下, 原 operation 可以被合法簡化?\boxed{ \text{在什麼判定域、什麼 target structure、 什麼 correction class 與什麼 reconstruction requirement 下, 原 operation 可以被合法簡化?} }

Paper 07 的跨領域壓力測試顯示:

這套語言可以同時描述:

  • 正實數乘除;
  • 有限 partial map;
  • exact discrete recovery;
  • 巨大組合數與乘積;
  • 物理 dimension;
  • complex phase;
  • Fourier cancellation;
  • Excel ordinary arithmetic;
  • noncommuting matrix obstruction。

也就是:

同一框架不只適用於單一 logarithm toy example。\boxed{ \text{同一框架不只適用於單一 logarithm toy example。} }

更重要的是,

測試中同時存在:

成功案例\boxed{ \text{成功案例} }

與:

預期失敗案例.\boxed{ \text{預期失敗案例}. }

這使 Series A 不只是:

找一堆能被 log 化的東西。

而是形成一套真正具有判別力的 classification framework。

至此,Series A 的七篇理論/驗證論文完成。

下一步不再增加論文數量。

後續工作為:

Operation Translation Runtime\boxed{ \text{Operation Translation Runtime} }

統一工程白皮書與 Runtime MVP,

把:

Domain ValidatorChart / Representation SelectorAlgebra ClassifierLinear CoreCorrection EngineInverse TransformError / Decision Validator\text{Domain Validator} \rightarrow \text{Chart / Representation Selector} \rightarrow \text{Algebra Classifier} \rightarrow \text{Linear Core} \rightarrow \text{Correction Engine} \rightarrow \text{Inverse Transform} \rightarrow \text{Error / Decision Validator}

正式做成可執行架構。


參考資料

  1. NIST Digital Library of Mathematical Functions, Chapter 4: Logarithm, Exponential, Powers.
  2. Python Documentation, decimal — Decimal fixed-point and floating-point arithmetic.
  3. SciPy Documentation, scipy.special.logsumexp.
  4. Microsoft Support, Excel specifications and limits — number precision: 15 digits.
  5. Microsoft Support, Excel precision / rounding documentation.
  6. Series A Papers 01–06 and accompanying computational artifacts.