code: public static void computeY_t_V( final int blockLength , final DSubmatrixD1 Y ,
                                       final int col , final double []temp )
      {
          final int widthB = Y.col1-Y.col0;
  
          for( int j = 0; j < col; j++ ) {
              temp[j] = innerProdCol(blockLength,Y,col,widthB,j,widthB);
          }
      }
docs: Computes Y<sup>T</sup>v<sup>(j)</sup>.  Where Y are the columns before 'col' and v is the column
  at 'col'.  The zeros and ones are taken in account.  The solution is a vector with 'col' elements.
  
  width of Y must be along the block of original matrix A
  
  @param temp Temporary storage of least length 'col'