numerical-computing

    5熱度

    1回答

    我正在開發一個需要浮點確定性的C應用程序。我也希望浮點運算速度相當快。這包括IEEE754沒有規定的標準超越函數,如正弦和對數。與硬件浮點相比,我考慮的軟件浮點實現相對較慢,所以我正在考慮簡單地將每個答案中的一個或兩個最低有效位四捨五入。精確度的損失對我的應用來說是一個妥協的妥協,但這足以確保跨平臺的確定性結果嗎?所有浮點值都將是雙精度值。 我意識到操作順序是在浮點結果差異的另一個潛在來源。我有辦

    4熱度

    1回答

    我有一個問題,它是由不完整的伽瑪函數和逆不完整的伽瑪函數構成的。回想一下,在Gamma函數是從0到無窮大的特定積分的情況下,不完全伽馬函數會添加另一個參數x,並從0到x(較低的不完整)以及從x到無窮大(上部不完整)積分。我知道Matlab的功能涵蓋了所有這些功能:完整的Gamma函數,上下不完整的Gamma函數以及逆函數。然而,Matlab將其不完整的Gamma函數實現爲正則化爲不完整的伽馬函數,

    1熱度

    3回答

    我需要將一個數字範圍劃分爲一些長度相同的段。但我無法決定哪種方式更準確。例如: double r1 = 100.0, r2 = 1000.0, r = r2 - r1; int n = 30; double[] position = new double[n]; for (int i = 0; i < n; i++) { position[i] = r1 + (double)i/

    1熱度

    1回答

    f(x) = (exp(x)-1)/x; g(x) = (exp(x)-1)/log(exp(x)) 從分析角度而言,f(x) = g(x)適用於所有的x。 當x接近0,無論f(x)和g(x)方法1. % Compute y against x for k = 1:15 x(k) = 10^(-k); f(k) =(exp(x(k))-1)/x(k); De(

    5熱度

    1回答

    我有兩個函數是純函數的組成部分。 第一個函數需要一個包裹,建立在其上的房屋,並拍攝一張照片來宣傳它在一本雜誌: let buildAndAdvertiseHouse parcel = parcel |> inspect |> buildWalls |> buildRoof |> takePhoto |> advertise 第二個功能

    0熱度

    1回答

    如果我編譯與命令g++ -std=c++11 Threads.cpp -lpthread -I../Libs/nr30,其中NR30是由http://www.nr.com/提供的庫將以下代碼,我得到沒有錯誤: #include <chrono> #include <thread> int main(void) { /* ... Sadly, No further code in h

    2熱度

    2回答

    如何計算Mathematica中數值的二重積分? Integrate[Exp[-0.099308 s] * Integrate[Exp[0.041657423 u] Exp[-3.1413 s + 3.12 u] * ((u/(s - u))^(1/2) BesselI[1,2 (u (s - u))^(1/2)] + 0.293 BesselI[0,2 (

    15熱度

    2回答

    我在Clojure中實現了一些基本的複數運算,並且發現它比大致相當的Java代碼慢10倍左右,即使是類型提示也是如此。 比較: (defn plus [[^double x1 ^double y1] [^double x2 ^double y2]] [(+ x1 x2) (+ y1 y2)]) (defn times [[^double x1 ^double y1] [^double

    1熱度

    2回答

    說我有一個巨大的numpy矩陣A佔用數十千兆字節。分配此內存需要花費不可忽視的時間。 假設我也有一個與numpy矩陣具有相同尺寸的scipy稀疏矩陣的集合。有時我想將這些稀疏矩陣之一轉換爲稠密矩陣來執行一些向量化操作。 我可以將這些稀疏矩陣之一加載到A而不是每次我想將稀疏矩陣轉換爲稠密矩陣時重新分配空間嗎?可用於scipy稀疏矩陣的.toarray()方法似乎沒有采用可選的密集數組參數,但也許有其

    4熱度

    2回答

    在設計特定於領域的數值計算庫時尋找適當的數據類型(例如IndexedSeq[Double])。對於這個問題,我將範圍限制爲使用1維數組Double。該庫將定義通常應用於1D陣列中每個元素的數字函數。 注意事項: 不想一成不變的數據類型,如Vector或IndexedSeq 要儘量減少數據轉換 在空間和時間上合理有效 友好使用其他人圖書館 高雅清潔的API 我應該使用更高一些的收集層次結構,如Seq