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