Correct setting of info in dr_block.m if eigenvalue is close to 0/0

Must be the same as in dyn_first_order_solver.m
time-shift
Johannes Pfeifer 2016-05-05 10:08:52 +02:00
parent bebcfffe04
commit c6f705a8b9
1 changed files with 8 additions and 2 deletions

View File

@ -446,8 +446,14 @@ for i = 1:Size;
disp(data(i).eigval);
end;
if info1
info(1) = 2;
info(2) = info1;
if info1 == -30
% one eigenvalue is close to 0/0
info(1) = 7;
else
info(1) = 2;
info(2) = info1;
info(3) = size(E,2);
end
return
end
nba = nd-sdim;