Bug in k_order_solver? ghx changes from order 1/2 to 3
Posted: Fri Aug 20, 2010 2:45 pm
Dear all,
I compared in two models the solutions using different solvers and approximation orders. Unfortunately, there are significant differences that should not occur from my point of view.
What I found is the following: The first order solution, i.e. in particular dr.ghx, is the same whether I use the k_order_solver option for order 1 or 2 or if I use the standard solver. However, independent of the solver I use, when I move from 2nd to 3rd order, the first order solution (i.e. dr.ghx) significantly changes. From my understanding this should not be the case as the first term of a Taylor approximation should be independent of higher order terms.
In contrast, the second order terms dr.ghxx and dr.ghxu are the same for 2nd and 3rd order approximations.
The above problem might also be the reason for my post about the strange error message " Caught Kord exception: NaN or Inf asserted in first order derivatives in FirstOrder::solve" (http://www.dynare.org/phpBB3/viewtopic.php?f=1&t=2720) as the problem obviously only occurs in the first order derivatives.
The problem occured with the unstable Dynare version of 08/20/2010 under Matlab 2010a 32bit with MSVC10.
I compared in two models the solutions using different solvers and approximation orders. Unfortunately, there are significant differences that should not occur from my point of view.
What I found is the following: The first order solution, i.e. in particular dr.ghx, is the same whether I use the k_order_solver option for order 1 or 2 or if I use the standard solver. However, independent of the solver I use, when I move from 2nd to 3rd order, the first order solution (i.e. dr.ghx) significantly changes. From my understanding this should not be the case as the first term of a Taylor approximation should be independent of higher order terms.
In contrast, the second order terms dr.ghxx and dr.ghxu are the same for 2nd and 3rd order approximations.
The above problem might also be the reason for my post about the strange error message " Caught Kord exception: NaN or Inf asserted in first order derivatives in FirstOrder::solve" (http://www.dynare.org/phpBB3/viewtopic.php?f=1&t=2720) as the problem obviously only occurs in the first order derivatives.
The problem occured with the unstable Dynare version of 08/20/2010 under Matlab 2010a 32bit with MSVC10.