docs: public static void rank1Update(double gamma,
                                     DMatrixRMaj A ,
                                     DMatrixRMaj u , DMatrixRMaj w ,
                                     DMatrixRMaj B )
      {
          int n = u.getNumElements();
  
          int matrixIndex = 0;
          for( int i = 0; i < n; i++ ) {
              double elementU = u.data[i];
  
              for( int j = 0; j < n; j++ , matrixIndex++) {
                  B.data[matrixIndex] = A.data[matrixIndex] + gamma*elementU*w.data[j];
              }
          }
      }