zlanhf.f N=1 NORM="M" weakness

[themos]

Hello

In zlanhf.f, first executable statements

Code: Select all

      IF( N.EQ.0 ) THEN
         ZLANHF = ZERO
         RETURN
      ELSE IF( N.EQ.1 ) THEN
         ZLANHF = ABS(A(0))
         RETURN
      END IF


When N=1, we should be ignoring any imaginary part (as we do for general diagonal elements later on) and so it should read

Code: Select all

      IF( N.EQ.0 ) THEN
         ZLANHF = ZERO
         RETURN
      ELSE IF( N.EQ.1 ) THEN
         ZLANHF = ABS(DBLE(A(0)))
         RETURN
      END IF


Themos Tsikas
NAG Ltd

[Julien Langou]

Hi Themos, thanks for the bug report, I have just fixed this (svn revision 1483). Cheers, Julien.