- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
基于@Piotr Skotnicki 的回答 here我定义了以下数据类型:
#ifndef TENSOR_HPP
#define TENSOR_HPP
#include <vector>
#include <array>
#include <boost/numeric/odeint.hpp>
namespace nex {
template <size_t...> struct seq {};
template <size_t N, size_t... Is> struct gen_seq : gen_seq<N-1, N-1, Is...> {};
template <size_t... Is> struct gen_seq<0, Is...> { using type = seq<Is...>; };
template<typename T, typename S>
class tensor;
template<typename T, size_t... Is>
class tensor<T, seq<Is...>>
{
typedef std::vector<T> vector;
public:
// constructor
tensor ();
tensor (decltype(Is)...);
// iterators
using iterator = typename vector::iterator;
using const_iterator = typename vector::const_iterator;
iterator begin()
{ return m_v.begin(); }
const_iterator begin() const
{ return m_v.begin(); }
iterator end()
{ return m_v.end(); }
const_iterator end() const
{ return m_v.end(); }
// resizing
size_t size() const;
void resize( const size_t );
// operators
const T& operator() (decltype(Is)...) const;
T& operator() (decltype(Is)...);
private:
// helper functions
template <size_t I>
auto index(const std::array<size_t, sizeof...(Is)>& a) const
-> typename std::enable_if<I == 0, size_t>::type
{
return a[I];
}
template <size_t I>
auto index(const std::array<size_t, sizeof...(Is)>& a) const
-> typename std::enable_if<I != 0, size_t>::type
{
return index<I-1>(a) * s[I] + a[I];
}
size_t mult(size_t N)
{
return N;
}
template <typename... S>
size_t mult(size_t N, S... Ns)
{
return N * mult(Ns...);
}
vector m_v;
const std::array<size_t, sizeof...(Is)> s;
}; // class tensor
template<typename T, size_t... Is
tensor<T, seq<Is...>> :: tensor ()
: m_v(), s{ { 0 } }
{}
template<typename T, size_t... Is>
tensor<T, seq<Is...>> :: tensor (decltype(Is)... size)
: m_v( mult( size... ) ), s{ { size... } }
{}
template<typename T, size_t... Is>
size_t tensor<T, seq<Is...>> :: size () const
{ return m_v.size(); }
template<typename T, size_t... Is>
void tensor<T, seq<Is...>> :: resize (const size_t n)
{ m_v.resize( n ); }
template<typename T, size_t... Is>
const T& tensor<T, seq<Is...>> :: operator() (decltype(Is)... n) const
{ return m_v.at(index<sizeof...(Is)-1>( { { n... } } ) ); }
template<typename T, size_t... Is>
T& tensor<T, seq<Is...>> :: operator() (decltype(Is)... n)
{ return m_v.at(index<sizeof...(Is)-1>( { { n... } } ) ); }
} // namespace nex
namespace boost { namespace numeric { namespace odeint {
template<typename T, size_t... Is>
struct is_resizeable< nex::tensor<T, nex::seq<Is...>> >
{
typedef boost::true_type type;
static const bool value = type::value;
};
} // namespace odeint
} // namespace numeric
} // namespace boost
template <typename T, std::size_t N>
using tensor = nex::tensor<T, typename nex::gen_seq<N>::type>;
#endif // TENSOR_HPP
然后我想在 boost::odeint 中使用这个数据类型来集成非常基本的问题
#include "tensor.hpp"
typedef tensor<double, 1> state_type;
void lorenz( const state_type &x , state_type &dxdt , const double t )
{
dxdt(0) = x(0);
dxdt(1) = x(1);
dxdt(2) = x(2);
}
void write_lorenz( const state_type &x , const double t )
{
std::cout << t << '\t' << x(0) << '\t' << x(1) << '\t' << x(2) << std::endl;
}
using namespace boost::numeric::odeint;
int main()
{
state_type x(3);
x(0) = 1.0 ; x(1) = 10.0 ; x(2) = 10.0;
integrate_const( runge_kutta4< state_type >() , lorenz , x , 0.0 , 1.0 , 0.1, write_lorenz );
}
即使对于不同的方程组(如洛伦兹吸引子),它也能正常工作。不幸的是,当我想使用张量的相同容量但定义不同时,一切都出错了
#include "tensor.hpp"
typedef tensor<double, 2> state_type;
void lorenz( const state_type &x , state_type &dxdt , const double t )
{
dxdt(0,0) = x(0,0);
dxdt(1,0) = x(1,0);
dxdt(2,0) = x(2,0);
}
void write_lorenz( const state_type &x , const double t )
{
std::cout << t << '\t' << x(0,0) << '\t' << x(1,0) << '\t' << x(2,0) << std::endl;
}
using namespace boost::numeric::odeint;
int main()
{
state_type x(3,1);
x(0,0) = 1.0 ; x(1,0) = 10.0 ; x(2,0) = 10.0;
integrate_const( runge_kutta4< state_type >() , lorenz , x , 0.0 , 1.0 , 0.1, write_lorenz );
}
您可以自己运行代码并检查 x(1,0) x(2,0) 是否根本没有改变。我认为,operator() 的适当定义足以让 odeint 正常工作。我对 boost::odeint 的了解相当少,所以也许这只是一个基本的错误。我指望你的帮助。
当我使用 std::array 将张量定义为固定大小的数组时,它工作正常。但这样我就不需要添加 boost::is_resizeable 了。我认为调整大小需要以不同的方式定义(我不太清楚)。
最佳答案
确实,问题出在调整大小上。目前,调整张量大小只会更改底层向量 m_v
的大小。虽然这确保了正确的内存分配,但这还不够。您还需要正确设置维度,即成员s
。您也可以这样想:仅给定一个整数,不可能调整多维张量的大小,您需要设置每个维度的大小。唯一的异常(exception)是一维张量,对于它来说一切都很好,对于任何更高的维度,这都必须失败。
以下tensor.hpp解决了这个问题,并且在我的快速测试中似乎运行良好:
#ifndef TENSOR_HPP
#define TENSOR_HPP
#include <vector>
#include <array>
#include <boost/numeric/odeint.hpp>
namespace nex {
template <size_t...> struct seq {};
template <size_t N, size_t... Is> struct gen_seq : gen_seq<N-1, N-1, Is...> {};
template <size_t... Is> struct gen_seq<0, Is...> { using type = seq<Is...>; };
template<typename T, typename S>
class tensor;
template<typename T, size_t... Is>
class tensor<T, seq<Is...>>
{
typedef std::vector<T> vector;
public:
// constructor
tensor ();
tensor (decltype(Is)...);
// iterators
using iterator = typename vector::iterator;
using const_iterator = typename vector::const_iterator;
iterator begin()
{ return m_v.begin(); }
const_iterator begin() const
{ return m_v.begin(); }
iterator end()
{ return m_v.end(); }
const_iterator end() const
{ return m_v.end(); }
// resizing
size_t size() const;
void resize( const tensor<T, seq<Is...>>& );
// operators
const T& operator() (decltype(Is)...) const;
T& operator() (decltype(Is)...);
private:
// helper functions
template <size_t I>
auto index(const std::array<size_t, sizeof...(Is)>& a) const
-> typename std::enable_if<I == 0, size_t>::type
{
return a[I];
}
template <size_t I>
auto index(const std::array<size_t, sizeof...(Is)>& a) const
-> typename std::enable_if<I != 0, size_t>::type
{
return index<I-1>(a) * s[I] + a[I];
}
size_t mult(size_t N)
{
return N;
}
template <typename... S>
size_t mult(size_t N, S... Ns)
{
return N * mult(Ns...);
}
vector m_v;
std::array<size_t, sizeof...(Is)> s;
}; // class tensor
template<typename T, size_t... Is>
tensor<T, seq<Is...>> :: tensor ()
: m_v(), s{ { 0 } }
{}
template<typename T, size_t... Is>
tensor<T, seq<Is...>> :: tensor (decltype(Is)... size)
: m_v( mult( size... ) ), s{ { size... } }
{}
template<typename T, size_t... Is>
size_t tensor<T, seq<Is...>> :: size () const
{ return m_v.size(); }
template<typename T, size_t... Is>
void tensor<T, seq<Is...>> :: resize (const tensor<T, seq<Is...>> &x)
{
m_v.resize( x.size() );
s = x.s;
}
template<typename T, size_t... Is>
const T& tensor<T, seq<Is...>> :: operator() (decltype(Is)... n) const
{ return m_v.at(index<sizeof...(Is)-1>( { { n... } } ) ); }
template<typename T, size_t... Is>
T& tensor<T, seq<Is...>> :: operator() (decltype(Is)... n)
{ return m_v.at(index<sizeof...(Is)-1>( { { n... } } ) ); }
} // namespace nex
namespace boost { namespace numeric { namespace odeint {
template<typename T, size_t... Is>
struct is_resizeable< nex::tensor<T, nex::seq<Is...>> >
{
typedef boost::true_type type;
static const bool value = type::value;
};
template<typename T, size_t... Is>
struct resize_impl< nex::tensor<T, nex::seq<Is...>>, nex::tensor<T, nex::seq<Is...>> >
{
typedef nex::tensor<T, nex::seq<Is...>> state_type;
static void resize( state_type &x1 , const state_type &x2 )
{
x1.resize(x2);
}
};
} // namespace odeint
} // namespace numeric
} // namespace boost
template <typename T, std::size_t N>
using tensor = nex::tensor<T, typename nex::gen_seq<N>::type>;
#endif // TENSOR_HPP
关于c++11 - 使用自定义张量数据结构的 boost::odeint 输出错误,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/33546040/
我已经使用 vue-cli 两个星期了,直到今天一切正常。我在本地建立这个项目。 https://drive.google.com/open?id=0BwGw1zyyKjW7S3RYWXRaX24tQ
您好,我正在尝试使用 python 库 pytesseract 从图像中提取文本。请找到代码: from PIL import Image from pytesseract import image_
我的错误 /usr/bin/ld: errno: TLS definition in /lib/libc.so.6 section .tbss mismatches non-TLS reference
我已经训练了一个模型,我正在尝试使用 predict函数但它返回以下错误。 Error in contrasts<-(*tmp*, value = contr.funs[1 + isOF[nn]])
根据Microsoft DataConnectors的信息我想通过 this ODBC driver 创建一个从 PowerBi 到 PostgreSQL 的连接器使用直接查询。我重用了 Micros
我已经为 SoundManagement 创建了一个包,其中有一个扩展 MediaPlayer 的类。我希望全局控制这个变量。这是我的代码: package soundmanagement; impo
我在Heroku上部署了一个应用程序。我正在使用免费服务。 我经常收到以下错误消息。 PG::Error: ERROR: out of memory 如果刷新浏览器,就可以了。但是随后,它又随机发生
我正在运行 LAMP 服务器,这个 .htaccess 给我一个 500 错误。其作用是过滤关键字并重定向到相应的域名。 Options +FollowSymLinks RewriteEngine
我有两个驱动器 A 和 B。使用 python 脚本,我在“A”驱动器中创建一些文件,并运行 powerscript,该脚本以 1 秒的间隔将驱动器 A 中的所有文件复制到驱动器 B。 我在 powe
下面的函数一直返回这个错误信息。我认为可能是 double_precision 字段类型导致了这种情况,我尝试使用 CAST,但要么不是这样,要么我没有做对...帮助? 这是错误: ERROR: i
这个问题已经有答案了: Syntax error due to using a reserved word as a table or column name in MySQL (1 个回答) 已关闭
我的数据库有这个小问题。 我创建了一个表“articoli”,其中包含商品的品牌、型号和价格。 每篇文章都由一个 id (ID_ARTICOLO)` 定义,它是一个自动递增字段。 好吧,现在当我尝试插
我是新来的。我目前正在 DeVry 在线学习中级 C++ 编程。我们正在使用 C++ Primer Plus 这本书,到目前为止我一直做得很好。我的老师最近向我们扔了一个曲线球。我目前的任务是这样的:
这个问题在这里已经有了答案: What is an undefined reference/unresolved external symbol error and how do I fix it?
我的网站中有一段代码有问题;此错误仅发生在 Internet Explorer 7 中。 我没有在这里发布我所有的 HTML/CSS 标记,而是发布了网站的一个版本 here . 如您所见,我在列中有
如果尝试在 USB 设备上构建 node.js 应用程序时在我的树莓派上使用 npm 时遇到一些问题。 package.json 看起来像这样: { "name" : "node-todo",
在 Python 中,您有 None单例,在某些情况下表现得很奇怪: >>> a = None >>> type(a) >>> isinstance(a,None) Traceback (most
这是我的 build.gradle (Module:app) 文件: apply plugin: 'com.android.application' android { compileSdkV
我是 android 的新手,我的项目刚才编译和运行正常,但在我尝试实现抽屉导航后,它给了我这个错误 FAILURE: Build failed with an exception. What wen
谁能解释一下?我想我正在做一些非常愚蠢的事情,并且急切地等待着启蒙。 我得到这个输出: phpversion() == 7.2.25-1+0~20191128.32+debian8~1.gbp108
我是一名优秀的程序员,十分优秀!