You are viewing a plain text version of this content. The canonical link for it is here.

Posted to dev@commons.apache.org by Greg Sterijevski <gs...@gmail.com> on 2011/08/22 23:32:14 UTC

[Math] TridiagonalTransformer

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