docs: public static void computeRowOfV( final int blockLength ,
                                        final DSubmatrixD1 A ,
                                        final DSubmatrixD1 V ,
                                        int row ,
                                        double gamma )
      {
          // val=(y^T*u)
          double val = BlockHouseHolder_DDRB.innerProdRow(blockLength,A,row,V,row,1);
  
          // take in account the one
          double before = A.get(row,row+1);
          A.set(row,row+1,1);
  
          // v = y - (1/2)gamma*val * u
          VectorOps_DDRB.add_row(blockLength,V,row,1,A,row,-0.5*gamma*val,V,row,row+1,A.col1-A.col0);
  
          A.set(row,row+1,before);
      }