From 9cf9c0591134c3b8738f9864abdb2b6951ba022b Mon Sep 17 00:00:00 2001 From: Johannes Pfeifer Date: Thu, 31 Aug 2017 11:26:41 +0200 Subject: [PATCH] lyapunov_solver: Eliminate confusing transpose The transpose was copied from `dsge_likelihood.m` in https://github.com/DynareTeam/dynare/commit/2c5b1fed2d43ffdeb60207295e0a03610864b487, while the calls in all other functions do not have it. As Q is symmetric, it does not matter --- matlab/lyapunov_solver.m | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/matlab/lyapunov_solver.m b/matlab/lyapunov_solver.m index 98d782cdf..2bafb065a 100644 --- a/matlab/lyapunov_solver.m +++ b/matlab/lyapunov_solver.m @@ -41,7 +41,7 @@ function P=lyapunov_solver(T,R,Q,DynareOptions) % --*-- Unitary tests --*-- % along with Dynare. If not, see . if DynareOptions.lyapunov_fp == 1 - P = lyapunov_symm(T,R*Q'*R',DynareOptions.lyapunov_fixed_point_tol,DynareOptions.qz_criterium,DynareOptions.lyapunov_complex_threshold, 3, DynareOptions.debug); + P = lyapunov_symm(T,R*Q*R',DynareOptions.lyapunov_fixed_point_tol,DynareOptions.qz_criterium,DynareOptions.lyapunov_complex_threshold, 3, DynareOptions.debug); elseif DynareOptions.lyapunov_db == 1 [P, errorflag] = disclyap_fast(T,R*Q*R',DynareOptions.lyapunov_doubling_tol); if errorflag %use Schur-based method @@ -183,4 +183,4 @@ end %$ end %$ %$ T = all(t); -%@eof:1 \ No newline at end of file +%@eof:1