- iOS/Objective-C 元类和类别
- objective-c - -1001 错误,当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
我打算在 MCMC 例程中使用 Boost odeint 库来估计 ODE 模型中的参数。由于这些 ODE 可能是刚性的,我需要能够将雅可比矩阵传递到具有导数的求解器中。我想创建一个类,它具有作为私有(private)成员的参数和初始值,然后是派生类、jacobian 和将参数更改为公共(public)方法的方法。我试图修改 odeint 网站上的僵硬示例,以使用一个包含两者的类,但在编译时收到错误“错误:没有匹配函数调用‘Fitzhugh::deriv()’”。我不是经验丰富的 C++ 程序员,所以这很可能是概念上的错误。这是代码。
/* Fitzhugh Nagumo Equation in odeint */
#include <iostream>
#include <fstream>
#include <utility>
#include <cmath>
#include <boost/numeric/odeint.hpp>
#include <boost/phoenix/core.hpp>
#include <boost/phoenix/operator.hpp>
using namespace std;
using namespace boost::numeric::odeint;
namespace phoenix = boost::phoenix;
//[ stiff_system_definition
typedef boost::numeric::ublas::vector< double > vector_type;
typedef boost::numeric::ublas::matrix< double > matrix_type;
class Fitzhugh
{
private:
double a;
double b;
double c;
public:
Fitzhugh( double a_, double b_, double c_ ) : a(a_), b(b_), c(c_) { }
void deriv ( const vector_type &x , vector_type &dxdt , double )
{
dxdt[ 0 ] = c*(x[0] - pow(x[0], 3.0)/3.0 + x[1]);
dxdt[ 1 ] = -(x[0] - a + b*x[1])/c;
}
void jac ( const vector_type &x, matrix_type &J , const double & , vector_type &dfdt )
{
J( 0 , 0 ) = c*(1 - pow(x[0], 2.0));
J( 0 , 1 ) = c;
J( 1 , 0 ) = -1/c;
J( 1 , 1 ) = -b/c;
dfdt[0] = 0.0;
dfdt[1] = 0.0;
}
};
int main( int argc , char **argv )
{
// typedef rosenbrock4< double > stepper_type;
// typedef rosenbrock4_controller< stepper_type > controlled_stepper_type;
// typedef rosenbrock4_dense_output< controlled_stepper_type > dense_output_type;
//[ integrate_stiff_system
vector_type x( 2 , 1.0 );
Fitzhugh fitzhugh(0.2, 0.2, 3.0);
size_t num_of_steps = integrate_const( make_dense_output< rosenbrock4< double > >( 1.0e-6 , 1.0e-6 ) ,
make_pair( fitzhugh.deriv() , fitzhugh.jac() ) ,
x , 0.0 , 50.0 , 0.01 ,
cout << phoenix::arg_names::arg2 << " " << phoenix::arg_names::arg1[0] << "\n" );
//]
clog << num_of_steps << endl;
return 0;
}
这是完整的输出
/home/chris/C++/examples/fitzhugh/main.cpp: In function ‘int main(int, char**)’:
/home/chris/C++/examples/fitzhugh/main.cpp:60:39: error: no matching function for call to ‘Fitzhugh::deriv()’
/home/chris/C++/examples/fitzhugh/main.cpp:60:39: note: candidate is:
/home/chris/C++/examples/fitzhugh/main.cpp:30:14: note: void Fitzhugh::deriv(const vector_type&, vector_type&, double)
/home/chris/C++/examples/fitzhugh/main.cpp:30:14: note: candidate expects 3 arguments, 0 provided
/home/chris/C++/examples/fitzhugh/main.cpp:60:56: error: no matching function for call to ‘Fitzhugh::jac()’
/home/chris/C++/examples/fitzhugh/main.cpp:60:56: note: candidate is:
/home/chris/C++/examples/fitzhugh/main.cpp:36:14: note: void Fitzhugh::jac(const vector_type&, matrix_type&, const double&, vector_type&)
/home/chris/C++/examples/fitzhugh/main.cpp:36:14: note: candidate expects 4 arguments, 0 provided
最佳答案
要传递成员函数,您需要使用绑定(bind)机制。如果你有一个 c++11 编译器,你可以使用 std::bind;包括 std::placeholders 命名空间:using namespace std::placeholders;
然后使用来自 <functional>
的 std::bind :
make_pair( bind( &Fitzhugh::deriv , &fitzhugh , _1 , _2 , _3 ) ,
bind( &Fitzhugh::jac , &fitzhugh , _1 , _2 , _3, _4 ) )
关于c++ - 使用导数和 jacobian boost odeint 类,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/18928733/
我正在通过 odeint 集成生化常微分方程(见下图),但是主输入函数似乎在调用时用奇怪的数字重新设置(或替换)输入参数。尽管参数“iu”(v 的感应率,应在整个过程中保持恒定)和“v”(结果之一,应
我一直在努力使用新版本的 boost。我正在使用多精度的 odeint。下面这段代码可以用boost version 1.67.0编译成功。但是,自 1.68.0 及更新版本以来,我无法再编译。在版本
我最近偶然发现了 boost.odeint 库,我对可能性和可配置性的数量感到惊讶。但是,在广泛使用 scipy.integrate.odeint(它本质上是 Fortran 中 ODEPACK 的包
我刚开始使用 Boost Odeint 来集成一个 ODE 系统。为方便起见,我想将它与 Armadillo 一起使用,因为两者都是具有方便 API 的现代 C++ 库。但是,如果我指定 arma::
我想使用 scipy 的 odeint 函数求解具有 15 个时间相关系数的 7 个常微分方程 (ODE) 系统。 我将系数存储在字典中,以便可以通过我定义为与 odeint() 一起使用的函数 (f
我尝试在 Mac OS X 10.9.2 g++ 5.1 上的 boost_1_55_0 中运行 [odeint 复杂状态类型示例代码。 下面的代码是网站上解决 Stuart-Landau 振荡器的拷
所以我试图求解一个包含三个 ODE 的系统,并开发了以下代码来使用 ODEint 求解它们。但是当我运行时,ODEint 在调用我的方程组函数时出现问题。 from scipy.integrate i
我有以下 odeint 程序: #include #include using namespace std; typedef boost::array state_type; void eqsys
我有一个使用“odeint”模拟种群动态的程序。我想设置一个 if 条件来禁止我的颂歌的结果为负。这是我的代码摘要: class Communities { public : type
我正在尝试求解一个简单的方程:dM/dr = r*p(r) 在 python 中。 我在 r 的某些值处有 p 的值: p(0)=1, p(1)=3, p(2)=5, p( 3)=7, p(4)=9,
我正在尝试使用 scipy 的 odeint 来求解一些常微分方程。唯一的问题是我只想定义一个参数,看来要组成一个元组,你至少需要两个值。 我的代码是这样的: def system(state, t,
我正在尝试用 odeint 求解微分方程。这里一些常量参数是固定的,一些在列表中。 from scipy.integrate import odeint import matplotl
使用Python 2.7.8。 我正在使用的微分方程是 x'=2-3*x。没那么难。正确的解是 y 截距为 2/3 的衰减指数。运动有三个初始条件。还必须在同一地 block 上有一个带有解决方案的斜
我想将 scipy 的 odeint 与一个函数一起使用 def func(y,t,a=123,b=456) 然后将其用作 odeint(func,y0,t) 如果我想使用 args 改变值 a 和
我对使用隐式方案使用 odeint 库求解 ODE 系统很感兴趣,但我很难实现一个简单的 implicit_euler 示例。 查看文档,我设法使工作显式步进器、自适应步进器以及 rosenbrock
我刚刚实现了一组耦合 ODE 的数值积分来自使用 odeint C++ 库的离散 PDE。它很好用并且快如闪电,但有一个问题: 我的 ODE 系统具有所谓的吸收边界条件:时间我的状态变量 n 的导数,
运行以下代码: #include #include using namespace std; using namespace boost::numeric::odeint; class CSyst
显然,getting a non-negative solution from an ODE solver is non-trivial .在 Matlab 中,有 NonNegative optio
我在使用受控错误步进器和复杂状态类型的 odeint 库时遇到了问题。我对具有复杂斯图尔特朗道方程的示例中的代码进行了修改,使其包含自适应积分器。代码现在看起来像这样: #include #incl
有一些使用 arbitrary precision 的例子和 matrices在 boost.odeint( boost 常微分方程求解器)中。 我想在不同类型的坐标(笛卡尔、极坐标或作用角)中使用
我是一名优秀的程序员,十分优秀!