Page 1 of 1

[CZ]GESVJ - OMPQ eval only when rotating

PostPosted: Wed Jan 27, 2016 11:35 am
by lawrence mulholland
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.

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

PostPosted: Thu Feb 11, 2016 12:24 am
by admin
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