commons-dev mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From Greg Sterijevski <gsterijev...@gmail.com>
Subject [Math] TridiagonalTransformer
Date Mon, 22 Aug 2011 21:32:14 GMT
Hello All,

In TriaDiagonalTransformer, I notice the following (commencing at line 104).


                final double[] hK = householderVectors[k - 1];
                final double inv = 1.0 / (secondary[k - 1] * hK[k]);
                cachedQt.setEntry(k, k, 1);
                if (hK[k] != 0.0) {


Shouldn't the line:
     final double inv = 1.0 / (secondary[k - 1] * hK[k]);

be after the test for hK[k] != 0?

-Greg

Mime
  • Unnamed multipart/alternative (inline, None, 0 bytes)
View raw message