f this array contain
          the cosines and sines of the symplectic Givens rotations
          applied from the right-hand side used to compute the
          reduction.

  TAUL    (output) DOUBLE PRECISION array, dimension (NB)
          On exit, the first NB elements of this array contain the
          scalar factors of some of the elementary reflectors
          applied form the left-hand side.

  TAUR    (output) DOUBLE PRECISION array, dimension (NB)
          On exit, the first NB-1 elements of this