The R Project SVN R

Rev

Rev 72877 | Rev 81909 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 72877 Rev 73873
Line 2682... Line 2682...
2682
*> \author Univ. of Tennessee
2682
*> \author Univ. of Tennessee
2683
*> \author Univ. of California Berkeley
2683
*> \author Univ. of California Berkeley
2684
*> \author Univ. of Colorado Denver
2684
*> \author Univ. of Colorado Denver
2685
*> \author NAG Ltd.
2685
*> \author NAG Ltd.
2686
*
2686
*
2687
*> \date June 2017
2687
*> \date November 2017
2688
*
2688
*
2689
*> \ingroup complex16GEcomputational
2689
*> \ingroup complex16GEcomputational
2690
*
2690
*
2691
*> \par Further Details:
2691
*> \par Further Details:
2692
*  =====================
2692
*  =====================
Line 2740... Line 2740...
2740
*>
2740
*>
2741
*  =====================================================================
2741
*  =====================================================================
2742
      SUBROUTINE ZGEBRD( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, LWORK,
2742
      SUBROUTINE ZGEBRD( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, LWORK,
2743
     $                   INFO )
2743
     $                   INFO )
2744
*
2744
*
2745
*  -- LAPACK computational routine (version 3.7.1) --
2745
*  -- LAPACK computational routine (version 3.8.0) --
2746
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
2746
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
2747
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
2747
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
2748
*     June 2017
2748
*     November 2017
2749
*
2749
*
2750
*     .. Scalar Arguments ..
2750
*     .. Scalar Arguments ..
2751
      INTEGER            INFO, LDA, LWORK, M, N
2751
      INTEGER            INFO, LDA, LWORK, M, N
2752
*     ..
2752
*     ..
2753
*     .. Array Arguments ..
2753
*     .. Array Arguments ..
Line 2762... Line 2762...
2762
      PARAMETER          ( ONE = ( 1.0D+0, 0.0D+0 ) )
2762
      PARAMETER          ( ONE = ( 1.0D+0, 0.0D+0 ) )
2763
*     ..
2763
*     ..
2764
*     .. Local Scalars ..
2764
*     .. Local Scalars ..
2765
      LOGICAL            LQUERY
2765
      LOGICAL            LQUERY
2766
      INTEGER            I, IINFO, J, LDWRKX, LDWRKY, LWKOPT, MINMN, NB,
2766
      INTEGER            I, IINFO, J, LDWRKX, LDWRKY, LWKOPT, MINMN, NB,
2767
     $                   NBMIN, NX
2767
     $                   NBMIN, NX, WS
2768
      DOUBLE PRECISION   WS
-
 
2769
*     ..
2768
*     ..
2770
*     .. External Subroutines ..
2769
*     .. External Subroutines ..
2771
      EXTERNAL           XERBLA, ZGEBD2, ZGEMM, ZLABRD
2770
      EXTERNAL           XERBLA, ZGEBD2, ZGEMM, ZLABRD
2772
*     ..
2771
*     ..
2773
*     .. Intrinsic Functions ..
2772
*     .. Intrinsic Functions ..
Line 7994... Line 7993...
7994
*> \author Univ. of Tennessee
7993
*> \author Univ. of Tennessee
7995
*> \author Univ. of California Berkeley
7994
*> \author Univ. of California Berkeley
7996
*> \author Univ. of Colorado Denver
7995
*> \author Univ. of Colorado Denver
7997
*> \author NAG Ltd.
7996
*> \author NAG Ltd.
7998
*
7997
*
7999
*> \date December 2016
7998
*> \date November 2017
8000
*
7999
*
8001
*> \ingroup complex16GEauxiliary
8000
*> \ingroup complex16GEauxiliary
8002
*
8001
*
8003
*> \par Contributors:
8002
*> \par Contributors:
8004
*  ==================
8003
*  ==================
Line 8007... Line 8006...
8007
*>     Umea University, S-901 87 Umea, Sweden.
8006
*>     Umea University, S-901 87 Umea, Sweden.
8008
*
8007
*
8009
*  =====================================================================
8008
*  =====================================================================
8010
      SUBROUTINE ZGESC2( N, A, LDA, RHS, IPIV, JPIV, SCALE )
