《冯洛伊曼物理常数库》

## 【前言：写给每一个被NaN熬过夜的开发者】 你是否经历过：跑了一夜的流体模拟，凌晨三点打开电脑，屏幕上铺满了NaN；调了一周的优化器，始终卡在局部最优一动不动；明明公式推导完全正确，代码跑出来却是热下沉冷上升的反物理世界。 计算机，本质上是数学的完美与物理的极限在硅片上的强行拼凑。数学允许除以零，允许无穷大，允许绝对连续；但冯·诺依曼架构的寄存器不行。NaN是瘟疫，0乘以它依然是NaN，传染整个数组；Inf是雪崩，瞬间摧毁所有边界。机器不懂你的推导，它只懂溢出。 比崩溃更让人掉头发的，是那些反人类的易错公式：扩散差分忘除dx²，Adam偏置校正写反，Boussinesq浮力少个负号……人类的直觉在这些繁复的数学细节面前不堪一击，一错就是灾难。 这个库，就是我花了五年时间，踩过所有你能想到和想不到的坑，攒出来的终极补丁集。把防爆、连续、反卡死、易错公式全焊死在底层。从此，数学归数学，物理归物理，寄存器的破事，不再忍受。 **库不再新增核心功能，只做后期维护，完全免费开源！** 无任何商业目的，无任何隐藏收费，仅为分享技术，避免同行重复踩坑。 网盘分享：silicon-physics 链接: https://pan.baidu.com/s/1saDFtatrAQ4SRrxSdpIl1w 提取码: 1949 --- ## 【我们敬畏的底线：写在所有代码之前】 这些常数不是玄学，是无数次崩溃换来的血的教训。动一个，你的程序就可能在某个凌晨悄无声息地死去。 - **EPSILON_MACHINE (1e-10)**：宇宙的真空下限，任何分母必须穿上这层防护服 - **CFL_SAFE_CONVECTIVE (0.45)**：对流因果律的锁，时间步长必须追上空间波速 - **CFL_SAFE_DIFFUSIVE (0.25)**：扩散因果律的锁，比对流更严！别混用 - **CUTOFF_STRESS (0.001)**：存在感的阈值，低于此值的扰动在宏观场不复存在 - **ENTROPY_FLOOR (1e-6)**：熵的地板，低于此值系统判定为热寂停滞，触发救援 ## 【本库四大铁律：违反一条，不如不用】 1. **常参是确定的**：不靠玄学调参，只查物理手册（水就是998，岩浆就是2600） 2. **连续即稳定**：用 tanh/softsign/smoothstep 替代一切截断，消灭数值爆炸 3. **无分支即并行**：一切逻辑转乘法掩码，GPU全员满载飞驰 4. **空间即因果**：降维映射防维度灾难，拓扑边传导定结构因果 --- ## 【十卷天书：解决你99%的底层计算问题】 不用再翻维基百科，不用再查论文附录，不用再对着Stack Overflow复制粘贴可能有错的代码。 1. **宇宙常参字典**：18大类数百个物理常数，业务代码绝不出现魔法数字 2. **无分支算子核武库**：13层安全防护，从根源上消灭NaN和Inf 3. **概率与随机过程**：精确实现所有常用分布和随机过程，含Itô修正项 4. **优化器与神经网络**：Adam偏置校正再也不会写反，激活函数全防爆 5. **信号、控制与几何**：滤波、PID、3D几何，所有易错公式已预编译 6. **高级插值与物理动力学**：动画、碰撞、积分，系数和公式一次写对 7. **金融工程**：Black-Scholes、风险度量，收益率方向再也不会搞反 8. **线性代数与傅里叶**：2x2/3x3矩阵运算、FFT，符号和索引再也不会错 9. **卡尔曼滤波与统计基础**：在线统计、传感器融合，内置反热寂机制 10. **卷积与图像算子**：常用卷积核一键调用，边界处理已优化 --- ## 【一行代码的差距：用与不用的区别】 ### 手写代码（90%的开发者会在这里犯错） ```rust // 可能的错误： // 1. 忘记加eps导致除零爆炸 // 2. 浮力方向搞反（热应该上升，这里写成了正y方向） // 3. 时间步长随便写0.001，可能导致CFL不稳定 // 4. 没有反停滞机制，系统会卡在局部最优 let buoyancy = G_EARTH * alpha * (t_local - t_ref); let dt = 0.001; ``` ### 使用本库（零错误，零思考，直接用） ```rust // 自动包含符号修正、除零保护、CFL安全步长和反停滞救援 let buoyancy = boussinesq_buoyancy(G_EARTH, water.alpha_thermal, 400.0, 300.0); let dt = safe_dt_convection(c_max, dx).min(safe_dt_diffusion(nu, dx)); let sigma = adaptive_sigma(variance, SIGMA_NOISE, STAGNATION_THRESHOLD, THERMAL_KICK_STRENGTH); ``` **别人写100行还可能有bug的代码，你3行搞定，而且绝对稳定。** 省下的时间，去喝咖啡，去陪家人，去做真正有创造性的工作。 --- ## 【V4.0 致命修复：这些坑我们已经替你踩过了】 - 修复了norm_inv函数掩码方向反转的致命错误，所有分位数计算结果完全正确 - 彻底解决了NaN穿透问题（0*NaN=NaN的IEEE754铁律），改用位运算实现真正的nan_to_zero - 修复了统计函数在count=0时的下溢panic，debug模式再也不会突然崩溃 - 为卡尔曼滤波加入了反热寂机制，Q=0时再也不会拒绝一切观测 - 全库统一使用safe_log和safe_exp，从源头杜绝Inf爆炸 --- ## 【零依赖，开箱即用】 - 纯Rust编写，零外部依赖，编译速度极快 - 支持所有主流平台：Windows/Linux/macOS - 所有接口经过严格测试，覆盖所有边界条件 - 完整的文档和示例，上手即用 ## 【最后】 这个库是我无数个熬夜调试的夜晚换来的。现在，我把它免费分享给你。希望它能帮你少掉一些头发，多睡几个好觉。 如果觉得有用，不妨分享给更多需要的人。技术的意义，在于互相照亮。 作者联系方式： 邮箱：jxjx24no@163.com QQ：2623233520 **物理法则，调用完毕。无分支，无爆炸，只有真理。** jxjx24no@163.com qq 2623233520