work1[pj] = sPar * diag[pj]; determineLMDirection(qy, work1, work2, internalData, solvedCols, work3, lmDir);
work1[pj] = sPar * diag[pj]; determineLMDirection(qy, work1, work2, internalData, solvedCols, work3, lmDir);
work1[pj] = sPar * diag[pj]; determineLMDirection(qy, work1, work2, internalData, solvedCols, work3, lmDir);