private void copyZeros(FSubmatrixD1 subU)
{
int N = Math.Min(A.blockLength, subU.col1 - subU.col0);
for (int i = 0; i < N; i++)
{
// save the zeros
for (int j = 0; j <= i; j++)
{
zerosM.unsafe_set(i, j, subU.get(i, j));
subU.set(i, j, 0);
}
// save the one
if (subU.col0 + i + 1 < subU.original.numCols)
{
zerosM.unsafe_set(i, i + 1, subU.get(i, i + 1));
subU.set(i, i + 1, 1);
}
}
}