Method qrUpdateR

JSAT/src/jsat/linear/GenericMatrix.java:748–766

(int k, int N, Matrix A, double[] vk, double TwoOverBeta, int M)

    }

    private void qrUpdateR(int k, int N, Matrix A, double[] vk, double TwoOverBeta, int M)
    {
        //First run of loop removed, as it will be setting zeros. More accurate to just set them ourselves
        if(k < N)
        {
            qrUpdateRInitalLoop(k, A, vk, TwoOverBeta, M);
        }
        //The rest of the normal look
        for(int j = k+1; j < N; j++)
        {
            double y = 0;//y = vk dot A_j
            for(int i = k; i < A.cols(); i++)
                y += vk[i]*A.get(j, i);
    
            y *= TwoOverBeta;
            for(int i = k; i < M; i++)
                A.increment(j, i, -y*vk[i]);
        }
    }

    private void qrUpdateRInitalLoop(int k, Matrix A, double[] vk, double TwoOverBeta, int M)
    {