armadillo

1热度

1回答

我试图在我的Ubuntu 14.04中安装Armadillo。我运行cmake .命令在终端的Armadillo文件所在的文件夹，然后运行make命令。我得到了以下错误。 /usr/bin/ld: /usr/local/lib/libsuperlu.a(util.c.o): relocation R_X86_64_32 against `.rodata' can not be used when

0热度

1回答

设置Armadillo的线程数

如何让Armadillo在OpenMP上使用OpenBLAS？我用--with-openmp选项安装了带有Homebrew的OpenBLAS，并且正在使用：/usr/local/Cellar/gcc/6.3.0_1/bin/g++-6 -DARMA_DONT_USE_WRAPPER -fopenmp -L/usr/local/Cellar/armadillo/7.600.2/lib -L/usr/

0热度

2回答

Spectra的SymEigsShiftSolver是否不返回特征向量？

我已经实现了SymEigsShiftSolver来计算大型稀疏矩阵的特征值，但是它并没有返回给我的特征向量。它尚未实施到目前为止？ void Eigens::computeEigenvectors(Matrices m) { SparseSymShiftSolve<double> op(m.Lpl); SymEigsShiftSolver< double, SMALLEST_MAGN, Sp

0热度

1回答

犰狳n_rows不起作用

我使用犰狳来获取矩阵，但似乎它给了我关于cols和行的错误数字。第一张图片我有一个叫做se的30x30矩阵。我打印行数和列数。然后我将se放入一个名为Acqui的函数中，并打印该函数中的行和列。 //define 'mat se' as a 30 by 30 matrix and give it certain value printf("%u \n",se.n_rows); printf(

0热度

1回答

如何使用犰狳

谁能告诉我什么是以下操作的犰狳线性代数包相当于 L到得到正规拉普拉斯矩阵= d^-0.5 * A * d^-0.5 一般如何来计算^ n或在犰狳a^-0.5其中A是一个方阵

0热度

1回答

RcppArmadillo样本上犰狳矢量类

我们一直在使用sample函数从RcppArmadillo随机抽样一个NumericVector对象。但是，我们注意到在犰狳类型上不可能使用相同的功能（vec或uvec）。我们已经看了sample.h文件中的函数定义，它看起来像一个模板化的函数，应该能够处理这些类型，但是我们还没有能够弄清楚如何使它与Armadillo类一起工作而没有做很多来自Rcpp库的NumericVector或Integer

2热度

1回答

Armadillo C++：如何使用另一个矩阵中的多个元素修改矩阵的多个数组元素，特别是在多维数据集结构中？

在MATLAB，如果我有两个3×3矩阵A和B，并希望B的最后两列复制到前两个A的，我会用下面的语法： A(:,1:2) = B(:,2:3) 我试图完成在C++中使用Armadillo的相同动作，专门用于多维数据集结构。在犰狳，如果我有两个立方体A和B具有九个片（每片是一个3×3矩阵），I假设我会使用以下方法来执行相同的列元素的更新： A(span(0,2),span(0,1),span(i)

0热度

2回答

Armadillo与Boost Odeint冲突：Odeint在整合期间将状态向量调整为零

我刚开始使用Boost Odeint来整合ODE系统。为了方便起见，我想将它与Armadillo一起使用，因为它们都是具有方便API的现代C++库。但是，如果我将arma::vec指定为状态类型，立即证明在集成的第一步中，integrate_adaptive()将状态向量调整为0x1。我在这里发布一个简单的例子： #include <iostream> #include <armadillo>

1热度

1回答

Armadillo中的SpMat <Type>的迭代器只访问非零条目吗？

我想知道如何在Armadillo中遍历所有非零条目的sp_umat（即SpMat<unsigned int>），并且遇到了这个相关的问题（link）。该帖子建议使用const_iterator来检索sp_mat中的非零位置和值。假设所有迭代器sp_mat（和犰狳中的其他相关类型的稀疏矩阵;在我的情况下为sp_umat）仅支持只访问非零条目的迭代器？我无法从文档中找出这个问题。另一个相关的问题也出现

2热度

1回答

如何根据行总数条件在Armadillo中删除SpMat <unsigned int>的行？

是否有一种有效的方法仅保留Armadillo稀疏矩阵的行，总和至少达到矩阵列的总计数的某个级别？例如，如果其值为>=C，其中C是某个选定值，我想要保留第i行。犰狳的文档说只有连续的子矩阵视图才被允许使用稀疏矩阵。所以我猜这是不容易通过子设置获得的。是否有一种替代方法可以简单地循环遍历元素，并创建一个新的稀疏矩阵，并使用符合所需条件的新位置，值和colPtr设置？谢谢！