public MatrixF getL() {
if (!isDecomposed()) {
decompose();
}
MatrixF l = MatrixF.empty(m, n);
for (int i = 0; i < m; i++) {
for (int j = 0; j < n; j++) {
if (i > j) {
l.setData(i, j, lu[i][j]);
} else if (i == j) {
l.setData(i, j, 1D);
} else {
l.setData(i, j, 0F);
}
}
}
return l;
}