LAPACK 3.12.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches
zerrps.f
Go to the documentation of this file.
1*> \brief \b ZERRPS
2*
3* =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6* http://www.netlib.org/lapack/explore-html/
7*
8* Definition:
9* ===========
10*
11* SUBROUTINE ZERRPS( PATH, NUNIT )
12*
13* .. Scalar Arguments ..
14* INTEGER NUNIT
15* CHARACTER*3 PATH
16* ..
17*
18*
19*> \par Purpose:
20* =============
21*>
22*> \verbatim
23*>
24*> ZERRPS tests the error exits for the COMPLEX routines
25*> for ZPSTRF.
26*> \endverbatim
27*
28* Arguments:
29* ==========
30*
31*> \param[in] PATH
32*> \verbatim
33*> PATH is CHARACTER*3
34*> The LAPACK path name for the routines to be tested.
35*> \endverbatim
36*>
37*> \param[in] NUNIT
38*> \verbatim
39*> NUNIT is INTEGER
40*> The unit number for output.
41*> \endverbatim
42*
43* Authors:
44* ========
45*
46*> \author Univ. of Tennessee
47*> \author Univ. of California Berkeley
48*> \author Univ. of Colorado Denver
49*> \author NAG Ltd.
50*
51*> \ingroup complex16_lin
52*
53* =====================================================================
54 SUBROUTINE zerrps( PATH, NUNIT )
55*
56* -- LAPACK test routine --
57* -- LAPACK is a software package provided by Univ. of Tennessee, --
58* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
59*
60* .. Scalar Arguments ..
61 INTEGER NUNIT
62 CHARACTER*3 PATH
63* ..
64*
65* =====================================================================
66*
67* .. Parameters ..
68 INTEGER NMAX
69 parameter( nmax = 4 )
70* ..
71* .. Local Scalars ..
72 INTEGER I, INFO, J, RANK
73* ..
74* .. Local Arrays ..
75 COMPLEX*16 A( NMAX, NMAX )
76 DOUBLE PRECISION RWORK( 2*NMAX )
77 INTEGER PIV( NMAX )
78* ..
79* .. External Subroutines ..
80 EXTERNAL alaesm, chkxer, zpstf2, zpstrf
81* ..
82* .. Scalars in Common ..
83 INTEGER INFOT, NOUT
84 LOGICAL LERR, OK
85 CHARACTER*32 SRNAMT
86* ..
87* .. Common blocks ..
88 COMMON / infoc / infot, nout, ok, lerr
89 COMMON / srnamc / srnamt
90* ..
91* .. Intrinsic Functions ..
92 INTRINSIC dble
93* ..
94* .. Executable Statements ..
95*
96 nout = nunit
97 WRITE( nout, fmt = * )
98*
99* Set the variables to innocuous values.
100*
101 DO 110 j = 1, nmax
102 DO 100 i = 1, nmax
103 a( i, j ) = 1.d0 / dble( i+j )
104*
105 100 CONTINUE
106 piv( j ) = j
107 rwork( j ) = 0.d0
108 rwork( nmax+j ) = 0.d0
109*
110 110 CONTINUE
111 ok = .true.
112*
113*
114* Test error exits of the routines that use the Cholesky
115* decomposition of an Hermitian positive semidefinite matrix.
116*
117* ZPSTRF
118*
119 srnamt = 'ZPSTRF'
120 infot = 1
121 CALL zpstrf( '/', 0, a, 1, piv, rank, -1.d0, rwork, info )
122 CALL chkxer( 'ZPSTRF', infot, nout, lerr, ok )
123 infot = 2
124 CALL zpstrf( 'U', -1, a, 1, piv, rank, -1.d0, rwork, info )
125 CALL chkxer( 'ZPSTRF', infot, nout, lerr, ok )
126 infot = 4
127 CALL zpstrf( 'U', 2, a, 1, piv, rank, -1.d0, rwork, info )
128 CALL chkxer( 'ZPSTRF', infot, nout, lerr, ok )
129*
130* ZPSTF2
131*
132 srnamt = 'ZPSTF2'
133 infot = 1
134 CALL zpstf2( '/', 0, a, 1, piv, rank, -1.d0, rwork, info )
135 CALL chkxer( 'ZPSTF2', infot, nout, lerr, ok )
136 infot = 2
137 CALL zpstf2( 'U', -1, a, 1, piv, rank, -1.d0, rwork, info )
138 CALL chkxer( 'ZPSTF2', infot, nout, lerr, ok )
139 infot = 4
140 CALL zpstf2( 'U', 2, a, 1, piv, rank, -1.d0, rwork, info )
141 CALL chkxer( 'ZPSTF2', infot, nout, lerr, ok )
142*
143*
144* Print a summary line.
145*
146 CALL alaesm( path, ok, nout )
147*
148 RETURN
149*
150* End of ZERRPS
151*
152 END
subroutine alaesm(path, ok, nout)
ALAESM
Definition alaesm.f:63
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3224
subroutine zpstf2(uplo, n, a, lda, piv, rank, tol, work, info)
ZPSTF2 computes the Cholesky factorization with complete pivoting of a complex Hermitian positive sem...
Definition zpstf2.f:142
subroutine zpstrf(uplo, n, a, lda, piv, rank, tol, work, info)
ZPSTRF computes the Cholesky factorization with complete pivoting of a complex Hermitian positive sem...
Definition zpstrf.f:142
subroutine zerrps(path, nunit)
ZERRPS
Definition zerrps.f:55