問題タブ [numerical-computing]

IndexedSeq[Double]ドメイン固有の数値計算ライブラリを設計するときに使用する適切なデータ型 ( など) を探します。この質問では、範囲を の 1 次元配列の操作に限定していますDouble。ライブラリは、1D 配列の各要素に通常適用される数値関数を定義します。

考慮事項:

• Vectorまたはなどの不変のデータ型を優先するIndexedSeq
• データ変換を最小限に抑えたい
• スペースと時間の合理的な効率
• 図書館を利用する人にやさしい
• エレガントでクリーンな API

など、コレクション階層の上位にあるものを使用する必要がありますSeqか?

それとも、単一要素の関数を定義して、マッピング/反復をエンド ユーザーに任せた方がよいでしょうか?

これは効率が悪いように見えますが (一部の計算は一連の呼び出しごとに 1 回実行される可能性があるため)、同時に、より柔軟な API です。これは、どのタイプのコレクションでも機能するためです。

推奨事項はありますか？

与えられたベクトルからゼロ平均ベクトルを取得することがわかっているので、このベクトルの各メンバーから与えられたベクトルの平均を減算する必要があります。たとえば、次の例を見てみましょう

我々が得る

次の操作で別のベクトル s を作成しましょう

この後、私たちは得る

次の式でsの平均を計算すると

ans =を得る

実際、私が確認したように、この数は非常に小さいです

したがって、ベクトルの平均はゼロであると考える必要があります。つまり、0 からのわずかな偏差は、丸め誤差によるものです。たとえば、ホワイト ノイズなどの無作為な無相関のデータ シーケンスを作成する場合、実際にはすでに想定されています。 0から遠く離れたこのような小さなデータでも、平均はゼロであり、この場合、たとえばこの場合

ほぼ当然

入力が数値であるかどうかを評価しようとしているだけisNaNで、最善の方法であると考えました。ただし、これは信頼できない結果を引き起こします。たとえば、次の方法を使用します。

これらの値について:

このフィドルに示されています：http://jsfiddle.net/4nm7r/1

なぜisNaNいつもうまくいかないのですか？

純粋関数の合成である 2 つの関数があります。最初の関数は区画を受け取り、その上に家を建て、雑誌で宣伝するために写真を撮ります。

2 番目の関数も区画を受け取り、その上に家を建て、仕上げを加えます。

2 つの関数も純粋関数の合成であることは明らかです。今、区画があり、niceParcel両方の機能をそれに適用したいとします。ただし、最初の 3 つのサブ関数は計算に時間がかかり、2 つの関数間で共有されるため、2 回計算されることは避けたいと思います。

コードをリファクタリングして、これらの不要な計算を回避し、明確な意味を持つこれらの素晴らしい純粋な関数を維持するにはどうすればよいですか?

コマンド を使用して次のコードをコンパイルすると、g++ -std=c++11 Threads.cpp -lpthread -I../Libs/nr30nr30 はhttp://www.nr.com/が提供するライブラリであり、エラーは発生しません。

#include "nr3.h"上記と同じコードにヘッダーを追加すると、次のエラーが発生します。

同じ問題に遭遇しましたか？アドバイスありがとうございます。

R : で非正方線形システムを解く方法はA X = B?

(系に解がないか、解が無限にある場合)

例 ：

Pythonで書いた数式を次のように適用する必要があります。

さて、これは正しいですが、明らかに非効率的です。逆に、次のようにします。

ただし、2 番目の例では、PHi2 のすべての要素で同じ数が得られるため、これは高速ですが、答えは間違っています。どうすればこれをより正確かつ効率的に行うことができますか?

注: eig1 と eig2 は同じ次元 d であり、theta と PHi2 は同じ次元 D ですが、d!=D です。