debug

bsvars · Aug 10, 2024 · 455b3d0 · 455b3d0
1 parent 08412a5
commit 455b3d0
Showing 1 changed file with 10 additions and 4 deletions.
diff --git a/src/restrictions_zero.cpp b/src/restrictions_zero.cpp
@@ -42,20 +42,26 @@ arma::colvec zero_restrictions(
 
   Rcout << "debug 2" << std::endl;
 
-  field<mat> ZF = ZIRF(Z, inv(A0.t()));
+  // field<mat> ZF = ZIRF(Z, inv(A0.t()));
 
   Rcout << "debug 3" << std::endl;
 
   vec z = vec();
 
   Rcout << "debug 4" << std::endl;
 
-  for (int j=0; j<ZF.n_elem; j++) {
-    z = join_vert(z, ZF(j).col(j));
+  mat irf_0 = inv(A0.t());
+
+  for (int j=0; j<Z.n_elem; j++) {
+    mat ZF = Z(j) * irf_0;
     Rcout << "debug 5" << std::endl;
+    vec ZF_j = ZF.col(j);
+    Rcout << "debug 6" << std::endl;
+    z = join_vert(z, ZF_j);
+    Rcout << "debug 7" << std::endl;
   }
 
-  Rcout << "debug 6" << std::endl;
+  Rcout << "debug 8" << std::endl;
 
   return z;
 }