8009
      SUBROUTINE ZGESC2( N, A, LDA, RHS, IPIV, JPIV, SCALE )
8011
*
8010
*
8012
*  -- LAPACK auxiliary routine (version 3.7.0) --
8011
*  -- LAPACK auxiliary routine (version 3.8.0) --
8013
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
8012
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
8014
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
8013
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
8015
*     December 2016
8014
*     November 2017
8016
*
8015
*
8017
*     .. Scalar Arguments ..
8016
*     .. Scalar Arguments ..
8018
      INTEGER            LDA, N
8017
      INTEGER            LDA, N
8019
      DOUBLE PRECISION   SCALE
8018
      DOUBLE PRECISION   SCALE
8020
*     ..
8019
*     ..
Line 8033... Line 8032...
8033
      INTEGER            I, J
8032
      INTEGER            I, J
8034
      DOUBLE PRECISION   BIGNUM, EPS, SMLNUM
8033
      DOUBLE PRECISION   BIGNUM, EPS, SMLNUM
8035
      COMPLEX*16         TEMP
8034
      COMPLEX*16         TEMP
8036
*     ..
8035
*     ..
8037
*     .. External Subroutines ..
8036
*     .. External Subroutines ..
8038
      EXTERNAL           ZLASWP, ZSCAL
8037
      EXTERNAL           ZLASWP, ZSCAL, DLABAD
8039
*     ..
8038
*     ..
8040
*     .. External Functions ..
8039
*     .. External Functions ..
8041
      INTEGER            IZAMAX
8040
      INTEGER            IZAMAX
8042
      DOUBLE PRECISION   DLAMCH
8041
      DOUBLE PRECISION   DLAMCH
8043
      EXTERNAL           IZAMAX, DLAMCH
8042
      EXTERNAL           IZAMAX, DLAMCH
Line 14914... Line 14913...
14914
*>     Umea University, S-901 87 Umea, Sweden.
14913
*>     Umea University, S-901 87 Umea, Sweden.
14915
*
14914
*
14916
*  =====================================================================
14915
*  =====================================================================
14917
      SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
14916
      SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
14918
*
14917
*
14919
*  -- LAPACK auxiliary routine (version 3.7.0) --
14918
*  -- LAPACK auxiliary routine (version 3.8.0) --
14920
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
14919
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
14921
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
14920
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
14922
*     June 2016
14921
*     June 2016
14923
*
14922
*
14924
*     .. Scalar Arguments ..
14923
*     .. Scalar Arguments ..
Line 14938... Line 14937...
14938
*     .. Local Scalars ..
14937
*     .. Local Scalars ..
14939
      INTEGER            I, IP, IPV, J, JP, JPV
14938
      INTEGER            I, IP, IPV, J, JP, JPV
14940
      DOUBLE PRECISION   BIGNUM, EPS, SMIN, SMLNUM, XMAX
14939
      DOUBLE PRECISION   BIGNUM, EPS, SMIN, SMLNUM, XMAX
14941
*     ..
14940
*     ..
14942
*     .. External Subroutines ..
14941
*     .. External Subroutines ..
14943
      EXTERNAL           ZGERU, ZSWAP
14942
      EXTERNAL           ZGERU, ZSWAP, DLABAD
14944
*     ..
14943
*     ..
14945
*     .. External Functions ..
14944
*     .. External Functions ..
14946
      DOUBLE PRECISION   DLAMCH
14945
      DOUBLE PRECISION   DLAMCH
14947
      EXTERNAL           DLAMCH
14946
      EXTERNAL           DLAMCH
