Testsuite / PAC: use sparse representation for the model

Ref. #1859

tests/pac/trend-component-14/example1.mod

@@ -87,7 +87,8 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = example1.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
-  [residual, g1] = example1.dynamic(y, x', M_.params, oo_.steady_state, 1);
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = example1.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = example1.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   save('example1.mat', 'residual', 'g1', 'TrueData');
-end;
+end;
\ No newline at end of file

tests/pac/trend-component-14/substitution.mod

@@ -78,9 +78,10 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = substitution.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = substitution.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = substitution.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   example1 = load('example1.mat');
-  [residual, g1] = substitution.dynamic(y, x', M_.params, oo_.steady_state, 1);
 end;
 
 if max(abs(example1.TrueData.data(:)-TrueData.data(:)))>1e-9

tests/pac/var-5/example1.mod

@@ -67,7 +67,8 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = example1.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
-  [residual, g1] = example1.dynamic(y, x', M_.params, oo_.steady_state, 1);
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = example1.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = example1.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   save('example1.mat', 'residual', 'g1', 'TrueData');
 end;

tests/pac/var-5/substitution.mod

@@ -59,9 +59,10 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = substitution.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = substitution.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = substitution.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   example1 = load('example1.mat');
-  [residual, g1] = substitution.dynamic(y, x', M_.params, oo_.steady_state, 1);
 end;
 
 if max(abs(example1.TrueData.data(:)-TrueData.data(:)))>1e-9

tests/pac/var-6/example1.mod

@@ -71,7 +71,8 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = example1.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
-  [residual, g1] = example1.dynamic(y, x', M_.params, oo_.steady_state, 1);
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = example1.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = example1.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   save('example1.mat', 'residual', 'g1', 'TrueData');
 end;

tests/pac/var-6/substitution.mod

@@ -64,9 +64,10 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = substitution.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = substitution.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = substitution.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   example1 = load('example1.mat');
-  [residual, g1] = substitution.dynamic(y, x', M_.params, oo_.steady_state, 1);
 
   if max(abs(example1.TrueData.data(:)-TrueData.data(:)))>1e-9
     error('Simulations do not match.')

tests/pac/var-7/example1.mod

@@ -73,7 +73,8 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = example1.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
-  [residual, g1] = example1.dynamic(y, x', M_.params, oo_.steady_state, 1);
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = example1.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = example1.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   save('example1.mat', 'residual', 'g1', 'TrueData');
 end;

tests/pac/var-7/substitution.mod

@@ -67,9 +67,10 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = substitution.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = substitution.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = substitution.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   example1 = load('example1.mat');
-  [residual, g1] = substitution.dynamic(y, x', M_.params, oo_.steady_state, 1);
 
   if max(abs(example1.TrueData.data(:)-TrueData.data(:)))>1e-9
     error('Simulations do not match.')

tests/pac/var-8/example1.mod

@@ -84,7 +84,8 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = example1.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
-  [residual, g1] = example1.dynamic(y, x', M_.params, oo_.steady_state, 1);
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = example1.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = example1.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   save('example1.mat', 'residual', 'g1', 'TrueData');
 end;

tests/pac/var-8/substitution.mod

@@ -64,9 +64,10 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = substitution.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = substitution.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = substitution.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   example1 = load('example1.mat');
-  [residual, g1] = substitution.dynamic(y, x', M_.params, oo_.steady_state, 1);
 
   if max(abs(example1.TrueData.data(:)-TrueData.data(:)))>1e-9
     error('Simulations do not match.')

tests/pac/var-9/example1.mod

@@ -88,7 +88,8 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = example1.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
-  [residual, g1] = example1.dynamic(y, x', M_.params, oo_.steady_state, 1);
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = example1.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = example1.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   save('example1.mat', 'residual', 'g1', 'TrueData');
 end;

tests/pac/var-9/substitution.mod

@@ -65,9 +65,10 @@ verbatim;
   y(1:M_.orig_endo_nbr) = rand(M_.orig_endo_nbr, 1);
   x = randn(M_.exo_nbr,1);
   y = substitution.set_auxiliary_variables(y, x, M_.params);
-  y = [y(find(M_.lead_lag_incidence(1,:))); y];
+  y3n = repmat(y, 3, 1);
+  [residual, T_order, T] = substitution.sparse.dynamic_resid(y3n, x', M_.params, oo_.steady_state);
+  g1 = substitution.sparse.dynamic_g1(y3n, x', M_.params, oo_.steady_state, M_.dynamic_g1_sparse_rowval, M_.dynamic_g1_sparse_colval, M_.dynamic_g1_sparse_colptr, T_order, T);
   example1 = load('example1.mat');
-  [residual, g1] = substitution.dynamic(y, x', M_.params, oo_.steady_state, 1);
 
   if max(abs(example1.TrueData.data(:)-TrueData.data(:)))>1e-9
     error('Simulations do not match.')