LAPACK Archives

[Lapack] Bug in lapack-3.2.1/TESTING/EIG/zdrvgbx.f


Dear Michael, thanks a lot. Attached the patch for ZDRVGBX. The "same" 
patch did apply to the three other precisions. The patch is committed in 
our svn repository and will be released with 3.2.2, and is referenced as 
bug 41 at http://www.netlib.org/lapack/Errata/ . Best wishes, Julien.




============================================================================
Index: LIN/zdrvgbx.f
===================================================================
--- LIN/zdrvgbx.f       (revision 695)
+++ LIN/zdrvgbx.f       (working copy)
@@ -828,11 +828,11 @@
                                IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
       $                           CALL ALADHD( NOUT, PATH )
                                IF( PREFAC ) THEN
-                                 WRITE( NOUT, FMT = 9997 )'ZGBSVXX',
+                                 WRITE( NOUT, FMT = 9995 )'ZGBSVXX',
       $                                FACT, TRANS, N, KL, KU, EQUED,
       $                                IMAT, K, RESULT( K )
                                ELSE
-                                 WRITE( NOUT, FMT = 9998 )'ZGBSVXX',
+                                 WRITE( NOUT, FMT = 9996 )'ZGBSVXX',
       $                                FACT, TRANS, N, KL, KU, IMAT, K,
       $                                RESULT( K )
                                END IF
@@ -846,11 +846,11 @@
                             IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
       $                        CALL ALADHD( NOUT, PATH )
                             IF( PREFAC ) THEN
-                              WRITE( NOUT, FMT = 9997 )'ZGBSVXX', FACT,
+                              WRITE( NOUT, FMT = 9995 )'ZGBSVXX', FACT,
       $                             TRANS, N, KL, KU, EQUED, IMAT, 1,
       $                             RESULT( 1 )
                             ELSE
-                              WRITE( NOUT, FMT = 9998 )'ZGBSVXX', FACT,
+                              WRITE( NOUT, FMT = 9996 )'ZGBSVXX', FACT,
       $                             TRANS, N, KL, KU, IMAT, 1,
       $                             RESULT( 1 )
                             END IF
@@ -861,11 +861,11 @@
                             IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
       $                        CALL ALADHD( NOUT, PATH )
                             IF( PREFAC ) THEN
-                              WRITE( NOUT, FMT = 9997 )'ZGBSVXX', FACT,
+                              WRITE( NOUT, FMT = 9995 )'ZGBSVXX', FACT,
       $                             TRANS, N, KL, KU, EQUED, IMAT, 6,
       $                             RESULT( 6 )
                             ELSE
-                              WRITE( NOUT, FMT = 9998 )'ZGBSVXX', FACT,
+                              WRITE( NOUT, FMT = 9996 )'ZGBSVXX', FACT,
       $                             TRANS, N, KL, KU, IMAT, 6,
       $                             RESULT( 6 )
                             END IF
@@ -876,11 +876,11 @@
                             IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
       $                        CALL ALADHD( NOUT, PATH )
                             IF( PREFAC ) THEN
-                              WRITE( NOUT, FMT = 9997 )'ZGBSVXX', FACT,
+                              WRITE( NOUT, FMT = 9995 )'ZGBSVXX', FACT,
       $                             TRANS, N, KL, KU, EQUED, IMAT, 7,
       $                             RESULT( 7 )
                             ELSE
-                              WRITE( NOUT, FMT = 9998 )'ZGBSVXX', FACT,
+                              WRITE( NOUT, FMT = 9996 )'ZGBSVXX', FACT,
       $                             TRANS, N, KL, KU, IMAT, 7,
       $                             RESULT( 7 )
                             END IF
============================================================================



On Tue, 29 Sep 2009, Chuvelev, Michael wrote:

Hi,
?
please pay attention to the printing format 9998 misuse in zdrvgbx.f:
?
????????????????????????????? WRITE( NOUT, FMT = 9998 )'ZGBSVXX', FACT,
???? $???????????????????????????? TRANS, N, KL, KU, IMAT, 7,
???? $???????????????????????????? RESULT( 7 )
...
?
9998 FORMAT( ' *** In ZDRVGB, LAFB=', I5, ' is too small for N=', I5,
???? $????? ', KU=', I5, ', KL=', I5, /
???? $????? ' ==> Increase LAFB to at least ', I5 )
?
This would possibly lead to a segmentation fault in case of errors in zgbsvxx.
?
The same issue concerns other precisions too.
?
Best regards,
?
Michael.
?
?

 --------------------------------------------------------------------
Closed Joint Stock Company Intel A/O
Registered legal address: Krylatsky Hills Business Park, 
17 Krylatskaya Str., Bldg 4, Moscow 121614, 
Russian Federation

This e-mail and any attachments may contain confidential material for
the sole use of the intended recipient(s). Any review or distribution
by others is strictly prohibited. If you are not the intended
recipient, please contact the sender and delete all copies.



<Prev in Thread] Current Thread [Next in Thread>


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or