14948
*     ..
14947
*     ..
Line 33954... Line 33953...
33954
*> \author Univ. of Tennessee
33953
*> \author Univ. of Tennessee
33955
*> \author Univ. of California Berkeley
33954
*> \author Univ. of California Berkeley
33956
*> \author Univ. of Colorado Denver
33955
*> \author Univ. of Colorado Denver
33957
*> \author NAG Ltd.
33956
*> \author NAG Ltd.
33958
*
33957
*
33959
*> \date December 2016
33958
*> \date November 2017
33960
*
33959
*
33961
*> \ingroup complex16OTHERauxiliary
33960
*> \ingroup complex16OTHERauxiliary
33962
*
33961
*
33963
*  =====================================================================
33962
*  =====================================================================
33964
      SUBROUTINE ZLARFG( N, ALPHA, X, INCX, TAU )
33963
      SUBROUTINE ZLARFG( N, ALPHA, X, INCX, TAU )
33965
*
33964
*
33966
*  -- LAPACK auxiliary routine (version 3.7.0) --
33965
*  -- LAPACK auxiliary routine (version 3.8.0) --
33967
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
33966
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
33968
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
33967
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
33969
*     December 2016
33968
*     November 2017
33970
*
33969
*
33971
*     .. Scalar Arguments ..
33970
*     .. Scalar Arguments ..
33972
      INTEGER            INCX, N
33971
      INTEGER            INCX, N
33973
      COMPLEX*16         ALPHA, TAU
33972
      COMPLEX*16         ALPHA, TAU
33974
*     ..
33973
*     ..
Line 34030... Line 34029...
34030
            KNT = KNT + 1
34029
            KNT = KNT + 1
34031
            CALL ZDSCAL( N-1, RSAFMN, X, INCX )
34030
            CALL ZDSCAL( N-1, RSAFMN, X, INCX )
34032
            BETA = BETA*RSAFMN
34031
            BETA = BETA*RSAFMN
34033
            ALPHI = ALPHI*RSAFMN
34032
            ALPHI = ALPHI*RSAFMN
34034
            ALPHR = ALPHR*RSAFMN
34033
            ALPHR = ALPHR*RSAFMN
34035
            IF( ABS( BETA ).LT.SAFMIN )
34034
            IF( (ABS( BETA ).LT.SAFMIN) .AND. (KNT .LT. 20) )
34036
     $         GO TO 10
34035
     $         GO TO 10
34037
*
34036
*
34038
*           New BETA is at most 1, at least SAFMIN
34037
*           New BETA is at most 1, at least SAFMIN
34039
*
34038
*
34040
            XNORM = DZNRM2( N-1, X, INCX )
34039
            XNORM = DZNRM2( N-1, X, INCX )
Line 38357... Line 38356...
38357
*> \author Univ. of Tennessee
38356
*> \author Univ. of Tennessee
38358
*> \author Univ. of California Berkeley
38357
*> \author Univ. of California Berkeley
38359
*> \author Univ. of Colorado Denver
38358
*> \author Univ. of Colorado Denver
38360
*> \author NAG Ltd.
38359
*> \author NAG Ltd.
38361
*
38360
*
38362
*> \date December 2016
38361
*> \date November 2017
38363
*
38362
*
38364
*> \ingroup complex16OTHERauxiliary
38363
*> \ingroup complex16OTHERauxiliary
38365
*
38364
*
38366
*> \par Further Details:
38365
*> \par Further Details:
38367
*  =====================
38366
*  =====================
Line 38437... Line 38436...
38437
*>
38436
*>
38438
*  =====================================================================
38437
*  =====================================================================
38439
      SUBROUTINE ZLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
38438
      SUBROUTINE ZLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
38440
     $                   CNORM, INFO )
38439
     $                   CNORM, INFO )
38441
*
38440
*
38442
*  -- LAPACK auxiliary routine (version 3.7.0) --
38441
*  -- LAPACK auxiliary routine (version 3.8.0) --
38443
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
38442
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
38444
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
38443
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
38445
*     December 2016
38444
*     November 2017
38446
*
38445
*
38447
*     .. Scalar Arguments ..
38446
*     .. Scalar Arguments ..
38448
      CHARACTER          DIAG, NORMIN, TRANS, UPLO
38447
      CHARACTER          DIAG, NORMIN, TRANS, UPLO
38449
      INTEGER            INFO, LDA, N
38448
      INTEGER            INFO, LDA, N
