問題タブ [numerical]

説明 (これは hwk の質問です):

ここからどこから始めればよいかわかりません。ラプラス展開を使用する予定ですが、nxn 行列に実装する方法がわかりません。どんな助けでも大歓迎です。

注: nxn 行列のランダム行列を生成する関数が既にあります。また、計算のタイミングも問題ありません。私が問題を抱えている唯一のことは、行列式を計算する方法です。

クラス ポリシーの質問の説明 b/c を削除する必要がありました。

10進数のxが与えられた場合、xが分母が9999以下の有理数の10^-12以内にあるかどうかをテストしたいと思います。明らかに、x、2x、3xなどを調べて、これらのいずれかが整数に十分に近いかどうかを確認することで、それを行うことができます。しかし、より効率的なアルゴリズムはありますか？

私は数値コードでBoostのuBLASを使用しており、「重い」ソルバーを配置しています。

http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion

コードはうまく機能しますが、非常に遅いです。いくつかの調査の結果、スパース行列ソルバーであるUMFPACKを見つけました（とりわけ）。私のコードは、非常に頻繁に反転する必要がある大きなスパース行列を生成します（より正確には、逆行列の値は関係ありません）。したがって、UMFPACkとBOOSTのSparse_Matrixクラスは幸せな組み合わせのようです。

UMFPACKは、エントリ数、行インデックス、およびエントリの3つのベクトルで指定されたスパース行列を要求します。（例を参照してください）。

私の質問は要約すると、BOOSTのスパース行列クラスからこれらの3つのベクトルを効率的に取得できますか？

C++ で重力の数値シミュレーションを行っていますが、1 ステップがカウントされるたびに結果をバックアップしたいと考えています。

ただし、現在の方法では、プログラムは常にファイルを上書きします。テキストを常に別のファイルまたは変数に保存することで解決できると思いますが、上書きしないようにテキストファイルを開く簡単な方法があるのではないかと思います。

現在の「バックアップ コード」は次のようになります。

関数を提供するオープンソースの数値 C ライブラリを知っている人はいlogsumexpますか?

このlogsumexp(a)関数は、数値オーバーフローを回避して、配列 a のコンポーネントの指数 log(e^{a_1}+...e^{a_n}) の合計を計算します。

CUDA プログラムで、最近、次を使用して無限大のテストから切り替えました。

math.hINFINITYから

異なる結果が得られたことに非常に驚きました。gnu.orgは、実際には同様に動作する必要があることを示唆しています。何か不足していますか？INFINITYCUDA カーネルでの使用は許可されていませんか?

編集： 私はちょうど発見isinfし、使用してチェックしていることに気づきました

INFINITY チェックと同じ結果になります。なぜではないの isfinite(x)==!isinf(x)ですか？

optim を使用して、BFGS (および L-BFGS-B) で非線形最小二乗問題を当てはめようとしています。分析勾配を指定すると、ラインサーチが異常終了し、最終解は常に開始点に非常に近くなります。ただし、勾配を指定しないと、うまく収束するようです。これは誰かに数値的な問題を示唆していますか? グラデーションが正しいと確信しています。それはスケーリングの問題でしょうか？助けてくれてありがとう。

Java / Scala で使用する科学計算用の数学ライブラリを探しています。特に、完全楕円積分と修正ベッセル関数が必要です。オープンソースであればいいのですが、そこにあるものは何でも取り入れなければならないと思います。scipy (科学計算用の python lib) の交換は素晴らしいでしょう :-)

プログラムの速度に問題があります。巨大な配列の 4 つの近傍の平均を計算したいと考えています。これが私のコードの一部です。最後の行を変更する方法はありますか? または、別の配列を使用する必要がありますか?

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

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

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

心から。ジョージ