stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Block decomposition: remove unused 2nd and 3rd derivatives in output of dynamic functions

time-shift
Sébastien Villemot 2020-06-03 16:59:59 +02:00
parent 0d11f8182a
commit 2ec5e4075a
No known key found for this signature in database
GPG Key ID: 2CECE9350ECEBE4A
3 changed files with 5 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
matlab/solve_one_boundary.m
View File

@ -72,8 +72,6 @@ function [y, T, oo_, info] = solve_one_boundary(fname, y, x, params, steady_stat
Blck_size=size(y_index_eq,2);
g2 = [];
g3 = [];
correcting_factor=0.01;
ilu_setup.droptol=1e-10;
max_resa=1e100;
@ -94,7 +92,7 @@ for it_=start:incr:finish
g1=spalloc( Blck_size, Blck_size, nze);
while ~(cvg==1 || iter>maxit_)
if is_dynamic
[r, y, T, g1, g2, g3] = feval(fname, y, x, params, steady_state, T, it_, false);
[r, y, T, g1] = feval(fname, y, x, params, steady_state, T, it_, false);
else
[r, y, T, g1] = feval(fname, y, x, params, T);
end
@ -264,8 +262,8 @@ for it_=start:incr:finish
y(y_index_eq) = phat;
end
if is_dynamic
[r, y, T, g1, g2, g3] = feval(fname, y, x, params, ...
steady_state, T, it_, false);
[r, y, T, g1] = feval(fname, y, x, params, ...
steady_state, T, it_, false);
else
[r, y, T, g1] = feval(fname, y, x, params, T);
end

4
matlab/solve_two_boundaries.m
View File

@ -69,8 +69,6 @@ verbose = options.verbosity;
cvg=0;
iter=0;
Per_u_=0;
g2 = [];
g3 = [];
Blck_size=size(y_index,2);
correcting_factor=0.01;
ilu_setup.droptol=1e-10;
@ -84,7 +82,7 @@ Jacobian_Size=Blck_size*(y_kmin+y_kmax_l +periods);
g1=spalloc( Blck_size*periods, Jacobian_Size, nze*periods);
reduced = 0;
while ~(cvg==1 || iter>maxit_)
[r, y, T, g1, g2, g3, b]=feval(fname, y, x, params, steady_state, T, periods, false, y_kmin, Blck_size,options.periods);
[r, y, T, g1, b]=feval(fname, y, x, params, steady_state, T, periods, false, y_kmin, Blck_size,options.periods);
preconditioner = 2;
g1a=g1(:, y_kmin*Blck_size+1:(periods+y_kmin)*Blck_size);
term1 = g1(:, 1:y_kmin_l*Blck_size)*reshape(y(1+y_kmin-y_kmin_l:y_kmin,y_index)',1,y_kmin_l*Blck_size)';

2
preprocessor

@ -1 +1 @@
Subproject commit f0c8a145518cb86953a44fca4ba2b0ff7c88cf0a
Subproject commit 050f1cdd74c92eb7c65c7151980ce5c411fed378