r sets TOL(4) &lt;= 0, then an
          implicitly computed, default tolerance, defined by
          TOLDEF2 = N*EPS, is used instead. The 1-norms of A and E
          are also taken into account.  TOL(4) &lt; 1.

</PRE>
<B>Workspace</B>
<PRE>
  IWORK   INTEGER array, dimension (LIWORK)
          LIWORK &gt;= 1, if MIN(N,P,M) = 0, or B = 0, or C = 0; else
          LIWORK &gt;= MAX(1,N), if DICO = 'C', JOBE = 'I', and
                    JOBD &lt;&gt; 'D';
          LIWORK &gt;= 2*N + M + P + R + 12, other