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

Fix bug in linsolve when second argument is a zero-column matrix 

Fixes crash of 1st-order purely forward models under Octave
time-shift
Sébastien Villemot 2013-08-15 12:19:51 +02:00
parent 4a0e674029
commit c34c90e8cb
1 changed files with 7 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
mex/sources/linsolve/linsolve.cc
View File

@ -10,7 +10,7 @@
*/
/*
* Copyright (C) 2012 Dynare Team
* Copyright (C) 2012-2013 Dynare Team
*
* This file is part of Dynare.
*
@ -114,7 +114,12 @@ If requested, @var{r} will contain the reciprocal condition number.\n\
retval(0) = A.solve(typA, B, info, rcond, NULL, true, transa ? blas_trans : blas_no_trans);
if (nargout == 2)
retval(1) = rcond;
{
if (dimB(1) > 0)
retval(1) = rcond;
else // If B has zero columns, A.solve() apparently does not compute A's rcond
retval(1) = A.rcond();
}
return retval;
}