LAPACK Archives

[Lapack] possible bugfix for floating point exception in LAPACK 3.1.0 te

Hello,

the lapack 3.1.0 testing breaks on a CRAY with a floating point 
exception. I suppose the problem is an uninitialized asignment in cunmrz.f.

So i am proposing the following patch, which seems to fix the problem:

kite:/ha/rull/lapack-3.1.0/SRC$ cat cunmrz.f.patch
*** cunmrz.f.orig       Tue Nov 14 01:32:02 2006
--- cunmrz.f    Thu Jan 25 16:39:41 2007
***************
*** 171,178 ****
        END IF
  *
        IF( INFO.EQ.0 ) THEN
           IF( M.EQ.0 .OR. N.EQ.0 ) THEN
!             LWKOPT = 1
  *
  *           Determine the block size.  NB may be at most NBMAX, where
  *           NBMAX is used to define the local array T.
--- 171,179 ----
        END IF
  *
        IF( INFO.EQ.0 ) THEN
+          LWKOPT = 1
           IF( M.EQ.0 .OR. N.EQ.0 ) THEN
! ******            LWKOPT = 1
  *
  *           Determine the block size.  NB may be at most NBMAX, where
  *           NBMAX is used to define the local array T.
kite:/ha/rull/lapack-3.1.0/SRC$

It's only a formal try, since i don't really understand this part, but 
the test is passed successfully. Perhaps it's helpful.

With regards
Ingo


<Prev in Thread] Current Thread [Next in Thread>
  • [Lapack] possible bugfix for floating point exception in LAPACK 3.1.0 test, Dr.-Ing. Ingo D. Rullhusen <=


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