38450
      DOUBLE PRECISION   SCALE
38449
      DOUBLE PRECISION   SCALE
Line 38475... Line 38474...
38475
      COMPLEX*16         ZDOTC, ZDOTU, ZLADIV
38474
      COMPLEX*16         ZDOTC, ZDOTU, ZLADIV
38476
      EXTERNAL           LSAME, IDAMAX, IZAMAX, DLAMCH, DZASUM, ZDOTC,
38475
      EXTERNAL           LSAME, IDAMAX, IZAMAX, DLAMCH, DZASUM, ZDOTC,
38477
     $                   ZDOTU, ZLADIV
38476
     $                   ZDOTU, ZLADIV
38478
*     ..
38477
*     ..
38479
*     .. External Subroutines ..
38478
*     .. External Subroutines ..
38480
      EXTERNAL           DSCAL, XERBLA, ZAXPY, ZDSCAL, ZTRSV
38479
      EXTERNAL           DSCAL, XERBLA, ZAXPY, ZDSCAL, ZTRSV, DLABAD
38481
*     ..
38480
*     ..
38482
*     .. Intrinsic Functions ..
38481
*     .. Intrinsic Functions ..
38483
      INTRINSIC          ABS, DBLE, DCMPLX, DCONJG, DIMAG, MAX, MIN
38482
      INTRINSIC          ABS, DBLE, DCMPLX, DCONJG, DIMAG, MAX, MIN
38484
*     ..
38483
*     ..
38485
*     .. Statement Functions ..
38484
*     .. Statement Functions ..
Line 47509... Line 47508...
47509
*> \author Univ. of Tennessee
47508
*> \author Univ. of Tennessee
47510
*> \author Univ. of California Berkeley
47509
*> \author Univ. of California Berkeley
47511
*> \author Univ. of Colorado Denver
47510
*> \author Univ. of Colorado Denver
47512
*> \author NAG Ltd.
47511
*> \author NAG Ltd.
47513
*
47512
*
47514
*> \date December 2016
47513
*> \date November 2017
47515
*
47514
*
47516
*> \ingroup complex16OTHERcomputational
47515
*> \ingroup complex16OTHERcomputational
47517
*
47516
*
47518
*> \par Further Details:
47517
*> \par Further Details:
47519
*  =====================
47518
*  =====================
Line 47531... Line 47530...
47531
*>
47530
*>
47532
*  =====================================================================
47531
*  =====================================================================
47533
      SUBROUTINE ZTREVC( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR,
47532
      SUBROUTINE ZTREVC( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR,
47534
     $                   LDVR, MM, M, WORK, RWORK, INFO )
47533
     $                   LDVR, MM, M, WORK, RWORK, INFO )
47535
*
47534
*
47536
*  -- LAPACK computational routine (version 3.7.0) --
47535
*  -- LAPACK computational routine (version 3.8.0) --
47537
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
47536
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
47538
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
47537
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
47539
*     December 2016
47538
*     November 2017
47540
*
47539
*
47541
*     .. Scalar Arguments ..
47540
*     .. Scalar Arguments ..
47542
      CHARACTER          HOWMNY, SIDE
47541
      CHARACTER          HOWMNY, SIDE
47543
      INTEGER            INFO, LDT, LDVL, LDVR, M, MM, N
47542
      INTEGER            INFO, LDT, LDVL, LDVR, M, MM, N
47544
*     ..
47543
*     ..
Line 47569... Line 47568...
47569
      INTEGER            IZAMAX
47568
      INTEGER            IZAMAX
47570
      DOUBLE PRECISION   DLAMCH, DZASUM
47569
      DOUBLE PRECISION   DLAMCH, DZASUM
47571
      EXTERNAL           LSAME, IZAMAX, DLAMCH, DZASUM
47570
      EXTERNAL           LSAME, IZAMAX, DLAMCH, DZASUM
47572
*     ..
47571
*     ..
47573
*     .. External Subroutines ..
47572
*     .. External Subroutines ..
47574
      EXTERNAL           XERBLA, ZCOPY, ZDSCAL, ZGEMV, ZLATRS
