odeint

0熱

1答えて

カスタム状態タイプを実装するためにvector_space_norm_infを適切にオーバーロードするためにコードを変更する必要がありますか？ template<typename T> class Traits{ public: typedef T value_t; typedef complex<value_t> scalar_t; typedef vector

-1熱

1答えて

Xcode 7.3.1でOdeIntからの暗黙ソルバを正しく使用できない

いくつかの硬い常微分方程式を解く必要があるので、私はOdeintが提供する暗黙の方法について知らせました。しかし、Odeint、 http://headmyshoulder.github.io/odeint-v2/doc/boost_numeric_odeint/tutorial/stiff_systems.html のドキュメントの例のXcodeをコピーしてコンパイルすると、の暗黙的なrose

0熱

2答えて

python - 補間軸に沿って異なる配列の長さですか？

与えられたxの値yを得るためにPythonの補間関数を使用しようとしていますが、エラー "raise ValueError（" xとy配列は補間軸に沿って同じ長さでなければなりません "私のコードで.shapeを使用するときに得られるものに応じて同じサイズと形状の両方があります。プログラミングには新しかったので、配列内で他に何が違うかチェックする方法はわかりません。 s = [] def sl

0熱

1答えて

odeintブーストライブラリでのdopri5の使用

system of equations こんにちは。私は0から10^16までの時間でそれらの方程式を進化させたいと思います。そして、最初の条件はx（0）= 10^8とy（0）= 0.5です。方程式が分母にxに依存するので、runge_kutta_dopri5でodeintを使うのは適応的なステップ制御のために良い選択です。事は私が実際にこれを行う方法を少し考えているので、私はC++とodeintで

0熱

1答えて

最終的なodeint出力のみを出力する

これは結構な疑問に思えるかもしれませんが、最終的な出力値だけを出力することが可能かどうかを尋ねる必要があります。実際には、無作為に生成された時間間隔の2つの結合微分方程式を解いて、各区間ごとに最終出力のみを出力しようとしています。

0熱

1答えて

odeintでopenmpを使用する

odeintとopenmpを使用してopenmpを使用してコードを並列化しようとしていますが、スレッドの数を変更すると並列性がうまくいかず、実行を確定する時間が改善されません。私が間違っていることは何ですか？ここでは、コードの重要な部分です： using namespace boost::numeric::odeint; using namespace std; typedef std::v

0熱

1答えて

C++

に時間を逆方向にして常微分方程式のシステムを解くことは残念ながら、私はこれは実際には何もしません typedef std::vector<double> state_type; void ode_function(const state_type &x, state_type &dxdt, const double { dxdt[0] = x[0]; } using namespa

0熱

1答えて

Python odeint-初期条件y0は1次元でなければならない

odeintを使って状態空間形式で動的システムをシミュレートしようとしています。私のA行列は12 * 12で、B行列は12 * 4（行* cols）なので、私の初期状態ベクトルは12 * 1です。私は状態ベクトルを定義すると考えることができ Traceback (most recent call last): File "sim.py", line 17, in <module>

1熱

2答えて

ルックアップテーブルからの時間可変パラメータをboost :: odeint、C++に組み込む方法

boost :: odeintを使用して非線形システムを数値的に統合しようとしています。システムには外部から生成される時変パラメータがあり、これをプログラムに組み込みたいと考えています。これはodeintで可能ですか？ Matlabでは、同様のことをする場合は、値が利用可能になると補間する必要があります。ありがとうございました！

1熱

1答えて

boost/odeintのstatetypeとしていくつかの固有行列を使用

私はODEに従ってパラメータを展開しなければならない物理問題に取り組んでいます。対角化などのルーチンで使用できるデータ型を持つように、時にはそれらを操作する必要があります。したがって、メンバーとして固有クラス::行列を実装し、統合を実行したいodeintと一緒に。単一の固有値::行列に対して、これはうまくいきました。、私が知っている（Macの場合はG ++、基本的に、これは私が「の手順をint型