?LAHQR: docs for WI

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

?LAHQR: docs for WI

Postby vchalupecky » Mon Jun 20, 2016 9:29 pm

The comment for the WI parameter says that "...if H(i:i+1,i:i+1) is a 2-by-2 diagonal block, WI(i) = sqrt(H(i+1,i)*H(i,i+1)) ...".
However, the product of the off-diagonal elements in the 2x2 block is negative, so it should say "...WI(i) = SQRT(ABS(H(i+1,i)*H(i,i+1)))..." instead.
Posts: 6
Joined: Fri Jun 10, 2016 2:54 am

Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest