openmp_range_algebra を使用して、boost::odeint で制御されたステッパーを構築しようとしています
typedef vector< complex< double > > state_type;
typedef runge_kutta_dopri5< state_type > error_stepper_type;
typedef controlled_runge_kutta< error_stepper_type > controlled_stepper_type;
controlled_stepper_type controlled_stepper(default_error_checker< double, openmp_range_algebra >;
ただし、odeint にはそのようなコンストラクターが存在しないため、コードはコンパイルされません。
私の質問: OpenMP で使用できるように、制御された runge_kutte-dopri5 ステッパーを作成するにはどうすればよいですか?
状態ベクトルが長い (長さ: 2^20) ため、プログラムの中で最も時間がかかる部分であるため、アダプティブ ステッパーを並列化したいのです。
お手伝いいただきありがとうございます