- Java 双重比较
- java - 比较器与 Apache BeanComparator
- Objective-C 完成 block 导致额外的方法调用?
- database - RESTful URI 是否应该公开数据库主键?
在下面的代码中,ExtractSubArray
函数是完全通用的,而 ExtractSubArrayCornerAndExtent
需要在编写代码时了解维度(以构造RangeType 参数)。有什么方法可以编写通用的 ExtractSubArrayCornerAndExtent
(没有 SFINAE 为每个维度使用不同的函数(这会很烦人,并且需要一组固定的可能维度)。
#include <boost/multi_array.hpp>
template <unsigned int Dimension>
boost::multi_array<double, Dimension> ExtractSubArray(const boost::multi_array<double, Dimension>& array, const typename boost::detail::multi_array::index_gen<Dimension, Dimension>& indices)
{
using ArrayType = boost::multi_array<double, Dimension>;
using IndexType = boost::array<double, Dimension>;
typename ArrayType::template const_array_view<Dimension>::type view = array[indices];
IndexType subArraySize;
for(size_t dim = 0 ; dim < Dimension; ++dim) {
subArraySize[dim] = indices.ranges_[dim].finish_ - indices.ranges_[dim].start_;
}
ArrayType subArray = view;
return subArray;
}
template <unsigned int Dimension>
boost::multi_array<double, Dimension> ExtractSubArrayCornerAndExtent(const boost::multi_array<double, Dimension>& array, const boost::array<double, Dimension>& corner,
const boost::array<double, Dimension>& subarraySize)
{
using ArrayType = boost::multi_array<double, Dimension>;
using RangeType = typename ArrayType::index_range;
// Here I have assumed Dimension=3 to produce the second argument. How do you construct this second argument when you don't know Dimension ahead of time.
return ExtractSubArray<Dimension>(array, boost::indices[RangeType(corner[0],subarraySize[0])][RangeType(corner[1],subarraySize[1])][RangeType(corner[2],subarraySize[2])]);
}
int main()
{
using ArrayType = boost::multi_array<double, 3>;
using IndexType = boost::array<double, 3>;
ArrayType myArray(IndexType({{3,3,3}}));
std::vector<double> data(9);
for(size_t i = 0; i < data.size(); ++i) {
data[i] = i;
}
boost::detail::multi_array::index_gen<3,3> indices = boost::indices[ArrayType::index_range(0,1)][ArrayType::index_range(0,1)][ArrayType::index_range(0,1)];
ArrayType subArray = ExtractSubArray<3>(myArray, indices);
IndexType corner = {0,0,0};
IndexType subarraySize = {1,1,1};
ArrayType subArray2 = ExtractSubArrayCornerAndExtent<3>(myArray, corner, subarraySize);
return 0;
}
最佳答案
编译时迭代通常必须用递归来实现。您可以使用递归模板构建索引对象。它看起来像这样:
#include <boost/multi_array.hpp>
template <typename T, size_t Dimension>
boost::multi_array<T, Dimension> ExtractSubArray(
const boost::multi_array<T, Dimension>& array,
const typename boost::detail::multi_array::index_gen<Dimension, Dimension>& indices)
{
using ArrayType = boost::multi_array<T, Dimension>;
using IndexType = boost::array<size_t, Dimension>;
typename ArrayType::template const_array_view<Dimension>::type view = array[indices];
IndexType subArraySize;
for(size_t dim = 0 ; dim < Dimension; ++dim) {
subArraySize[dim] = indices.ranges_[dim].finish_ - indices.ranges_[dim].start_;
}
ArrayType subArray = view;
return subArray;
}
// Helper functor to build indices.
template<typename RangeArrayType, size_t Dimension>
struct IndicesBuilder {
// Recursively invoke the functor for the next lowest dimension and
// add the next range.
static auto build(const RangeArrayType& range)
-> decltype(IndicesBuilder<RangeArrayType, Dimension - 1>::build(range)[range[Dimension - 1]]) {
return IndicesBuilder<RangeArrayType, Dimension - 1>::build(range)[range[Dimension - 1]];
}
};
// Helper functor specialization to terminate recursion.
template<typename RangeArrayType>
struct IndicesBuilder<RangeArrayType, 1> {
static auto build(const RangeArrayType& range)
-> decltype(boost::indices[range[0]]) {
return boost::indices[range[0]];
}
};
template <typename T, size_t Dimension>
boost::multi_array<T, Dimension> ExtractSubArrayCornerAndExtent(
const boost::multi_array<T, Dimension>& array,
const boost::array<size_t, Dimension>& corner,
const boost::array<size_t, Dimension>& subarraySize)
{
using ArrayType = boost::multi_array<T, Dimension>;
using RangeType = typename ArrayType::index_range;
// Build a random-access container with the ranges.
std::vector<RangeType> range;
for (size_t i = 0; i < Dimension; ++i)
range.push_back(RangeType(corner[i], corner[i] + subarraySize[i]));
// Use the helper functor to build the index object.
const auto index = IndicesBuilder<decltype(range), Dimension>::build(range);
return ExtractSubArray<T, Dimension>(array, index);
}
int main() {
using ArrayType = boost::multi_array<double, 3>;
using IndexType = boost::array<size_t, 3>;
ArrayType myArray(IndexType({{3,3,3}}));
IndexType corner = {0,0,0};
IndexType subarraySize = {1,1,1};
ArrayType subArray2 = ExtractSubArrayCornerAndExtent(myArray, corner, subarraySize);
return 0;
};
添加到您的代码中的是 IndicesBuilder
模板结构,它具有用于构建索引对象的一维的静态函数。该静态函数调用自身(使用不同的模板参数)来构建所有维度。递归由返回第一个维度的 IndicesBuilder
的特化终止。尾随返回类型万岁!
我编译并运行了这段代码,没有出现 fatal error ,但我没有做任何进一步的测试,因此可能存在一些小错误。不过,一般方法应该可行。
关于c++ - 如何在运行时构建 multi_array::index_gen,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36501811/
我在使用 gradle 构建一个特定应用程序时遇到问题。该应用程序可以用 eclipse 编译和构建,它在平板电脑上运行良好。当我尝试使用 Gradle 构建它时,“compileDebugJava”
我有一个 C 程序,是一位离开的开发人员留给我的。我试图弄清楚他到底在做什么,并将软件重新安排成更合乎逻辑的东西,这样我就可以更轻松地构建它。我正在使用 CMake 构建,而他使用的是 Make。 有
我刚开始阅读“Pro Spring MVC with web flow”,它附带了一个我想遵循的代码示例。 我要什么 - 我想像书中那样构建应用程序,使用 Gradle 有什么问题 - 我没用过 Gr
我希望有人已经这样做了。我正在尝试为我的一个 angular 2 项目在 teamcity 中建立一个连续的构建。在做了一些研究之后,我按照以下步骤操作: 构建步骤 1:为 teamcity 安装 j
我有一个旧的 ASP.Net 网站解决方案,看起来像: 当我在 Visual Studio 中构建解决方案时,我得到以下输出: ------ Build started: Project: C:\..
我使用 gulp-usref、gulp-if、gulp-uglify、gulp-csso 和 gulp-file-include 来构建我的应用程序。除了 HTML 保持原样外,构建中的一切都运行良好
我正在使用 ionic2 开发内部移动应用程序。我可以通过以下方式成功构建 ios: ionic build ios and ionic build ios --prod 但当我这样做时,它一直失败
我是一位经验丰富的 .NET/C# 开发人员,但对这里的几乎所有技术/库(包括 SQL/DB 工作)都是新手。 我正在开发一个具有 Azure/Entity Framework .NET 后端和可移植
我正在使用 VS 2008。我可以使用 IDE 成功编译我的解决方案。但是,当我尝试使用 devenv.com 构建它时,它失败并提示“错误:找不到项目输出组'(无法确定名称)的输出”。该组、其配置或
版本: ember.js 2.7,ember-data 2.7 ember-cli 2.9.1//同样适用于 ember-cli 2.7 node 6.9.1, npm 3.10.9//也适用于 no
我第一次修补 AzureDevops,设置一些 CI 任务。 我有一个公共(public)存储库(开源)和一个包含 3 个 F# 项目的解决方案(.sln)。该解决方案在 Windows/Mac/Li
目前 5.1.5 版本或 STLPort CVS 存储库似乎仍不支持 VS2008。如果有人已经完成了这项工作,那么如果可能的话,分享会很有用:) 同样,了解 VS2005 或 2008 x64 构建
我有一个 Python 2.7 项目,到目前为止一直使用 gfortran 和 MinGW 来构建扩展。我使用 MinGW,因为它似乎支持 Fortran 代码中的写入语句和可分配数组,而 MSVC
关闭。这个问题是off-topic .它目前不接受答案。 想改进这个问题? Update the question所以它是on-topic对于堆栈溢出。 9年前关闭。 Improve this que
我想知道为什么在 Zimbra Wiki 中只列出了构建过程的特定平台。这意味着不可能在其他 Linux 发行版上构建 Zimbra? Zimbra 社区选择一个特殊的 Linux 发行版来构建 Zi
我将在 Swift 中构建一个 CLI 工具。我用这个命令创建了项目 swift package init --type executable当我构建我的项目并解析 时读取别名 Xcode 中的参数并
我想为添加到 docker 镜像的文件设置文件权限。我有这个简单的 Dockerfile: FROM ubuntu:utopic WORKDIR /app RUN groupadd -g 1000 b
当我使用 clBuildProgram在我的 OpenCl 代码中,它失败并显示错误代码 -11,没有任何日志信息。 这是我的代码的样子: ret = clBuildProgram(program
我有一个底部导航栏,它有一个列表页面,该页面使用状态块。 class _MainPageState extends State { int _index = 0; @override Wi
我在本地计算机上使用Jenkins(Jenkins URL未通过Internet公开,但该计算机上已启用Internet。) 我进行了以下配置更改: 在Jenkins工具上安装了Git和Github插
我是一名优秀的程序员,十分优秀!