47573
      EXTERNAL           XERBLA, ZCOPY, ZDSCAL, ZGEMV, ZLATRS, DLABAD
47575
*     ..
47574
*     ..
47576
*     .. Intrinsic Functions ..
47575
*     .. Intrinsic Functions ..
47577
      INTRINSIC          ABS, DBLE, DCMPLX, DCONJG, DIMAG, MAX
47576
      INTRINSIC          ABS, DBLE, DCMPLX, DCONJG, DIMAG, MAX
47578
*     ..
47577
*     ..
47579
*     .. Statement Functions ..
47578
*     .. Statement Functions ..
Line 48021... Line 48020...
48021
*> \author Univ. of Tennessee
48020
*> \author Univ. of Tennessee
48022
*> \author Univ. of California Berkeley
48021
*> \author Univ. of California Berkeley
48023
*> \author Univ. of Colorado Denver
48022
*> \author Univ. of Colorado Denver
48024
*> \author NAG Ltd.
48023
*> \author NAG Ltd.
48025
*
48024
*
48026
*> \date June 2017
48025
*> \date November 2017
48027
*
48026
*
48028
*  @precisions fortran z -> c
48027
*  @precisions fortran z -> c
48029
*
48028
*
48030
*> \ingroup complex16OTHERcomputational
48029
*> \ingroup complex16OTHERcomputational
48031
*
48030
*
Line 48046... Line 48045...
48046
*  =====================================================================
48045
*  =====================================================================
48047
      SUBROUTINE ZTREVC3( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR,
48046
      SUBROUTINE ZTREVC3( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR,
48048
     $                    LDVR, MM, M, WORK, LWORK, RWORK, LRWORK, INFO)
48047
     $                    LDVR, MM, M, WORK, LWORK, RWORK, LRWORK, INFO)
48049
      IMPLICIT NONE
48048
      IMPLICIT NONE
48050
*
48049
*
48051
*  -- LAPACK computational routine (version 3.7.1) --
48050
*  -- LAPACK computational routine (version 3.8.0) --
48052
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
48051
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
48053
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
48052
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
48054
*     June 2017
48053
*     November 2017
48055
*
48054
*
48056
*     .. Scalar Arguments ..
48055
*     .. Scalar Arguments ..
48057
      CHARACTER          HOWMNY, SIDE
48056
      CHARACTER          HOWMNY, SIDE
48058
      INTEGER            INFO, LDT, LDVL, LDVR, LWORK, LRWORK, M, MM, N
48057
      INTEGER            INFO, LDT, LDVL, LDVR, LWORK, LRWORK, M, MM, N
48059
*     ..
48058
*     ..
Line 48087... Line 48086...
48087
      DOUBLE PRECISION   DLAMCH, DZASUM
48086
      DOUBLE PRECISION   DLAMCH, DZASUM
48088
      EXTERNAL           LSAME, ILAENV, IZAMAX, DLAMCH, DZASUM
48087
      EXTERNAL           LSAME, ILAENV, IZAMAX, DLAMCH, DZASUM
48089
*     ..
48088
*     ..
48090
*     .. External Subroutines ..
48089
*     .. External Subroutines ..
48091
      EXTERNAL           XERBLA, ZCOPY, ZDSCAL, ZGEMV, ZLATRS,
48090
      EXTERNAL           XERBLA, ZCOPY, ZDSCAL, ZGEMV, ZLATRS,
48092
     $                   ZGEMM, DLABAD, ZLASET
48091
     $                   ZGEMM, DLABAD, ZLASET, ZLACPY
48093
*     ..
48092
*     ..
48094
*     .. Intrinsic Functions ..
48093
*     .. Intrinsic Functions ..
48095
      INTRINSIC          ABS, DBLE, DCMPLX, CONJG, AIMAG, MAX
48094
      INTRINSIC          ABS, DBLE, DCMPLX, CONJG, AIMAG, MAX
48096
*     ..
48095
*     ..
48097
*     .. Statement Functions ..
48096
*     .. Statement Functions ..