diff options
| author | james <james@8a072113-8704-0410-8d35-dd094bca7971> | 2012-04-11 01:50:07 +0000 |
|---|---|---|
| committer | james <james@8a072113-8704-0410-8d35-dd094bca7971> | 2012-04-11 01:50:07 +0000 |
| commit | c254495b4de7667073f941bdfe89f8325c66687b (patch) | |
| tree | 3f7ceda8c608e0a9f2cf235a8b2059719b3ca40b /SRC/clanhf.f | |
| parent | 1b5afc4229b27adcfb5a1f761e308710ff4c7d56 (diff) | |
added code to handle N=1 case to fix divide by zero problem that occurs
when the single (diagonal) element is negative
see: http://icl.cs.utk.edu/lapack-forum/archives/lapack/msg01195.html
Diffstat (limited to 'SRC/clanhf.f')
| -rw-r--r-- | SRC/clanhf.f | 3 |
1 files changed, 3 insertions, 0 deletions
diff --git a/SRC/clanhf.f b/SRC/clanhf.f index af5c16f8..988cd708 100644 --- a/SRC/clanhf.f +++ b/SRC/clanhf.f @@ -286,6 +286,9 @@ IF( N.EQ.0 ) THEN CLANHF = ZERO RETURN + ELSE IF( N.EQ.1 ) THEN + CLANHF = ABS( A(0) ) + RETURN END IF * * set noe = 1 if n is odd. if n is even set noe=0 |