问题标签 [numerical]

我正在开发一款游戏，我有必要检查给定数字高度图中的闭合路径：服务器和客户端使用此高度图来设置正确的坐标以移动等......现在，当用户继续前进时一个“特殊”的瓷砖，它会亮起......我的问题是：当用户在这些瓷砖上行走时会创建一条包含空瓷砖的封闭路径，服务器应该自动填充这条路径中的瓷砖......

它应该这样做： http ://www.youtube.com/watch?v= kAVUNE2NTUQ - 1:32

我确定我必须在这里或那里使用一些数学，但我不知道如何......我可以做一个“for”循环，但它会太长，问题是服务器需要做这个循环每次用户走路时...提前感谢您的回答，希望有人可以帮助我。

PS：我正在使用 C#

编辑：当用户在瓷砖上行走时，服务器会自动将 heightmap[X, Y] 替换为表示用户颜色的整数

描述（这是一个 hwk 问题）：

我不知道从哪里开始。我计划使用拉普拉斯扩展，但我不确定如何为 nxn 矩阵实现它。任何帮助，将不胜感激。

注意：我已经有了一个为 nxn 矩阵生成随机矩阵的函数。计算的时间也不成问题。我唯一遇到的问题是如何计算行列式。

不得不删除我的班级政策的问题描述 b/c。

给定一个小数 x，我想测试 x 是否在分母 9999 或更小的有理数的 10^-12 以内。显然，我可以通过查看 x、2x、3x 等，并查看其中任何一个是否足够接近整数来做到这一点。但是有没有更有效的算法？

我在数字代码中使用 Boost 的 uBLAS 并有一个“重型”求解器：

http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion

该代码运行良好，但是速度非常慢。经过一番研究，我发现了 UMFPACK，它是一个稀疏矩阵求解器（除其他外）。我的代码生成大型稀疏矩阵，我需要非常频繁地求逆（更正确地解决，逆矩阵的值无关紧要），因此 UMFPAck 和 BOOST 的 Sparse_Matrix 类似乎是幸福的婚姻。

UMFPACK 要求由三个向量指定的稀疏矩阵：条目计数、行索引和条目。（见示例）。

我的问题归结为，我可以从 BOOST 的稀疏矩阵类中有效地获取这三个向量吗？

我正在用 C++ 对重力进行数值模拟，我想在每次计算一个步骤时备份我的结果。

但是，按照我现在的方式，程序总是会覆盖文件。我想我可以通过始终将文本保存在不同的文件或变量中来解决，但我想知道是否有更简单的方法可以打开文本文件，这样我就不会覆盖它。

我当前的“备份代码”如下所示：

logsumexp有人知道提供-function的开源数字 C 库吗？

该logsumexp(a)函数计算数组 a 的分量的指数 log(e^{a_1}+...e^{a_n}) 之和，避免数值溢出。

在一个 CUDA 程序中，我最近使用

INFINITY从 math.h 到哪里

并且对得到不同的结果感到非常惊讶。gnu.org建议它们实际上应该表现类似。我错过了什么吗？是否不允许INFINITY在 CUDA 内核中使用？

编辑： 我刚刚发现isinf并注意到使用

给出与 INFINITY 检查相同的结果。为什么不是 isfinite(x)==!isinf(x)？

我正在尝试使用 optim 拟合 BFGS（和 L-BFGS-B）的非线性最小二乘问题。当我提供分析梯度时，线搜索异常终止，最终解决方案总是非常接近起点。但是，当我不提供渐变时，它似乎收敛得很好。这是否向任何人暗示了任何数字问题？我很确定渐变是正确的。这可能是一个缩放问题吗？谢谢你的帮助。

我正在寻找一个用于科学计算的数学库以在 Java / Scala 中使用。特别是我需要完整的椭圆积分和修正的贝塞尔函数。如果它是开源的，我会很好，但我想我将不得不采用那里的任何东西。一个 scipy（用于科学计算的 python 库）替换会很棒:-)

我的程序速度有问题。我想计算一个巨大数组中四个邻居的平均值。这是我的代码的一部分。你有什么想法如何改变最后一行？或者我应该使用另一个数组？