問題タブ [numerical]

この問題を本当に解決する必要があるので、助けてください。

何よりも、お時間をいただきありがとうございます...

私の問題: 行列 (826x826 double) があり、この行列を (826x1 double) のベクトルに関して統合したいのですが、これの機能はありません。ベクトルに関して行列の積分を取るコマンドまたはアルゴリズムはありますか? 本当に助けが必要です。私はmatlabの初心者です。

心から。ジョージ

「C の数値レシピ」という本を読み始めました...プログラムをより効率的かつ高速にしようとしています...つまり、C または C++ を使用することは同じことですか? C++ が C のスーパーセットであることは知っていますが、「cmath」ライブラリ (c++) と「math.h」ライブラリ (c) には違いがありますか? 実行速度に介入しますか?...そのために... Cプログラムと何の違いもなくC++を使用できますか?

matlab を使用して 3 体の問題をプログラムしようとしています。宇宙での回転座標系における月の軌道の式を教えてもらいました。基本的には ydotdot、xdotdot=GM/(x^2+y^2)^3/2 式です。式が何であるかは、それほど重要ではありません。

私が直面している問題は、月の軌道方程式を数値的に解くプログラムをコード化することになっているということです。私の目標はODE45と同じ結果を得ることなので、ODE45を使用して比較しています。私の究極の問題は、時間を日数で反復したいので、tspan = [0 365]. 薄いのは、重力定数を秒に変換してから tspace= [0 365] を実行すると、まったく異なる結果が得られます。[0 365*3600*34] を実行すると、1 年の秒と G= 6.67e を表します。 -11. 私のユニットは非常に奇妙に思えます。

ODE 45 を使用しているときになぜこのようなことが起こるのか、誰か説明できるかどうか疑問に思っていました。ODE45 を使用して秒を日に明確に変換できないのはなぜですか? 私がしなければならない余分なステップはありますか？私の問題の唯一の他の変数は、半径、距離、および 3 つの物体の質量です。

どうもありがとう。私はこれに非常に長い間取り組んできました。どんな助けでも大歓迎です。

次のアルゴリズムに従って gauss-legendre コードを作成しようとしています:

n ポイント

つまり、2n 方程式系が作成されます (次数 2n-1 の多項式に対して正確である必要がある場合、

ti は次数 n のLegendre 多項式の根です。Legendre 多項式は次のように与えられます。

と wi :

私のコードは次のとおりです。

私の問題は次のとおりです。

1) コンパイラは私に :error: 型 'double*' および 'double' の無効なオペランドをバイナリ 'operator*' に与えます --> at result=w*f(*t);

2）すべてを正しく行ったかどうかはわかりません。つまり、すべてを組み合わせて、アルゴリズムを正しく実装したかどうかです。

D プログラミング言語の並列処理ライブラリをデバッグしようとしています。タスクを使用して実行される一部の浮動小数点演算の下位ビットが実行間で非決定論的であることを示すバグレポートが最近提出されました。(レポートを読む場合は、決定論的な方法でタスクを作成することにより、内部で並列リデュースが機能していることに注意してください。)

丸めモードを手動で設定しようとしたため、これは丸めモードの問題ではないようです。また、これは並行性のバグではないと確信しています。ライブラリは十分にテストされており ( Jinxストレス テストに合格することを含む)、問題は常に下位ビットに限定されており、低レベルのメモリ モデルの問題がそれほど問題にならないシングルコア マシンでも発生します。操作がスケジュールされているスレッドによって浮動小数点の結果が異なる可能性がある他の理由は何ですか?

編集: ここでいくつかの printf デバッグを行っていますが、個々のタスクの結果が実行ごとに異なる場合があるようです。

編集 # 2: 次のコードは、この問題をより簡単な方法で再現します。メインスレッドで配列の項を合計し、新しいスレッドを開始してまったく同じ関数を実行します。このコードは私のライブラリを使用していないため、問題は間違いなく私のライブラリのバグではありません。

