Page 1 of 1

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

PostPosted: Mon Jun 16, 2014 7:34 am
by 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

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

PostPosted: Mon Jun 16, 2014 1:41 pm
by Julien Langou
Hi Themos, thanks for the bug report, I have just fixed this (svn revision 1483). Cheers, Julien.