跨領域壓力測試與可計算性驗證
——從純數學、有限離散域、物理公式到複數相位、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 所組成的理論框架。本文不再新增母理論,而是回答一個更直接的問題:
本文整理並重現本系列先前完成的七組壓力測試:
- 高動態範圍實數乘除與 residual-log coordinate;
- 有限 Collatz-like partial map 與 exact discrete recovery;
- 維度感知物理公式;
- 複數相位、干涉與 Fourier / DFT;
- 純數學巨大乘積、組合數、Vandermonde、Euler product 與多項式乘積;
- commuting / noncommuting matrix 對照與 BCH correction hierarchy;
- Microsoft Excel 原生函數 backend 驗證。
本文特別建立四層證據制度:
其目的在於避免把「數學上已證明」「高精度數值吻合」「普通軟體跑得動」與「某次實驗沒有出錯」混成同一種證據。
在 80 decimal-digit Python Decimal 實驗中,表示
成功處理包括
等跨越數千數量級的輸入;500 組隨機乘除壓力測試中,所使用 precision 下重建相對誤差為零。零值乘法可由特殊 zero state 處理,而除以零被 Domain Validator 正確拒絕。
有限 partial-map 測試取
共檢查 66,667 個合法狀態,其中 50,000 個 even branch 與 16,667 個 odd branch;33,333 個狀態因超出 partial-map 判定域而被標記為 undefined。所有合法狀態的 exact integer recovery failure 為:
最大 decode-target 絕對差約為:
遠低於該有限判定域的 uniform log-error margin:
純數學壓力測試涵蓋一個 601 位十進位數的中央二項式係數:
284 位數的 factorial ratio:
190 個 multiplicative factors 所構成的 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 執行,包含:
blackbody flux、relativistic energy 與極端 Lorentz-factor cancellation。轉譯路徑與直接高精度參考值的相對差大多位於:
量級。更重要的是,dimension-aware validator 成功拒絕「meter + kelvin」等數學形式可寫、物理維度不合法的加法,顯示判定域不應只檢查數值 domain,也可以包含型別/量綱規則。
複數與 Fourier 測試則展示另一種情況:multiplication 可以被 magnitude-log 與 phase addition 自然處理,但 addition / destructive interference 必須進入 correction / cancellation path。對
高精度模型得到幅值:
且相對誤差約:
128 點純 tone DFT 的 target bin 精確得到幅值 128,而其他 bins 僅留下約 的 numerical residues。Gaussian Fourier integral 的 coordinate path 與同一 numerical quadrature 的 relative difference 約:
而與 analytic full-line answer 的約 差距來自有限積分窗口與 quadrature approximation,不是 operation translation 本身。
非交換矩陣測試則構成 negative control。Commuting diagonal matrices 有:
且:
只留下 ordinary floating backend 尺度的誤差。對 noncommuting matrices,naive additive approximation error 為約:
加入:
後降至約:
再加入三階 BCH commutators 後降至:
這提供了 Series A 最重要的結構失敗對照:不是所有乘法都能 faithful exact 地變成 ordinary addition;當 source operation 非交換時,correction hierarchy 是結構必要項。
最後,本文使用 Microsoft Excel 原生公式建立 ordinary-backend workbook。Excel 的數值精度規格為 15 digits;因此該測試刻意不要求 symbolic exactness,而只檢查這套方法是否能在一般辦公軟體中被實際執行。結果顯示 multiplication、division、powers、 、 、Stefan–Boltzmann、有限 partial map、DFT 等皆可用標準 LN、EXP、POWER、PI、SQRT、COS、SIN、SUM 等公式重現;其殘餘誤差落在 ordinary double-precision 所預期的尺度。
本文的最終結論不是「operation translation 比傳統運算永遠更快」,也不是「所有運算都能被加法化」,而是:
其中 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 不再追求新的抽象結構。
它是一篇:
2. 驗證不是證明的替代品
必須先固定一個研究原則:
例如:
不需要靠 500 次 random test 證明。
它本身已是 theorem-level identity。
Random test 的用途是確認:
- implementation 是否正確;
- chosen chart 是否處理 edge cases;
- numerical backend 是否符合預期;
- Runtime 判定域是否正確。
3. 四層 Evidence Model
:Theorem / Identity
例如:
或:
這是數學結構證據。
:High-Precision Computational Validation
使用 arbitrary / adjustable precision backend,
檢查:
與高精度 reference 是否一致。
:Ordinary Backend Reproducibility
例如:
- binary64;
- Microsoft Excel;
- standard spreadsheet functions。
它檢查:
不靠特殊 symbolic engine,方法是否仍可實際執行?
:Negative Control / Structural Failure
刻意選一個理論預測「不能 exact additive」的案例。
例如:
若 naive:
真的失敗,
而 BCH correction 改善,
則支持我們的分類機制確實有區別能力。
4. Backend 與母理論分離
本系列反覆強調:
Python decimal 支援可調 precision 的 decimal arithmetic。
Microsoft Excel 則是典型 ordinary finite-precision spreadsheet backend。
所以同一:
可以在不同 backend 實現。
Paper 07 的目的之一就是確認這一點。
5. Stress Test A:高動態範圍實數乘除
實驗 representation:
其中:
這是 Paper 02 的 engineering chart。
6. 測試尺度
實驗包含:
其乘積:
而商:
也包含:
此時:
7. 實數高精度結果
在 80 decimal-digit precision 下,
上述案例的 reconstructed multiplication / division relative error 在輸出 precision 上為:
另外進行:
組 random stress tests,
得到:
以及:
在該 Decimal context 下成立。
注意這裡的「0」表示:
在所設定 decimal precision 與所比較的 Decimal reference 下沒有可見差異。
它不應被誤寫成:
有限精度機器已經變成數學實數。
8. Zero State 與 Domain Validator
實驗:
得到:
而:
被拒絕。
這驗證 Paper 02 的結論:
以及 Paper 01:
9. Stress Test B:有限 Collatz-like Partial Map
本文不研究原始 Collatz conjecture。
只借用其兩種 operation branch:
與:
限制:
定義 partial map:
10. 有限 Partial-Map 統計
取:
總合法 states:
其中:
個 even branch,
個 odd branch。
另外:
個 state 為 boundary / undefined。
這些 undefined 狀態不是 failure。
它們正是 partial operation domain 的一部分。
11. Exact Recovery 結果
使用 80 decimal-digit precision,
全部合法 states:
個中,
exact integer recovery failure:
最大:
約:
最大 log update error 約:
12. Decision Margin
由 Paper 04:
若:
則:
對所有:
成立。
在:
時:
實際誤差約:
因此距離 decision boundary 尚有約:
個十進位數量級以上的安全差距。
13. 這個測試驗證的是 Decision Exactness
它不是說:
成為 exact real equality。
它說:
所以:
這正是 Paper 04 所定義的 Type III exactness。
14. Stress Test C:純數學巨大乘積
純數學測試刻意選擇:
- 大階乘比;
- 多因子乘積;
- prime products;
- 高次 polynomial product;
避免只測小型 toy arithmetic。
15. 中央二項式係數
測試:
它具有:
位十進位數。
使用:
再經:
與 nearest integer decoder,
得到:
16. Multinomial / Factorial Ratio
測試:
其有:
位十進位數。
轉成:
nearest integer recovery:
17. Vandermonde Determinant
取:
Vandermonde determinant:
共有:
個 multiplicative factors。
log-coordinate:
determinant 有:
位數,
並成功 exact integer recovery。
18. Euler Product
測試:
其中 prime 數量:
直接有限乘積:
log-coordinate reconstruction 得到相同顯示值。
兩條有限乘積路徑的 relative difference 約:
19. Truncation Error 與 Translation Error 分離
該 finite Euler product 與:
仍有約:
relative difference。
原因是只取:
因此:
這是 Paper 07 的核心證據分類之一。
20. Degree-80 Root-Product Polynomial
測試:
在:
處。
直接值約:
representation:
relative difference 約:
21. Stress Test D:維度感知物理公式
物理測試不是為了證明物理常數。
而是測:
多個 、冪次、根號、尺度與物理 dimension metadata 混合後,operation translation 是否仍可執行?
22. Physical Quantity Representation
使用:
其中:
是 SI base dimension exponent vector,
而:
乘法:
除法:
23. 測試公式
包括:
以及:
24. 120-Digit 結果
在 120 decimal-digit precision 下,
代表性的 relative differences:
sphere:
phase cell:
Stefan–Boltzmann:
blackbody flux:
relativistic energy:
這些數值不是物理精度聲明,
而是 transformed path 與同 backend 高精度 direct reference 的一致性測試。
25. Dimension Illegality Test
測試:
系統拒絕:
這說明判定域可以包括:
而不限於:
等純數值限制。
26. Extreme Cancellation
測:
在 high-precision backend 中:
但 ordinary binary64:
會 round 回:
因此:
若完全依賴 binary64 原始 subtraction,
會直接進錯誤 domain。
high-precision log-diff path 則得到:
這是一個非常清楚的 backend-vs-theory 分界案例。
27. Stress Test E:複數相位與干涉
非零複數表示:
multiplication:
所以 complex multiplication 仍然高度 additive。
但 addition 需要 phasor sum。
28. 近完全相消
測試:
其幅值約:
120-digit computation 得:
relative difference 約:
這顯示 corrected addition 能處理非常強的 destructive interference,
前提是 backend precision 足夠。
29. Symbolic Phase Test 的正確定位
先前 V2 prototype 使用:
保存 rational- tag,
所以:
可被 symbolic rule 判為 exact zero。
但 Paper 07 不把這視為母理論必要條件。
它只是:
是否支援 cyclotomic / root-of-unity exact algebra,
屬於 symbolic-engine scope。
30. 128-Point DFT
純 tone:
測試:
target bin:
代表性 off-target bins 僅留下:
量級 numerical residual。
因此 DFT 本身不需要被重新發明;
測試只說明 complex phase + summation correction path 可以重建 Fourier cancellation structure。
31. Gaussian Fourier Integral
log-coordinate quadrature:
direct same-quadrature:
coordinate relative difference:
analytic full-line answer 與有限 numerical quadrature 差約:
所以:
32. Stress Test F:Noncommutative Negative Control
這是整篇最重要的 negative control。
如果所有案例都成功,
我們無法知道分類是不是只是:
只挑適合案例。
所以特意加入 noncommuting matrices。
33. Commuting Matrix Case
取 commuting diagonal matrices。
得到:
並有:
這和 ordinary floating backend noise 同尺度。
34. Noncommuting Matrix Case
取:
得到:
35. Naive Addition Failure
只取:
與:
比較:
這不是 floating-point 細微殘差。
它是明顯 structural mismatch。
36. 加入 BCH 二階
取:
誤差:
相對 naive model 改善約:
倍。
37. 加入 BCH 三階
加入:
誤差降為:
又改善約:
倍。
這正符合 Paper 06:
38. 這個實驗不證明 BCH 全域收斂
它只驗證:
對此局部矩陣案例,加入預測的 commutator terms 會逐步恢復 non-Abelian structure。
是否收斂仍需 BCH convergence theory 與 domain metadata。
因此:
不能冒充:
39. Stress Test G:Microsoft Excel
為了測試:
是否需要特殊 research backend 才能執行這些 transformed formulas?
建立:
excel_log_coordinate_function_stress_test.xlsx
使用原生函數:
LNEXPPOWERPISQRTSIGNIFMODCOSSINSUMSUMPRODUCT
等。
40. Excel 是刻意選擇的低門檻 Backend
Microsoft Excel 官方規格中,
number precision 為:
所以 Excel 測試的目的不是:
證明 symbolic exactness。
而是:
在一般 office spreadsheet 的有限精度下,transformed formulas 是否仍能正確工作到 backend 所允許的程度?
41. Excel Real-Log Results
Real-log multiplication / division transformed formulas 的代表性 difference 落在:
量級。
這符合 ordinary finite-precision environment 的預期。
42. Excel 與物理公式
使用:
roots,
Stefan–Boltzmann 型公式,
transformed path 均能以 ordinary formulas 計算。
因此:
本身不是 operation translation 的核心障礙。
在 log-space:
43. Excel Finite Partial Map
有限:
與:
branch 測試可用 ordinary worksheet formulas 執行,
並可由 rounding 回復正確 finite state。
這再次展示:
不要求 symbolic CAS。
44. Excel DFT
對:
target bin:
得到:
off-target bins 留下:
量級 residual。
這正是:
ordinary-backend evidence。
45. Excel Exact-Zero Residual
例如:
在普通 finite floating representation 中可能留下:
Paper 07 對此不視為 framework failure。
因為:
是 CAS / exact phase backend 的責任,
而 ordinary Excel 的工作是 finite-precision numerical realization。
46. Validation Matrix
| 測試 | Source operation | Transform type | 分類 | 證據 |
|---|---|---|---|---|
| Positive real multiply/divide | log / residual log | Exact | ||
| Finite Collatz-like | partial affine branches | log + correction | Corrected + Decision Exact | |
| Binomial / factorial ratio | products / ratios | log-gamma additive | Exact + Decision Exact | |
| Vandermonde | large product | sum of logs | Exact + Decision Exact | |
| Euler product | large product | sum of logs | Exact finite product | |
| Factored polynomial | product | sign + log magnitude | Exact | |
| Physics multiplicative formulas | products / powers | dimension + log | Exact on valid dimension domain | |
| Positive sum | log-sum-exp | Corrected | ||
| Positive difference | log-diff-exp | Corrected, singular at zero | ||
| Complex multiplication | log magnitude + phase | Exact modulo phase topology | ||
| DFT / phasor sum | complex correction | Corrected | ||
| Commuting matrices | multiplication | local matrix log | Exact local additive | |
| Noncommuting matrices | multiplication | log + BCH | Corrected | |
| Excel workbook | mixed | native formulas | Backend reproducible |
47. 通過了什麼?
這一輪測試支持以下命題。
第一:
確實能在多種巨大乘積結構上工作。
第二:
不是只對 有用,
也涵蓋 log-sum-exp、log-diff-exp 與 BCH。
第三:
可以把 approximate coordinate 重新投影為 exact state。
第四:
可以包含:
- 除零;
- positivity;
- partial map boundary;
- physical dimensions;
- logarithm branch;
- noncommutative structure。
48. 沒有證明什麼?
本文沒有證明:
沒有證明:
也沒有證明:
這些都需要指定 cost model。
49. Factorization Cost 仍然存在
Paper 04 valuation model:
數學上完全 exact。
但若輸入是一個 arbitrary large integer,
取得:
可能需要 factorization。
因此:
不代表:
Paper 07 保留這個重要負面結果。
50. Addition Cost 也只是重新分配
Logarithmic coordinate 讓:
變:
但 ordinary:
變成:
所以:
不能只比較單一 primitive。
51. Backend Failure 不等於 Theory Failure
例如 binary64:
的情況。
這是有限 representation 無法分辨兩數。
不是:
在數學上真的等於:
同理 Excel 的:
tiny residual,
也不能推翻:
52. Theory Failure 也不能推給 Backend
反過來,
noncommuting matrix:
時:
不是提高到 1000 digits 就會消失。
它是:
這兩種 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。
因此:
55. 但「可計算」不等於「值得替換現有硬體」
是否值得建立:
- dedicated LNS hardware;
- GPU kernels;
- SIMD;
- cache table;
- exact valuation engine;
屬於 Engineering Whitepaper 的問題。
Paper 07 只證明:
56. Runtime Validation Requirements
由本篇結果,
Runtime 至少需要六類 validator:
- Domain Validator
- Chart / Transform Validator
- Algebra / Commutativity Classifier
- Correction Validity Checker
- Precision / Error Validator
- Decision Recovery Validator
如果缺一層,
某些壓力測試就會失敗。
57. Reproducibility Artifacts
本系列已建立:
topological_log_coordinate_mvp.pyfinite_collatz_additive_coordinates.pydimension_aware_log_physics_stress.pycomplex_phase_log_coordinate_stress_v2.pypure_math_additive_coordinate_stress.pyexcel_log_coordinate_function_stress_test.xlsx
以及各自的結果與 README。
這些 artifact 將在 Series A 最終工程 bundle 中統一歸檔。
58. Series A 的實證狀態
至 Paper 07,
各主要類別都已有代表案例:
Exact
valuation multiplication,
Vandermonde product,
finite Euler product。
Corrected
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 從一個非常簡單的問題開始:
經過七篇論文,
問題已變成:
Paper 07 的跨領域壓力測試顯示:
這套語言可以同時描述:
- 正實數乘除;
- 有限 partial map;
- exact discrete recovery;
- 巨大組合數與乘積;
- 物理 dimension;
- complex phase;
- Fourier cancellation;
- Excel ordinary arithmetic;
- noncommuting matrix obstruction。
也就是:
更重要的是,
測試中同時存在:
與:
這使 Series A 不只是:
找一堆能被 log 化的東西。
而是形成一套真正具有判別力的 classification framework。
至此,Series A 的七篇理論/驗證論文完成。
下一步不再增加論文數量。
後續工作為:
統一工程白皮書與 Runtime MVP,
把:
正式做成可執行架構。
參考資料
- NIST Digital Library of Mathematical Functions, Chapter 4: Logarithm, Exponential, Powers.
- Python Documentation,
decimal— Decimal fixed-point and floating-point arithmetic. - SciPy Documentation,
scipy.special.logsumexp. - Microsoft Support, Excel specifications and limits — number precision: 15 digits.
- Microsoft Support, Excel precision / rounding documentation.
- Series A Papers 01–06 and accompanying computational artifacts.