[CZ]GESVJ - OMPQ eval only when rotating

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

[CZ]GESVJ - OMPQ eval only when rotating

Postby lawrence mulholland » Wed Jan 27, 2016 11:35 am

SRC/[cz]gesvj.f

lines 908 and 1129:
OMPQ = AAPQ / ABS(AAPQ)

need to move inside IF( ROTOK ) THEN block
that is, after lines 926 and 1142 respectively.

The variable OMPQ is only used inside these blocks and can, when rotation is unsafe,
overflow when evaluated.
lawrence mulholland
 
Posts: 25
Joined: Mon Jun 11, 2012 6:33 am
Location: NAG Ltd, Oxford, UK

Re: [CZ]GESVJ - OMPQ eval only when rotating

Postby admin » Thu Feb 11, 2016 12:24 am

Dear Lawrence,
Thank you for reporting the issue
This issue has been listed as bug 155 and corrected in revision 1688 in our svn repository.
Julie
admin
Site Admin
 
Posts: 608
Joined: Wed Dec 08, 2004 7:07 pm


Return to Bug report

Who is online

Users browsing this forum: No registered users and 2 guests