From 455b3d0a75fc4c696f98da3afaa1596507f93cff Mon Sep 17 00:00:00 2001
From: adamwang15 <adamwang15@gmail.com>
Date: Sat, 10 Aug 2024 13:02:53 +1000
Subject: [PATCH] debug

---
 src/restrictions_zero.cpp | 14 ++++++++++----
 1 file changed, 10 insertions(+), 4 deletions(-)

diff --git a/src/restrictions_zero.cpp b/src/restrictions_zero.cpp
index d553359..2370545 100644
--- a/src/restrictions_zero.cpp
+++ b/src/restrictions_zero.cpp
@@ -42,7 +42,7 @@ 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;
   
@@ -50,12 +50,18 @@ arma::colvec zero_restrictions(
   
   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;
 }