[Math] TridiagonalTransformer

[Greg Sterijevski]

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