* the algorithm fails.
*/
private IDoubleArray computeDecomposition(IDoubleArray input) {
int n = input.rows();
IDoubleArray result = Doubles.create.array(n, n);
result.zero();
for (int i = 0; i < n; i++) {
for (int j = 0; j <= i; j++) {
double sum = 0.0;
for (int k = 0; k < j; k++) {