[BACK]Return to zgelsd.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/zgelsd.f between versions 1.4 and 1.5

version 1.4, 2010/08/06 15:32:38 version 1.5, 2010/08/07 13:18:08
Line 114 Line 114
 *  *
 *  RWORK   (workspace) DOUBLE PRECISION array, dimension (MAX(1,LRWORK))  *  RWORK   (workspace) DOUBLE PRECISION array, dimension (MAX(1,LRWORK))
 *          LRWORK >=  *          LRWORK >=
 *              10*N + 2*N*SMLSIZ + 8*N*NLVL + 3*SMLSIZ*NRHS +  *             10*N + 2*N*SMLSIZ + 8*N*NLVL + 3*SMLSIZ*NRHS +
 *             (SMLSIZ+1)**2  *             MAX( (SMLSIZ+1)**2, N*(1+NRHS) + 2*NRHS )
 *          if M is greater than or equal to N or  *          if M is greater than or equal to N or
 *             10*M + 2*M*SMLSIZ + 8*M*NLVL + 3*SMLSIZ*NRHS +  *             10*M + 2*M*SMLSIZ + 8*M*NLVL + 3*SMLSIZ*NRHS +
 *             (SMLSIZ+1)**2  *             MAX( (SMLSIZ+1)**2, N*(1+NRHS) + 2*NRHS )
 *          if M is less than N, the code will execute correctly.  *          if M is less than N, the code will execute correctly.
 *          SMLSIZ is returned by ILAENV and is equal to the maximum  *          SMLSIZ is returned by ILAENV and is equal to the maximum
 *          size of the subproblems at the bottom of the computation  *          size of the subproblems at the bottom of the computation
Line 229 Line 229
 *              Path 1 - overdetermined or exactly determined.  *              Path 1 - overdetermined or exactly determined.
 *  *
                LRWORK = 10*N + 2*N*SMLSIZ + 8*N*NLVL + 3*SMLSIZ*NRHS +                 LRWORK = 10*N + 2*N*SMLSIZ + 8*N*NLVL + 3*SMLSIZ*NRHS +
      $                  ( SMLSIZ + 1 )**2       $                  MAX( (SMLSIZ+1)**2, N*(1+NRHS) + 2*NRHS )
                MAXWRK = MAX( MAXWRK, 2*N + ( MM + N )*ILAENV( 1,                 MAXWRK = MAX( MAXWRK, 2*N + ( MM + N )*ILAENV( 1,
      $                       'ZGEBRD', ' ', MM, N, -1, -1 ) )       $                       'ZGEBRD', ' ', MM, N, -1, -1 ) )
                MAXWRK = MAX( MAXWRK, 2*N + NRHS*ILAENV( 1, 'ZUNMBR',                 MAXWRK = MAX( MAXWRK, 2*N + NRHS*ILAENV( 1, 'ZUNMBR',
Line 241 Line 241
             END IF              END IF
             IF( N.GT.M ) THEN              IF( N.GT.M ) THEN
                LRWORK = 10*M + 2*M*SMLSIZ + 8*M*NLVL + 3*SMLSIZ*NRHS +                 LRWORK = 10*M + 2*M*SMLSIZ + 8*M*NLVL + 3*SMLSIZ*NRHS +
      $                  ( SMLSIZ + 1 )**2       $                  MAX( (SMLSIZ+1)**2, N*(1+NRHS) + 2*NRHS )
                IF( N.GE.MNTHR ) THEN                 IF( N.GE.MNTHR ) THEN
 *  *
 *                 Path 2a - underdetermined, with many more columns  *                 Path 2a - underdetermined, with many more columns
Removed from v.1.4  
changed lines
  Added in v.1.5

CVSweb interface <joel.bertrand@systella.fr>