mathlib-lapack-arpack.c
1 //! @file mathlib-lapack-arpack.c
2 //
3 //! @section copyright
4 //
5 // This file is part of VIF - vintage fortran compiler.
6 // Copyright 2020-2025 J. Marcel van der Veer <algol68g@xs4all.nl>.
7 //
8 //! @section license
9 //
10 // This program is free software; you can redistribute it and/or modify it
11 // under the terms of the gnu general public license as published by the
12 // free software foundation; either version 3 of the license, or
13 // (at your option) any later version.
14 //
15 // This program is distributed in the hope that it will be useful, but
16 // without any warranty; without even the implied warranty of merchantability
17 // or fitness for a particular purpose. See the GNU general public license for
18 // more details. you should have received a copy of the GNU general public
19 // license along with this program. If not, see <http://www.gnu.org/licenses/>.
20 //
21 //! @Synopsis
22 //
23 //! 144 subprograms from LAPACK-ARPACK
24
25 // Compiled from Fortran source code by VIF.
26 // Selected subprograms are:
27 //
28 // CGBTF2 CLARNV DGBTRF DLAPY3 ICMAX1 SLAEV2 SLANV2 SLASR ZGBTF2 ZLARFG
29 // CGBTRF CLARTG DGBTRS DLARAN IZMAX1 SLAEXC SLAPTM SLASRT ZGBTRF ZLARNV
30 // CGBTRS CLASCL DGTTRF DLARFX SCSUM1 SLAGTM SLAPY2 SLASSQ ZGBTRS ZLARTG
31 // CGEQR2 CLASET DGTTRS DLARND SGBTF2 SLAHQR SLAPY3 SLASWP ZGEQR2 ZLASCL
32 // CGTTRF CLASSQ DLACON DLARNV SGBTRF SLALN2 SLARAN SLASY2 ZGTTRF ZLASET
33 // CGTTRS CLASWP DLADIV DLARUV SGBTRS SLAMC1 SLARF SORM2R ZGTTRS ZLASSQ
34 // CLACGV CLATRS DLAE2 DLASY2 SGEQR2 SLAMC2 SLARFG SPTTRF ZLACGV ZLASWP
35 // CLACON CMACH DLAEV2 DPTTRF SGTTRF SLAMC3 SLARFX SPTTRS ZLACON ZLATRS
36 // CLACPY CROT DLAEXC DPTTRS SGTTRS SLAMC4 SLARND SSTEQR ZLACPY ZROT
37 // CLADIV CTREVC DLAGTM DSTEQR SLABAD SLAMC5 SLARNV STREVC ZLADIV ZTREVC
38 // CLAHQR CTREXC DLAHQR DTREVC SLACON SLAMCH SLARTG STREXC ZLAHQR ZTREXC
39 // CLANGE CTRSEN DLALN2 DTREXC SLACPY SLANGE SLARUV STRSEN ZLANGE ZTRSEN
40 // CLANHS CTRSYL DLANHS DTRSEN SLADIV SLANHS SLASCL STRSYL ZLANHS ZTRSYL
41 // CLARF CUNM2R DLANV2 DTRSYL SLAE2 SLANST SLASET XLAENV ZLARF ZUNM2R
42 // CLARFG DGBTF2 DLAPTM DZSUM1
43
44 // LAPACK routines to support ARPACK.
45 // These LAPACK routines MUST be from the public release.
46
47 // The license for LAPACK Fortran source code is:
48 //
49 // BSD Software License
50 //
51 // Pertains to ARPACK and P_ARPACK
52 //
53 // Copyright (c) 1996-2008 Rice University.
54 // Developed by D.C. Sorensen, R.B. Lehoucq, C. Yang, and K. Maschhoff.
55 // All rights reserved.
56 //
57 // Redistribution and use in source and binary forms, with or without
58 // modification, are permitted provided that the following conditions are
59 // met:
60 //
61 // - Redistributions of source code must retain the above copyright
62 // notice, this list of conditions and the following disclaimer.
63 //
64 // - Redistributions in binary form must reproduce the above copyright
65 // notice, this list of conditions and the following disclaimer listed
66 // in this license in the documentation and/or other materials
67 // provided with the distribution.
68 //
69 // - Neither the name of the copyright holders nor the names of its
70 // contributors may be used to endorse or promote products derived from
71 // this software without specific prior written permission.
72 //
73 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
74 // "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
75 // LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
76 // A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
77 // OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
78 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
79 // LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
80 // DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
81 // THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
82 // (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
83 // OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
84
85 /*
86 Generated by VIF - experimental VIntage Fortran compiler.
87 VIF release 1.2.10
88 */
196 static CALLS __calls[__ncalls] = {
197 {"cgbtf2", 0}, // subroutine
198 {"cgbtrf", 0}, // subroutine
199 {"cgbtrs", 0}, // subroutine
200 {"cgeqr2", 0}, // subroutine
201 {"cgttrf", 0}, // subroutine
202 {"cgttrs", 0}, // subroutine
203 {"clacgv", 0}, // subroutine
204 {"clacon", 0}, // subroutine
205 {"clacpy", 0}, // subroutine
206 {"cladiv", 0}, // complex*8 function
207 {"clahqr", 0}, // subroutine
208 {"clange", 0}, // real*4 function
209 {"clanhs", 0}, // real*4 function
210 {"clarf", 0}, // subroutine
211 {"clarfg", 0}, // subroutine
212 {"clarnv", 0}, // subroutine
213 {"clartg", 0}, // subroutine
214 {"clascl", 0}, // subroutine
215 {"claset", 0}, // subroutine
216 {"classq", 0}, // subroutine
217 {"claswp", 0}, // subroutine
218 {"clatrs", 0}, // subroutine
219 {"cmach", 0}, // real*4 function
220 {"crot", 0}, // subroutine
221 {"ctrevc", 0}, // subroutine
222 {"ctrexc", 0}, // subroutine
223 {"ctrsen", 0}, // subroutine
224 {"ctrsyl", 0}, // subroutine
225 {"cunm2r", 0}, // subroutine
226 {"dgbtf2", 0}, // subroutine
227 {"dgbtrf", 0}, // subroutine
228 {"dgbtrs", 0}, // subroutine
229 {"dgttrf", 0}, // subroutine
230 {"dgttrs", 0}, // subroutine
231 {"dlacon", 0}, // subroutine
232 {"dladiv", 0}, // subroutine
233 {"dlae2", 0}, // subroutine
234 {"dlaev2", 0}, // subroutine
235 {"dlaexc", 0}, // subroutine
236 {"dlagtm", 0}, // subroutine
237 {"dlahqr", 0}, // subroutine
238 {"dlaln2", 0}, // subroutine
239 {"dlanhs", 0}, // real*8 function
240 {"dlanv2", 0}, // subroutine
241 {"dlaptm", 0}, // subroutine
242 {"dlapy3", 0}, // real*8 function
243 {"dlaran", 0}, // real*8 function
244 {"dlarfx", 0}, // subroutine
245 {"dlarnd", 0}, // real*8 function
246 {"dlarnv", 0}, // subroutine
247 {"dlaruv", 0}, // subroutine
248 {"dlasy2", 0}, // subroutine
249 {"dpttrf", 0}, // subroutine
250 {"dpttrs", 0}, // subroutine
251 {"dsteqr", 0}, // subroutine
252 {"dtrevc", 0}, // subroutine
253 {"dtrexc", 0}, // subroutine
254 {"dtrsen", 0}, // subroutine
255 {"dtrsyl", 0}, // subroutine
256 {"dzsum1", 0}, // real*8 function
257 {"icmax1", 0}, // integer*4 function
258 {"izmax1", 0}, // integer*4 function
259 {"scsum1", 0}, // real*4 function
260 {"sgbtf2", 0}, // subroutine
261 {"sgbtrf", 0}, // subroutine
262 {"sgbtrs", 0}, // subroutine
263 {"sgeqr2", 0}, // subroutine
264 {"sgttrf", 0}, // subroutine
265 {"sgttrs", 0}, // subroutine
266 {"slabad", 0}, // subroutine
267 {"slacon", 0}, // subroutine
268 {"slacpy", 0}, // subroutine
269 {"sladiv", 0}, // subroutine
270 {"slae2", 0}, // subroutine
271 {"slaev2", 0}, // subroutine
272 {"slaexc", 0}, // subroutine
273 {"slagtm", 0}, // subroutine
274 {"slahqr", 0}, // subroutine
275 {"slaln2", 0}, // subroutine
276 {"slamch", 0}, // real*4 function
277 {"slamc1", 0}, // subroutine
278 {"slamc2", 0}, // subroutine
279 {"slamc3", 0}, // real*4 function
280 {"slamc4", 0}, // subroutine
281 {"slamc5", 0}, // subroutine
282 {"slange", 0}, // real*4 function
283 {"slanhs", 0}, // real*4 function
284 {"slanst", 0}, // real*4 function
285 {"slanv2", 0}, // subroutine
286 {"slaptm", 0}, // subroutine
287 {"slapy2", 0}, // real*4 function
288 {"slapy3", 0}, // real*4 function
289 {"slaran", 0}, // real*4 function
290 {"slarf", 0}, // subroutine
291 {"slarfg", 0}, // subroutine
292 {"slarfx", 0}, // subroutine
293 {"slarnd", 0}, // real*4 function
294 {"slarnv", 0}, // subroutine
295 {"slartg", 0}, // subroutine
296 {"slaruv", 0}, // subroutine
297 {"slascl", 0}, // subroutine
298 {"slaset", 0}, // subroutine
299 {"slasr", 0}, // subroutine
300 {"slasrt", 0}, // subroutine
301 {"slassq", 0}, // subroutine
302 {"slaswp", 0}, // subroutine
303 {"slasy2", 0}, // subroutine
304 {"sorm2r", 0}, // subroutine
305 {"spttrf", 0}, // subroutine
306 {"spttrs", 0}, // subroutine
307 {"ssteqr", 0}, // subroutine
308 {"strevc", 0}, // subroutine
309 {"strexc", 0}, // subroutine
310 {"strsen", 0}, // subroutine
311 {"strsyl", 0}, // subroutine
312 {"xlaenv", 0}, // subroutine
313 {"zgbtf2", 0}, // subroutine
314 {"zgbtrf", 0}, // subroutine
315 {"zgbtrs", 0}, // subroutine
316 {"zgeqr2", 0}, // subroutine
317 {"zgttrf", 0}, // subroutine
318 {"zgttrs", 0}, // subroutine
319 {"zlacgv", 0}, // subroutine
320 {"zlacon", 0}, // subroutine
321 {"zlacpy", 0}, // subroutine
322 {"zladiv", 0}, // complex*16 function
323 {"zlahqr", 0}, // subroutine
324 {"zlange", 0}, // real*8 function
325 {"zlanhs", 0}, // real*8 function
326 {"zlarf", 0}, // subroutine
327 {"zlarfg", 0}, // subroutine
328 {"zlarnv", 0}, // subroutine
329 {"zlartg", 0}, // subroutine
330 {"zlascl", 0}, // subroutine
331 {"zlaset", 0}, // subroutine
332 {"zlassq", 0}, // subroutine
333 {"zlaswp", 0}, // subroutine
334 {"zlatrs", 0}, // subroutine
335 {"zrot", 0}, // subroutine
336 {"ztrevc", 0}, // subroutine
337 {"ztrexc", 0}, // subroutine
338 {"ztrsen", 0}, // subroutine
339 {"ztrsyl", 0}, // subroutine
340 {"zunm2r", 0}, // subroutine
341 {NULL, 0}
342 };
343 static struct {
344 int_4 iparms_[100];
345 } claenv_;
347 int_4 _p_ ipiv_, int_4 _p_ info_);
349 int_4 _p_ ipiv_, int_4 _p_ info_);
351 ab_, int_4 _p_ ldab_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
353 _p_ work_, int_4 _p_ info_);
355 int_4 _p_ ipiv_, int_4 _p_ info_);
357 _p_ du_, complex_8 _p_ du2_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
361 , int_4 _p_ ldb_);
364 complex_8 _p_ h_, int_4 _p_ ldh_, complex_8 _p_ w_, int_4 _p_ iloz_, int_4 _p_ ihiz_, complex_8 _p_ z_, int_4 _p_ ldz_,
365 int_4 _p_ info_);
367 work_);
370 tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_);
375 , int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
377 complex_8 _p_ a_, int_4 _p_ lda_);
380 int_4 _p_ incx_);
382 a_, int_4 _p_ lda_, complex_8 _p_ x_, real_4 _p_ scale_, real_4 _p_ cnorm_, int_4 _p_ info_);
385 c_, complex_8 _p_ s_);
387 _p_ ldt_, complex_8 _p_ vl_, int_4 _p_ ldvl_, complex_8 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_,
388 complex_8 _p_ work_, real_4 _p_ rwork_, int_4 _p_ info_);
390 ldq_, int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_);
392 _p_ ldt_, complex_8 _p_ q_, int_4 _p_ ldq_, complex_8 _p_ w_, int_4 _p_ m_, real_4 _p_ s_, real_4 _p_ sep_, complex_8
393 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
395 a_, int_4 _p_ lda_, complex_8 _p_ b_, int_4 _p_ ldb_, complex_8 _p_ c_, int_4 _p_ ldc_, real_4 _p_ scale_, int_4 _p_
396 info_);
398 int_4 _p_ lda_, complex_8 _p_ tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_, int_4 _p_ info_);
400 int_4 _p_ ipiv_, int_4 _p_ info_);
402 int_4 _p_ ipiv_, int_4 _p_ info_);
404 int_4 _p_ ldab_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
406 int_4 _p_ info_);
408 real_8 _p_ du2_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
410 ;
414 , real_8 _p_ sn1_);
416 ldq_, int_4 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ work_, int_4 _p_ info_);
418 d_, real_8 _p_ du_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_);
420 real_8 _p_ h_, int_4 _p_ ldh_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ iloz_, int_4 _p_ ihiz_, real_8 _p_ z_, int_4
421 _p_ ldz_, int_4 _p_ info_);
423 _p_ a_, int_4 _p_ lda_, real_8 _p_ d1_, real_8 _p_ d2_, real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ wr_, real_8 _p_ wi_,
424 real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ scale_, real_8 _p_ xnorm_, int_4 _p_ info_);
427 , real_8 _p_ rt2r_, real_8 _p_ rt2i_, real_8 _p_ cs_, real_8 _p_ sn_);
429 int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_);
433 int_4 _p_ ldc_, real_8 _p_ work_);
438 real_8 _p_ tl_, int_4 _p_ ldtl_, real_8 _p_ tr_, int_4 _p_ ldtr_, real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ scale_,
439 real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ xnorm_, int_4 _p_ info_);
442 int_4 _p_ info_);
444 real_8 _p_ work_, int_4 _p_ info_);
446 _p_ ldt_, real_8 _p_ vl_, int_4 _p_ ldvl_, real_8 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_
447 work_, int_4 _p_ info_);
449 int_4 _p_ ifst_, int_4 _p_ ilst_, real_8 _p_ work_, int_4 _p_ info_);
451 ldt_, real_8 _p_ q_, int_4 _p_ ldq_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ m_, real_8 _p_ s_, real_8 _p_ sep_,
452 real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_);
454 int_4 _p_ lda_, real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ scale_, int_4 _p_ info_);
460 int_4 _p_ ipiv_, int_4 _p_ info_);
462 int_4 _p_ ipiv_, int_4 _p_ info_);
464 int_4 _p_ ldab_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
466 int_4 _p_ info_);
468 int_4 _p_ info_);
470 real_4 _p_ du2_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
473 ;
475 int_4 _p_ ldb_);
479 , real_4 _p_ sn1_);
481 ldq_, int_4 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_ work_, int_4 _p_ info_);
483 d_, real_4 _p_ du_, real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_);
485 real_4 _p_ h_, int_4 _p_ ldh_, real_4 _p_ wr_, real_4 _p_ wi_, int_4 _p_ iloz_, int_4 _p_ ihiz_, real_4 _p_ z_, int_4
486 _p_ ldz_, int_4 _p_ info_);
488 _p_ a_, int_4 _p_ lda_, real_4 _p_ d1_, real_4 _p_ d2_, real_4 _p_ b_, int_4 _p_ ldb_, real_4 _p_ wr_, real_4 _p_ wi_,
489 real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ scale_, real_4 _p_ xnorm_, int_4 _p_ info_);
493 _p_ rmin_, int_4 _p_ emax_, real_4 _p_ rmax_);
497 _p_ rmax_);
502 , real_4 _p_ rt2r_, real_4 _p_ rt2i_, real_4 _p_ cs_, real_4 _p_ sn_);
504 int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_);
509 real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_);
512 int_4 _p_ ldc_, real_4 _p_ work_);
518 , int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
520 , int_4 _p_ lda_);
522 real_4 _p_ s_, real_4 _p_ a_, int_4 _p_ lda_);
526 int_4 _p_ incx_);
528 real_4 _p_ tl_, int_4 _p_ ldtl_, real_4 _p_ tr_, int_4 _p_ ldtr_, real_4 _p_ b_, int_4 _p_ ldb_, real_4 _p_ scale_,
529 real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ xnorm_, int_4 _p_ info_);
531 int_4 _p_ lda_, real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_, int_4 _p_ info_);
534 int_4 _p_ info_);
536 real_4 _p_ work_, int_4 _p_ info_);
538 _p_ ldt_, real_4 _p_ vl_, int_4 _p_ ldvl_, real_4 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_, real_4 _p_
539 work_, int_4 _p_ info_);
541 int_4 _p_ ifst_, int_4 _p_ ilst_, real_4 _p_ work_, int_4 _p_ info_);
543 ldt_, real_4 _p_ q_, int_4 _p_ ldq_, real_4 _p_ wr_, real_4 _p_ wi_, int_4 _p_ m_, real_4 _p_ s_, real_4 _p_ sep_,
544 real_4 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_);
546 int_4 _p_ lda_, real_4 _p_ b_, int_4 _p_ ldb_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ scale_, int_4 _p_ info_);
549 int_4 _p_ ipiv_, int_4 _p_ info_);
551 int_4 _p_ ipiv_, int_4 _p_ info_);
553 ab_, int_4 _p_ ldab_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
555 _p_ work_, int_4 _p_ info_);
557 int_4 _p_ ipiv_, int_4 _p_ info_);
559 complex_16 _p_ du_, complex_16 _p_ du2_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
563 b_, int_4 _p_ ldb_);
566 complex_16 _p_ h_, int_4 _p_ ldh_, complex_16 _p_ w_, int_4 _p_ iloz_, int_4 _p_ ihiz_, complex_16 _p_ z_, int_4 _p_
567 ldz_, int_4 _p_ info_);
569 work_);
572 tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_);
577 , int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
579 complex_16 _p_ a_, int_4 _p_ lda_);
582 , int_4 _p_ incx_);
584 _p_ a_, int_4 _p_ lda_, complex_16 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_);
586 _p_ c_, complex_16 _p_ s_);
588 int_4 _p_ ldt_, complex_16 _p_ vl_, int_4 _p_ ldvl_, complex_16 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_,
589 complex_16 _p_ work_, real_8 _p_ rwork_, int_4 _p_ info_);
591 ldq_, int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_);
593 _p_ ldt_, complex_16 _p_ q_, int_4 _p_ ldq_, complex_16 _p_ w_, int_4 _p_ m_, real_8 _p_ s_, real_8 _p_ sep_,
594 complex_16 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
596 a_, int_4 _p_ lda_, complex_16 _p_ b_, int_4 _p_ ldb_, complex_16 _p_ c_, int_4 _p_ ldc_, real_8 _p_ scale_, int_4 _p_
597 info_);
599 int_4 _p_ lda_, complex_16 _p_ tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_, int_4 _p_ info_);
600 int_4 _cgbtf2 (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, complex_8 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
601 ipiv_, int_4 _p_ info_)
602 { // ** body not listed **
675 }
676 int_4 _cgbtrf (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, complex_8 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
677 ipiv_, int_4 _p_ info_)
678 { // ** body not listed **
875 }
876 int_4 _cgbtrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, int_4 _p_ nrhs_, complex_8 _p_ ab_, int_4
877 _p_ ldab_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
878 { // ** body not listed **
979 }
980 int_4 _cgeqr2 (int_4 _p_ m_, int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, complex_8 _p_ tau_, complex_8 _p_ work_,
981 int_4 _p_ info_)
982 { // ** body not listed **
1017 }
1018 int_4 _cgttrf (int_4 _p_ n_, complex_8 _p_ dl_, complex_8 _p_ d_, complex_8 _p_ du_, complex_8 _p_ du2_, int_4 _p_
1019 ipiv_, int_4 _p_ info_)
1020 { // ** body not listed **
1068 }
1069 if (d_[(*n_) - 1] == CMPLXF (0.0, 0.0) && (*info_) == 0) {
1070 (*info_) = (*n_);
1071 goto _l0;
1072 }
1073 _l0:;
1074 __calls[4].calls++;
1075 return 0;
1076 }
1077 int_4 _cgttrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ nrhs_, complex_8 _p_ dl_, complex_8 _p_ d_, complex_8 _p_ du_,
1078 complex_8 _p_ du2_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
1079 { // ** body not listed **
1181 }
1182 int_4 _clacgv (int_4 _p_ n_, complex_8 _p_ x_, int_4 _p_ incx_)
1183 { // ** body not listed **
1202 }
1203 int_4 _clacon (int_4 _p_ n_, complex_8 _p_ v_, complex_8 _p_ x_, real_4 _p_ est_, int_4 _p_ kase_)
1204 { // ** body not listed **
1299 }
1300 int_4 _clacpy (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, complex_8 _p_ b_, int_4
1301 _p_ ldb_)
1302 { // ** body not listed **
1328 }
1329 complex_8 _cladiv (complex_8 _p_ x_, complex_8 _p_ y_)
1330 { // ** body not listed **
1339 }
1340 int_4 _clahqr (logical_4 _p_ wantt_, logical_4 _p_ wantz_, int_4 _p_ n_, int_4 _p_ ilo_, int_4 _p_ ihi_, complex_8 _p_
1341 h_, int_4 _p_ ldh_, complex_8 _p_ w_, int_4 _p_ iloz_, int_4 _p_ ihiz_, complex_8 _p_ z_, int_4 _p_ ldz_, int_4 _p_
1342 info_)
1343 { // ** body not listed **
1515 }
1516 real_4 _clange (char _p_ norm_, int_4 _p_ m_, int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
1517 { // ** body not listed **
1568 }
1569 real_4 _clanhs (char _p_ norm_, int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
1570 { // ** body not listed **
1622 }
1623 int_4 _clarf (char _p_ side_, int_4 _p_ m_, int_4 _p_ n_, complex_8 _p_ v_, int_4 _p_ incv_, complex_8 _p_ tau_,
1624 complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_)
1625 { // ** body not listed **
1646 }
1647 int_4 _clarfg (int_4 _p_ n_, complex_8 _p_ alpha_, complex_8 _p_ x_, int_4 _p_ incx_, complex_8 _p_ tau_)
1648 { // ** body not listed **
1702 }
1703 int_4 _clarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, complex_8 _p_ x_)
1704 { // ** body not listed **
1739 }
1740 int_4 _clartg (complex_8 _p_ f_, complex_8 _p_ g_, real_4 _p_ cs_, complex_8 _p_ sn_, complex_8 _p_ r_)
1741 { // ** body not listed **
1785 }
1786 int_4 _clascl (char _p_ type_, int_4 _p_ kl_, int_4 _p_ ku_, real_4 _p_ cfrom_, real_4 _p_ cto_, int_4 _p_ m_, int_4
1787 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
1788 { // ** body not listed **
1938 }
1939 int_4 _claset (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, complex_8 _p_ alpha_, complex_8 _p_ beta_, complex_8 _p_ a_,
1940 int_4 _p_ lda_)
1941 { // ** body not listed **
1976 }
1977 int_4 _classq (int_4 _p_ n_, complex_8 _p_ x_, int_4 _p_ incx_, real_4 _p_ scale_, real_4 _p_ sumsq_)
1978 { // ** body not listed **
2009 }
2010 int_4 _claswp (int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ k1_, int_4 _p_ k2_, int_4 _p_ ipiv_, int_4 _p_
2011 incx_)
2012 { // ** body not listed **
2053 }
2054 int_4 _clatrs (char _p_ uplo_, char _p_ trans_, char _p_ diag_, char _p_ normin_, int_4 _p_ n_, complex_8 _p_ a_, int_4
2055 _p_ lda_, complex_8 _p_ x_, real_4 _p_ scale_, real_4 _p_ cnorm_, int_4 _p_ info_)
2056 { // ** body not listed **
2525 }
2526 real_4 _cmach (int_4 _p_ job_)
2527 { // ** body not listed **
2567 }
2568 int_4 _crot (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_, complex_8 _p_ cy_, int_4 _p_ incy_, real_4 _p_ c_,
2569 complex_8 _p_ s_)
2570 { // ** body not listed **
2604 }
2605 int_4 _ctrevc (char _p_ side_, char _p_ howmny_, logical_4 _p_ select_, int_4 _p_ n_, complex_8 _p_ t_, int_4 _p_ ldt_,
2606 complex_8 _p_ vl_, int_4 _p_ ldvl_, complex_8 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_, complex_8 _p_
2607 work_, real_4 _p_ rwork_, int_4 _p_ info_)
2608 { // ** body not listed **
2731 }
2732 if (leftv_) {
2733 is_ = 1;
2734 for (ki_ = 1; ki_ <= (*n_); (ki_)++) {
2735 if (somev_) {
2736 if (! (select_[ki_ - 1])) {
2737 goto _l130;
2738 }
2739 }
2740 smin_ = _max (ulp_ * (({complex_8 cdum_1_ = t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))]; _abs ((real_4) (crealf (cdum_1_))) +
2741 _abs (cimagf (cdum_1_));})), smlnum_);
2742 work_[(*n_) - 1] = CMPLXF (1.0, 0.0);
2743 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2744 work_[k_ - 1] = -conjf (t_[ki_ - 1 + ((*ldt_)) * ((k_ - 1))]);
2745 }
2746 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2747 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] - t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))];
2748 if (({complex_8 cdum_1_ = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))]; _abs ((real_4) (crealf (cdum_1_))) + _abs (cimagf
2749 (cdum_1_));}) < smin_) {
2750 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = smin_;
2751 }
2752 }
2753 if (ki_ < (*n_)) {
2754 (void) _clatrs (_dc_11, _dc_14, _dc_12, _dc_35, ($6_ = (*n_) - ki_, &$6_), &t_[ki_ + 1 - 1 + ((*ldt_)) * ((ki_ + 1 -
2755 1))], ldt_, &work_[ki_ + 1 - 1], &scale_, rwork_, info_);
2756 work_[ki_ - 1] = scale_;
2757 }
2758 if (! (over_)) {
2759 (void) _ccopy (($7_ = (*n_) - ki_ + 1, &$7_), &work_[ki_ - 1], &_k1, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
2760 ii_ = _icamax (($8_ = (*n_) - ki_ + 1, &$8_), &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1) + ki_ - 1;
2761 remax_ = 1.0 / ({complex_8 cdum_1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((is_ - 1))]; _abs ((real_4) (crealf (cdum_1_))) + _abs
2762 (cimagf (cdum_1_));});
2763 (void) _csscal (($9_ = (*n_) - ki_ + 1, &$9_), &remax_, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
2764 for (k_ = 1; k_ <= ki_ - 1; (k_)++) {
2765 vl_[k_ - 1 + ((*ldvl_)) * ((is_ - 1))] = CMPLXF (0.0, 0.0);
2766 }
2767 }
2768 else {
2769 if (ki_ < (*n_)) {
2770 (void) _cgemv (_dc_7, n_, ($10_ = (*n_) - ki_, &$10_), ($11_ = CMPLXF (1.0, 0.0), &$11_), &vl_[0 + ((*ldvl_)) * ((ki_ +
2771 1 - 1))], ldvl_, &work_[ki_ + 1 - 1], &_k1, ($12_ = (complex_8) (scale_), &$12_), &vl_[0 + ((*ldvl_)) * ((ki_ - 1))],
2772 &_k1);
2773 }
2774 ii_ = _icamax (n_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
2775 remax_ = 1.0 / ({complex_8 cdum_1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((ki_ - 1))]; _abs ((real_4) (crealf (cdum_1_))) + _abs
2776 (cimagf (cdum_1_));});
2777 (void) _csscal (n_, &remax_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
2778 }
2779 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2780 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = work_[k_ + (*n_) - 1];
2781 }
2782 is_ = is_ + 1;
2783 _l130:;
2784 }
2785 }
2786 _l0:;
2787 __calls[24].calls++;
2788 return 0;
2789 }
2790 int_4 _ctrexc (char _p_ compq_, int_4 _p_ n_, complex_8 _p_ t_, int_4 _p_ ldt_, complex_8 _p_ q_, int_4 _p_ ldq_, int_4
2791 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_)
2792 { // ** body not listed **
2858 }
2859 int_4 _ctrsen (char _p_ job_, char _p_ compq_, logical_4 _p_ select_, int_4 _p_ n_, complex_8 _p_ t_, int_4 _p_ ldt_,
2860 complex_8 _p_ q_, int_4 _p_ ldq_, complex_8 _p_ w_, int_4 _p_ m_, real_4 _p_ s_, real_4 _p_ sep_, complex_8 _p_ work_,
2861 int_4 _p_ lwork_, int_4 _p_ info_)
2862 { // ** body not listed **
2960 }
2961 int_4 _ctrsyl (char _p_ trana_, char _p_ tranb_, int_4 _p_ isgn_, int_4 _p_ m_, int_4 _p_ n_, complex_8 _p_ a_, int_4
2962 _p_ lda_, complex_8 _p_ b_, int_4 _p_ ldb_, complex_8 _p_ c_, int_4 _p_ ldc_, real_4 _p_ scale_, int_4 _p_ info_)
2963 { // ** body not listed **
3150 }
3151 int_4 _cunm2r (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, complex_8 _p_ a_, int_4 _p_
3152 lda_, complex_8 _p_ tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_, int_4 _p_ info_)
3153 { // ** body not listed **
3239 }
3240 int_4 _dgbtf2 (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, real_8 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
3241 ipiv_, int_4 _p_ info_)
3242 { // ** body not listed **
3315 }
3316 int_4 _dgbtrf (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, real_8 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
3317 ipiv_, int_4 _p_ info_)
3318 { // ** body not listed **
3514 }
3515 int_4 _dgbtrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, int_4 _p_ nrhs_, real_8 _p_ ab_, int_4 _p_
3516 ldab_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
3517 { // ** body not listed **
3597 }
3598 int_4 _dgttrf (int_4 _p_ n_, real_8 _p_ dl_, real_8 _p_ d_, real_8 _p_ du_, real_8 _p_ du2_, int_4 _p_ ipiv_, int_4 _p_
3599 info_)
3600 { // ** body not listed **
3655 }
3656 int_4 _dgttrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ dl_, real_8 _p_ d_, real_8 _p_ du_, real_8
3657 _p_ du2_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
3658 { // ** body not listed **
3736 }
3737 int_4 _dlacon (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ x_, int_4 _p_ isgn_, real_8 _p_ est_, int_4 _p_ kase_)
3738 { // ** body not listed **
3833 }
3834 int_4 _dladiv (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ d_, real_8 _p_ p_, real_8 _p_ q_)
3835 { // ** body not listed **
3851 }
3852 int_4 _dlae2 (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ rt1_, real_8 _p_ rt2_)
3853 { // ** body not listed **
3892 }
3893 int_4 _dlaev2 (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ rt1_, real_8 _p_ rt2_, real_8 _p_ cs1_, real_8
3894 _p_ sn1_)
3895 { // ** body not listed **
3968 }
3969 int_4 _dlaexc (logical_4 _p_ wantq_, int_4 _p_ n_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_ q_, int_4 _p_ ldq_, int_4
3970 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ work_, int_4 _p_ info_)
3971 { // ** body not listed **
4141 }
4142 int_4 _dlagtm (char _p_ trans_, int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ alpha_, real_8 _p_ dl_, real_8 _p_ d_, real_8
4143 _p_ du_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_)
4144 { // ** body not listed **
4243 }
4244 int_4 _dlahqr (logical_4 _p_ wantt_, logical_4 _p_ wantz_, int_4 _p_ n_, int_4 _p_ ilo_, int_4 _p_ ihi_, real_8 _p_ h_,
4245 int_4 _p_ ldh_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ iloz_, int_4 _p_ ihiz_, real_8 _p_ z_, int_4 _p_ ldz_, int_4
4246 _p_ info_)
4247 { // ** body not listed **
4442 }
4443 int_4 _dlaln2 (logical_4 _p_ ltrans_, int_4 _p_ na_, int_4 _p_ nw_, real_8 _p_ smin_, real_8 _p_ ca_, real_8 _p_ a_,
4444 int_4 _p_ lda_, real_8 _p_ d1_, real_8 _p_ d2_, real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ wr_, real_8 _p_ wi_, real_8
4445 _p_ x_, int_4 _p_ ldx_, real_8 _p_ scale_, real_8 _p_ xnorm_, int_4 _p_ info_)
4446 { // ** body not listed **
4771 }
4772 real_8 _dlanhs (char _p_ norm_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
4773 { // ** body not listed **
4826 }
4827 int_4 _dlanv2 (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ d_, real_8 _p_ rt1r_, real_8 _p_ rt1i_, real_8
4828 _p_ rt2r_, real_8 _p_ rt2i_, real_8 _p_ cs_, real_8 _p_ sn_)
4829 { // ** body not listed **
4911 }
4912 int_4 _dlaptm (int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ alpha_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ x_, int_4 _p_
4913 ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_)
4914 { // ** body not listed **
4972 }
4973 real_8 _dlapy3 (real_8 _p_ x_, real_8 _p_ y_, real_8 _p_ z_)
4974 { // ** body not listed **
4991 }
4992 real_8 _dlaran (int_4 _p_ iseed_)
4993 { // ** body not listed **
5015 }
5016 int_4 _dlarfx (char _p_ side_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_
5017 ldc_, real_8 _p_ work_)
5018 { // ** body not listed **
5525 }
5526 real_8 _dlarnd (int_4 _p_ idist_, int_4 _p_ iseed_)
5527 { // ** body not listed **
5543 }
5544 int_4 _dlarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, real_8 _p_ x_)
5545 { // ** body not listed **
5575 }
5576 int_4 _dlaruv (int_4 _p_ iseed_, int_4 _p_ n_, real_8 _p_ x_)
5577 { // ** body not listed **
7666 }
7667 int_4 _dlasy2 (logical_4 _p_ ltranl_, logical_4 _p_ ltranr_, int_4 _p_ isgn_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_
7668 tl_, int_4 _p_ ldtl_, real_8 _p_ tr_, int_4 _p_ ldtr_, real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ scale_, real_8 _p_ x_,
7669 int_4 _p_ ldx_, real_8 _p_ xnorm_, int_4 _p_ info_)
7670 { // ** body not listed **
7980 }
7981 int_4 _dpttrf (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ info_)
7982 { // ** body not listed **
8015 }
8016 int_4 _dpttrs (int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_
8017 info_)
8018 { // ** body not listed **
8053 }
8054 int_4 _dsteqr (char _p_ compz_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ z_, int_4 _p_ ldz_, real_8 _p_
8055 work_, int_4 _p_ info_)
8056 { // ** body not listed **
8385 }
8386 int_4 _dtrevc (char _p_ side_, char _p_ howmny_, logical_4 _p_ select_, int_4 _p_ n_, real_8 _p_ t_, int_4 _p_ ldt_,
8387 real_8 _p_ vl_, int_4 _p_ ldvl_, real_8 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ work_, int_4
8388 _p_ info_)
8389 { // ** body not listed **
9002 }
9003 int_4 _dtrexc (char _p_ compq_, int_4 _p_ n_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_ q_, int_4 _p_ ldq_, int_4 _p_
9004 ifst_, int_4 _p_ ilst_, real_8 _p_ work_, int_4 _p_ info_)
9005 { // ** body not listed **
9195 }
9196 int_4 _dtrsen (char _p_ job_, char _p_ compq_, logical_4 _p_ select_, int_4 _p_ n_, real_8 _p_ t_, int_4 _p_ ldt_,
9197 real_8 _p_ q_, int_4 _p_ ldq_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ m_, real_8 _p_ s_, real_8 _p_ sep_, real_8 _p_
9198 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
9199 { // ** body not listed **
9361 }
9362 int_4 _dtrsyl (char _p_ trana_, char _p_ tranb_, int_4 _p_ isgn_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_
9363 lda_, real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ scale_, int_4 _p_ info_)
9364 { // ** body not listed **
10102 }
10103 real_8 _dzsum1 (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_)
10104 { // ** body not listed **
10130 }
10131 int_4 _icmax1 (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_)
10132 { // ** body not listed **
10159 }
10160 goto _l0;
10161 _l30:;
10162 smax_ = ({complex_8 zdum_1_ = cx_[0]; _abs ((real_4) (crealf (zdum_1_)));});
10163 for (i_ = 2; i_ <= (*n_); (i_)++) {
10164 if (({complex_8 zdum_1_ = cx_[i_ - 1]; _abs ((real_4) (crealf (zdum_1_)));}) <= smax_) {
10165 goto _l40;
10166 }
10167 icmax1_ = i_;
10168 smax_ = ({complex_8 zdum_1_ = cx_[i_ - 1]; _abs ((real_4) (crealf (zdum_1_)));});
10169 _l40:;
10170 }
10171 _l0:;
10172 __calls[60].calls++;
10173 return icmax1_;
10174 }
10175 int_4 _izmax1 (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_)
10176 { // ** body not listed **
10203 }
10204 goto _l0;
10205 _l30:;
10206 smax_ = ({complex_16 zdum_1_ = cx_[0]; _abs ((real_8) (creal (zdum_1_)));});
10207 for (i_ = 2; i_ <= (*n_); (i_)++) {
10208 if (({complex_16 zdum_1_ = cx_[i_ - 1]; _abs ((real_8) (creal (zdum_1_)));}) <= smax_) {
10209 goto _l40;
10210 }
10211 izmax1_ = i_;
10212 smax_ = ({complex_16 zdum_1_ = cx_[i_ - 1]; _abs ((real_8) (creal (zdum_1_)));});
10213 _l40:;
10214 }
10215 _l0:;
10216 __calls[61].calls++;
10217 return izmax1_;
10218 }
10219 real_4 _scsum1 (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_)
10220 { // ** body not listed **
10246 }
10247 int_4 _sgbtf2 (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, real_4 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
10248 ipiv_, int_4 _p_ info_)
10249 { // ** body not listed **
10322 }
10323 int_4 _sgbtrf (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, real_4 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
10324 ipiv_, int_4 _p_ info_)
10325 { // ** body not listed **
10520 }
10521 int_4 _sgbtrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, int_4 _p_ nrhs_, real_4 _p_ ab_, int_4 _p_
10522 ldab_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
10523 { // ** body not listed **
10603 }
10604 int_4 _sgeqr2 (int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, real_4 _p_ tau_, real_4 _p_ work_, int_4 _p_
10605 info_)
10606 { // ** body not listed **
10640 }
10641 int_4 _sgttrf (int_4 _p_ n_, real_4 _p_ dl_, real_4 _p_ d_, real_4 _p_ du_, real_4 _p_ du2_, int_4 _p_ ipiv_, int_4 _p_
10642 info_)
10643 { // ** body not listed **
10698 }
10699 int_4 _sgttrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ nrhs_, real_4 _p_ dl_, real_4 _p_ d_, real_4 _p_ du_, real_4
10700 _p_ du2_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
10701 { // ** body not listed **
10779 }
10780 int_4 _slabad (real_4 _p_ small_, real_4 _p_ large_)
10781 { // ** body not listed **
10788 }
10789 int_4 _slacon (int_4 _p_ n_, real_4 _p_ v_, real_4 _p_ x_, int_4 _p_ isgn_, real_4 _p_ est_, int_4 _p_ kase_)
10790 { // ** body not listed **
10885 }
10886 int_4 _slacpy (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, real_4 _p_ b_, int_4 _p_ ldb_
10887 )
10888 { // ** body not listed **
10914 }
10915 int_4 _sladiv (real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ d_, real_4 _p_ p_, real_4 _p_ q_)
10916 { // ** body not listed **
10932 }
10933 int_4 _slae2 (real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ rt1_, real_4 _p_ rt2_)
10934 { // ** body not listed **
10973 }
10974 int_4 _slaev2 (real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ rt1_, real_4 _p_ rt2_, real_4 _p_ cs1_, real_4
10975 _p_ sn1_)
10976 { // ** body not listed **
11049 }
11050 int_4 _slaexc (logical_4 _p_ wantq_, int_4 _p_ n_, real_4 _p_ t_, int_4 _p_ ldt_, real_4 _p_ q_, int_4 _p_ ldq_, int_4
11051 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_ work_, int_4 _p_ info_)
11052 { // ** body not listed **
11217 }
11218 int_4 _slagtm (char _p_ trans_, int_4 _p_ n_, int_4 _p_ nrhs_, real_4 _p_ alpha_, real_4 _p_ dl_, real_4 _p_ d_, real_4
11219 _p_ du_, real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_)
11220 { // ** body not listed **
11319 }
11320 int_4 _slahqr (logical_4 _p_ wantt_, logical_4 _p_ wantz_, int_4 _p_ n_, int_4 _p_ ilo_, int_4 _p_ ihi_, real_4 _p_ h_,
11321 int_4 _p_ ldh_, real_4 _p_ wr_, real_4 _p_ wi_, int_4 _p_ iloz_, int_4 _p_ ihiz_, real_4 _p_ z_, int_4 _p_ ldz_, int_4
11322 _p_ info_)
11323 { // ** body not listed **
11515 }
11516 int_4 _slaln2 (logical_4 _p_ ltrans_, int_4 _p_ na_, int_4 _p_ nw_, real_4 _p_ smin_, real_4 _p_ ca_, real_4 _p_ a_,
11517 int_4 _p_ lda_, real_4 _p_ d1_, real_4 _p_ d2_, real_4 _p_ b_, int_4 _p_ ldb_, real_4 _p_ wr_, real_4 _p_ wi_, real_4
11518 _p_ x_, int_4 _p_ ldx_, real_4 _p_ scale_, real_4 _p_ xnorm_, int_4 _p_ info_)
11519 { // ** body not listed **
11843 }
11844 real_4 _slamch (char _p_ cmach_)
11845 { // ** body not listed **
11918 }
11919 int_4 _slamc1 (int_4 _p_ beta_, int_4 _p_ t_, logical_4 _p_ rnd_, logical_4 _p_ ieee1_)
11920 { // ** body not listed **
11996 }
11997 int_4 _slamc2 (int_4 _p_ beta_, int_4 _p_ t_, logical_4 _p_ rnd_, real_4 _p_ eps_, int_4 _p_ emin_, real_4 _p_ rmin_,
11998 int_4 _p_ emax_, real_4 _p_ rmax_)
11999 { // ** body not listed **
12172 }
12173 real_4 _slamc3 (real_4 _p_ a_, real_4 _p_ b_)
12174 { // ** body not listed **
12179 }
12180 int_4 _slamc4 (int_4 _p_ emin_, real_4 _p_ start_, int_4 _p_ base_)
12181 { // ** body not listed **
12215 }
12216 int_4 _slamc5 (int_4 _p_ beta_, int_4 _p_ p_, int_4 _p_ emin_, logical_4 _p_ ieee_, int_4 _p_ emax_, real_4 _p_ rmax_)
12217 { // ** body not listed **
12270 }
12271 real_4 _slange (char _p_ norm_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
12272 { // ** body not listed **
12323 }
12324 real_4 _slanhs (char _p_ norm_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
12325 { // ** body not listed **
12377 }
12378 real_4 _slanst (char _p_ norm_, int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_)
12379 { // ** body not listed **
12419 }
12420 int_4 _slanv2 (real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ d_, real_4 _p_ rt1r_, real_4 _p_ rt1i_, real_4
12421 _p_ rt2r_, real_4 _p_ rt2i_, real_4 _p_ cs_, real_4 _p_ sn_)
12422 { // ** body not listed **
12503 }
12504 int_4 _slaptm (int_4 _p_ n_, int_4 _p_ nrhs_, real_4 _p_ alpha_, real_4 _p_ d_, real_4 _p_ e_, real_4 _p_ x_, int_4 _p_
12505 ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_)
12506 { // ** body not listed **
12564 }
12565 real_4 _slapy2 (real_4 _p_ x_, real_4 _p_ y_)
12566 { // ** body not listed **
12582 }
12583 real_4 _slapy3 (real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ z_)
12584 { // ** body not listed **
12601 }
12602 real_4 _slaran (int_4 _p_ iseed_)
12603 { // ** body not listed **
12625 }
12626 int_4 _slarf (char _p_ side_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ v_, int_4 _p_ incv_, real_4 _p_ tau_, real_4 _p_
12627 c_, int_4 _p_ ldc_, real_4 _p_ work_)
12628 { // ** body not listed **
12647 }
12648 int_4 _slarfg (int_4 _p_ n_, real_4 _p_ alpha_, real_4 _p_ x_, int_4 _p_ incx_, real_4 _p_ tau_)
12649 { // ** body not listed **
12696 }
12697 int_4 _slarfx (char _p_ side_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ v_, real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_
12698 ldc_, real_4 _p_ work_)
12699 { // ** body not listed **
13206 }
13207 real_4 _slarnd (int_4 _p_ idist_, int_4 _p_ iseed_)
13208 { // ** body not listed **
13224 }
13225 int_4 _slarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, real_4 _p_ x_)
13226 { // ** body not listed **
13256 }
13257 int_4 _slartg (real_4 _p_ f_, real_4 _p_ g_, real_4 _p_ cs_, real_4 _p_ sn_, real_4 _p_ r_)
13258 { // ** body not listed **
13342 }
13343 int_4 _slaruv (int_4 _p_ iseed_, int_4 _p_ n_, real_4 _p_ x_)
13344 { // ** body not listed **
15434 }
15435 int_4 _slascl (char _p_ type_, int_4 _p_ kl_, int_4 _p_ ku_, real_4 _p_ cfrom_, real_4 _p_ cto_, int_4 _p_ m_, int_4
15436 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
15437 { // ** body not listed **
15587 }
15588 int_4 _slaset (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ alpha_, real_4 _p_ beta_, real_4 _p_ a_, int_4
15589 _p_ lda_)
15590 { // ** body not listed **
15619 }
15620 int_4 _slasr (char _p_ side_, char _p_ pivot_, char _p_ direct_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ c_, real_4 _p_
15621 s_, real_4 _p_ a_, int_4 _p_ lda_)
15622 { // ** body not listed **
15828 }
15829 int_4 _slasrt (char _p_ id_, int_4 _p_ n_, real_4 _p_ d_, int_4 _p_ info_)
15830 { // ** body not listed **
16004 }
16005 int_4 _slassq (int_4 _p_ n_, real_4 _p_ x_, int_4 _p_ incx_, real_4 _p_ scale_, real_4 _p_ sumsq_)
16006 { // ** body not listed **
16027 }
16028 int_4 _slaswp (int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ k1_, int_4 _p_ k2_, int_4 _p_ ipiv_, int_4 _p_
16029 incx_)
16030 { // ** body not listed **
16071 }
16072 int_4 _slasy2 (logical_4 _p_ ltranl_, logical_4 _p_ ltranr_, int_4 _p_ isgn_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_
16073 tl_, int_4 _p_ ldtl_, real_4 _p_ tr_, int_4 _p_ ldtr_, real_4 _p_ b_, int_4 _p_ ldb_, real_4 _p_ scale_, real_4 _p_ x_,
16074 int_4 _p_ ldx_, real_4 _p_ xnorm_, int_4 _p_ info_)
16075 { // ** body not listed **
16384 }
16385 int_4 _sorm2r (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_4 _p_ a_, int_4 _p_ lda_
16386 , real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_, int_4 _p_ info_)
16387 { // ** body not listed **
16467 }
16468 int_4 _spttrf (int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_, int_4 _p_ info_)
16469 { // ** body not listed **
16502 }
16503 int_4 _spttrs (int_4 _p_ n_, int_4 _p_ nrhs_, real_4 _p_ d_, real_4 _p_ e_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_
16504 info_)
16505 { // ** body not listed **
16540 }
16541 int_4 _ssteqr (char _p_ compz_, int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_, real_4 _p_ z_, int_4 _p_ ldz_, real_4 _p_
16542 work_, int_4 _p_ info_)
16543 { // ** body not listed **
16864 }
16865 int_4 _strevc (char _p_ side_, char _p_ howmny_, logical_4 _p_ select_, int_4 _p_ n_, real_4 _p_ t_, int_4 _p_ ldt_,
16866 real_4 _p_ vl_, int_4 _p_ ldvl_, real_4 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_, real_4 _p_ work_, int_4
16867 _p_ info_)
16868 { // ** body not listed **
17479 }
17480 int_4 _strexc (char _p_ compq_, int_4 _p_ n_, real_4 _p_ t_, int_4 _p_ ldt_, real_4 _p_ q_, int_4 _p_ ldq_, int_4 _p_
17481 ifst_, int_4 _p_ ilst_, real_4 _p_ work_, int_4 _p_ info_)
17482 { // ** body not listed **
17672 }
17673 int_4 _strsen (char _p_ job_, char _p_ compq_, logical_4 _p_ select_, int_4 _p_ n_, real_4 _p_ t_, int_4 _p_ ldt_,
17674 real_4 _p_ q_, int_4 _p_ ldq_, real_4 _p_ wr_, real_4 _p_ wi_, int_4 _p_ m_, real_4 _p_ s_, real_4 _p_ sep_, real_4 _p_
17675 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
17676 { // ** body not listed **
17836 }
17837 int_4 _strsyl (char _p_ trana_, char _p_ tranb_, int_4 _p_ isgn_, int_4 _p_ m_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_
17838 lda_, real_4 _p_ b_, int_4 _p_ ldb_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ scale_, int_4 _p_ info_)
17839 { // ** body not listed **
18574 }
18575 int_4 _xlaenv (int_4 _p_ ispec_, int_4 _p_ nvalue_)
18576 { // ** body not listed **
18582 }
18583 int_4 _zgbtf2 (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, complex_16 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
18584 ipiv_, int_4 _p_ info_)
18585 { // ** body not listed **
18658 }
18659 int_4 _zgbtrf (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, complex_16 _p_ ab_, int_4 _p_ ldab_, int_4 _p_
18660 ipiv_, int_4 _p_ info_)
18661 { // ** body not listed **
18858 }
18859 int_4 _zgbtrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ kl_, int_4 _p_ ku_, int_4 _p_ nrhs_, complex_16 _p_ ab_, int_4
18860 _p_ ldab_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
18861 { // ** body not listed **
18962 }
18963 int_4 _zgeqr2 (int_4 _p_ m_, int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, complex_16 _p_ tau_, complex_16 _p_ work_
18964 , int_4 _p_ info_)
18965 { // ** body not listed **
19002 }
19003 int_4 _zgttrf (int_4 _p_ n_, complex_16 _p_ dl_, complex_16 _p_ d_, complex_16 _p_ du_, complex_16 _p_ du2_, int_4 _p_
19004 ipiv_, int_4 _p_ info_)
19005 { // ** body not listed **
19053 }
19054 if (d_[(*n_) - 1] == CMPLX (0.0, 0.0) && (*info_) == 0) {
19055 (*info_) = (*n_);
19056 goto _l0;
19057 }
19058 _l0:;
19059 __calls[120].calls++;
19060 return 0;
19061 }
19062 int_4 _zgttrs (char _p_ trans_, int_4 _p_ n_, int_4 _p_ nrhs_, complex_16 _p_ dl_, complex_16 _p_ d_, complex_16 _p_
19063 du_, complex_16 _p_ du2_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
19064 { // ** body not listed **
19169 }
19170 int_4 _zlacgv (int_4 _p_ n_, complex_16 _p_ x_, int_4 _p_ incx_)
19171 { // ** body not listed **
19192 }
19193 int_4 _zlacon (int_4 _p_ n_, complex_16 _p_ v_, complex_16 _p_ x_, real_8 _p_ est_, int_4 _p_ kase_)
19194 { // ** body not listed **
19290 }
19291 int_4 _zlacpy (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, complex_16 _p_ b_, int_4
19292 _p_ ldb_)
19293 { // ** body not listed **
19319 }
19320 complex_16 _zladiv (complex_16 _p_ x_, complex_16 _p_ y_)
19321 { // ** body not listed **
19330 }
19331 int_4 _zlahqr (logical_4 _p_ wantt_, logical_4 _p_ wantz_, int_4 _p_ n_, int_4 _p_ ilo_, int_4 _p_ ihi_, complex_16 _p_
19332 h_, int_4 _p_ ldh_, complex_16 _p_ w_, int_4 _p_ iloz_, int_4 _p_ ihiz_, complex_16 _p_ z_, int_4 _p_ ldz_, int_4 _p_
19333 info_)
19334 { // ** body not listed **
19510 }
19511 real_8 _zlange (char _p_ norm_, int_4 _p_ m_, int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
19512 { // ** body not listed **
19563 }
19564 real_8 _zlanhs (char _p_ norm_, int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
19565 { // ** body not listed **
19617 }
19618 int_4 _zlarf (char _p_ side_, int_4 _p_ m_, int_4 _p_ n_, complex_16 _p_ v_, int_4 _p_ incv_, complex_16 _p_ tau_,
19619 complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_)
19620 { // ** body not listed **
19641 }
19642 int_4 _zlarfg (int_4 _p_ n_, complex_16 _p_ alpha_, complex_16 _p_ x_, int_4 _p_ incx_, complex_16 _p_ tau_)
19643 { // ** body not listed **
19698 }
19699 int_4 _zlarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, complex_16 _p_ x_)
19700 { // ** body not listed **
19736 }
19737 int_4 _zlartg (complex_16 _p_ f_, complex_16 _p_ g_, real_8 _p_ cs_, complex_16 _p_ sn_, complex_16 _p_ r_)
19738 { // ** body not listed **
19784 }
19785 int_4 _zlascl (char _p_ type_, int_4 _p_ kl_, int_4 _p_ ku_, real_8 _p_ cfrom_, real_8 _p_ cto_, int_4 _p_ m_, int_4
19786 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
19787 { // ** body not listed **
19938 }
19939 int_4 _zlaset (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, complex_16 _p_ alpha_, complex_16 _p_ beta_, complex_16 _p_
19940 a_, int_4 _p_ lda_)
19941 { // ** body not listed **
19976 }
19977 int_4 _zlassq (int_4 _p_ n_, complex_16 _p_ x_, int_4 _p_ incx_, real_8 _p_ scale_, real_8 _p_ sumsq_)
19978 { // ** body not listed **
20009 }
20010 int_4 _zlaswp (int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ k1_, int_4 _p_ k2_, int_4 _p_ ipiv_, int_4
20011 _p_ incx_)
20012 { // ** body not listed **
20053 }
20054 int_4 _zlatrs (char _p_ uplo_, char _p_ trans_, char _p_ diag_, char _p_ normin_, int_4 _p_ n_, complex_16 _p_ a_,
20055 int_4 _p_ lda_, complex_16 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_)
20056 { // ** body not listed **
20529 }
20530 int_4 _zrot (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_, complex_16 _p_ cy_, int_4 _p_ incy_, real_8 _p_ c_,
20531 complex_16 _p_ s_)
20532 { // ** body not listed **
20568 }
20569 int_4 _ztrevc (char _p_ side_, char _p_ howmny_, logical_4 _p_ select_, int_4 _p_ n_, complex_16 _p_ t_, int_4 _p_ ldt_
20570 , complex_16 _p_ vl_, int_4 _p_ ldvl_, complex_16 _p_ vr_, int_4 _p_ ldvr_, int_4 _p_ mm_, int_4 _p_ m_, complex_16 _p_
20571 work_, real_8 _p_ rwork_, int_4 _p_ info_)
20572 { // ** body not listed **
20699 }
20700 if (leftv_) {
20701 is_ = 1;
20702 for (ki_ = 1; ki_ <= (*n_); (ki_)++) {
20703 if (somev_) {
20704 if (! (select_[ki_ - 1])) {
20705 goto _l130;
20706 }
20707 }
20708 smin_ = _max (ulp_ * (({complex_16 cdum_1_ = t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))]; _abs ((real_8) (creal (cdum_1_))) +
20709 _abs (cimag (cdum_1_));})), smlnum_);
20710 work_[(*n_) - 1] = CMPLX (1.0, 0.0);
20711 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20712 work_[k_ - 1] = -__dconjg (&t_[ki_ - 1 + ((*ldt_)) * ((k_ - 1))]);
20713 }
20714 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20715 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] - t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))];
20716 if (({complex_16 cdum_1_ = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))]; _abs ((real_8) (creal (cdum_1_))) + _abs (cimag
20717 (cdum_1_));}) < smin_) {
20718 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = smin_;
20719 }
20720 }
20721 if (ki_ < (*n_)) {
20722 (void) _zlatrs (_dc_11, _dc_14, _dc_12, _dc_35, ($6_ = (*n_) - ki_, &$6_), &t_[ki_ + 1 - 1 + ((*ldt_)) * ((ki_ + 1 -
20723 1))], ldt_, &work_[ki_ + 1 - 1], &scale_, rwork_, info_);
20724 work_[ki_ - 1] = scale_;
20725 }
20726 if (! (over_)) {
20727 (void) _zcopy (($7_ = (*n_) - ki_ + 1, &$7_), &work_[ki_ - 1], &_k1, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
20728 ii_ = _izamax (($8_ = (*n_) - ki_ + 1, &$8_), &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1) + ki_ - 1;
20729 remax_ = 1.0 / ({complex_16 cdum_1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((is_ - 1))]; _abs ((real_8) (creal (cdum_1_))) + _abs
20730 (cimag (cdum_1_));});
20731 (void) _zdscal (($9_ = (*n_) - ki_ + 1, &$9_), &remax_, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
20732 for (k_ = 1; k_ <= ki_ - 1; (k_)++) {
20733 vl_[k_ - 1 + ((*ldvl_)) * ((is_ - 1))] = CMPLX (0.0, 0.0);
20734 }
20735 }
20736 else {
20737 if (ki_ < (*n_)) {
20738 (void) _zgemv (_dc_7, n_, ($10_ = (*n_) - ki_, &$10_), ($11_ = CMPLX (1.0, 0.0), &$11_), &vl_[0 + ((*ldvl_)) * ((ki_ +
20739 1 - 1))], ldvl_, &work_[ki_ + 1 - 1], &_k1, ($12_ = (complex_16) (scale_), &$12_), &vl_[0 + ((*ldvl_)) * ((ki_ - 1))],
20740 &_k1);
20741 }
20742 ii_ = _izamax (n_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
20743 remax_ = 1.0 / ({complex_16 cdum_1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((ki_ - 1))]; _abs ((real_8) (creal (cdum_1_))) + _abs
20744 (cimag (cdum_1_));});
20745 (void) _zdscal (n_, &remax_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
20746 }
20747 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20748 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = work_[k_ + (*n_) - 1];
20749 }
20750 is_ = is_ + 1;
20751 _l130:;
20752 }
20753 }
20754 _l0:;
20755 __calls[139].calls++;
20756 return 0;
20757 }
20758 int_4 _ztrexc (char _p_ compq_, int_4 _p_ n_, complex_16 _p_ t_, int_4 _p_ ldt_, complex_16 _p_ q_, int_4 _p_ ldq_,
20759 int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_)
20760 { // ** body not listed **
20828 }
20829 int_4 _ztrsen (char _p_ job_, char _p_ compq_, logical_4 _p_ select_, int_4 _p_ n_, complex_16 _p_ t_, int_4 _p_ ldt_,
20830 complex_16 _p_ q_, int_4 _p_ ldq_, complex_16 _p_ w_, int_4 _p_ m_, real_8 _p_ s_, real_8 _p_ sep_, complex_16 _p_
20831 work_, int_4 _p_ lwork_, int_4 _p_ info_)
20832 { // ** body not listed **
20930 }
20931 int_4 _ztrsyl (char _p_ trana_, char _p_ tranb_, int_4 _p_ isgn_, int_4 _p_ m_, int_4 _p_ n_, complex_16 _p_ a_, int_4
20932 _p_ lda_, complex_16 _p_ b_, int_4 _p_ ldb_, complex_16 _p_ c_, int_4 _p_ ldc_, real_8 _p_ scale_, int_4 _p_ info_)
20933 { // ** body not listed **
21125 }
21126 int_4 _zunm2r (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, complex_16 _p_ a_, int_4 _p_
21127 lda_, complex_16 _p_ tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_, int_4 _p_ info_)
21128 { // ** body not listed **
21216 }
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
