使用 matlab 编程三体问题
我正在尝试使用 matlab 对三体问题进行编程。我得到了月球在太空旋转框架中的轨迹公式。它基本上是 ydotdot, xdotdot=GM/(x^2+y^2)^3/2 公式。公式是…
optim 的 lnsrch 问题 - R、BFGS
我正在尝试使用 optim 来拟合 BFGS （和 L-BFGS-B）的非线性最小二乘问题。当我提供解析梯度时，线性搜索异常终止，并且最终解总是非常接近起点。然而…
在 CUDA 中测试无穷大
在 CUDA 程序中，我最近从使用 return x==INFINITY || x==-INFINITY 来自 math.h 的 INFINITY 来测试无穷大，转换为使用来自 math.h 的无穷大测试， r…
logsumexp 在 C 中的实现？
有人知道提供 logsumexp 函数的开源数字 C 库吗？ logsumexp(a) 函数计算数组 a 各分量的指数 log(e^{a_1}+...e^{a_n}) 之和，避免数值溢出。…
UMFPACK 和 BOOST 的 uBLAS 稀疏矩阵
我在数字代码中使用 Boost 的 uBLAS，并有一个“重”解算器： http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion…
使用 Python 求解行列式而不使用 scipy.linalg.det 的代码
描述（这是一个 hwk 问题）： 我不知道从哪里开始。我计划使用拉普拉斯展开式，但我不确定如何将其应用于 nxn 矩阵。任何帮助将不胜感激。 注意：我已…
将双精度数与整数进行同一性比较时，什么会转换成什么？
好的，我知道您通常不应该比较两个浮点数是否相等。然而，在 William Kahan 的浮点计算中对舍入的无意识评估是多么徒劳？ 他展示了以下代码（我相信是…
“e”是“e”吗？用于表示图形数据中无效数字的字符？
我有一个图形程序，它导出带有数字的数据文件，例如：-1.33227e-015和-4.02456e-016。 我长期以来对“e-”符号感到困惑。它是用来表示无效数字吗？我…
