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.
|