public inbox for gcc-bugs@sourceware.org
help / color / mirror / Atom feed
* [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f
@ 2005-01-29 13:48 billingd at gcc dot gnu dot org
2005-01-29 13:51 ` [Bug fortran/19693] " pinskia at gcc dot gnu dot org
` (9 more replies)
0 siblings, 10 replies; 11+ messages in thread
From: billingd at gcc dot gnu dot org @ 2005-01-29 13:48 UTC (permalink / raw)
To: gcc-bugs
With gfortran 20050129 I get a testsuite failure at -O2 that I don't see at -
O0. The failure is in the CGT tests in ctest.in:
CGT: General tridiagonal
Matrix types (1-6 have specified condition numbers):
1. Diagonal 7. Random, unspecified CNDNUM
2. Random, CNDNUM = 2 8. First column zero
3. Random, CNDNUM = sqrt(0.1/EPS) 9. Last column zero
4. Random, CNDNUM = 0.1/EPS 10. Last n/2 columns zero
5. Scaled near underflow 11. Scaled near underflow
6. Scaled near overflow 12. Scaled near overflow
Test ratios:
1: norm( L * U - A ) / ( N * norm(A) * EPS )
2: norm( B - A * X ) / ( norm(A) * norm(X) * EPS )
3: norm( X - XACT ) / ( norm(XACT) * CNDNUM * EPS )
4: norm( X - XACT ) / ( norm(XACT) * CNDNUM * EPS ), refined
5: norm( X - XACT ) / ( norm(XACT) * (error bound) )
6: (backward error) / EPS
7: RCOND * CNDNUM - 1.0
Messages:
NORM ='O', N = 1, type 5, test( 7) = 0.16777E+08
NORM ='I', N = 1, type 5, test( 7) = 0.16777E+08
TRANS='N', N = 1, NRHS= 1, type 5, test( 3) =
I have reduced it to the miscompilation of routine cgtts2.f at -O2. I have
isolated a single failing case, but if I try to change the routine the problem
evaporates.
Attached is:
- driv.f: a small driver program
- cgtts2.f: the LAPACK routine with the problem.
To reproduce:
gfortran -o pass.exe -O0 driv.f cgtts2.f
gfortran -o fail.exe -O0 driv.f cgtts2.f
./pass
B(1,1) = ( 4.0564819E+31, 0.000000 )
./fail
At line 15 of file driv.f
Internal Error: printf is broken
B(1,1) = (
Note: The internal error is PR 19363
--
Summary: optimization problem with LAPACK routine cgtts2.f
Product: gcc
Version: 4.0.0
Status: UNCONFIRMED
Severity: normal
Priority: P2
Component: fortran
AssignedTo: unassigned at gcc dot gnu dot org
ReportedBy: billingd at gcc dot gnu dot org
CC: gcc-bugs at gcc dot gnu dot org
GCC build triplet: i686-pc-cygwin
GCC host triplet: i686-pc-cygwin
GCC target triplet: i686-pc-cygwin
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
@ 2005-01-29 13:51 ` pinskia at gcc dot gnu dot org
2005-01-29 13:51 ` billingd at gcc dot gnu dot org
` (8 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: pinskia at gcc dot gnu dot org @ 2005-01-29 13:51 UTC (permalink / raw)
To: gcc-bugs
--
What |Removed |Added
----------------------------------------------------------------------------
Keywords| |wrong-code
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
2005-01-29 13:51 ` [Bug fortran/19693] " pinskia at gcc dot gnu dot org
@ 2005-01-29 13:51 ` billingd at gcc dot gnu dot org
2005-01-29 13:52 ` billingd at gcc dot gnu dot org
` (7 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: billingd at gcc dot gnu dot org @ 2005-01-29 13:51 UTC (permalink / raw)
To: gcc-bugs
------- Additional Comments From billingd at gcc dot gnu dot org 2005-01-29 13:51 -------
Created an attachment (id=8099)
--> (http://gcc.gnu.org/bugzilla/attachment.cgi?id=8099&action=view)
driv.f: test driver
driv.f is driver routine
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
2005-01-29 13:51 ` [Bug fortran/19693] " pinskia at gcc dot gnu dot org
2005-01-29 13:51 ` billingd at gcc dot gnu dot org
@ 2005-01-29 13:52 ` billingd at gcc dot gnu dot org
2005-01-29 13:55 ` pinskia at gcc dot gnu dot org
` (6 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: billingd at gcc dot gnu dot org @ 2005-01-29 13:52 UTC (permalink / raw)
To: gcc-bugs
------- Additional Comments From billingd at gcc dot gnu dot org 2005-01-29 13:52 -------
Created an attachment (id=8100)
--> (http://gcc.gnu.org/bugzilla/attachment.cgi?id=8100&action=view)
cgtts2.f is LAPACK routine that is miscompiled
cgtts2.f is LAPACK routine that is miscompiled
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (2 preceding siblings ...)
2005-01-29 13:52 ` billingd at gcc dot gnu dot org
@ 2005-01-29 13:55 ` pinskia at gcc dot gnu dot org
2005-01-30 12:15 ` billingd at gcc dot gnu dot org
` (5 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: pinskia at gcc dot gnu dot org @ 2005-01-29 13:55 UTC (permalink / raw)
To: gcc-bugs
------- Additional Comments From pinskia at gcc dot gnu dot org 2005-01-29 13:55 -------
At -O0 on powerpc-darwin I get:
B(1,1) = ( +Infinity, NaN)
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (3 preceding siblings ...)
2005-01-29 13:55 ` pinskia at gcc dot gnu dot org
@ 2005-01-30 12:15 ` billingd at gcc dot gnu dot org
2005-01-30 13:57 ` pinskia at gcc dot gnu dot org
` (4 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: billingd at gcc dot gnu dot org @ 2005-01-30 12:15 UTC (permalink / raw)
To: gcc-bugs
------- Additional Comments From billingd at gcc dot gnu dot org 2005-01-30 12:15 -------
This goes away when we set flag_complex_divide_method = 1 in toplev.c -
overflow in complex division - which is consistent with the magnitude of the
inputs.
PR 18902 again.
--
What |Removed |Added
----------------------------------------------------------------------------
BugsThisDependsOn| |18902
OtherBugsDependingO| |5900
nThis| |
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (4 preceding siblings ...)
2005-01-30 12:15 ` billingd at gcc dot gnu dot org
@ 2005-01-30 13:57 ` pinskia at gcc dot gnu dot org
2005-02-24 14:28 ` rth at gcc dot gnu dot org
` (3 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: pinskia at gcc dot gnu dot org @ 2005-01-30 13:57 UTC (permalink / raw)
To: gcc-bugs
--
What |Removed |Added
----------------------------------------------------------------------------
Status|UNCONFIRMED |NEW
Ever Confirmed| |1
Last reconfirmed|0000-00-00 00:00:00 |2005-01-30 13:57:14
date| |
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (5 preceding siblings ...)
2005-01-30 13:57 ` pinskia at gcc dot gnu dot org
@ 2005-02-24 14:28 ` rth at gcc dot gnu dot org
2005-03-01 20:22 ` tobi at gcc dot gnu dot org
` (2 subsequent siblings)
9 siblings, 0 replies; 11+ messages in thread
From: rth at gcc dot gnu dot org @ 2005-02-24 14:28 UTC (permalink / raw)
To: gcc-bugs
--
Bug 19693 depends on bug 18902, which changed state.
Bug 18902 Summary: Naive (default) complex division algorithm
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=18902
What |Old Value |New Value
----------------------------------------------------------------------------
Status|ASSIGNED |RESOLVED
Resolution| |FIXED
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (6 preceding siblings ...)
2005-02-24 14:28 ` rth at gcc dot gnu dot org
@ 2005-03-01 20:22 ` tobi at gcc dot gnu dot org
2005-03-02 0:09 ` billingd at gcc dot gnu dot org
2005-04-20 2:36 ` pinskia at gcc dot gnu dot org
9 siblings, 0 replies; 11+ messages in thread
From: tobi at gcc dot gnu dot org @ 2005-03-01 20:22 UTC (permalink / raw)
To: gcc-bugs
------- Additional Comments From tobi at gcc dot gnu dot org 2005-03-01 20:22 -------
Can this be closed with the recent change of default complex divide algorithm?
--
What |Removed |Added
----------------------------------------------------------------------------
CC| |tobi at gcc dot gnu dot org
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (7 preceding siblings ...)
2005-03-01 20:22 ` tobi at gcc dot gnu dot org
@ 2005-03-02 0:09 ` billingd at gcc dot gnu dot org
2005-04-20 2:36 ` pinskia at gcc dot gnu dot org
9 siblings, 0 replies; 11+ messages in thread
From: billingd at gcc dot gnu dot org @ 2005-03-02 0:09 UTC (permalink / raw)
To: gcc-bugs
------- Additional Comments From billingd at gcc dot gnu dot org 2005-03-02 00:07 -------
Fixed by changes to complex divide algorithm.
See comment #46 on PR5900 for LAPACK testsuite results at -O2.
--
What |Removed |Added
----------------------------------------------------------------------------
Status|NEW |RESOLVED
Resolution| |FIXED
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
* [Bug fortran/19693] optimization problem with LAPACK routine cgtts2.f
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
` (8 preceding siblings ...)
2005-03-02 0:09 ` billingd at gcc dot gnu dot org
@ 2005-04-20 2:36 ` pinskia at gcc dot gnu dot org
9 siblings, 0 replies; 11+ messages in thread
From: pinskia at gcc dot gnu dot org @ 2005-04-20 2:36 UTC (permalink / raw)
To: gcc-bugs
--
What |Removed |Added
----------------------------------------------------------------------------
Target Milestone|--- |4.0.0
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=19693
^ permalink raw reply [flat|nested] 11+ messages in thread
end of thread, other threads:[~2005-04-20 2:35 UTC | newest]
Thread overview: 11+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2005-01-29 13:48 [Bug fortran/19693] New: optimization problem with LAPACK routine cgtts2.f billingd at gcc dot gnu dot org
2005-01-29 13:51 ` [Bug fortran/19693] " pinskia at gcc dot gnu dot org
2005-01-29 13:51 ` billingd at gcc dot gnu dot org
2005-01-29 13:52 ` billingd at gcc dot gnu dot org
2005-01-29 13:55 ` pinskia at gcc dot gnu dot org
2005-01-30 12:15 ` billingd at gcc dot gnu dot org
2005-01-30 13:57 ` pinskia at gcc dot gnu dot org
2005-02-24 14:28 ` rth at gcc dot gnu dot org
2005-03-01 20:22 ` tobi at gcc dot gnu dot org
2005-03-02 0:09 ` billingd at gcc dot gnu dot org
2005-04-20 2:36 ` pinskia at gcc dot gnu dot org
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).