boost odeint中的受控误差步进器支持复杂数据类型吗？
在boost odeint中，受控误差步进器（controlled error stepper）是一种用于求解常微分方程（ODE）的数值方法。它可以根据给定的误差容限自适应地调整步长，以保证数值解的精度。
受控误差步进器支持复杂数据类型。在boost odeint中，可以使用自定义的数据类型作为ODE的状态变量。这意味着您可以使用任何复杂的数据类型，包括自定义类、结构体或复合类型作为ODE的状态变量。
使用复杂数据类型的优势在于可以更好地表示和处理实际问题中的数据。例如，对于涉及多个物理量的问题，可以使用自定义类来表示这些物理量的集合，并定义相应的运算规则和数值方法。这样可以更方便地进行问题建模和求解。
boost odeint提供了丰富的功能和灵活的接口，以支持各种类型的ODE求解。您可以根据具体的问题需求选择合适的受控误差步进器和数据类型，并根据需要进行定制和扩展。
腾讯云提供了一系列与云计算相关的产品和服务，包括云服务器、云数据库、云存储等。这些产品可以帮助用户快速搭建和部署云计算环境，并提供高可用性、高性能和高安全性的解决方案。您可以访问腾讯云官方网站（https://cloud.tencent.com/）了解更多关于腾讯云的产品和服务信息。
相关·内容
我在odeint库中遇到了一个问题，它使用了一个受控的错误步进器和一个复杂的状态类型。示例中的代码带有复杂的stuart landau方程，我对其进行了修改，使其包含一个自适应积分器。代码现在看起来像这样： #include <iostream>
#include <complex>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
//[ s
浏览 16提问于2019-10-09得票数 4
回答已采纳
我想使用boost::sandbox中的odeint模块，但我不确定如何将从svn检索到的沙箱代码集成到我的boost安装中。我已经将boost安装到$BOOST_ROOT=/usr/local/boost_1_45_0。将odeint的内容合并到我的boost路径中是否合适，该路径会将头文件放在$BOOST_ROOT/boost/numeric/odeint.hpp中。
这是正确的方法吗?还是有一种标准的方法来获取和安装boost::sandbox代码？
浏览 0提问于2010-12-29得票数 4
回答已采纳
1回答
我是C++世界的新手，在使用boost库时遇到了一些麻烦。在我的问题中，我想要解决一个有5个方程的常微分方程系统。这不是一个僵硬的问题。作为迭代方法，我同时使用了integreate(rhs, x0, t0, tf, size_step, write_output)和integreate_adaptive(stepper, sys, x0, t0, tf, size_step, write_output)。这两种方法实际上都集成了方程，但给出了毫无意义的结果，几乎随机地将步长从0.001更改为5。方程和数据是正确的。我能做些什么来解决这个问题？代码如下： #include <iostr
浏览 42提问于2020-11-04得票数 0
回答已采纳
我最近偶然发现了boost.odeint库，我对它的可能性和可配置性感到惊讶。但是，由于广泛使用了scipy.integrate.odeint (它本质上是fortran中ODEPACK的包装器)，我想知道它们的性能如何比较。我知道boost.odeint还附带了并行化，这在scipy中是不可能的(据我所知)，这将大大提高性能，但我要求的是单一核心情况。然而，由于在这种情况下，我必须将boost.odeint (使用cython或boost.python)包装到python中，也许你们中已经有人这样做了？这将是一个伟大的成就，因为所有的分析可能性在python中都要高级得多。
浏览 0提问于2016-02-04得票数 5
1回答
我尝试在Mac 10.9.2 boost_1_55_0 5.1上运行[odeint复杂状态类型示例代码]。
下面的代码是解决斯图亚特-兰道振荡器的网站上的副本。
#include <iostream>
#include <complex>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
//[ stuart_landau_system_func
浏览 0提问于2014-04-04得票数 1
回答已采纳
1回答
我想用boost::numeric::odeint::runge_kutta4.求解一个耦合微分方程组这是一个三维的晶格系统，所以我很自然(也很方便)使用3D阵列。runge_kutta4是否有一种方法可以使用用户定义的数据结构或增强multi_array？
浏览 2提问于2015-10-03得票数 3
1回答
我正在尝试创建一个简单的程序，以熟悉推进器的GPU计算能力和odeint的ODE解决能力。我希望能够使用GPU上的Runge方法来求解简单的ODE(即dy/dx = 3x^2y)，希望以后能够继续处理更复杂的问题。我能够相对轻松地使用odeint来完成这一任务：
#include <boost/lambda/lambda.hpp>
#include <boost/numeric/odeint.hpp>
#include <iostream>
#include <iterator>
#include <algorithm>
using
浏览 3提问于2016-04-10得票数 1
我正在使用boost odeint来进化一个微分方程。时间和值是由包装类围绕一个普通数(进行高斯误差传播)给出的。
using Time = Number<double>;
using Value = Number<double>;
状态及其导数由带有声明的自定义类给出。
class State : boost::addable<State>, boost::multipliable<State, Time>, boost::dividable<State>
{
public:
void Add(Group const&
浏览 2提问于2017-06-15得票数 0
回答已采纳
我正在为新版本的boost而苦苦挣扎。我使用多精度的odeint。以下代码可以使用boost版本1.67.0成功编译。但是，从1.68.0版本开始，我就不能再编译了。在1.68.0版本中引入了复杂的多精度数字，我希望在我的代码中也使用这种数据类型。有人能告诉我代码中的错误在哪里吗？ #include <iostream>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
#include <boost/multiprecision/gmp.hpp>
using na
浏览 40提问于2019-10-10得票数 2
回答已采纳
我正在尝试安装odeint，但是每次编译代码时，我都会得到这个错误。
c:...odeint\headmyshoulder-odeint-v2-f496df3\include\boost\numeric\odeint\config.hpp(44): fatal error C1083: Cannot open include file: 'boost/config.hpp': No such file or directory
我已经包括了包含路径。我检查了这个文件boost/config.hpp，这不是<config.hpp>的完整路径，所以我将它更改为<b
浏览 9提问于2014-10-31得票数 1
回答已采纳
我正在尝试将boost：：combine：：odeint与我自己的System类的实现结合起来(参见System.hpp)。
在BatchFilter类方法中使用(模板)系统对象，如下所示：
# BatchFilter.cpp
# template <typename state_type> BatchFilter class {...}
System<state_type> dynamics(this -> args,
N_true,
this -> true_dynamics_fun );
浏览 2提问于2017-12-16得票数 0
回答已采纳
3回答
调用另一个函数重载
、、、、
我要从odeint库中了解c++ boost，我需要知道哪个部分做什么。在boost/numeric/odeint/integrate/integrate_adaptive.hpp中，有一个叫做integrate_adaptive的函数。这个函数有一些过载。使用我的几个操作简化的文件如下：
integrate_adaptive_minimal.hpp
#define BOOST_NUMERIC_ODEINT_INTEGRATE_INTEGRATE_ADAPTIVE_HPP_INCLUDED
#include <boost/type_traits/is_same.hpp>
#in
浏览 4提问于2015-01-13得票数 4
回答已采纳
2回答
我在编程方面真的是经验不足。我在安装odeint (和boost)来解微分方程系统时遇到了麻烦。我已经看了很多关于安装boost的问题，我尝试了很多，但我还不能使用odeint。我正在使用Mac Yosemite 10.10.5，下面是我的问题，非常感谢你的帮助：
1，如果我可以成功安装boost，那么当我想使用odeint时，我是否只需要将#include<boost/odeint/optionalSubdirectories>作为我的文件头？我弄糊涂了，因为在上，有可供下载的文件，其中包括示例。我想知道在我下载boost之后，我是否必须下载这些文件才能使odeint工作？如果是
浏览 1提问于2016-06-23得票数 1
我正在运行Pop!_OS (ubuntu导数)和apt从默认存储库安装的boost (libboost-all dev)。我知道它安装得很好，因为我可以用GCC编译并运行简单的boost odeint示例。
然而，当我试图在jupyter笔记本中使用zeus-cling内核运行相同的示例时，我在包含odeint标题时出错了。我可以在执行以下代码时重新创建错误：
#pragma cling add_include_path("/usr/include")
#include <boost/numeric/odeint.hpp>
我得到的错误消息是：
In file in
浏览 7提问于2020-04-14得票数 1
回答已采纳
我正在尝试创建一个boost odeint解析器。我可以从方程中得到双数的解，但我需要一个在(3D)点上工作的解(我用它作为开普勒轨道解算器的学习曲线)，我得到：
没有可行的重载'-=‘
编译器中没有可行的重载错误指向operators.hpp中的第261行，即：
BOOST_BINARY_OPERATOR_NON_COMMUTATIVE( subtractable, - )
导致此错误的最小代码如下。注意，Point3D取自odeind中的lorenz_point代码。
#include <boost/operators.hpp>
#include <boost/
浏览 7提问于2020-08-14得票数 0
1回答
我试着用boost中的odeint解算器，但是我一开始就被一个愚蠢的错误困住了，我似乎无法修复。我将boost文件夹复制到我的项目文件夹中，并使用
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
从一个例子-代码，已经不工作(包括-文件找不到)。将< and >更改为引号(")有效(visual studio查找文件)。但是“使用命名空间的提升：：数
浏览 6提问于2022-01-24得票数 0
我使用brew install boost安装boost是为了使用odeint库( 说: odeint是一个仅用于头文件的库，不需要链接到预编译的代码)。
我在Mac Yosemite 10.10.5上。现在，当我cd到/usr/local/include时，我可以在那里看到boost目录。在boost (/usr/local/include/boost)中，有我的项目所需的所有头文件，以及我的特定用途所需的数字/odeint目录。同时，当我用cd连接到usr\local\lib时，我可以看到很多libboost_* .dylib和.a文件。
但是，当我尝试编译一个临时保存在~/Downloa
浏览 0提问于2016-06-23得票数 0
3回答
我正在尝试在C++中解决一组耦合的ODE，我已经决定尝试odeint ()，因为我已经读到了很多关于它的好东西。然而，我正在编译包中包含的示例。情况是这样的：
到目前为止，我只下载并解压了这个包。我在Windows XP上运行Visual C++ 6，并且正在编译示例"simple1d“。我收到错误消息：“无法打开包含文件：'boost/config.hpp'：没有这样的文件或目录”。
我不清楚我错过了什么。我需要完整的boost库才能工作吗？我的印象是ODEINT是一个“独立”的包。
浏览 2提问于2012-10-15得票数 3
回答已采纳
我使用的是scipy.integrate包中的odeint函数：
r0 = np.array([1,2,3,4])
t=np.linspace(0,1,20)
def drdt(r,t):
return r # or whatever else
r = odeint(drdt,r0,t)
r0是一个numpy数组，它包含一定数量的点的初始位置。在脚本的结尾，不出所料，我得到了20个时间步长上的点的位置。
现在，我想在满足r上的条件时停止odeint求解器。特别是，当其中两个点接近某个阈值时，我希望停止odeint，对r向量进行一些更改，并使用新的初始位置继续odeint求解器。有没有办
浏览 1提问于2014-12-03得票数 1
在R中使用Rcpp编译我的.cpp文件时，会出现以下错误消息：
对`boost：：system：：generic_category()的未定义引用
但是当我删除// [[Rcpp::plugins(cpp11)]]行时，就没有任何错误了。为什么？
下面是我的最小可重现性例子。
// include Rcpp, it takes care of most other headers you need
#include <Rcpp.h>
#include <boost/array.hpp>
// include Boost's odeint
#includ
浏览 0提问于2018-05-22得票数 0
回答已采纳
我正在尝试使用boost::odeint对一个非线性系统进行数值积分。该系统具有外部生成的时变参数，我希望将其合并到我的程序中。使用odeint可以做到这一点吗？在Matlab中，如果要执行类似的操作，则需要在值可用时对其进行插值。
提前感谢您的帮助！
浏览 16提问于2016-07-29得票数 1
回答已采纳
我想在runge_kutta4中使用方法。我已经用Matlab解决了这个问题。下面是我在Matlab中用初值x'' = -x - g*x' ( x1 = 1，x2 = 0 )求解的代码如下
main.m
clear all
clc
t = 0:0.1:10;
x0 = [1; 0];
[t, x] = ode45('ODESolver', t, x0);
plot(t, x(:,1));
title('Position');
xlabel('time (sec)');
ylabel('x(t)');
ODE
浏览 12提问于2014-11-05得票数 4
回答已采纳
2回答
为了了解C++是如何工作的，为了解决微分方程，我尝试了在互联网教程中找到的以下脚本。应该管用的。
我不直接使用C++。我使用R和Rcpp库来编译/运行脚本。下面是脚本：
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace Rcpp;
using namespace std;
using namespace boost::numeric::odeint;
/* we solve the simple ODE x' = 3/(2t^2) + x/(2t)
* with
浏览 0提问于2018-04-24得票数 0
回答已采纳
如果我使用一个类来定义导数，将一个系统集成到odeint模块中，那么这个类的析构函数就会经常被调用。
是为了这个行为吗？为什么是这样?如果我想在这个类中动态分配数组，我应该做什么？
例如，这段代码：
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
class foo
{
public:
virtual ~foo() {
std::cout << &
浏览 3提问于2019-10-24得票数 3
回答已采纳
我想用boost::odeint并行地求解不同变量集合的微分方程--比如说std::vector类型。当然，一种解决方案是将所有变量组合成一个大的向量，然后将其用作状态变量。
然而，我更喜欢一种更优雅的解决方案，比如使用boost::fusion作为状态类型，然后保存不同的向量。据我所知，从和的实现来看，存在-in原则-这样做没有障碍。我只遗漏了一些关于具体实现的提示--特别是关于正确的
代数、运算和调整大小
例如，创建错误步进器时需要。哪些现有的实现-例如odeint::fusion_algebra-可以直接使用，在这种情况下还需要做些什么？
浏览 0提问于2013-03-11得票数 2
回答已采纳
我正在尝试最小化两个有界函数之间的均方误差，curve_fit做得很好，但当两个函数之间的均方误差小于0.1时，我想停止计算。下面是一个简单的示例代码 import numpy as np
from scipy import optimize, integrate
def sir_model(y, x, beta, gamma):
sus = -beta * y[0] * y[1] / N
rec = gamma * y[1]
inf = -(sus + rec)
return sus, inf, rec
def fit_odeint(x, beta, g
浏览 65提问于2021-07-06得票数 1
1回答
我尝试在Debian挤压上运行boost_1_54_0中的g++4.4。
Lorenz系统工作良好，但简单的一维ode
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
/* we solve the simple ODE x' = 3/(2t^2) + x/(2t)
* with initial condition x(1) = 0.
* Analytic solut
浏览 3提问于2013-09-15得票数 3
回答已采纳
我正在测试Eigen::VectorXd作为boost::odeint的state_type。我使用以下代码：
#include <Eigen/Eigen>
#include <boost/numeric/odeint.hpp>
#include <boost/numeric/odeint/external/eigen/eigen_algebra.hpp>
#include <iostream>
#include <vector>
template<class T>
struct push_back_state_and
浏览 0提问于2020-06-06得票数 0
2回答
我要解几个矩阵微分方程，形式是d/dt (X) = F(X)，其中X是一个很大的复矩阵，F表示它的一些函数。我尝试使用Boost的odeint和state_type作为Armadillo的cx_mat。但它对受控步进类型会产生编译错误。我的示例代码如下
#include <armadillo>
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace arma;
using namespace boost::numeric::od
浏览 1提问于2017-07-08得票数 1
回答已采纳
下面概述了我正在使用的类步骤的类型签名：
它可以实例化如下：
boost::numeric::odeint::runge_kutta_dopri5< state_type_ > stepper;
到目前一切尚好。它起作用了。
我计划将我的程序移植到cuda (使用推力)，然后再移植到openmp。我将声明改为：
boost::numeric::odeint::runge_kutta_dopri5< state_type_
, double
, state_type_
, double
, boost::num
浏览 6提问于2016-04-12得票数 1
回答已采纳
我打算使用MCMC例程中的Boost odeint库来估计ODE模型中的参数。由于这些颂歌可能是僵硬的，所以我需要用导数将雅可比传递到求解器中。我想做一个类，它的参数和初始值是私有成员，然后是导数，jacobian，以及将参数更改为公共方法的方法。我试图修改odeint网站中的刚性示例，使其使用一个包含两者的类，但是在进行调用时，会收到一个错误' error : no匹配函数--调用'Fitzhugh::deriv()‘。我不是一个经验丰富的C++程序员，所以这可能是一个概念上的错误。这是密码。
/* Fitzhugh Nagumo Equation in odeint */
浏览 3提问于2013-09-21得票数 2
回答已采纳
2回答
带odeint的两点边界
、、、、
我正在尝试用odeint来解决两点边界问题。我的方程式的形式是
y'' + a*y' + b*y + c = 0
当我有y(x_1) = y_1 , y'(x_2) = y_2的边界条件时，这是相当微不足道的，但是当边界条件为y(x_1) = y_1 , y(x_2) = y_2时，我就迷失了。有没有人知道用odeint或其他科学图书馆解决这样的问题的方法？
浏览 0提问于2013-07-19得票数 3
回答已采纳
我正在使用Ubuntu12.04&已经在/usr/include中使用了一些增强的fies。我做了一个
sudo apt-get install libboost-all-dev
也安装了很多文件。我不想删除这个boost并从源代码中安装，因为其他几个包依赖于ubuntu的版本。这是我想运行的示例代码:-
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
typedef v
浏览 4提问于2013-11-09得票数 0
回答已采纳
为了评估R和C++在求解ODE方面的速度差异，我在R中构建了以下ODE系统： modelsir_cpp =function(t,x){
S = x[1]
I1 = x[2]
I2 = x[3]
N=S+I1+I2
with(as.list(parm), {
dS=B*I1-mu*S-beta*(S*(I1+I2)/N)
dI1=beta*(S*(I1+I2)/N)-B*I1-lambda12*I1
dI2=lambda12*I1
res=c(dS,dI1,dI2)
return(res)
})
} 为了解决这个问题，我使用了deSolve包。 times =
浏览 42提问于2018-04-26得票数 1
回答已采纳
1回答
假设我有以下boost::odeint代码：
#include <iostream>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
const double sigma = 10.0;
const double R = 28.0;
const double b = 8.0 / 3.0;
typedef boost::array< double , 3
浏览 0提问于2013-01-22得票数 4
1回答
我尝试使用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
浏览 5提问于2014-12-28得票数 4
我想检查位精度如何影响DNN的准确性。你知道有什么C/C++/Python库不需要大量的返工来支持量化类型的推断吗？例如，我想将所有操作数(权重、激活、误差和梯度)设置为8/16/32位精度，并在测试后检查每个精度的准确性。
浏览 20提问于2019-11-14得票数 1
回答已采纳
如果您有任何使用本机Android的经验，您就会知道每个EditText本身都保存了它的价值。我的意思是，当用户输入任何内容时，您必须更新state，但这不是很糟糕吗？
在android中，您可以访问Edittext，然后调用一些方法，如getText()或某些东西，但是在react本机中，您必须读取状态中的值！这意味着更多的耦合！因为您正在为某个组件使用状态，例如，在另一个组件中。
出于什么原因，本机的反应是使用这种方法，这似乎是旧的！在较老的C#中，你有类似的方法。
浏览 2提问于2019-05-04得票数 1
回答已采纳
1回答
我使用boost odeint来积分2n耦合方程。我的基本问题很笼统。如何使用观测器不仅返回每个时间步骤的状态x_i，而且还返回导数d(x_i)/dt？我知道如何在给定的时间返回状态或状态的某些函数，但我不知道如何返回某个给定时间和时间内状态的函数。
更确切地说，这里是我的代码的简化版本。
#include <boost/numeric/odeint.hpp>
#include <boost/numeric/ublas/matrix.hpp>
typedef boost::numeric::ublas::matrix<double> matrixd;
usi
浏览 1提问于2021-02-19得票数 1
回答已采纳
1回答
我是个学生..。我试图将查询保存到文档中的一个文件中，但我得到了一条错误消息。
找不到c:\user\my_account\documents\script.sql文件
请检查文件名，然后重试。
其他问题：
生成脚本失败
无法从word_db获取对象列表
结果误差
浏览 1提问于2018-08-09得票数 2
我正在尝试使用boost odeint集成一个非常简单的ODE。在某些情况下，值与Python的函数相同(或非常类似)。但对于其他初始条件，数值则大不相同。
函数是: d( uhat ) / dt =- alpha^2 * kappa^2 *uhat，其中alpha为1.0，kappa是一个常量，视情况而定(见下面的值)。
我尝试过几种不同的boost算法，但似乎没有一种方法有效。
更新:下面的代码正在工作。
在下面的代码中，第一种情况给出了几乎相同的结果，第二种情况有点琐碎(但令人放心)，第三种情况在C++版本中给出了错误的答案。
以下是C++版本：
#include <boost/n
浏览 4提问于2021-10-13得票数 2
回答已采纳
2回答
如何安装第三方库
、、、
我对C++很陌生，我决定尝试使用odeint进行一些模拟，因为python对于我的需要来说太慢了。
我找到了，我想和它一起玩。我只是不太清楚如何安装这些库，或者把这些库放在哪里。C++有类似于python的pip install的东西吗？
附带注意:我正在尝试使用Eclipse Kepler来完成这个任务，但我并没有接受这个想法。
浏览 8提问于2014-02-21得票数 7
回答已采纳
我用来解微分方程。在visual studio 2010中，没有出现错误，但当我使用visual studio 2013时，会得到此错误。
xutility(2132): error C4996: 'std::_Copy_impl': Function call with parameters that may be unsafe - this call relies on the caller to check that the passed values are correct. To disable this warning, use -D_SCL_SECURE_NO_W
浏览 1提问于2015-05-27得票数 0
回答已采纳
我有DB Vertica 8.0。创建模式仓库，表按类型包含字段LONG VARCHAR。现在我尝试执行SELECT，例如
SELECT * FROM WAREHOUSE.ALL_EVENTS a
WHERE
a.original_data like '%d963%'
返回误差
SQL Error [4286] [42883]: [Vertica][VJDBC](4286) ERROR: Operator does not exist: long varchar ~~ unknown
[Vertica][VJDBC](4286) ERROR: Operator does
浏览 0提问于2017-10-20得票数 2
回答已采纳
6回答
我是VxWorks的新手，我正在VxWorks平台上使用C++开发一个软件。我想知道VxWorks编译器是否支持C++ 11标准。我之所以问这个问题，是因为没有可用的shrink_to_fit() std：：向量函数(这个函数是在c++ 11标准中引入的)。所以我想知道有什么方法可以用C++ 11标准在VxWorks中编译代码吗？
浏览 20提问于2016-03-30得票数 7
回答已采纳
我不知道如何让scoped_ptr或scoped_array使用自定义删除器。也许还有另一个类似于shared_ptr的允许受控删除的实现。
顺便问一下，为什么shared_ptr允许自定义删除，而scoped_ptr不允许？只是好奇而已。
浏览 13提问于2011-04-15得票数 7
回答已采纳
我在使用'dense output‘的state_type ublas::vector<std::complex<double>>时收到以下错误，但只使用runge_kutta_dopri5 stepper进行编译。
C2440 'return'：无法在boost\numeric\odeint\stepper\controlled_runge_kutta.hpp 89中将‘std：：C2440’转换为'double‘
#include <iostream>
#include <complex>
#include &
浏览 1提问于2017-07-28得票数 0
回答已采纳
我正在研究使用最优控制方法(轨迹生成和优化)的无限维优化算法。我想要做这个轨迹优化的系统是非线性的，形式是$\dot{x}(t) = f(x(t)，u(t)，t)$。换句话说，我有时变的非线性输入。
是否有可能使用boost::odeint解决这样的ODE问题？我在文档中没有发现任何提示，但我可能没有看到它。
浏览 3提问于2014-04-09得票数 1
回答已采纳
1回答
不久前，我在这里发布了一个模型，我正在尝试使用pycuda构建一个模型，它可以解决大约9000个耦合ODE。然而，我的模型太慢了，一位SO成员建议，内存从主机传输到GPU可能是罪魁祸首。
目前，cuda仅用于计算我正在处理的9000个物种中每个物种的变化率。由于我从主机向GPU传递了一个数组来执行此计算，并从GPU返回一个数组以集成到主机上，因此我可以看到这将如何减慢速度。
boost能解决我的问题吗？据我所知，boost允许c++和python之间的互操作性。它还包括我读到的c++ odeint，它与推力相结合，可以在GPU上实现快速缩减和集成。我的理解正确吗？
谢谢你，卡斯汀
浏览 1提问于2013-05-23得票数 0
回答已采纳
我有一个由4个耦合方程组成的系统要求解，还有一个参数Gammai要迭代。因为我是C++的新手，所以我的代码非常初级。如果它在某些部分看起来复杂而优雅，那是因为我改编了odeint作者的代码。:)
这个问题与()有关，但并不完全相同。请不要删除此内容。:(
在代码行之间插入了问题。
#include <iostream>
#include <iterator>
#include <algorithm>
#include <boost/numeric/odeint.hpp>
#include <cmath>
#include <ve
浏览 2提问于2012-08-28得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
