我正在使用Rcpp和RcppArmadillo编写函数,我收到以下错误:
error: copy into submatrix: incompatible matrix dimensions: 1x40 and 40x1
哪些RcppArmadillo操作或功能可能导致此错误消息?我很难跟踪产生此错误的内容。
我基本上只是用矩阵行和列将矩阵与矩阵和矩阵相乘。
修改
该功能在https://github.com/jdanielnd/bfm/blob/master/src/bfm_gibbs.cpp
处可用x数据矩阵由https://github.com/jdanielnd/bfm/blob/master/R/generate_data.r
处的函数生成library(Rcpp)
source("generate_data.r")
sourceCpp("bfm_gibbs.cpp")
x <- sample_data()
res <- bfmu_c(x, ite = 100)