double[] factorWork;
lapack.dgeqp3(m, n, Afact.getData(), Matrices.ld(m), jpvt, tau,
factorWorkOptimalSize, -1, info);
factorWork = new double[(int) factorWorkOptimalSize[0]];
lapack.dgeqp3(m, n, Afact.getData(), Matrices.ld(m), jpvt, tau,
factorWork, factorWork.length, info);