numerical-computing

0热度

2回答

如何在导入org.apache.commons.math3.analysis.integration.SimpsonIntegrator时使用标准的SimpsonIntegrator;

public double evalute(double distance){ /** * equation (3.2) */ this.from = 0; this.to = distance; this.n = 2; return - 10 * Math.log10(Math.exp(-IntSimpson(this.fr

0热度

1回答

加速到位移

我正在做一个实验程序。我在汽车仪表板上安装了一个智能手机，在旅行期间，我从应用程序中读取了一些惯性传感器的价值。准确地说，我读了每秒注册的加速度计数据时间（以固定间隔）。所以，现在我想从垂直加速度转移到垂直位移，为此，我应该进行双重积分。我试过欧拉法。从最初的条件，如： v0=0.v0=0这是零时间的初始速度。 x0=0.x0=0这是零时刻的初始位置。限定， deltaT=registra

-4热度

1回答

如何在Python

解决指数方程现在我有一个方程求解： EXP（X * A）-exp（X * B）= c，其中a，b和c是已知的常数。我试图sympy和scipy.optimize.fsolve，甚至brenth和牛顿。没什么好的。我对python很陌生，就像2周。所以请帮我解决这个问题。谢谢！

3热度

1回答

在计算图中避免虚函数调用

我正在使用DAG（有向无环图）来表示和计算表达式;每个节点代表一个操作（+,-, /,*，积累等等），并且通过按照拓扑排序顺序依次评估每个节点来实现对整个表达式的评估。根据它代表的运算符，每个节点继承基类RefNode并实现虚函数评估。 Node类在代表操作符的仿函数上进行模板化。节点评估顺序保存在vector<RefNode*>中，调用每个元素的->evaluate()。某些快速剖析显示虚拟e

0热度

1回答

Romberg积分算法

我想解决在给定区间[a，b]上的基本积分，对于Romberg积分中的R（n，m）中的任意n和m值。我从梯形法则中推导出布尔的规则，所以我知道如何在纸上做到这一点。我甚至画出了一张显示所有依赖关系的流程图。这不是帮助我编码。我有一种感觉，这需要某种递归。我正在用Java编程。编辑：我没有问任何人为我代码。上述所有信息，包括其他人编辑的信息，都是为了说明我的相关经验和对技术词汇的理解以及我应用

0热度

2回答

控制C代码中的操作顺序

我读到Fortran对评估表达式的顺序有严格的规则。对于某些数值算法，这非常重要。数值C程序如何控制浮点操作的顺序，并防止编译器“优化”到不合需要的操作顺序，例如将(a*b)*c更改为a*(b*c)。

17热度

5回答

为什么Elixir在Ruby中最慢并且在解决Project Euler＃5方面走得很慢？

更新：药剂不慢，我的算法是。我的算法甚至不是苹果来比较苹果。有关Ruby和Go等效算法，请参阅下面的罗马答案。还要感谢José，只需在MIX_ENV = prod的前面加上我的慢速算法即可。我已经更新了问题中的统计数据。原始问题：我正在研究多种语言的欧拉问题，只是为了了解语言的多产和多快。在problem #5中，我们被要求找到能被1到20中的所有数字均匀划分的最小正数。我以多种语言实现了解

2热度

1回答

防止下溢性计算正常样品落在一定间隔在给定两个实数<code>x</code>和<code>y</code>蟒

的概率的对数时，我想计算在python下面的函数： log Pr [ x <= t <= y ], 其中t被采样从正态分布。一个天真的实现是使用scipy.stats.norm。 np.log(scipy.stats.norm.cdf(y) - scipy.stats.norm.cdf(x)) 不幸的是，这会导致下溢时x和y远离0。如何避免这样的数字错误？

2热度

1回答

QuickHull最坏情况

QHull（也许其他好的实现QuickHull）在许多情况下工作得非常好，速度很快。然而，我们从理论上知道它的最坏情况可能是O（n^2）。在实践中，我还没有看到QHull工作很差的许多维度（即20或100）的数值例子。你知道一个数值例子，其中QHull工作得不好，或者给出了错误的结果，或者任何表明它不能在这里应用的数值例子。

0热度

1回答

使用动态数组的C++正向差分表

我已经使用静态数组创建了前向差分表，但无法通过使用动态数组解决它，请帮我解决此表。这是静态数组的代码，但我需要它与动态数组 #include<iostream> int main() { inti,noOfDifferenceColumns; //j is no of columnds. Initially it is 2. One for x, and other f