1回答
我正在尝试模仿R函数，它允许基于索引向量运行列和行矩阵排列。sample(c(0:9),5*5,T),ncol=5,nrow=5)rand=sample(c(1:5))我在c++中尝试了以下代码：
Library(Rcpp
浏览 0提问于2018-02-27得票数 1
回答已采纳
1回答
C++特征LU分解-行列式的符号
、、、、
我使用C++ 14和Eigen，我想使用LU分解来计算方阵的行列式。()，但我有一些问题: A是主矩阵(n大小)，rez是A的LU形式。; det_r1 = det_r1 * r1(i,i);cout<<A.determinant();
r1的行列式是主对角线元素的乘积
浏览 33提问于2020-07-20得票数 0
1回答
我是Rcpp的新手。我试图使用R包RcppEigen来得到矩阵的行列式。下面的代码保存在一个文件中，我使用sourceCpp来使用它。使用sourceCpp时没有编译错误。在R中使用getDeterminant(A)时，A是一个矩阵。它总是抱怨以下错误。#include <RcppEigen.h>
using Eigen::Map; // 'm
浏览 4提问于2015-02-02得票数 5
回答已采纳
我被要求使用值0:9找出2乘2矩阵的所有可能排列的行列式。然而，我似乎找不到所有可能的矩阵。我尝试过使用combn函数，但它只会返回具有唯一元素的矩阵(换句话说，它不会返回具有重复元素的矩阵，例如1,1,0,2)。有人能给我指个方向吗？ 我已经尝试过类似的东西。x <- 0:9
combn(x, 4) 最终结果应该是使用数字0:9的所有可能的2乘2矩阵的行列式的向量。
浏览 9提问于2019-09-07得票数 3
我正在尝试使用boost c++库来计算行列式。我找到了下面复制的函数InvertMatrix()的代码。每次我计算这个倒数时，我也想要行列式。我有一个很好的想法，通过将LU分解得到的U矩阵的对角线向下乘以来计算。有一个问题，我能够正确地计算行列式，除了符号。根据旋转的不同，有一半的时间我得到的符号是不正确的。= 0 ) return false;
这就是我在计算行列式时的最佳尝试。
浏览 0提问于2009-09-14得票数 5
回答已采纳
我正在尝试使用RcppArmadillo编写一个MCMC过程，它涉及到计算大约30,000 x 30,000个稀疏矩阵的对数行列式。看起来Armadillo中的log_det()目前不支持sp_mat，所以我这样做：// [[Rcpp::depends(RcppEigen::as<Eigen::SparseMatrix<double>>(Rcpp::wrap(arma_
浏览 0提问于2018-08-01得票数 1
2回答
我有一个问题，我被要求证明矩阵B的行列式等于0。矩阵B定义为： import numpy as np [ -7. -3.j -4. -9.j -3. -2.j]
[ 11. -3.j -16.-12.j 6. -5.j]] 这个行列式使用我用https://www.symbolab.com/
浏览 25提问于2020-01-30得票数 0
回答已采纳
1回答
Excel -确定排列的奇偶
、、、、
我正在处理一个Excel表格，在这里我需要确定一个大小为N的垂直数组的奇偶性。该数组包含每个从1到N的每个数字，每一次正好一次。例如，数组{3;1;2;4}具有偶数奇偶，因为它需要两个掉期(至少)才能转换为{1;2;3;4}，但总是需要偶数的掉期。见下文。1 --> 3 --> 24 4 4
另一个例子：{2;1;4;5;3}具有奇偶性，因为它需要
浏览 3提问于2017-12-20得票数 2
回答已采纳
1回答
矩阵的列序
、、、、
我想知道如何通过给定向量的索引来构造一个按列排列的矩阵。我知道如何在R中做这件事，但在RcppArmadillo中它不起作用。,] 1 3我使用RcppArmadillo尝试了以下几种方法：using namespace Rcpp;
Li
浏览 2提问于2018-01-16得票数 0
回答已采纳
我试着按以下方式进行转换，但这是行不通的：{
Eigen::Matrix3d m3 = m.cast();
浏览 0提问于2018-08-27得票数 0
回答已采纳
例如，2 1 4 3 5是恒等置换1 2 3 4 5的偶数排列，因为它需要偶数的交换才能从一个向量得到另一向量。这似乎是人们经常想要做的事情，而且可能有内置的功能。如果我能让Matlab给我一个置换矩阵，那么我就可以得到这个矩阵的行列式，但是我还没有想出如何去做，我也不知道是否有一个快速的方法。
浏览 2提问于2013-05-03得票数 3
回答已采纳
1回答
但是为什么我可以得到一个行列式(用det命令)，比如9.143243e-16，而我只处理只有整数的3x3矩阵提前谢谢你。
浏览 2提问于2013-10-18得票数 0
回答已采纳
我希望通过存储二维矩阵的内存来优化内存命中。我打算将2D矩阵折叠成一个一维连续块，并想知道按行或按列将数据存储为连续块是否更有意义。我所考虑的运算类型是较昂贵的运算，例如乘法和奇异值分解。考虑一个3x3矩阵[a21 a22 a23]逐行存储矩阵是否更有意义？对于一个3x3矩阵来说，这可能不是一个问题，因为这两种情况都需要加载三个块才能完成计算，而且这三个矩阵一次都可以很容易地装入我们的缓存中，但是当我们有非常大的矩阵</
浏览 2提问于2014-12-06得票数 1
以下Rcpp代码是生成相同编译错误的大得多的代码的最小可重现示例。似乎我不能将一个数值矩阵赋值给一个列表，然后将该列表赋值给另一个矩阵。#include <Rcpp.h>List return_a(NumericMatrix a, NumericMatrixonly returns the input matrix a re
浏览 9提问于2019-02-20得票数 3
回答已采纳
在我的RcppArmadillo项目中，我有一些矩阵(例如mat A，B，C;)和一个S4对象，例如D(来自R中的外部函数)。由于我需要在这些矩阵和D之间进行一些计算，所以我希望在RcppArmadillo中将"D“转换为合适的数据类型，例如arma::mat D。有可能吗？做这件事最好的方法是什么？这是一个类似的代码：// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp</e
浏览 3提问于2017-06-12得票数 1
回答已采纳
1回答
RcppEigen列标签
、、、
如何从Rcpp：：特征矩阵中提取列标签？例如，这里有一些R矩阵colnames(mat) <- c("col1","col2")void fun_NM(const Rcpp::NumericMatrixXd& Mat) Rcpp::CharacterVector c
浏览 4提问于2020-12-11得票数 1
回答已采纳
1回答
我有一个关于从矩阵到向量减法的问题。用户可以显式地给出索引矩阵(它是一个与M相同大小的矩阵，如果不需要条目，则为0，如果必须提取条目，则为1)。如果提供了索引矩阵，那么我们只需对其进行子集，如果没有提供索引矩阵(索引矩阵= NULL)，则使用type1 (取true或false)构建索引矩阵。只有两种类型的索引矩阵是可能的。我使用了中提供的细分技术// [[
浏览 4提问于2016-01-24得票数 1
我想知道如何使5_5矩阵行列式。我已经将1_1矩阵行列式转化为4_4矩阵行列式，但不知道如何用C语言构造5_5矩阵行列式。谢谢你对我的照顾。
浏览 4提问于2021-11-11得票数 0
++){ *((Mat+i*COL_SIZE) + j) = value;} 但是，当我使用Rcpp// [[Rcpp::export]] NumericMatrix xx(4, 5);# [3,] 102 106 110 114 118
# [4,] 103 107 111
浏览 11提问于2019-02-25得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
