(int k, Matrix A, double[] vk, double TwoOverBeta, int M)
| 766 | } |
| 767 | |
| 768 | private void qrUpdateRInitalLoop(int k, Matrix A, double[] vk, double TwoOverBeta, int M) |
| 769 | { |
| 770 | double y = 0;//y = vk dot A_j |
| 771 | for(int i = k; i < A.cols(); i++) |
| 772 | y += vk[i]*A.get(k, i); |
| 773 | |
| 774 | y *= TwoOverBeta; |
| 775 | A.increment(k, k, -y*vk[k]); |
| 776 | |
| 777 | for(int i = k+1; i < M; i++) |
| 778 | A.set(k, i, 0.0); |
| 779 | } |
| 780 | |
| 781 | private void qrUpdateQ(Matrix Q, int k, double[] vk, double TwoOverBeta) |
| 782 | { |