fixed the bug with as_scalar #1

bsvars · Jan 18, 2024 · 710e88a · 710e88a
1 parent ad66aab
commit 710e88a
Show file tree

Hide file tree

Showing 3 changed files with 2 additions and 1 deletion.
diff --git a/src/bsvarSIGNs.so b/src/bsvarSIGNs.so
diff --git a/src/bsvars_in_sign.cpp b/src/bsvars_in_sign.cpp
@@ -1,6 +1,7 @@
 
 #include <RcppArmadillo.h>
 #include "progress.hpp"
+#include "Rcpp/Rmath.h"
 
 using namespace Rcpp;
 using namespace arma;
@@ -91,7 +92,7 @@ void sample_hyperparameters (
     shape_tmp         = prior_hyper_a_B + prior_hyper_nu_B / 2;
     aux_hyper(N + n, 0) = R::rgamma(shape_tmp, scale_tmp);
 
-    scale_tmp         = aux_hyper(N + n, 0) + as_scalar(aux_B * prior_B_V_inv * aux_B.t());
+    scale_tmp         = aux_hyper(N + n, 0) + as_scalar(aux_B.row(n) * prior_B_V_inv * aux_B.row(n).t());
     shape_tmp         = prior_hyper_nu_B + rn;
     aux_hyper(n, 0)   = scale_tmp / R::rchisq(shape_tmp);
 

diff --git a/src/bsvars_in_sign.o b/src/bsvars_in_sign.o