throwing exception if dgges reports an error

time-shift
Michel Juillard 2016-07-17 14:40:15 +02:00
parent 5d38280f27
commit e7a1988223
1 changed files with 4 additions and 0 deletions

View File

@ -198,6 +198,10 @@ the same. The difference is only numerical error.
matD.getData().base(), &n, &sdim2, alphar.base(), alphai.base(),
beta.base(), vsl.getData().base(), &n, vsr.getData().base(), &n,
work.base(), &lwork, bwork, &info);
if(info){
throw KordException(__FILE__,__LINE__,
"DGGES returns an error in FirstOrder::solve");
}
sdim = sdim2;
bk_cond = (sdim == ypart.nys());
delete[] bwork;