v4.1: adding penalty when Jacobian is complex

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2522 ac1d8469-bf42-47a9-8791-bf33cf982152
2009-03-26 09:03:08 +00:00 · 2009-03-26 09:03:08 +00:00 · ff473fedfd
parent e2a0ecba86
commit ff473fedfd
1 changed files with 1 additions and 0 deletions
--- a/matlab/dr1.m
+++ b/matlab/dr1.m
@ -210,6 +210,7 @@ function [dr,info,M_,options_,oo_] = dr1(dr,task,M_,options_,oo_)
            jacobia_ = real(jacobia_);
        else
            info(1) = 6;
+            info(2) = sum(sum(imag(jacobia_).^2));
            return
        end
    end