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

Bytecode: fix simulation of blocks “Solve forward/backward complete” with stack_solve_algo={0,1,4} 

The problem would also manifest at steady state computation with solve_algo=6.

It could crash or give incorrect results.
trustregion
Sébastien Villemot 2022-02-11 11:43:14 +01:00
parent 9a27735f47
commit 188df7711d
No known key found for this signature in database
GPG Key ID: 2CECE9350ECEBE4A
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
mex/sources/bytecode/SparseMatrix.cc
View File

@ -3908,7 +3908,7 @@ dynSparseMatrix::Simulate_One_Boundary(int block_num, int y_size, int y_kmin, in
else if ((solve_algo == 8 && steady_state) || (stack_solve_algo == 3 && !steady_state))
Solve_Matlab_BiCGStab(A_m, b_m, size, slowc, block_num, false, it_, x0_m, preconditioner);
else if ((solve_algo == 6 && steady_state) || ((stack_solve_algo == 0 || stack_solve_algo == 1 || stack_solve_algo == 4) && !steady_state))
Solve_LU_UMFPack(Ap, Ai, Ax, b, size, size, slowc, true, 0);
Solve_LU_UMFPack(Ap, Ai, Ax, b, size, size, slowc, false, it_);
}
return singular_system;
}