[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #63785] float compared to 0 in __ichol__.cc
From: |
Dmitri A. Sergatskov |
Subject: |
[Octave-bug-tracker] [bug #63785] float compared to 0 in __ichol__.cc |
Date: |
Sun, 12 Feb 2023 15:23:53 -0500 (EST) |
Follow-up Comment #2, bug #63785 (project octave):
May be this is a gcc bug. I played with code:
for ii = 1:100 a = randn(N) + i*randn(N); b = R*sparse(a'*a); ii, ichol(b);
endfor
(with different N and R)
and situation gets progressively worth for the larger N.
In clang compiled code (with the same optimizations)
Im part is strictly 0.
E.g. gcc (N =4, R=1):
ii = 29
Real: 5.37065 Im: 0
Real: 5.05964 Im: 8.51589e-19
Real: 7.94997 Im: 5.10252e-17
Real: 0.0602546 Im: 4.3797e-16
error: ichol: non-real pivot encountered. The matrix must be Hermitian.
error: called from
ichol at line 246 column 9
Dmitri.
--
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?63785>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/