We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 4dac1b2 commit 3630e02Copy full SHA for 3630e02
src/high-level/matrix-functions/hermitian-eig.lisp
@@ -10,7 +10,9 @@
10
(let* ((ajk2 (* ajk (conjugate ajk)))
11
(delta (/ (- ajj akk) 2))
12
(d2 (* delta delta))
13
- (sgnd (if (minusp delta) -1d0 1d0)) ; we do NOT want SIGNUM here
+ ;; XXX: DELTA should be real here? It's possibly of type
14
+ ;; COMPLEX.
15
+ (sgnd (if (minusp (realpart delta)) -1d0 1d0)) ; we do NOT want SIGNUM here
16
)
17
(- akk
18
(/ ajk2
0 commit comments