From f5cdededb32ba1e5cbff702d1b028d76605d9e51 Mon Sep 17 00:00:00 2001 From: james Date: Mon, 23 Apr 2012 23:16:24 +0000 Subject: modified exceptional shift for complex versions - solves problem of nan's appearing on some platforms, most notably IBM/XLF, causing an infinite loop in xLARTG (which also should be fixed by enforcing a max iteration count when computing SCALE) --- SRC/zhgeqz.f | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) (limited to 'SRC/zhgeqz.f') diff --git a/SRC/zhgeqz.f b/SRC/zhgeqz.f index 2fb01caf..3f28131b 100644 --- a/SRC/zhgeqz.f +++ b/SRC/zhgeqz.f @@ -744,7 +744,8 @@ * * Exceptional shift. Chosen for no particularly good reason. * - ESHIFT = ESHIFT + H(ILAST,ILAST-1)/T(ILAST-1,ILAST-1) + ESHIFT = ESHIFT + (ASCALE*H(ILAST,ILAST-1))/ + $ (BSCALE*T(ILAST-1,ILAST-1)) SHIFT = ESHIFT END IF * -- cgit v1.2.3