[CZ]GESVJ - OMPQ eval only when rotating
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.
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.