Annotation of rpl/lapack/lapack/ilaenv.f, revision 1.9

1.9     ! bertrand    1: *> \brief \b ILAENV
        !             2: *
        !             3: *  =========== DOCUMENTATION ===========
1.1       bertrand    4: *
1.9     ! bertrand    5: * Online html documentation available at 
        !             6: *            http://www.netlib.org/lapack/explore-html/ 
1.1       bertrand    7: *
1.9     ! bertrand    8: *> \htmlonly
        !             9: *> Download ILAENV + dependencies 
        !            10: *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/ilaenv.f"> 
        !            11: *> [TGZ]</a> 
        !            12: *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/ilaenv.f"> 
        !            13: *> [ZIP]</a> 
        !            14: *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/ilaenv.f"> 
        !            15: *> [TXT]</a>
        !            16: *> \endhtmlonly 
        !            17: *
        !            18: *  Definition:
        !            19: *  ===========
        !            20: *
        !            21: *       INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
        !            22: * 
        !            23: *       .. Scalar Arguments ..
        !            24: *       CHARACTER*( * )    NAME, OPTS
        !            25: *       INTEGER            ISPEC, N1, N2, N3, N4
        !            26: *       ..
        !            27: *  
        !            28: *
        !            29: *> \par Purpose:
        !            30: *  =============
        !            31: *>
        !            32: *> \verbatim
        !            33: *>
        !            34: *> ILAENV is called from the LAPACK routines to choose problem-dependent
        !            35: *> parameters for the local environment.  See ISPEC for a description of
        !            36: *> the parameters.
        !            37: *>
        !            38: *> ILAENV returns an INTEGER
        !            39: *> if ILAENV >= 0: ILAENV returns the value of the parameter specified by ISPEC
        !            40: *> if ILAENV < 0:  if ILAENV = -k, the k-th argument had an illegal value.
        !            41: *>
        !            42: *> This version provides a set of parameters which should give good,
        !            43: *> but not optimal, performance on many of the currently available
        !            44: *> computers.  Users are encouraged to modify this subroutine to set
        !            45: *> the tuning parameters for their particular machine using the option
        !            46: *> and problem size information in the arguments.
        !            47: *>
        !            48: *> This routine will not function correctly if it is converted to all
        !            49: *> lower case.  Converting it to all upper case is allowed.
        !            50: *> \endverbatim
        !            51: *
        !            52: *  Arguments:
        !            53: *  ==========
        !            54: *
        !            55: *> \param[in] ISPEC
        !            56: *> \verbatim
        !            57: *>          ISPEC is INTEGER
        !            58: *>          Specifies the parameter to be returned as the value of
        !            59: *>          ILAENV.
        !            60: *>          = 1: the optimal blocksize; if this value is 1, an unblocked
        !            61: *>               algorithm will give the best performance.
        !            62: *>          = 2: the minimum block size for which the block routine
        !            63: *>               should be used; if the usable block size is less than
        !            64: *>               this value, an unblocked routine should be used.
        !            65: *>          = 3: the crossover point (in a block routine, for N less
        !            66: *>               than this value, an unblocked routine should be used)
        !            67: *>          = 4: the number of shifts, used in the nonsymmetric
        !            68: *>               eigenvalue routines (DEPRECATED)
        !            69: *>          = 5: the minimum column dimension for blocking to be used;
        !            70: *>               rectangular blocks must have dimension at least k by m,
        !            71: *>               where k is given by ILAENV(2,...) and m by ILAENV(5,...)
        !            72: *>          = 6: the crossover point for the SVD (when reducing an m by n
        !            73: *>               matrix to bidiagonal form, if max(m,n)/min(m,n) exceeds
        !            74: *>               this value, a QR factorization is used first to reduce
        !            75: *>               the matrix to a triangular form.)
        !            76: *>          = 7: the number of processors
        !            77: *>          = 8: the crossover point for the multishift QR method
        !            78: *>               for nonsymmetric eigenvalue problems (DEPRECATED)
        !            79: *>          = 9: maximum size of the subproblems at the bottom of the
        !            80: *>               computation tree in the divide-and-conquer algorithm
        !            81: *>               (used by xGELSD and xGESDD)
        !            82: *>          =10: ieee NaN arithmetic can be trusted not to trap
        !            83: *>          =11: infinity arithmetic can be trusted not to trap
        !            84: *>          12 <= ISPEC <= 16:
        !            85: *>               xHSEQR or one of its subroutines,
        !            86: *>               see IPARMQ for detailed explanation
        !            87: *> \endverbatim
        !            88: *>
        !            89: *> \param[in] NAME
        !            90: *> \verbatim
        !            91: *>          NAME is CHARACTER*(*)
        !            92: *>          The name of the calling subroutine, in either upper case or
        !            93: *>          lower case.
        !            94: *> \endverbatim
        !            95: *>
        !            96: *> \param[in] OPTS
        !            97: *> \verbatim
        !            98: *>          OPTS is CHARACTER*(*)
        !            99: *>          The character options to the subroutine NAME, concatenated
        !           100: *>          into a single character string.  For example, UPLO = 'U',
        !           101: *>          TRANS = 'T', and DIAG = 'N' for a triangular routine would
        !           102: *>          be specified as OPTS = 'UTN'.
        !           103: *> \endverbatim
        !           104: *>
        !           105: *> \param[in] N1
        !           106: *> \verbatim
        !           107: *>          N1 is INTEGER
        !           108: *> \endverbatim
        !           109: *>
        !           110: *> \param[in] N2
        !           111: *> \verbatim
        !           112: *>          N2 is INTEGER
        !           113: *> \endverbatim
        !           114: *>
        !           115: *> \param[in] N3
        !           116: *> \verbatim
        !           117: *>          N3 is INTEGER
        !           118: *> \endverbatim
        !           119: *>
        !           120: *> \param[in] N4
        !           121: *> \verbatim
        !           122: *>          N4 is INTEGER
        !           123: *>          Problem dimensions for the subroutine NAME; these may not all
        !           124: *>          be required.
        !           125: *> \endverbatim
        !           126: *
        !           127: *  Authors:
        !           128: *  ========
        !           129: *
        !           130: *> \author Univ. of Tennessee 
        !           131: *> \author Univ. of California Berkeley 
        !           132: *> \author Univ. of Colorado Denver 
        !           133: *> \author NAG Ltd. 
        !           134: *
        !           135: *> \date November 2011
        !           136: *
        !           137: *> \ingroup auxOTHERauxiliary
        !           138: *
        !           139: *> \par Further Details:
        !           140: *  =====================
        !           141: *>
        !           142: *> \verbatim
        !           143: *>
        !           144: *>  The following conventions have been used when calling ILAENV from the
        !           145: *>  LAPACK routines:
        !           146: *>  1)  OPTS is a concatenation of all of the character options to
        !           147: *>      subroutine NAME, in the same order that they appear in the
        !           148: *>      argument list for NAME, even if they are not used in determining
        !           149: *>      the value of the parameter specified by ISPEC.
        !           150: *>  2)  The problem dimensions N1, N2, N3, N4 are specified in the order
        !           151: *>      that they appear in the argument list for NAME.  N1 is used
        !           152: *>      first, N2 second, and so on, and unused problem dimensions are
        !           153: *>      passed a value of -1.
        !           154: *>  3)  The parameter value returned by ILAENV is checked for validity in
        !           155: *>      the calling subroutine.  For example, ILAENV is used to retrieve
        !           156: *>      the optimal blocksize for STRTRI as follows:
        !           157: *>
        !           158: *>      NB = ILAENV( 1, 'STRTRI', UPLO // DIAG, N, -1, -1, -1 )
        !           159: *>      IF( NB.LE.1 ) NB = MAX( 1, N )
        !           160: *> \endverbatim
        !           161: *>
        !           162: *  =====================================================================
        !           163:       INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
1.1       bertrand  164: *
1.9     ! bertrand  165: *  -- LAPACK auxiliary routine (version 3.4.0) --
1.1       bertrand  166: *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
                    167: *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
1.9     ! bertrand  168: *     November 2011
1.1       bertrand  169: *
                    170: *     .. Scalar Arguments ..
                    171:       CHARACTER*( * )    NAME, OPTS
                    172:       INTEGER            ISPEC, N1, N2, N3, N4
                    173: *     ..
                    174: *
                    175: *  =====================================================================
                    176: *
                    177: *     .. Local Scalars ..
                    178:       INTEGER            I, IC, IZ, NB, NBMIN, NX
                    179:       LOGICAL            CNAME, SNAME
                    180:       CHARACTER          C1*1, C2*2, C4*2, C3*3, SUBNAM*6
                    181: *     ..
                    182: *     .. Intrinsic Functions ..
                    183:       INTRINSIC          CHAR, ICHAR, INT, MIN, REAL
                    184: *     ..
                    185: *     .. External Functions ..
                    186:       INTEGER            IEEECK, IPARMQ
                    187:       EXTERNAL           IEEECK, IPARMQ
                    188: *     ..
                    189: *     .. Executable Statements ..
                    190: *
                    191:       GO TO ( 10, 10, 10, 80, 90, 100, 110, 120,
                    192:      $        130, 140, 150, 160, 160, 160, 160, 160 )ISPEC
                    193: *
                    194: *     Invalid value for ISPEC
                    195: *
                    196:       ILAENV = -1
                    197:       RETURN
                    198: *
                    199:    10 CONTINUE
                    200: *
                    201: *     Convert NAME to upper case if the first character is lower case.
                    202: *
                    203:       ILAENV = 1
                    204:       SUBNAM = NAME
                    205:       IC = ICHAR( SUBNAM( 1: 1 ) )
                    206:       IZ = ICHAR( 'Z' )
                    207:       IF( IZ.EQ.90 .OR. IZ.EQ.122 ) THEN
                    208: *
                    209: *        ASCII character set
                    210: *
                    211:          IF( IC.GE.97 .AND. IC.LE.122 ) THEN
                    212:             SUBNAM( 1: 1 ) = CHAR( IC-32 )
                    213:             DO 20 I = 2, 6
                    214:                IC = ICHAR( SUBNAM( I: I ) )
                    215:                IF( IC.GE.97 .AND. IC.LE.122 )
                    216:      $            SUBNAM( I: I ) = CHAR( IC-32 )
                    217:    20       CONTINUE
                    218:          END IF
                    219: *
                    220:       ELSE IF( IZ.EQ.233 .OR. IZ.EQ.169 ) THEN
                    221: *
                    222: *        EBCDIC character set
                    223: *
                    224:          IF( ( IC.GE.129 .AND. IC.LE.137 ) .OR.
                    225:      $       ( IC.GE.145 .AND. IC.LE.153 ) .OR.
                    226:      $       ( IC.GE.162 .AND. IC.LE.169 ) ) THEN
                    227:             SUBNAM( 1: 1 ) = CHAR( IC+64 )
                    228:             DO 30 I = 2, 6
                    229:                IC = ICHAR( SUBNAM( I: I ) )
                    230:                IF( ( IC.GE.129 .AND. IC.LE.137 ) .OR.
                    231:      $             ( IC.GE.145 .AND. IC.LE.153 ) .OR.
                    232:      $             ( IC.GE.162 .AND. IC.LE.169 ) )SUBNAM( I:
                    233:      $             I ) = CHAR( IC+64 )
                    234:    30       CONTINUE
                    235:          END IF
                    236: *
                    237:       ELSE IF( IZ.EQ.218 .OR. IZ.EQ.250 ) THEN
                    238: *
                    239: *        Prime machines:  ASCII+128
                    240: *
                    241:          IF( IC.GE.225 .AND. IC.LE.250 ) THEN
                    242:             SUBNAM( 1: 1 ) = CHAR( IC-32 )
                    243:             DO 40 I = 2, 6
                    244:                IC = ICHAR( SUBNAM( I: I ) )
                    245:                IF( IC.GE.225 .AND. IC.LE.250 )
                    246:      $            SUBNAM( I: I ) = CHAR( IC-32 )
                    247:    40       CONTINUE
                    248:          END IF
                    249:       END IF
                    250: *
                    251:       C1 = SUBNAM( 1: 1 )
                    252:       SNAME = C1.EQ.'S' .OR. C1.EQ.'D'
                    253:       CNAME = C1.EQ.'C' .OR. C1.EQ.'Z'
                    254:       IF( .NOT.( CNAME .OR. SNAME ) )
                    255:      $   RETURN
                    256:       C2 = SUBNAM( 2: 3 )
                    257:       C3 = SUBNAM( 4: 6 )
                    258:       C4 = C3( 2: 3 )
                    259: *
                    260:       GO TO ( 50, 60, 70 )ISPEC
                    261: *
                    262:    50 CONTINUE
                    263: *
                    264: *     ISPEC = 1:  block size
                    265: *
                    266: *     In these examples, separate code is provided for setting NB for
                    267: *     real and complex.  We assume that NB will take the same value in
                    268: *     single or double precision.
                    269: *
                    270:       NB = 1
                    271: *
                    272:       IF( C2.EQ.'GE' ) THEN
                    273:          IF( C3.EQ.'TRF' ) THEN
                    274:             IF( SNAME ) THEN
                    275:                NB = 64
                    276:             ELSE
                    277:                NB = 64
                    278:             END IF
                    279:          ELSE IF( C3.EQ.'QRF' .OR. C3.EQ.'RQF' .OR. C3.EQ.'LQF' .OR.
                    280:      $            C3.EQ.'QLF' ) THEN
                    281:             IF( SNAME ) THEN
                    282:                NB = 32
                    283:             ELSE
                    284:                NB = 32
                    285:             END IF
                    286:          ELSE IF( C3.EQ.'HRD' ) THEN
                    287:             IF( SNAME ) THEN
                    288:                NB = 32
                    289:             ELSE
                    290:                NB = 32
                    291:             END IF
                    292:          ELSE IF( C3.EQ.'BRD' ) THEN
                    293:             IF( SNAME ) THEN
                    294:                NB = 32
                    295:             ELSE
                    296:                NB = 32
                    297:             END IF
                    298:          ELSE IF( C3.EQ.'TRI' ) THEN
                    299:             IF( SNAME ) THEN
                    300:                NB = 64
                    301:             ELSE
                    302:                NB = 64
                    303:             END IF
                    304:          END IF
                    305:       ELSE IF( C2.EQ.'PO' ) THEN
                    306:          IF( C3.EQ.'TRF' ) THEN
                    307:             IF( SNAME ) THEN
                    308:                NB = 64
                    309:             ELSE
                    310:                NB = 64
                    311:             END IF
                    312:          END IF
                    313:       ELSE IF( C2.EQ.'SY' ) THEN
                    314:          IF( C3.EQ.'TRF' ) THEN
                    315:             IF( SNAME ) THEN
                    316:                NB = 64
                    317:             ELSE
                    318:                NB = 64
                    319:             END IF
                    320:          ELSE IF( SNAME .AND. C3.EQ.'TRD' ) THEN
                    321:             NB = 32
                    322:          ELSE IF( SNAME .AND. C3.EQ.'GST' ) THEN
                    323:             NB = 64
                    324:          END IF
                    325:       ELSE IF( CNAME .AND. C2.EQ.'HE' ) THEN
                    326:          IF( C3.EQ.'TRF' ) THEN
                    327:             NB = 64
                    328:          ELSE IF( C3.EQ.'TRD' ) THEN
                    329:             NB = 32
                    330:          ELSE IF( C3.EQ.'GST' ) THEN
                    331:             NB = 64
                    332:          END IF
                    333:       ELSE IF( SNAME .AND. C2.EQ.'OR' ) THEN
                    334:          IF( C3( 1: 1 ).EQ.'G' ) THEN
                    335:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    336:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    337:      $           THEN
                    338:                NB = 32
                    339:             END IF
                    340:          ELSE IF( C3( 1: 1 ).EQ.'M' ) THEN
                    341:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    342:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    343:      $           THEN
                    344:                NB = 32
                    345:             END IF
                    346:          END IF
                    347:       ELSE IF( CNAME .AND. C2.EQ.'UN' ) THEN
                    348:          IF( C3( 1: 1 ).EQ.'G' ) THEN
                    349:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    350:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    351:      $           THEN
                    352:                NB = 32
                    353:             END IF
                    354:          ELSE IF( C3( 1: 1 ).EQ.'M' ) THEN
                    355:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    356:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    357:      $           THEN
                    358:                NB = 32
                    359:             END IF
                    360:          END IF
                    361:       ELSE IF( C2.EQ.'GB' ) THEN
                    362:          IF( C3.EQ.'TRF' ) THEN
                    363:             IF( SNAME ) THEN
                    364:                IF( N4.LE.64 ) THEN
                    365:                   NB = 1
                    366:                ELSE
                    367:                   NB = 32
                    368:                END IF
                    369:             ELSE
                    370:                IF( N4.LE.64 ) THEN
                    371:                   NB = 1
                    372:                ELSE
                    373:                   NB = 32
                    374:                END IF
                    375:             END IF
                    376:          END IF
                    377:       ELSE IF( C2.EQ.'PB' ) THEN
                    378:          IF( C3.EQ.'TRF' ) THEN
                    379:             IF( SNAME ) THEN
                    380:                IF( N2.LE.64 ) THEN
                    381:                   NB = 1
                    382:                ELSE
                    383:                   NB = 32
                    384:                END IF
                    385:             ELSE
                    386:                IF( N2.LE.64 ) THEN
                    387:                   NB = 1
                    388:                ELSE
                    389:                   NB = 32
                    390:                END IF
                    391:             END IF
                    392:          END IF
                    393:       ELSE IF( C2.EQ.'TR' ) THEN
                    394:          IF( C3.EQ.'TRI' ) THEN
                    395:             IF( SNAME ) THEN
                    396:                NB = 64
                    397:             ELSE
                    398:                NB = 64
                    399:             END IF
                    400:          END IF
                    401:       ELSE IF( C2.EQ.'LA' ) THEN
                    402:          IF( C3.EQ.'UUM' ) THEN
                    403:             IF( SNAME ) THEN
                    404:                NB = 64
                    405:             ELSE
                    406:                NB = 64
                    407:             END IF
                    408:          END IF
                    409:       ELSE IF( SNAME .AND. C2.EQ.'ST' ) THEN
                    410:          IF( C3.EQ.'EBZ' ) THEN
                    411:             NB = 1
                    412:          END IF
                    413:       END IF
                    414:       ILAENV = NB
                    415:       RETURN
                    416: *
                    417:    60 CONTINUE
                    418: *
                    419: *     ISPEC = 2:  minimum block size
                    420: *
                    421:       NBMIN = 2
                    422:       IF( C2.EQ.'GE' ) THEN
                    423:          IF( C3.EQ.'QRF' .OR. C3.EQ.'RQF' .OR. C3.EQ.'LQF' .OR. C3.EQ.
                    424:      $       'QLF' ) THEN
                    425:             IF( SNAME ) THEN
                    426:                NBMIN = 2
                    427:             ELSE
                    428:                NBMIN = 2
                    429:             END IF
                    430:          ELSE IF( C3.EQ.'HRD' ) THEN
                    431:             IF( SNAME ) THEN
                    432:                NBMIN = 2
                    433:             ELSE
                    434:                NBMIN = 2
                    435:             END IF
                    436:          ELSE IF( C3.EQ.'BRD' ) THEN
                    437:             IF( SNAME ) THEN
                    438:                NBMIN = 2
                    439:             ELSE
                    440:                NBMIN = 2
                    441:             END IF
                    442:          ELSE IF( C3.EQ.'TRI' ) THEN
                    443:             IF( SNAME ) THEN
                    444:                NBMIN = 2
                    445:             ELSE
                    446:                NBMIN = 2
                    447:             END IF
                    448:          END IF
                    449:       ELSE IF( C2.EQ.'SY' ) THEN
                    450:          IF( C3.EQ.'TRF' ) THEN
                    451:             IF( SNAME ) THEN
                    452:                NBMIN = 8
                    453:             ELSE
                    454:                NBMIN = 8
                    455:             END IF
                    456:          ELSE IF( SNAME .AND. C3.EQ.'TRD' ) THEN
                    457:             NBMIN = 2
                    458:          END IF
                    459:       ELSE IF( CNAME .AND. C2.EQ.'HE' ) THEN
                    460:          IF( C3.EQ.'TRD' ) THEN
                    461:             NBMIN = 2
                    462:          END IF
                    463:       ELSE IF( SNAME .AND. C2.EQ.'OR' ) THEN
                    464:          IF( C3( 1: 1 ).EQ.'G' ) THEN
                    465:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    466:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    467:      $           THEN
                    468:                NBMIN = 2
                    469:             END IF
                    470:          ELSE IF( C3( 1: 1 ).EQ.'M' ) THEN
                    471:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    472:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    473:      $           THEN
                    474:                NBMIN = 2
                    475:             END IF
                    476:          END IF
                    477:       ELSE IF( CNAME .AND. C2.EQ.'UN' ) THEN
                    478:          IF( C3( 1: 1 ).EQ.'G' ) THEN
                    479:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    480:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    481:      $           THEN
                    482:                NBMIN = 2
                    483:             END IF
                    484:          ELSE IF( C3( 1: 1 ).EQ.'M' ) THEN
                    485:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    486:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    487:      $           THEN
                    488:                NBMIN = 2
                    489:             END IF
                    490:          END IF
                    491:       END IF
                    492:       ILAENV = NBMIN
                    493:       RETURN
                    494: *
                    495:    70 CONTINUE
                    496: *
                    497: *     ISPEC = 3:  crossover point
                    498: *
                    499:       NX = 0
                    500:       IF( C2.EQ.'GE' ) THEN
                    501:          IF( C3.EQ.'QRF' .OR. C3.EQ.'RQF' .OR. C3.EQ.'LQF' .OR. C3.EQ.
                    502:      $       'QLF' ) THEN
                    503:             IF( SNAME ) THEN
                    504:                NX = 128
                    505:             ELSE
                    506:                NX = 128
                    507:             END IF
                    508:          ELSE IF( C3.EQ.'HRD' ) THEN
                    509:             IF( SNAME ) THEN
                    510:                NX = 128
                    511:             ELSE
                    512:                NX = 128
                    513:             END IF
                    514:          ELSE IF( C3.EQ.'BRD' ) THEN
                    515:             IF( SNAME ) THEN
                    516:                NX = 128
                    517:             ELSE
                    518:                NX = 128
                    519:             END IF
                    520:          END IF
                    521:       ELSE IF( C2.EQ.'SY' ) THEN
                    522:          IF( SNAME .AND. C3.EQ.'TRD' ) THEN
                    523:             NX = 32
                    524:          END IF
                    525:       ELSE IF( CNAME .AND. C2.EQ.'HE' ) THEN
                    526:          IF( C3.EQ.'TRD' ) THEN
                    527:             NX = 32
                    528:          END IF
                    529:       ELSE IF( SNAME .AND. C2.EQ.'OR' ) THEN
                    530:          IF( C3( 1: 1 ).EQ.'G' ) THEN
                    531:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    532:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    533:      $           THEN
                    534:                NX = 128
                    535:             END IF
                    536:          END IF
                    537:       ELSE IF( CNAME .AND. C2.EQ.'UN' ) THEN
                    538:          IF( C3( 1: 1 ).EQ.'G' ) THEN
                    539:             IF( C4.EQ.'QR' .OR. C4.EQ.'RQ' .OR. C4.EQ.'LQ' .OR. C4.EQ.
                    540:      $          'QL' .OR. C4.EQ.'HR' .OR. C4.EQ.'TR' .OR. C4.EQ.'BR' )
                    541:      $           THEN
                    542:                NX = 128
                    543:             END IF
                    544:          END IF
                    545:       END IF
                    546:       ILAENV = NX
                    547:       RETURN
                    548: *
                    549:    80 CONTINUE
                    550: *
                    551: *     ISPEC = 4:  number of shifts (used by xHSEQR)
                    552: *
                    553:       ILAENV = 6
                    554:       RETURN
                    555: *
                    556:    90 CONTINUE
                    557: *
                    558: *     ISPEC = 5:  minimum column dimension (not used)
                    559: *
                    560:       ILAENV = 2
                    561:       RETURN
                    562: *
                    563:   100 CONTINUE
                    564: *
                    565: *     ISPEC = 6:  crossover point for SVD (used by xGELSS and xGESVD)
                    566: *
                    567:       ILAENV = INT( REAL( MIN( N1, N2 ) )*1.6E0 )
                    568:       RETURN
                    569: *
                    570:   110 CONTINUE
                    571: *
                    572: *     ISPEC = 7:  number of processors (not used)
                    573: *
                    574:       ILAENV = 1
                    575:       RETURN
                    576: *
                    577:   120 CONTINUE
                    578: *
                    579: *     ISPEC = 8:  crossover point for multishift (used by xHSEQR)
                    580: *
                    581:       ILAENV = 50
                    582:       RETURN
                    583: *
                    584:   130 CONTINUE
                    585: *
                    586: *     ISPEC = 9:  maximum size of the subproblems at the bottom of the
                    587: *                 computation tree in the divide-and-conquer algorithm
                    588: *                 (used by xGELSD and xGESDD)
                    589: *
                    590:       ILAENV = 25
                    591:       RETURN
                    592: *
                    593:   140 CONTINUE
                    594: *
                    595: *     ISPEC = 10: ieee NaN arithmetic can be trusted not to trap
                    596: *
                    597: *     ILAENV = 0
                    598:       ILAENV = 1
                    599:       IF( ILAENV.EQ.1 ) THEN
                    600:          ILAENV = IEEECK( 1, 0.0, 1.0 )
                    601:       END IF
                    602:       RETURN
                    603: *
                    604:   150 CONTINUE
                    605: *
                    606: *     ISPEC = 11: infinity arithmetic can be trusted not to trap
                    607: *
                    608: *     ILAENV = 0
                    609:       ILAENV = 1
                    610:       IF( ILAENV.EQ.1 ) THEN
                    611:          ILAENV = IEEECK( 0, 0.0, 1.0 )
                    612:       END IF
                    613:       RETURN
                    614: *
                    615:   160 CONTINUE
                    616: *
                    617: *     12 <= ISPEC <= 16: xHSEQR or one of its subroutines. 
                    618: *
                    619:       ILAENV = IPARMQ( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
                    620:       RETURN
                    621: *
                    622: *     End of ILAENV
                    623: *
                    624:       END

CVSweb interface <joel.bertrand@systella.fr>