Commit 8f738f10 authored by Sébastien Villemot's avatar Sébastien Villemot

Dynare++: fix check for QZ eigenvalue criterium

The criterium was previously incorrectly applied to the square absolute value
of eigenvalues. Rather apply it to the absolute value itself (as now done in
Dynare).

Ref #1632

(manually cherry picked from commit 8ebcd7ca)
parent 0f8c5ef2
Pipeline #658 failed with stages
in 332 minutes and 21 seconds
...@@ -22,7 +22,7 @@ to select (return true) the pairs for which $\alpha<\beta$. ...@@ -22,7 +22,7 @@ to select (return true) the pairs for which $\alpha<\beta$.
@<|order_eigs| function code@>= @<|order_eigs| function code@>=
lapack_int order_eigs(const double* alphar, const double* alphai, const double* beta) lapack_int order_eigs(const double* alphar, const double* alphai, const double* beta)
{ {
return (*alphar * *alphar + *alphai * *alphai < *beta * *beta * qz_criterium); return (*alphar * *alphar + *alphai * *alphai < *beta * *beta * qz_criterium * qz_criterium);
} }
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment