--- rpl/lapack/lapack/dorcsd.f	2010/12/21 13:53:34	1.2
+++ rpl/lapack/lapack/dorcsd.f	2011/07/22 07:38:08	1.3
@@ -5,7 +5,7 @@
      $                             LDV2T, WORK, LWORK, IWORK, INFO )
       IMPLICIT NONE
 *
-*  -- LAPACK routine (version 3.3.0) --
+*  -- LAPACK routine (version 3.3.1) --
 *
 *  -- Contributed by Brian Sutton of the Randolph-Macon College --
 *  -- November 2010
@@ -13,6 +13,8 @@
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--     
 *
+* @precisions normal d -> s
+*
 *     .. Scalar Arguments ..
       CHARACTER          JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS
       INTEGER            INFO, LDU1, LDU2, LDV1T, LDV2T, LDX11, LDX12,
@@ -144,7 +146,7 @@
 *          this value as the first entry of the work array, and no error
 *          message related to LWORK is issued by XERBLA.
 *
-*  IWORK   (workspace) INTEGER array, dimension (M-Q)
+*  IWORK   (workspace) INTEGER array, dimension (M-MIN(P, M-P, Q, M-Q))
 *
 *  INFO    (output) INTEGER
 *          = 0:  successful exit.
@@ -302,7 +304,7 @@
      $              IORBDB + LORBDBWORKOPT, IBBCSD + LBBCSDWORKOPT ) - 1
          LWORKMIN = MAX( IORGQR + LORGQRWORKMIN, IORGLQ + LORGLQWORKMIN,
      $              IORBDB + LORBDBWORKOPT, IBBCSD + LBBCSDWORKMIN ) - 1
-         WORK(1) = LWORKOPT
+         WORK(1) = MAX(LWORKOPT,LWORKMIN)
 *
          IF( LWORK .LT. LWORKMIN .AND. .NOT. LQUERY ) THEN
             INFO = -22