move rungibbs inside getspillover

ygeunkim · Jan 19, 2025 · eb8d72a · eb8d72a
1 parent bd9f332
commit eb8d72a
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/inst/include/bvharspillover.h b/inst/include/bvharspillover.h
@@ -410,6 +410,7 @@ class DynamicLdltSpillover {
 		model[window][chain].reset();
 	}
 	void getSpillover(int window, int chain) {
+		runGibbs(window, chain);
 		spillover[window][chain]->computeSpillover();
 		to_sp[window][chain] = spillover[window][chain]->returnTo();
 		from_sp[window][chain] = spillover[window][chain]->returnFrom();
@@ -423,7 +424,6 @@ class DynamicLdltSpillover {
 			#pragma omp parallel for num_threads(nthreads)
 		#endif
 			for (int window = 0; window < num_horizon; ++window) {
-				runGibbs(window, 0);
 				getSpillover(window, 0);
 			}
 		} else {
@@ -432,7 +432,6 @@ class DynamicLdltSpillover {
 		#endif
 			for (int window = 0; window < num_horizon; ++window) {
 				for (int chain = 0; chain < num_chains; ++chain) {
-					runGibbs(window, chain);
 					getSpillover(window, chain);
 				}
 			}