出力：

丸めモード: 0

0.7853986633972191094

丸めモード: 0

0.7853986633972437348

別の編集

代わりに16進数で印刷したときの出力は次のとおりです。

丸めモード: 0

0x1.921fc60b39f1331cp-1

丸めモード: 0

0x1.921fc60b39ff1p-1

また、これは Windows でのみ発生するようです。このコードを Linux VM で実行すると、両方のスレッドで同じ結果が得られます。

ANSWER : 根本的な原因は、浮動小数点の状態が D の Windows の他のスレッドとは異なる方法でメイン スレッドで初期化されることであることが判明しました。私が提出したばかりのバグ レポートを参照してください。

インプレース LU 行列分解のために、C の Numerical Recipes の提供されたソース コードから文字どおりコピー アンド ペーストしましたが、問題は機能していません。

私は愚かなことをしていると確信していますが、誰かが私を正しい方向に向けることができれば幸いです。私は一日中取り組んできましたが、何が間違っているのかわかりません。

回答後の更新: プロジェクトは終了し、機能しています。皆様のご指導に感謝いたします。

WolframAlphaは、答えは

私は得ています：

これまでのところ、「同じ」アルゴリズムの少なくとも 3 つの異なるバージョンを見つけたので、完全に混乱しています。

PSはい、これを行うには少なくとも12の異なるライブラリがあることは知っていますが、正しい答えよりも自分が間違っていることを理解することに興味があります。

PPS LU 分解では、下の結果の行列は 1 であり、Crouts アルゴリズムを (私が思うに) 実装されているように使用すると、配列インデックスへのアクセスは依然として安全であり、L と U の両方をインプレースで重ね合わせることができます。したがって、これに対する単一の結果行列です。

私は、一般的な .NET アプリケーションを開発するという観点から Sho を検討していますが、.NET プラットフォーム用の標準の数学ライブラリがないということは、私にとって過去に難しかったことです。商用実装とは別に、現在 2 つの .NET ライブラリがその深さと幅の広さで際立っています。

1) Sho (クローズド ソース?、より新しいが、MS .NET フレームワークのコアに含まれるようになる予定)

2) Math.NET 数値 (オープンソースですが、「ベータ版」では、それ自体がかなり成熟したプロジェクトであり、MSoft の他のメンバーがプロジェクトに取り組んでいます)

類似点が多すぎるように思えるので、実際に協力を開始する必要があるのではないかと考えています。

どちらもMKL、ACML、LAPACKなどのネイティブ高速ライブラリへの接続を許可するため、パフォーマンスを考慮して上記の2つに選択を制限しました...

他の誰かが同様のことを調べていますか?

ありがとう

サム

私はC＃と.NETを使用していますが、次のような正規表現があります

これはすべての非数値文字に一致しますが、小数点（。）と負の符号（-）に一致させたくありません。正規表現でそれを行うにはどうすればよいですか？

だから私はクリスを試してみました、そしてそれが機能するようにいくつかの調整をしました：

（「オリジナル」という名前のTextBoxがあります）

これOriginal.SelectionStart = Original.TextLength;は、置き換えられるたびに選択範囲が最初に配置され、ユーザーには少し奇妙に見えるためです...

I have differential equations derived from epidemic spreading. I want to obtain the numerical solutions. Here's the equations,

t is a independent variable and ranges from [0,100]. The initial value is

At first, I planned to deal these with ode45 function in matlab, however, I don't know how to express the series and the combination. So I'm asking for help here.

**

The problem is how to express the right side of the equations as the odefun, which is a parameter in the ode45 function.

**

このことの呼び方も説明の仕方も知らないので、タイトルが少し誤解を招くかもしれません。

最初に添付されたグラフは pyplot で作成されました。現在使用している 3 つの赤い点の代わりに、すべてのグラフを通る直線を描きたいと思います。pyplotで可能ですか？2番目の画像は私が探しているものです。