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.17
88 */
89 #if defined (__GNUC__)
90 #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
91 #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
92 #if (__GNUC__ >= 14)
93 #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
94 #pragma GCC diagnostic ignored "-Wimplicit-int"
95 #pragma GCC diagnostic ignored "-Wint-conversion"
96 #pragma GCC diagnostic ignored "-Wreturn-mismatch"
97 #endif
98 #else
99 #error VIF requires GCC
100 #endif
208 static CALLS __calls[__ncalls] = {
209 {"cgbtf2", 0}, // subroutine
210 {"cgbtrf", 0}, // subroutine
211 {"cgbtrs", 0}, // subroutine
212 {"cgeqr2", 0}, // subroutine
213 {"cgttrf", 0}, // subroutine
214 {"cgttrs", 0}, // subroutine
215 {"clacgv", 0}, // subroutine
216 {"clacon", 0}, // subroutine
217 {"clacpy", 0}, // subroutine
218 {"cladiv", 0}, // complex*8 function
219 {"clahqr", 0}, // subroutine
220 {"clange", 0}, // real*4 function
221 {"clanhs", 0}, // real*4 function
222 {"clarf", 0}, // subroutine
223 {"clarfg", 0}, // subroutine
224 {"clarnv", 0}, // subroutine
225 {"clartg", 0}, // subroutine
226 {"clascl", 0}, // subroutine
227 {"claset", 0}, // subroutine
228 {"classq", 0}, // subroutine
229 {"claswp", 0}, // subroutine
230 {"clatrs", 0}, // subroutine
231 {"cmach", 0}, // real*4 function
232 {"crot", 0}, // subroutine
233 {"ctrevc", 0}, // subroutine
234 {"ctrexc", 0}, // subroutine
235 {"ctrsen", 0}, // subroutine
236 {"ctrsyl", 0}, // subroutine
237 {"cunm2r", 0}, // subroutine
238 {"dgbtf2", 0}, // subroutine
239 {"dgbtrf", 0}, // subroutine
240 {"dgbtrs", 0}, // subroutine
241 {"dgttrf", 0}, // subroutine
242 {"dgttrs", 0}, // subroutine
243 {"dlacon", 0}, // subroutine
244 {"dladiv", 0}, // subroutine
245 {"dlae2", 0}, // subroutine
246 {"dlaev2", 0}, // subroutine
247 {"dlaexc", 0}, // subroutine
248 {"dlagtm", 0}, // subroutine
249 {"dlahqr", 0}, // subroutine
250 {"dlaln2", 0}, // subroutine
251 {"dlanhs", 0}, // real*8 function
252 {"dlanv2", 0}, // subroutine
253 {"dlaptm", 0}, // subroutine
254 {"dlapy3", 0}, // real*8 function
255 {"dlaran", 0}, // real*8 function
256 {"dlarfx", 0}, // subroutine
257 {"dlarnd", 0}, // real*8 function
258 {"dlarnv", 0}, // subroutine
259 {"dlaruv", 0}, // subroutine
260 {"dlasy2", 0}, // subroutine
261 {"dpttrf", 0}, // subroutine
262 {"dpttrs", 0}, // subroutine
263 {"dsteqr", 0}, // subroutine
264 {"dtrevc", 0}, // subroutine
265 {"dtrexc", 0}, // subroutine
266 {"dtrsen", 0}, // subroutine
267 {"dtrsyl", 0}, // subroutine
268 {"dzsum1", 0}, // real*8 function
269 {"icmax1", 0}, // integer*4 function
270 {"izmax1", 0}, // integer*4 function
271 {"scsum1", 0}, // real*4 function
272 {"sgbtf2", 0}, // subroutine
273 {"sgbtrf", 0}, // subroutine
274 {"sgbtrs", 0}, // subroutine
275 {"sgeqr2", 0}, // subroutine
276 {"sgttrf", 0}, // subroutine
277 {"sgttrs", 0}, // subroutine
278 {"slabad", 0}, // subroutine
279 {"slacon", 0}, // subroutine
280 {"slacpy", 0}, // subroutine
281 {"sladiv", 0}, // subroutine
282 {"slae2", 0}, // subroutine
283 {"slaev2", 0}, // subroutine
284 {"slaexc", 0}, // subroutine
285 {"slagtm", 0}, // subroutine
286 {"slahqr", 0}, // subroutine
287 {"slaln2", 0}, // subroutine
288 {"slamch", 0}, // real*4 function
289 {"slamc1", 0}, // subroutine
290 {"slamc2", 0}, // subroutine
291 {"slamc3", 0}, // real*4 function
292 {"slamc4", 0}, // subroutine
293 {"slamc5", 0}, // subroutine
294 {"slange", 0}, // real*4 function
295 {"slanhs", 0}, // real*4 function
296 {"slanst", 0}, // real*4 function
297 {"slanv2", 0}, // subroutine
298 {"slaptm", 0}, // subroutine
299 {"slapy2", 0}, // real*4 function
300 {"slapy3", 0}, // real*4 function
301 {"slaran", 0}, // real*4 function
302 {"slarf", 0}, // subroutine
303 {"slarfg", 0}, // subroutine
304 {"slarfx", 0}, // subroutine
305 {"slarnd", 0}, // real*4 function
306 {"slarnv", 0}, // subroutine
307 {"slartg", 0}, // subroutine
308 {"slaruv", 0}, // subroutine
309 {"slascl", 0}, // subroutine
310 {"slaset", 0}, // subroutine
311 {"slasr", 0}, // subroutine
312 {"slasrt", 0}, // subroutine
313 {"slassq", 0}, // subroutine
314 {"slaswp", 0}, // subroutine
315 {"slasy2", 0}, // subroutine
316 {"sorm2r", 0}, // subroutine
317 {"spttrf", 0}, // subroutine
318 {"spttrs", 0}, // subroutine
319 {"ssteqr", 0}, // subroutine
320 {"strevc", 0}, // subroutine
321 {"strexc", 0}, // subroutine
322 {"strsen", 0}, // subroutine
323 {"strsyl", 0}, // subroutine
324 {"xlaenv", 0}, // subroutine
325 {"zgbtf2", 0}, // subroutine
326 {"zgbtrf", 0}, // subroutine
327 {"zgbtrs", 0}, // subroutine
328 {"zgeqr2", 0}, // subroutine
329 {"zgttrf", 0}, // subroutine
330 {"zgttrs", 0}, // subroutine
331 {"zlacgv", 0}, // subroutine
332 {"zlacon", 0}, // subroutine
333 {"zlacpy", 0}, // subroutine
334 {"zladiv", 0}, // complex*16 function
335 {"zlahqr", 0}, // subroutine
336 {"zlange", 0}, // real*8 function
337 {"zlanhs", 0}, // real*8 function
338 {"zlarf", 0}, // subroutine
339 {"zlarfg", 0}, // subroutine
340 {"zlarnv", 0}, // subroutine
341 {"zlartg", 0}, // subroutine
342 {"zlascl", 0}, // subroutine
343 {"zlaset", 0}, // subroutine
344 {"zlassq", 0}, // subroutine
345 {"zlaswp", 0}, // subroutine
346 {"zlatrs", 0}, // subroutine
347 {"zrot", 0}, // subroutine
348 {"ztrevc", 0}, // subroutine
349 {"ztrexc", 0}, // subroutine
350 {"ztrsen", 0}, // subroutine
351 {"ztrsyl", 0}, // subroutine
352 {"zunm2r", 0}, // subroutine
353 {NULL, 0}
354 };
355
356 static struct {
357 int_4 iparms_[100];
358 } claenv_;
360 int_4 _p_ ipiv_, int_4 _p_ info_);
362 int_4 _p_ ipiv_, int_4 _p_ info_);
364 ab_, int_4 _p_ ldab_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
366 _p_ work_, int_4 _p_ info_);
368 int_4 _p_ ipiv_, int_4 _p_ info_);
370 _p_ du_, complex_8 _p_ du2_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
374 , int_4 _p_ ldb_);
377 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_,
378 int_4 _p_ info_);
380 work_);
383 tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_);
388 , int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
390 complex_8 _p_ a_, int_4 _p_ lda_);
393 int_4 _p_ incx_);
395 a_, int_4 _p_ lda_, complex_8 _p_ x_, real_4 _p_ scale_, real_4 _p_ cnorm_, int_4 _p_ info_);
398 c_, complex_8 _p_ s_);
400 _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_,
401 complex_8 _p_ work_, real_4 _p_ rwork_, int_4 _p_ info_);
403 ldq_, int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_);
405 _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
406 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
408 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_
409 info_);
411 int_4 _p_ lda_, complex_8 _p_ tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_, int_4 _p_ info_);
413 int_4 _p_ ipiv_, int_4 _p_ info_);
415 int_4 _p_ ipiv_, int_4 _p_ info_);
417 int_4 _p_ ldab_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
419 int_4 _p_ info_);
421 real_8 _p_ du2_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
423 ;
427 , real_8 _p_ sn1_);
429 ldq_, int_4 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ work_, int_4 _p_ info_);
431 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_);
433 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
434 _p_ ldz_, int_4 _p_ info_);
436 _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_,
437 real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ scale_, real_8 _p_ xnorm_, int_4 _p_ info_);
440 , real_8 _p_ rt2r_, real_8 _p_ rt2i_, real_8 _p_ cs_, real_8 _p_ sn_);
442 int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_);
446 int_4 _p_ ldc_, real_8 _p_ work_);
451 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_,
452 real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ xnorm_, int_4 _p_ info_);
455 int_4 _p_ info_);
457 real_8 _p_ work_, int_4 _p_ info_);
459 _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_
460 work_, int_4 _p_ info_);
462 int_4 _p_ ifst_, int_4 _p_ ilst_, real_8 _p_ work_, int_4 _p_ info_);
464 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_,
465 real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_);
467 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_);
473 int_4 _p_ ipiv_, int_4 _p_ info_);
475 int_4 _p_ ipiv_, int_4 _p_ info_);
477 int_4 _p_ ldab_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
479 int_4 _p_ info_);
481 int_4 _p_ info_);
483 real_4 _p_ du2_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
486 ;
488 int_4 _p_ ldb_);
492 , real_4 _p_ sn1_);
494 ldq_, int_4 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_ work_, int_4 _p_ info_);
496 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_);
498 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
499 _p_ ldz_, int_4 _p_ info_);
501 _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_,
502 real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ scale_, real_4 _p_ xnorm_, int_4 _p_ info_);
506 _p_ rmin_, int_4 _p_ emax_, real_4 _p_ rmax_);
510 _p_ rmax_);
515 , real_4 _p_ rt2r_, real_4 _p_ rt2i_, real_4 _p_ cs_, real_4 _p_ sn_);
517 int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_);
522 real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_);
525 int_4 _p_ ldc_, real_4 _p_ work_);
531 , int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
533 , int_4 _p_ lda_);
535 real_4 _p_ s_, real_4 _p_ a_, int_4 _p_ lda_);
539 int_4 _p_ incx_);
541 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_,
542 real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ xnorm_, int_4 _p_ info_);
544 int_4 _p_ lda_, real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_, int_4 _p_ info_);
547 int_4 _p_ info_);
549 real_4 _p_ work_, int_4 _p_ info_);
551 _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_
552 work_, int_4 _p_ info_);
554 int_4 _p_ ifst_, int_4 _p_ ilst_, real_4 _p_ work_, int_4 _p_ info_);
556 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_,
557 real_4 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_);
559 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_);
562 int_4 _p_ ipiv_, int_4 _p_ info_);
564 int_4 _p_ ipiv_, int_4 _p_ info_);
566 ab_, int_4 _p_ ldab_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
568 _p_ work_, int_4 _p_ info_);
570 int_4 _p_ ipiv_, int_4 _p_ info_);
572 complex_16 _p_ du_, complex_16 _p_ du2_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
576 b_, int_4 _p_ ldb_);
579 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_
580 ldz_, int_4 _p_ info_);
582 work_);
585 tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_);
590 , int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
592 complex_16 _p_ a_, int_4 _p_ lda_);
595 , int_4 _p_ incx_);
597 _p_ a_, int_4 _p_ lda_, complex_16 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_);
599 _p_ c_, complex_16 _p_ s_);
601 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_,
602 complex_16 _p_ work_, real_8 _p_ rwork_, int_4 _p_ info_);
604 ldq_, int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_);
606 _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_,
607 complex_16 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
609 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_
610 info_);
612 int_4 _p_ lda_, complex_16 _p_ tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_, int_4 _p_ info_);
613 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_
614 ipiv_, int_4 _p_ info_)
615 { // ** body not listed **
688 }
689
690 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_
691 ipiv_, int_4 _p_ info_)
692 { // ** body not listed **
889 }
890
891 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
892 _p_ ldab_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
893 { // ** body not listed **
994 }
995
996 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_,
997 int_4 _p_ info_)
998 { // ** body not listed **
1033 }
1034
1035 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_
1036 ipiv_, int_4 _p_ info_)
1037 { // ** body not listed **
1085 }
1086
1087 if (d_[(*n_) - 1] == CMPLXF (0.0, 0.0) && (*info_) == 0) {
1088 (*info_) = (*n_);
1089 goto _l0;
1090 }
1091
1092 _l0:;
1093 __calls[4].calls++;
1094 return 0;
1095 }
1096
1097 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_,
1098 complex_8 _p_ du2_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
1099 { // ** body not listed **
1201 }
1202
1203 int_4 _clacgv (int_4 _p_ n_, complex_8 _p_ x_, int_4 _p_ incx_)
1204 { // ** body not listed **
1223 }
1224
1225 int_4 _clacon (int_4 _p_ n_, complex_8 _p_ v_, complex_8 _p_ x_, real_4 _p_ est_, int_4 _p_ kase_)
1226 { // ** body not listed **
1321 }
1322 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
1323 _p_ ldb_)
1324 { // ** body not listed **
1350 }
1351
1352 complex_8 _cladiv (complex_8 _p_ x_, complex_8 _p_ y_)
1353 { // ** body not listed **
1362 }
1363
1364 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_
1365 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_
1366 info_)
1367 { // ** body not listed **
1539 }
1540
1541 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_)
1542 { // ** body not listed **
1593 }
1594
1595 real_4 _clanhs (char _p_ norm_, int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
1596 { // ** body not listed **
1648 }
1649
1650 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_,
1651 complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_)
1652 { // ** body not listed **
1673 }
1674
1675 int_4 _clarfg (int_4 _p_ n_, complex_8 _p_ alpha_, complex_8 _p_ x_, int_4 _p_ incx_, complex_8 _p_ tau_)
1676 { // ** body not listed **
1730 }
1731
1732 int_4 _clarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, complex_8 _p_ x_)
1733 { // ** body not listed **
1768 }
1769
1770 int_4 _clartg (complex_8 _p_ f_, complex_8 _p_ g_, real_4 _p_ cs_, complex_8 _p_ sn_, complex_8 _p_ r_)
1771 { // ** body not listed **
1815 }
1816
1817 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
1818 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
1819 { // ** body not listed **
1969 }
1970
1971 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_,
1972 int_4 _p_ lda_)
1973 { // ** body not listed **
2008 }
2009
2010 int_4 _classq (int_4 _p_ n_, complex_8 _p_ x_, int_4 _p_ incx_, real_4 _p_ scale_, real_4 _p_ sumsq_)
2011 { // ** body not listed **
2042 }
2043
2044 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_
2045 incx_)
2046 { // ** body not listed **
2087 }
2088
2089 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
2090 _p_ lda_, complex_8 _p_ x_, real_4 _p_ scale_, real_4 _p_ cnorm_, int_4 _p_ info_)
2091 { // ** body not listed **
2560 }
2561
2562 real_4 _cmach (int_4 _p_ job_)
2563 { // ** body not listed **
2603 }
2604
2605 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_,
2606 complex_8 _p_ s_)
2607 { // ** body not listed **
2641 }
2642
2643 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_,
2644 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_
2645 work_, real_4 _p_ rwork_, int_4 _p_ info_)
2646 { // ** body not listed **
2769 }
2770
2771 if (leftv_) {
2772 is_ = 1;
2773 for (ki_ = 1; ki_ <= (*n_); (ki_)++) {
2774 if (somev_) {
2775 if (! (select_[ki_ - 1])) {
2776 goto _l130;
2777 }
2778 }
2779 smin_ = _max (ulp_ * (({complex_8 cdum__1_ = t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))]; _abs ((real_4) (crealf (cdum__1_)))
2780 + _abs (cimagf (cdum__1_));})), smlnum_);
2781 work_[(*n_) - 1] = CMPLXF (1.0, 0.0);
2782 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2783 work_[k_ - 1] = -conjf (t_[ki_ - 1 + ((*ldt_)) * ((k_ - 1))]);
2784 }
2785 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2786 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] - t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))];
2787 if (({complex_8 cdum__1_ = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))]; _abs ((real_4) (crealf (cdum__1_))) + _abs (cimagf
2788 (cdum__1_));}) < smin_) {
2789 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = smin_;
2790 }
2791 }
2792 if (ki_ < (*n_)) {
2793 (void) _clatrs (_dc_11, _dc_14, _dc_12, _dc_35, ($6_ = (*n_) - ki_, &$6_), &t_[ki_ + 1 - 1 + ((*ldt_)) * ((ki_ + 1 -
2794 1))], ldt_, &work_[ki_ + 1 - 1], &scale_, rwork_, info_);
2795 work_[ki_ - 1] = scale_;
2796 }
2797 if (! (over_)) {
2798 (void) _ccopy (($7_ = (*n_) - ki_ + 1, &$7_), &work_[ki_ - 1], &_k1, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
2799 ii_ = _icamax (($8_ = (*n_) - ki_ + 1, &$8_), &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1) + ki_ - 1;
2800 remax_ = 1.0 / ({complex_8 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((is_ - 1))]; _abs ((real_4) (crealf (cdum__1_))) +
2801 _abs (cimagf (cdum__1_));});
2802 (void) _csscal (($9_ = (*n_) - ki_ + 1, &$9_), &remax_, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
2803 for (k_ = 1; k_ <= ki_ - 1; (k_)++) {
2804 vl_[k_ - 1 + ((*ldvl_)) * ((is_ - 1))] = CMPLXF (0.0, 0.0);
2805 }
2806 }
2807 else {
2808 if (ki_ < (*n_)) {
2809 (void) _cgemv (_dc_7, n_, ($10_ = (*n_) - ki_, &$10_), ($11_ = CMPLXF (1.0, 0.0), &$11_), &vl_[0 + ((*ldvl_)) * ((ki_ +
2810 1 - 1))], ldvl_, &work_[ki_ + 1 - 1], &_k1, ($12_ = (complex_8) (scale_), &$12_), &vl_[0 + ((*ldvl_)) * ((ki_ - 1))],
2811 &_k1);
2812 }
2813 ii_ = _icamax (n_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
2814 remax_ = 1.0 / ({complex_8 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((ki_ - 1))]; _abs ((real_4) (crealf (cdum__1_))) +
2815 _abs (cimagf (cdum__1_));});
2816 (void) _csscal (n_, &remax_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
2817 }
2818 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2819 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = work_[k_ + (*n_) - 1];
2820 }
2821 is_ = is_ + 1;
2822 _l130:;
2823 }
2824
2825 }
2826
2827 _l0:;
2828 __calls[24].calls++;
2829 return 0;
2830 }
2831
2832 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
2833 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_)
2834 { // ** body not listed **
2900 }
2901
2902 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_,
2903 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_,
2904 int_4 _p_ lwork_, int_4 _p_ info_)
2905 { // ** body not listed **
3003 }
3004
3005 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
3006 _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_)
3007 { // ** body not listed **
3194 }
3195
3196 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_
3197 lda_, complex_8 _p_ tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_, int_4 _p_ info_)
3198 { // ** body not listed **
3284 }
3285
3286 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_
3287 ipiv_, int_4 _p_ info_)
3288 { // ** body not listed **
3361 }
3362
3363 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_
3364 ipiv_, int_4 _p_ info_)
3365 { // ** body not listed **
3561 }
3562
3563 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_
3564 ldab_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
3565 { // ** body not listed **
3645 }
3646
3647 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_
3648 info_)
3649 { // ** body not listed **
3704 }
3705
3706 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
3707 _p_ du2_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
3708 { // ** body not listed **
3786 }
3787
3788 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_)
3789 { // ** body not listed **
3884 }
3885
3886 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_)
3887 { // ** body not listed **
3903 }
3904
3905 int_4 _dlae2 (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ rt1_, real_8 _p_ rt2_)
3906 { // ** body not listed **
3945 }
3946
3947 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
3948 _p_ sn1_)
3949 { // ** body not listed **
4022 }
4023
4024 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
4025 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ work_, int_4 _p_ info_)
4026 { // ** body not listed **
4196 }
4197
4198 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
4199 _p_ du_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_)
4200 { // ** body not listed **
4299 }
4300
4301 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_,
4302 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
4303 _p_ info_)
4304 { // ** body not listed **
4499 }
4500
4501 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_,
4502 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
4503 _p_ x_, int_4 _p_ ldx_, real_8 _p_ scale_, real_8 _p_ xnorm_, int_4 _p_ info_)
4504 { // ** body not listed **
4829 }
4830
4831 real_8 _dlanhs (char _p_ norm_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
4832 { // ** body not listed **
4885 }
4886
4887 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
4888 _p_ rt2r_, real_8 _p_ rt2i_, real_8 _p_ cs_, real_8 _p_ sn_)
4889 { // ** body not listed **
4971 }
4972
4973 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_
4974 ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_)
4975 { // ** body not listed **
5033 }
5034
5035 real_8 _dlapy3 (real_8 _p_ x_, real_8 _p_ y_, real_8 _p_ z_)
5036 { // ** body not listed **
5053 }
5054
5055 real_8 _dlaran (int_4 _p_ iseed_)
5056 { // ** body not listed **
5078 }
5079
5080 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_
5081 ldc_, real_8 _p_ work_)
5082 { // ** body not listed **
5589 }
5590
5591 real_8 _dlarnd (int_4 _p_ idist_, int_4 _p_ iseed_)
5592 { // ** body not listed **
5608 }
5609
5610 int_4 _dlarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, real_8 _p_ x_)
5611 { // ** body not listed **
5641 }
5642
5643 int_4 _dlaruv (int_4 _p_ iseed_, int_4 _p_ n_, real_8 _p_ x_)
5644 { // ** body not listed **
7733 }
7734
7735 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_
7736 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_,
7737 int_4 _p_ ldx_, real_8 _p_ xnorm_, int_4 _p_ info_)
7738 { // ** body not listed **
8048 }
8049
8050 int_4 _dpttrf (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ info_)
8051 { // ** body not listed **
8084 }
8085
8086 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_
8087 info_)
8088 { // ** body not listed **
8123 }
8124
8125 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_
8126 work_, int_4 _p_ info_)
8127 { // ** body not listed **
8456 }
8457
8458 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_,
8459 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
8460 _p_ info_)
8461 { // ** body not listed **
9074 }
9075
9076 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_
9077 ifst_, int_4 _p_ ilst_, real_8 _p_ work_, int_4 _p_ info_)
9078 { // ** body not listed **
9268 }
9269
9270 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_,
9271 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_
9272 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
9273 { // ** body not listed **
9435 }
9436
9437 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_
9438 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_)
9439 { // ** body not listed **
10177 }
10178
10179 real_8 _dzsum1 (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_)
10180 { // ** body not listed **
10206 }
10207
10208 int_4 _icmax1 (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_)
10209 { // ** body not listed **
10236 }
10237
10238 goto _l0;
10239 _l30:;
10240 smax_ = ({complex_8 zdum__1_ = cx_[0]; _abs ((real_4) (crealf (zdum__1_)));});
10241 for (i_ = 2; i_ <= (*n_); (i_)++) {
10242 if (({complex_8 zdum__1_ = cx_[i_ - 1]; _abs ((real_4) (crealf (zdum__1_)));}) <= smax_) {
10243 goto _l40;
10244 }
10245
10246 icmax1_ = i_;
10247 smax_ = ({complex_8 zdum__1_ = cx_[i_ - 1]; _abs ((real_4) (crealf (zdum__1_)));});
10248 _l40:;
10249 }
10250
10251 _l0:;
10252 __calls[60].calls++;
10253 return icmax1_;
10254 }
10255
10256 int_4 _izmax1 (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_)
10257 { // ** body not listed **
10284 }
10285
10286 goto _l0;
10287 _l30:;
10288 smax_ = ({complex_16 zdum__1_ = cx_[0]; _abs ((real_8) (creal (zdum__1_)));});
10289 for (i_ = 2; i_ <= (*n_); (i_)++) {
10290 if (({complex_16 zdum__1_ = cx_[i_ - 1]; _abs ((real_8) (creal (zdum__1_)));}) <= smax_) {
10291 goto _l40;
10292 }
10293
10294 izmax1_ = i_;
10295 smax_ = ({complex_16 zdum__1_ = cx_[i_ - 1]; _abs ((real_8) (creal (zdum__1_)));});
10296 _l40:;
10297 }
10298
10299 _l0:;
10300 __calls[61].calls++;
10301 return izmax1_;
10302 }
10303
10304 real_4 _scsum1 (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_)
10305 { // ** body not listed **
10331 }
10332
10333 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_
10334 ipiv_, int_4 _p_ info_)
10335 { // ** body not listed **
10408 }
10409
10410 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_
10411 ipiv_, int_4 _p_ info_)
10412 { // ** body not listed **
10607 }
10608 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_
10609 ldab_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
10610 { // ** body not listed **
10690 }
10691
10692 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_
10693 info_)
10694 { // ** body not listed **
10728 }
10729
10730 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_
10731 info_)
10732 { // ** body not listed **
10787 }
10788
10789 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
10790 _p_ du2_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
10791 { // ** body not listed **
10869 }
10870
10871 int_4 _slabad (real_4 _p_ small_, real_4 _p_ large_)
10872 { // ** body not listed **
10879 }
10880
10881 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_)
10882 { // ** body not listed **
10977 }
10978
10979 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_
10980 )
10981 { // ** body not listed **
11007 }
11008
11009 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_)
11010 { // ** body not listed **
11026 }
11027
11028 int_4 _slae2 (real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ rt1_, real_4 _p_ rt2_)
11029 { // ** body not listed **
11068 }
11069
11070 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
11071 _p_ sn1_)
11072 { // ** body not listed **
11145 }
11146
11147 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
11148 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_ work_, int_4 _p_ info_)
11149 { // ** body not listed **
11314 }
11315
11316 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
11317 _p_ du_, real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_)
11318 { // ** body not listed **
11417 }
11418
11419 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_,
11420 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
11421 _p_ info_)
11422 { // ** body not listed **
11614 }
11615
11616 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_,
11617 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
11618 _p_ x_, int_4 _p_ ldx_, real_4 _p_ scale_, real_4 _p_ xnorm_, int_4 _p_ info_)
11619 { // ** body not listed **
11943 }
11944
11945 real_4 _slamch (char _p_ cmach_)
11946 { // ** body not listed **
12019 }
12020
12021 int_4 _slamc1 (int_4 _p_ beta_, int_4 _p_ t_, logical_4 _p_ rnd_, logical_4 _p_ ieee1_)
12022 { // ** body not listed **
12098 }
12099
12100 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_,
12101 int_4 _p_ emax_, real_4 _p_ rmax_)
12102 { // ** body not listed **
12275 }
12276
12277 real_4 _slamc3 (real_4 _p_ a_, real_4 _p_ b_)
12278 { // ** body not listed **
12283 }
12284
12285 int_4 _slamc4 (int_4 _p_ emin_, real_4 _p_ start_, int_4 _p_ base_)
12286 { // ** body not listed **
12320 }
12321
12322 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_)
12323 { // ** body not listed **
12376 }
12377
12378 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_)
12379 { // ** body not listed **
12430 }
12431
12432 real_4 _slanhs (char _p_ norm_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
12433 { // ** body not listed **
12485 }
12486
12487 real_4 _slanst (char _p_ norm_, int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_)
12488 { // ** body not listed **
12528 }
12529
12530 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
12531 _p_ rt2r_, real_4 _p_ rt2i_, real_4 _p_ cs_, real_4 _p_ sn_)
12532 { // ** body not listed **
12613 }
12614
12615 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_
12616 ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_)
12617 { // ** body not listed **
12675 }
12676
12677 real_4 _slapy2 (real_4 _p_ x_, real_4 _p_ y_)
12678 { // ** body not listed **
12694 }
12695
12696 real_4 _slapy3 (real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ z_)
12697 { // ** body not listed **
12714 }
12715
12716 real_4 _slaran (int_4 _p_ iseed_)
12717 { // ** body not listed **
12739 }
12740
12741 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_
12742 c_, int_4 _p_ ldc_, real_4 _p_ work_)
12743 { // ** body not listed **
12762 }
12763
12764 int_4 _slarfg (int_4 _p_ n_, real_4 _p_ alpha_, real_4 _p_ x_, int_4 _p_ incx_, real_4 _p_ tau_)
12765 { // ** body not listed **
12812 }
12813
12814 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_
12815 ldc_, real_4 _p_ work_)
12816 { // ** body not listed **
13323 }
13324
13325 real_4 _slarnd (int_4 _p_ idist_, int_4 _p_ iseed_)
13326 { // ** body not listed **
13342 }
13343
13344 int_4 _slarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, real_4 _p_ x_)
13345 { // ** body not listed **
13375 }
13376
13377 int_4 _slartg (real_4 _p_ f_, real_4 _p_ g_, real_4 _p_ cs_, real_4 _p_ sn_, real_4 _p_ r_)
13378 { // ** body not listed **
13462 }
13463
13464 int_4 _slaruv (int_4 _p_ iseed_, int_4 _p_ n_, real_4 _p_ x_)
13465 { // ** body not listed **
15555 }
15556
15557 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
15558 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
15559 { // ** body not listed **
15709 }
15710
15711 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
15712 _p_ lda_)
15713 { // ** body not listed **
15742 }
15743
15744 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_
15745 s_, real_4 _p_ a_, int_4 _p_ lda_)
15746 { // ** body not listed **
15952 }
15953
15954 int_4 _slasrt (char _p_ id_, int_4 _p_ n_, real_4 _p_ d_, int_4 _p_ info_)
15955 { // ** body not listed **
16129 }
16130
16131 int_4 _slassq (int_4 _p_ n_, real_4 _p_ x_, int_4 _p_ incx_, real_4 _p_ scale_, real_4 _p_ sumsq_)
16132 { // ** body not listed **
16153 }
16154
16155 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_
16156 incx_)
16157 { // ** body not listed **
16198 }
16199
16200 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_
16201 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_,
16202 int_4 _p_ ldx_, real_4 _p_ xnorm_, int_4 _p_ info_)
16203 { // ** body not listed **
16512 }
16513
16514 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_
16515 , real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_, int_4 _p_ info_)
16516 { // ** body not listed **
16596 }
16597
16598 int_4 _spttrf (int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_, int_4 _p_ info_)
16599 { // ** body not listed **
16632 }
16633
16634 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_
16635 info_)
16636 { // ** body not listed **
16671 }
16672
16673 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_
16674 work_, int_4 _p_ info_)
16675 { // ** body not listed **
16996 }
16997
16998 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_,
16999 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
17000 _p_ info_)
17001 { // ** body not listed **
17612 }
17613
17614 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_
17615 ifst_, int_4 _p_ ilst_, real_4 _p_ work_, int_4 _p_ info_)
17616 { // ** body not listed **
17806 }
17807
17808 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_,
17809 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_
17810 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
17811 { // ** body not listed **
17971 }
17972
17973 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_
17974 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_)
17975 { // ** body not listed **
18710 }
18711
18712 int_4 _xlaenv (int_4 _p_ ispec_, int_4 _p_ nvalue_)
18713 { // ** body not listed **
18719 }
18720
18721 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_
18722 ipiv_, int_4 _p_ info_)
18723 { // ** body not listed **
18796 }
18797
18798 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_
18799 ipiv_, int_4 _p_ info_)
18800 { // ** body not listed **
18997 }
18998
18999 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
19000 _p_ ldab_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
19001 { // ** body not listed **
19102 }
19103
19104 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_
19105 , int_4 _p_ info_)
19106 { // ** body not listed **
19143 }
19144
19145 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_
19146 ipiv_, int_4 _p_ info_)
19147 { // ** body not listed **
19195 }
19196
19197 if (d_[(*n_) - 1] == CMPLX (0.0, 0.0) && (*info_) == 0) {
19198 (*info_) = (*n_);
19199 goto _l0;
19200 }
19201
19202 _l0:;
19203 __calls[120].calls++;
19204 return 0;
19205 }
19206
19207 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_
19208 du_, complex_16 _p_ du2_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
19209 { // ** body not listed **
19314 }
19315
19316 int_4 _zlacgv (int_4 _p_ n_, complex_16 _p_ x_, int_4 _p_ incx_)
19317 { // ** body not listed **
19338 }
19339
19340 int_4 _zlacon (int_4 _p_ n_, complex_16 _p_ v_, complex_16 _p_ x_, real_8 _p_ est_, int_4 _p_ kase_)
19341 { // ** body not listed **
19437 }
19438
19439 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
19440 _p_ ldb_)
19441 { // ** body not listed **
19467 }
19468
19469 complex_16 _zladiv (complex_16 _p_ x_, complex_16 _p_ y_)
19470 { // ** body not listed **
19479 }
19480
19481 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_
19482 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_
19483 info_)
19484 { // ** body not listed **
19660 }
19661
19662 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_)
19663 { // ** body not listed **
19714 }
19715
19716 real_8 _zlanhs (char _p_ norm_, int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
19717 { // ** body not listed **
19769 }
19770
19771 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_,
19772 complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_)
19773 { // ** body not listed **
19794 }
19795
19796 int_4 _zlarfg (int_4 _p_ n_, complex_16 _p_ alpha_, complex_16 _p_ x_, int_4 _p_ incx_, complex_16 _p_ tau_)
19797 { // ** body not listed **
19852 }
19853
19854 int_4 _zlarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, complex_16 _p_ x_)
19855 { // ** body not listed **
19891 }
19892
19893 int_4 _zlartg (complex_16 _p_ f_, complex_16 _p_ g_, real_8 _p_ cs_, complex_16 _p_ sn_, complex_16 _p_ r_)
19894 { // ** body not listed **
19940 }
19941
19942 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
19943 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
19944 { // ** body not listed **
20095 }
20096
20097 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_
20098 a_, int_4 _p_ lda_)
20099 { // ** body not listed **
20134 }
20135
20136 int_4 _zlassq (int_4 _p_ n_, complex_16 _p_ x_, int_4 _p_ incx_, real_8 _p_ scale_, real_8 _p_ sumsq_)
20137 { // ** body not listed **
20168 }
20169
20170 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
20171 _p_ incx_)
20172 { // ** body not listed **
20213 }
20214
20215 int_4 _zlatrs (char _p_ uplo_, char _p_ trans_, char _p_ diag_, char _p_ normin_, int_4 _p_ n_, complex_16 _p_ a_,
20216 int_4 _p_ lda_, complex_16 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_)
20217 { // ** body not listed **
20690 }
20691
20692 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_,
20693 complex_16 _p_ s_)
20694 { // ** body not listed **
20730 }
20731
20732 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_
20733 , 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_
20734 work_, real_8 _p_ rwork_, int_4 _p_ info_)
20735 { // ** body not listed **
20862 }
20863
20864 if (leftv_) {
20865 is_ = 1;
20866 for (ki_ = 1; ki_ <= (*n_); (ki_)++) {
20867 if (somev_) {
20868 if (! (select_[ki_ - 1])) {
20869 goto _l130;
20870 }
20871 }
20872 smin_ = _max (ulp_ * (({complex_16 cdum__1_ = t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))]; _abs ((real_8) (creal (cdum__1_)))
20873 + _abs (cimag (cdum__1_));})), smlnum_);
20874 work_[(*n_) - 1] = CMPLX (1.0, 0.0);
20875 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20876 work_[k_ - 1] = -__dconjg (&t_[ki_ - 1 + ((*ldt_)) * ((k_ - 1))]);
20877 }
20878 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20879 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] - t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))];
20880 if (({complex_16 cdum__1_ = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))]; _abs ((real_8) (creal (cdum__1_))) + _abs (cimag
20881 (cdum__1_));}) < smin_) {
20882 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = smin_;
20883 }
20884 }
20885 if (ki_ < (*n_)) {
20886 (void) _zlatrs (_dc_11, _dc_14, _dc_12, _dc_35, ($6_ = (*n_) - ki_, &$6_), &t_[ki_ + 1 - 1 + ((*ldt_)) * ((ki_ + 1 -
20887 1))], ldt_, &work_[ki_ + 1 - 1], &scale_, rwork_, info_);
20888 work_[ki_ - 1] = scale_;
20889 }
20890 if (! (over_)) {
20891 (void) _zcopy (($7_ = (*n_) - ki_ + 1, &$7_), &work_[ki_ - 1], &_k1, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
20892 ii_ = _izamax (($8_ = (*n_) - ki_ + 1, &$8_), &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1) + ki_ - 1;
20893 remax_ = 1.0 / ({complex_16 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((is_ - 1))]; _abs ((real_8) (creal (cdum__1_))) +
20894 _abs (cimag (cdum__1_));});
20895 (void) _zdscal (($9_ = (*n_) - ki_ + 1, &$9_), &remax_, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
20896 for (k_ = 1; k_ <= ki_ - 1; (k_)++) {
20897 vl_[k_ - 1 + ((*ldvl_)) * ((is_ - 1))] = CMPLX (0.0, 0.0);
20898 }
20899 }
20900 else {
20901 if (ki_ < (*n_)) {
20902 (void) _zgemv (_dc_7, n_, ($10_ = (*n_) - ki_, &$10_), ($11_ = CMPLX (1.0, 0.0), &$11_), &vl_[0 + ((*ldvl_)) * ((ki_ +
20903 1 - 1))], ldvl_, &work_[ki_ + 1 - 1], &_k1, ($12_ = (complex_16) (scale_), &$12_), &vl_[0 + ((*ldvl_)) * ((ki_ - 1))],
20904 &_k1);
20905 }
20906 ii_ = _izamax (n_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
20907 remax_ = 1.0 / ({complex_16 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((ki_ - 1))]; _abs ((real_8) (creal (cdum__1_))) +
20908 _abs (cimag (cdum__1_));});
20909 (void) _zdscal (n_, &remax_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
20910 }
20911 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20912 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = work_[k_ + (*n_) - 1];
20913 }
20914 is_ = is_ + 1;
20915 _l130:;
20916 }
20917
20918 }
20919
20920 _l0:;
20921 __calls[139].calls++;
20922 return 0;
20923 }
20924
20925 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_,
20926 int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_)
20927 { // ** body not listed **
20995 }
20996
20997 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_,
20998 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_
20999 work_, int_4 _p_ lwork_, int_4 _p_ info_)
21000 { // ** body not listed **
21098 }
21099
21100 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
21101 _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_)
21102 { // ** body not listed **
21294 }
21295
21296 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_
21297 lda_, complex_16 _p_ tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_, int_4 _p_ info_)
21298 { // ** body not listed **
21386 }
21387
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|