76 COMPLEX*16 A( NMAX, NMAX ), T( NMAX, NMAX ), W( NMAX ),
89 COMMON / infoc / infot, nout, ok, lerr
90 COMMON / srnamc / srnamt
93 INTRINSIC dble, dcmplx
98 WRITE( nout, fmt = * )
104 a( i, j ) = 1.d0 / dcmplx( dble( i+j ), 0.d0 )
105 c( i, j ) = 1.d0 / dcmplx( dble( i+j ), 0.d0 )
106 t( i, j ) = 1.d0 / dcmplx( dble( i+j ), 0.d0 )
118 CALL zgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
119 CALL chkxer(
'ZGEQRT', infot, nout, lerr, ok )
121 CALL zgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'ZGEQRT', infot, nout, lerr, ok )
124 CALL zgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
125 CALL chkxer(
'ZGEQRT', infot, nout, lerr, ok )
127 CALL zgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
128 CALL chkxer(
'ZGEQRT', infot, nout, lerr, ok )
130 CALL zgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
131 CALL chkxer(
'ZGEQRT', infot, nout, lerr, ok )
137 CALL zgeqrt2( -1, 0, a, 1, t, 1, info )
138 CALL chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
140 CALL zgeqrt2( 0, -1, a, 1, t, 1, info )
141 CALL chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
143 CALL zgeqrt2( 2, 1, a, 1, t, 1, info )
144 CALL chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
146 CALL zgeqrt2( 2, 2, a, 2, t, 1, info )
147 CALL chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
153 CALL zgeqrt3( -1, 0, a, 1, t, 1, info )
154 CALL chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
156 CALL zgeqrt3( 0, -1, a, 1, t, 1, info )
157 CALL chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
159 CALL zgeqrt3( 2, 1, a, 1, t, 1, info )
160 CALL chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
162 CALL zgeqrt3( 2, 2, a, 2, t, 1, info )
163 CALL chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
169 CALL zgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
170 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
172 CALL zgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
175 CALL zgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
178 CALL zgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
181 CALL zgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
182 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
184 CALL zgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
187 CALL zgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
188 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
190 CALL zgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
191 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
193 CALL zgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
196 CALL zgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
197 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
199 CALL zgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
200 CALL chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
204 CALL alaesm( path, ok, nout )