Method run

JSAT/src/jsat/linear/DenseMatrix.java:914–952

()

        }
        
        public void run()
        {
            //Computing Q
            {
                //We are computing Q' in what we are treating as the column major order, which represents Q in row major order, which is what we want!
                for(int j = 0+threadID; j < Q.cols(); j+=LogicalCores)
                {
                    double[] Q_j = Q.matrix[j];
                    double y = 0;//y = vk dot A_j
                    for (int i = k; i < Q.cols(); i++)
                        y += vk[i] * Q_j[i];

                    y *= TwoOverBeta;
                    for (int i = k; i < Q.rows(); i++)
                    {
                        Q_j[i] -= y*vk[i];
                    }
                }
            }
            
            //First run of loop removed, as it will be setting zeros. More accurate to just set them ourselves
            if(k < N && threadID == 0)
            {
                qrUpdateRFirstIteration(A, k, vk, TwoOverBeta, M);
            }
            //The rest of the normal look
            for(int j = k+1+threadID; j < N; j+=LogicalCores)
            {
                double[] A_j = A.matrix[j];
                double y = 0;//y = vk dot A_j
                for(int i = k; i < A.cols(); i++)
                    y += vk[i]*A_j[i];
        
                y *= TwoOverBeta;
                for(int i = k; i < M; i++)
                    A_j[i] -= y*vk[i];
            }
            latch.countDown();
        }
        
    }