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

Post here if you want to report a bug to the LAPACK team

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

Postby themos » Mon Jun 16, 2014 7:34 am

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
themos
 
Posts: 6
Joined: Thu Mar 23, 2006 7:54 am

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

Postby Julien Langou » Mon Jun 16, 2014 1:41 pm

Hi Themos, thanks for the bug report, I have just fixed this (svn revision 1483). Cheers, Julien.
Julien Langou
 
Posts: 734
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA


Return to Bug report

Who is online

Users browsing this forum: Google [Bot] and 1 guest