Displaying 1 result from an estimated 1 matches for "pack_rh".
Did you mean:
pack_rhs
2024 Jun 10
1
changes in R-devel and zero-extent objects in Rcpp
...ctly where
it happens: line 100 of
RcppEigen/inst/include/Eigen/src/Core/products/GeneralMatrixMatrixTriangular.h:
for(Index k2=0; k2<depth; k2+=kc)
{
const Index actual_kc = (std::min)(k2+kc,depth)-k2;
// note that the actual rhs is the transpose/adjoint of mat
pack_rhs(blockB, rhs.getSubMapper(k2,0), actual_kc, size);
^^^^^^^^^^^^^^^^^^^^^^
where 'rhs' is an object wrapping the pointer with a method getSubMapper(i, j)
for accessing the data like a matrix. In the first loop iteration, you access
rhs[0]; there is no defensive test...