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.12
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
344 static struct {
345 int_4 iparms_[100];
346 } claenv_;
348 int_4 _p_ ipiv_, int_4 _p_ info_);
350 int_4 _p_ ipiv_, int_4 _p_ info_);
352 ab_, int_4 _p_ ldab_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
354 _p_ work_, int_4 _p_ info_);
356 int_4 _p_ ipiv_, int_4 _p_ info_);
358 _p_ du_, complex_8 _p_ du2_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
362 , int_4 _p_ ldb_);
365 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_,
366 int_4 _p_ info_);
368 work_);
371 tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_);
376 , int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
378 complex_8 _p_ a_, int_4 _p_ lda_);
381 int_4 _p_ incx_);
383 a_, int_4 _p_ lda_, complex_8 _p_ x_, real_4 _p_ scale_, real_4 _p_ cnorm_, int_4 _p_ info_);
386 c_, complex_8 _p_ s_);
388 _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_,
389 complex_8 _p_ work_, real_4 _p_ rwork_, int_4 _p_ info_);
391 ldq_, int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_);
393 _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
394 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
396 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_
397 info_);
399 int_4 _p_ lda_, complex_8 _p_ tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_, int_4 _p_ info_);
401 int_4 _p_ ipiv_, int_4 _p_ info_);
403 int_4 _p_ ipiv_, int_4 _p_ info_);
405 int_4 _p_ ldab_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
407 int_4 _p_ info_);
409 real_8 _p_ du2_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
411 ;
415 , real_8 _p_ sn1_);
417 ldq_, int_4 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ work_, int_4 _p_ info_);
419 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_);
421 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
422 _p_ ldz_, int_4 _p_ info_);
424 _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_,
425 real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ scale_, real_8 _p_ xnorm_, int_4 _p_ info_);
428 , real_8 _p_ rt2r_, real_8 _p_ rt2i_, real_8 _p_ cs_, real_8 _p_ sn_);
430 int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_);
434 int_4 _p_ ldc_, real_8 _p_ work_);
439 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_,
440 real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ xnorm_, int_4 _p_ info_);
443 int_4 _p_ info_);
445 real_8 _p_ work_, int_4 _p_ info_);
447 _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_
448 work_, int_4 _p_ info_);
450 int_4 _p_ ifst_, int_4 _p_ ilst_, real_8 _p_ work_, int_4 _p_ info_);
452 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_,
453 real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_);
455 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_);
461 int_4 _p_ ipiv_, int_4 _p_ info_);
463 int_4 _p_ ipiv_, int_4 _p_ info_);
465 int_4 _p_ ldab_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
467 int_4 _p_ info_);
469 int_4 _p_ info_);
471 real_4 _p_ du2_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
474 ;
476 int_4 _p_ ldb_);
480 , real_4 _p_ sn1_);
482 ldq_, int_4 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_ work_, int_4 _p_ info_);
484 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_);
486 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
487 _p_ ldz_, int_4 _p_ info_);
489 _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_,
490 real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ scale_, real_4 _p_ xnorm_, int_4 _p_ info_);
494 _p_ rmin_, int_4 _p_ emax_, real_4 _p_ rmax_);
498 _p_ rmax_);
503 , real_4 _p_ rt2r_, real_4 _p_ rt2i_, real_4 _p_ cs_, real_4 _p_ sn_);
505 int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_);
510 real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_);
513 int_4 _p_ ldc_, real_4 _p_ work_);
519 , int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
521 , int_4 _p_ lda_);
523 real_4 _p_ s_, real_4 _p_ a_, int_4 _p_ lda_);
527 int_4 _p_ incx_);
529 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_,
530 real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ xnorm_, int_4 _p_ info_);
532 int_4 _p_ lda_, real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_, int_4 _p_ info_);
535 int_4 _p_ info_);
537 real_4 _p_ work_, int_4 _p_ info_);
539 _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_
540 work_, int_4 _p_ info_);
542 int_4 _p_ ifst_, int_4 _p_ ilst_, real_4 _p_ work_, int_4 _p_ info_);
544 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_,
545 real_4 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_);
547 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_);
550 int_4 _p_ ipiv_, int_4 _p_ info_);
552 int_4 _p_ ipiv_, int_4 _p_ info_);
554 ab_, int_4 _p_ ldab_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
556 _p_ work_, int_4 _p_ info_);
558 int_4 _p_ ipiv_, int_4 _p_ info_);
560 complex_16 _p_ du_, complex_16 _p_ du2_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_);
564 b_, int_4 _p_ ldb_);
567 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_
568 ldz_, int_4 _p_ info_);
570 work_);
573 tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_);
578 , int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
580 complex_16 _p_ a_, int_4 _p_ lda_);
583 , int_4 _p_ incx_);
585 _p_ a_, int_4 _p_ lda_, complex_16 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_);
587 _p_ c_, complex_16 _p_ s_);
589 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_,
590 complex_16 _p_ work_, real_8 _p_ rwork_, int_4 _p_ info_);
592 ldq_, int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_);
594 _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_,
595 complex_16 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
597 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_
598 info_);
600 int_4 _p_ lda_, complex_16 _p_ tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_, int_4 _p_ info_);
601 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_
602 ipiv_, int_4 _p_ info_)
603 { // ** body not listed **
676 }
677
678 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_
679 ipiv_, int_4 _p_ info_)
680 { // ** body not listed **
877 }
878
879 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
880 _p_ ldab_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
881 { // ** body not listed **
982 }
983
984 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_,
985 int_4 _p_ info_)
986 { // ** body not listed **
1021 }
1022
1023 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_
1024 ipiv_, int_4 _p_ info_)
1025 { // ** body not listed **
1073 }
1074
1075 if (d_[(*n_) - 1] == CMPLXF (0.0, 0.0) && (*info_) == 0) {
1076 (*info_) = (*n_);
1077 goto _l0;
1078 }
1079
1080 _l0:;
1081 __calls[4].calls++;
1082 return 0;
1083 }
1084
1085 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_,
1086 complex_8 _p_ du2_, int_4 _p_ ipiv_, complex_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
1087 { // ** body not listed **
1189 }
1190
1191 int_4 _clacgv (int_4 _p_ n_, complex_8 _p_ x_, int_4 _p_ incx_)
1192 { // ** body not listed **
1211 }
1212
1213 int_4 _clacon (int_4 _p_ n_, complex_8 _p_ v_, complex_8 _p_ x_, real_4 _p_ est_, int_4 _p_ kase_)
1214 { // ** body not listed **
1309 }
1310 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
1311 _p_ ldb_)
1312 { // ** body not listed **
1338 }
1339
1340 complex_8 _cladiv (complex_8 _p_ x_, complex_8 _p_ y_)
1341 { // ** body not listed **
1350 }
1351
1352 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_
1353 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_
1354 info_)
1355 { // ** body not listed **
1527 }
1528
1529 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_)
1530 { // ** body not listed **
1581 }
1582
1583 real_4 _clanhs (char _p_ norm_, int_4 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
1584 { // ** body not listed **
1636 }
1637
1638 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_,
1639 complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_)
1640 { // ** body not listed **
1661 }
1662
1663 int_4 _clarfg (int_4 _p_ n_, complex_8 _p_ alpha_, complex_8 _p_ x_, int_4 _p_ incx_, complex_8 _p_ tau_)
1664 { // ** body not listed **
1718 }
1719
1720 int_4 _clarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, complex_8 _p_ x_)
1721 { // ** body not listed **
1756 }
1757
1758 int_4 _clartg (complex_8 _p_ f_, complex_8 _p_ g_, real_4 _p_ cs_, complex_8 _p_ sn_, complex_8 _p_ r_)
1759 { // ** body not listed **
1803 }
1804
1805 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
1806 _p_ n_, complex_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
1807 { // ** body not listed **
1957 }
1958
1959 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_,
1960 int_4 _p_ lda_)
1961 { // ** body not listed **
1996 }
1997
1998 int_4 _classq (int_4 _p_ n_, complex_8 _p_ x_, int_4 _p_ incx_, real_4 _p_ scale_, real_4 _p_ sumsq_)
1999 { // ** body not listed **
2030 }
2031
2032 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_
2033 incx_)
2034 { // ** body not listed **
2075 }
2076
2077 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
2078 _p_ lda_, complex_8 _p_ x_, real_4 _p_ scale_, real_4 _p_ cnorm_, int_4 _p_ info_)
2079 { // ** body not listed **
2548 }
2549
2550 real_4 _cmach (int_4 _p_ job_)
2551 { // ** body not listed **
2591 }
2592
2593 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_,
2594 complex_8 _p_ s_)
2595 { // ** body not listed **
2629 }
2630
2631 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_,
2632 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_
2633 work_, real_4 _p_ rwork_, int_4 _p_ info_)
2634 { // ** body not listed **
2757 }
2758
2759 if (leftv_) {
2760 is_ = 1;
2761 for (ki_ = 1; ki_ <= (*n_); (ki_)++) {
2762 if (somev_) {
2763 if (! (select_[ki_ - 1])) {
2764 goto _l130;
2765 }
2766 }
2767 smin_ = _max (ulp_ * (({complex_8 cdum__1_ = t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))]; _abs ((real_4) (crealf (cdum__1_)))
2768 + _abs (cimagf (cdum__1_));})), smlnum_);
2769 work_[(*n_) - 1] = CMPLXF (1.0, 0.0);
2770 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2771 work_[k_ - 1] = -conjf (t_[ki_ - 1 + ((*ldt_)) * ((k_ - 1))]);
2772 }
2773 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2774 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] - t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))];
2775 if (({complex_8 cdum__1_ = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))]; _abs ((real_4) (crealf (cdum__1_))) + _abs (cimagf
2776 (cdum__1_));}) < smin_) {
2777 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = smin_;
2778 }
2779 }
2780 if (ki_ < (*n_)) {
2781 (void) _clatrs (_dc_11, _dc_14, _dc_12, _dc_35, ($6_ = (*n_) - ki_, &$6_), &t_[ki_ + 1 - 1 + ((*ldt_)) * ((ki_ + 1 -
2782 1))], ldt_, &work_[ki_ + 1 - 1], &scale_, rwork_, info_);
2783 work_[ki_ - 1] = scale_;
2784 }
2785 if (! (over_)) {
2786 (void) _ccopy (($7_ = (*n_) - ki_ + 1, &$7_), &work_[ki_ - 1], &_k1, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
2787 ii_ = _icamax (($8_ = (*n_) - ki_ + 1, &$8_), &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1) + ki_ - 1;
2788 remax_ = 1.0 / ({complex_8 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((is_ - 1))]; _abs ((real_4) (crealf (cdum__1_))) +
2789 _abs (cimagf (cdum__1_));});
2790 (void) _csscal (($9_ = (*n_) - ki_ + 1, &$9_), &remax_, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
2791 for (k_ = 1; k_ <= ki_ - 1; (k_)++) {
2792 vl_[k_ - 1 + ((*ldvl_)) * ((is_ - 1))] = CMPLXF (0.0, 0.0);
2793 }
2794 }
2795 else {
2796 if (ki_ < (*n_)) {
2797 (void) _cgemv (_dc_7, n_, ($10_ = (*n_) - ki_, &$10_), ($11_ = CMPLXF (1.0, 0.0), &$11_), &vl_[0 + ((*ldvl_)) * ((ki_ +
2798 1 - 1))], ldvl_, &work_[ki_ + 1 - 1], &_k1, ($12_ = (complex_8) (scale_), &$12_), &vl_[0 + ((*ldvl_)) * ((ki_ - 1))],
2799 &_k1);
2800 }
2801 ii_ = _icamax (n_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
2802 remax_ = 1.0 / ({complex_8 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((ki_ - 1))]; _abs ((real_4) (crealf (cdum__1_))) +
2803 _abs (cimagf (cdum__1_));});
2804 (void) _csscal (n_, &remax_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
2805 }
2806 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
2807 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = work_[k_ + (*n_) - 1];
2808 }
2809 is_ = is_ + 1;
2810 _l130:;
2811 }
2812
2813 }
2814
2815 _l0:;
2816 __calls[24].calls++;
2817 return 0;
2818 }
2819
2820 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
2821 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_)
2822 { // ** body not listed **
2888 }
2889
2890 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_,
2891 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_,
2892 int_4 _p_ lwork_, int_4 _p_ info_)
2893 { // ** body not listed **
2991 }
2992
2993 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
2994 _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_)
2995 { // ** body not listed **
3182 }
3183
3184 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_
3185 lda_, complex_8 _p_ tau_, complex_8 _p_ c_, int_4 _p_ ldc_, complex_8 _p_ work_, int_4 _p_ info_)
3186 { // ** body not listed **
3272 }
3273
3274 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_
3275 ipiv_, int_4 _p_ info_)
3276 { // ** body not listed **
3349 }
3350
3351 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_
3352 ipiv_, int_4 _p_ info_)
3353 { // ** body not listed **
3549 }
3550
3551 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_
3552 ldab_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
3553 { // ** body not listed **
3633 }
3634
3635 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_
3636 info_)
3637 { // ** body not listed **
3692 }
3693
3694 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
3695 _p_ du2_, int_4 _p_ ipiv_, real_8 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
3696 { // ** body not listed **
3774 }
3775
3776 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_)
3777 { // ** body not listed **
3872 }
3873
3874 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_)
3875 { // ** body not listed **
3891 }
3892
3893 int_4 _dlae2 (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ rt1_, real_8 _p_ rt2_)
3894 { // ** body not listed **
3933 }
3934
3935 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
3936 _p_ sn1_)
3937 { // ** body not listed **
4010 }
4011
4012 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
4013 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ work_, int_4 _p_ info_)
4014 { // ** body not listed **
4184 }
4185
4186 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
4187 _p_ du_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_)
4188 { // ** body not listed **
4287 }
4288
4289 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_,
4290 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
4291 _p_ info_)
4292 { // ** body not listed **
4487 }
4488
4489 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_,
4490 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
4491 _p_ x_, int_4 _p_ ldx_, real_8 _p_ scale_, real_8 _p_ xnorm_, int_4 _p_ info_)
4492 { // ** body not listed **
4817 }
4818
4819 real_8 _dlanhs (char _p_ norm_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
4820 { // ** body not listed **
4873 }
4874
4875 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
4876 _p_ rt2r_, real_8 _p_ rt2i_, real_8 _p_ cs_, real_8 _p_ sn_)
4877 { // ** body not listed **
4959 }
4960
4961 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_
4962 ldx_, real_8 _p_ beta_, real_8 _p_ b_, int_4 _p_ ldb_)
4963 { // ** body not listed **
5021 }
5022
5023 real_8 _dlapy3 (real_8 _p_ x_, real_8 _p_ y_, real_8 _p_ z_)
5024 { // ** body not listed **
5041 }
5042
5043 real_8 _dlaran (int_4 _p_ iseed_)
5044 { // ** body not listed **
5066 }
5067
5068 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_
5069 ldc_, real_8 _p_ work_)
5070 { // ** body not listed **
5577 }
5578
5579 real_8 _dlarnd (int_4 _p_ idist_, int_4 _p_ iseed_)
5580 { // ** body not listed **
5596 }
5597
5598 int_4 _dlarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, real_8 _p_ x_)
5599 { // ** body not listed **
5629 }
5630
5631 int_4 _dlaruv (int_4 _p_ iseed_, int_4 _p_ n_, real_8 _p_ x_)
5632 { // ** body not listed **
7721 }
7722
7723 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_
7724 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_,
7725 int_4 _p_ ldx_, real_8 _p_ xnorm_, int_4 _p_ info_)
7726 { // ** body not listed **
8036 }
8037
8038 int_4 _dpttrf (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ info_)
8039 { // ** body not listed **
8072 }
8073
8074 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_
8075 info_)
8076 { // ** body not listed **
8111 }
8112
8113 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_
8114 work_, int_4 _p_ info_)
8115 { // ** body not listed **
8444 }
8445
8446 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_,
8447 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
8448 _p_ info_)
8449 { // ** body not listed **
9062 }
9063
9064 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_
9065 ifst_, int_4 _p_ ilst_, real_8 _p_ work_, int_4 _p_ info_)
9066 { // ** body not listed **
9256 }
9257
9258 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_,
9259 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_
9260 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
9261 { // ** body not listed **
9423 }
9424
9425 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_
9426 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_)
9427 { // ** body not listed **
10165 }
10166
10167 real_8 _dzsum1 (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_)
10168 { // ** body not listed **
10194 }
10195
10196 int_4 _icmax1 (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_)
10197 { // ** body not listed **
10224 }
10225
10226 goto _l0;
10227 _l30:;
10228 smax_ = ({complex_8 zdum__1_ = cx_[0]; _abs ((real_4) (crealf (zdum__1_)));});
10229 for (i_ = 2; i_ <= (*n_); (i_)++) {
10230 if (({complex_8 zdum__1_ = cx_[i_ - 1]; _abs ((real_4) (crealf (zdum__1_)));}) <= smax_) {
10231 goto _l40;
10232 }
10233
10234 icmax1_ = i_;
10235 smax_ = ({complex_8 zdum__1_ = cx_[i_ - 1]; _abs ((real_4) (crealf (zdum__1_)));});
10236 _l40:;
10237 }
10238
10239 _l0:;
10240 __calls[60].calls++;
10241 return icmax1_;
10242 }
10243
10244 int_4 _izmax1 (int_4 _p_ n_, complex_16 _p_ cx_, int_4 _p_ incx_)
10245 { // ** body not listed **
10272 }
10273
10274 goto _l0;
10275 _l30:;
10276 smax_ = ({complex_16 zdum__1_ = cx_[0]; _abs ((real_8) (creal (zdum__1_)));});
10277 for (i_ = 2; i_ <= (*n_); (i_)++) {
10278 if (({complex_16 zdum__1_ = cx_[i_ - 1]; _abs ((real_8) (creal (zdum__1_)));}) <= smax_) {
10279 goto _l40;
10280 }
10281
10282 izmax1_ = i_;
10283 smax_ = ({complex_16 zdum__1_ = cx_[i_ - 1]; _abs ((real_8) (creal (zdum__1_)));});
10284 _l40:;
10285 }
10286
10287 _l0:;
10288 __calls[61].calls++;
10289 return izmax1_;
10290 }
10291
10292 real_4 _scsum1 (int_4 _p_ n_, complex_8 _p_ cx_, int_4 _p_ incx_)
10293 { // ** body not listed **
10319 }
10320
10321 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_
10322 ipiv_, int_4 _p_ info_)
10323 { // ** body not listed **
10396 }
10397
10398 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_
10399 ipiv_, int_4 _p_ info_)
10400 { // ** body not listed **
10595 }
10596 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_
10597 ldab_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
10598 { // ** body not listed **
10678 }
10679
10680 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_
10681 info_)
10682 { // ** body not listed **
10716 }
10717
10718 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_
10719 info_)
10720 { // ** body not listed **
10775 }
10776
10777 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
10778 _p_ du2_, int_4 _p_ ipiv_, real_4 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
10779 { // ** body not listed **
10857 }
10858
10859 int_4 _slabad (real_4 _p_ small_, real_4 _p_ large_)
10860 { // ** body not listed **
10867 }
10868
10869 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_)
10870 { // ** body not listed **
10965 }
10966
10967 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_
10968 )
10969 { // ** body not listed **
10995 }
10996
10997 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_)
10998 { // ** body not listed **
11014 }
11015
11016 int_4 _slae2 (real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ rt1_, real_4 _p_ rt2_)
11017 { // ** body not listed **
11056 }
11057
11058 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
11059 _p_ sn1_)
11060 { // ** body not listed **
11133 }
11134
11135 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
11136 _p_ j1_, int_4 _p_ n1_, int_4 _p_ n2_, real_4 _p_ work_, int_4 _p_ info_)
11137 { // ** body not listed **
11302 }
11303
11304 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
11305 _p_ du_, real_4 _p_ x_, int_4 _p_ ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_)
11306 { // ** body not listed **
11405 }
11406
11407 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_,
11408 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
11409 _p_ info_)
11410 { // ** body not listed **
11602 }
11603
11604 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_,
11605 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
11606 _p_ x_, int_4 _p_ ldx_, real_4 _p_ scale_, real_4 _p_ xnorm_, int_4 _p_ info_)
11607 { // ** body not listed **
11931 }
11932
11933 real_4 _slamch (char _p_ cmach_)
11934 { // ** body not listed **
12007 }
12008
12009 int_4 _slamc1 (int_4 _p_ beta_, int_4 _p_ t_, logical_4 _p_ rnd_, logical_4 _p_ ieee1_)
12010 { // ** body not listed **
12086 }
12087
12088 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_,
12089 int_4 _p_ emax_, real_4 _p_ rmax_)
12090 { // ** body not listed **
12263 }
12264
12265 real_4 _slamc3 (real_4 _p_ a_, real_4 _p_ b_)
12266 { // ** body not listed **
12271 }
12272
12273 int_4 _slamc4 (int_4 _p_ emin_, real_4 _p_ start_, int_4 _p_ base_)
12274 { // ** body not listed **
12308 }
12309
12310 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_)
12311 { // ** body not listed **
12364 }
12365
12366 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_)
12367 { // ** body not listed **
12418 }
12419
12420 real_4 _slanhs (char _p_ norm_, int_4 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, real_4 _p_ work_)
12421 { // ** body not listed **
12473 }
12474
12475 real_4 _slanst (char _p_ norm_, int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_)
12476 { // ** body not listed **
12516 }
12517
12518 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
12519 _p_ rt2r_, real_4 _p_ rt2i_, real_4 _p_ cs_, real_4 _p_ sn_)
12520 { // ** body not listed **
12601 }
12602
12603 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_
12604 ldx_, real_4 _p_ beta_, real_4 _p_ b_, int_4 _p_ ldb_)
12605 { // ** body not listed **
12663 }
12664
12665 real_4 _slapy2 (real_4 _p_ x_, real_4 _p_ y_)
12666 { // ** body not listed **
12682 }
12683
12684 real_4 _slapy3 (real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ z_)
12685 { // ** body not listed **
12702 }
12703
12704 real_4 _slaran (int_4 _p_ iseed_)
12705 { // ** body not listed **
12727 }
12728
12729 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_
12730 c_, int_4 _p_ ldc_, real_4 _p_ work_)
12731 { // ** body not listed **
12750 }
12751
12752 int_4 _slarfg (int_4 _p_ n_, real_4 _p_ alpha_, real_4 _p_ x_, int_4 _p_ incx_, real_4 _p_ tau_)
12753 { // ** body not listed **
12800 }
12801
12802 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_
12803 ldc_, real_4 _p_ work_)
12804 { // ** body not listed **
13311 }
13312
13313 real_4 _slarnd (int_4 _p_ idist_, int_4 _p_ iseed_)
13314 { // ** body not listed **
13330 }
13331
13332 int_4 _slarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, real_4 _p_ x_)
13333 { // ** body not listed **
13363 }
13364
13365 int_4 _slartg (real_4 _p_ f_, real_4 _p_ g_, real_4 _p_ cs_, real_4 _p_ sn_, real_4 _p_ r_)
13366 { // ** body not listed **
13450 }
13451
13452 int_4 _slaruv (int_4 _p_ iseed_, int_4 _p_ n_, real_4 _p_ x_)
13453 { // ** body not listed **
15543 }
15544
15545 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
15546 _p_ n_, real_4 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
15547 { // ** body not listed **
15697 }
15698
15699 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
15700 _p_ lda_)
15701 { // ** body not listed **
15730 }
15731
15732 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_
15733 s_, real_4 _p_ a_, int_4 _p_ lda_)
15734 { // ** body not listed **
15940 }
15941
15942 int_4 _slasrt (char _p_ id_, int_4 _p_ n_, real_4 _p_ d_, int_4 _p_ info_)
15943 { // ** body not listed **
16117 }
16118
16119 int_4 _slassq (int_4 _p_ n_, real_4 _p_ x_, int_4 _p_ incx_, real_4 _p_ scale_, real_4 _p_ sumsq_)
16120 { // ** body not listed **
16141 }
16142
16143 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_
16144 incx_)
16145 { // ** body not listed **
16186 }
16187
16188 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_
16189 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_,
16190 int_4 _p_ ldx_, real_4 _p_ xnorm_, int_4 _p_ info_)
16191 { // ** body not listed **
16500 }
16501
16502 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_
16503 , real_4 _p_ tau_, real_4 _p_ c_, int_4 _p_ ldc_, real_4 _p_ work_, int_4 _p_ info_)
16504 { // ** body not listed **
16584 }
16585
16586 int_4 _spttrf (int_4 _p_ n_, real_4 _p_ d_, real_4 _p_ e_, int_4 _p_ info_)
16587 { // ** body not listed **
16620 }
16621
16622 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_
16623 info_)
16624 { // ** body not listed **
16659 }
16660
16661 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_
16662 work_, int_4 _p_ info_)
16663 { // ** body not listed **
16984 }
16985
16986 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_,
16987 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
16988 _p_ info_)
16989 { // ** body not listed **
17600 }
17601
17602 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_
17603 ifst_, int_4 _p_ ilst_, real_4 _p_ work_, int_4 _p_ info_)
17604 { // ** body not listed **
17794 }
17795
17796 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_,
17797 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_
17798 work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ liwork_, int_4 _p_ info_)
17799 { // ** body not listed **
17959 }
17960
17961 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_
17962 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_)
17963 { // ** body not listed **
18698 }
18699
18700 int_4 _xlaenv (int_4 _p_ ispec_, int_4 _p_ nvalue_)
18701 { // ** body not listed **
18707 }
18708
18709 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_
18710 ipiv_, int_4 _p_ info_)
18711 { // ** body not listed **
18784 }
18785
18786 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_
18787 ipiv_, int_4 _p_ info_)
18788 { // ** body not listed **
18985 }
18986
18987 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
18988 _p_ ldab_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
18989 { // ** body not listed **
19090 }
19091
19092 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_
19093 , int_4 _p_ info_)
19094 { // ** body not listed **
19131 }
19132
19133 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_
19134 ipiv_, int_4 _p_ info_)
19135 { // ** body not listed **
19183 }
19184
19185 if (d_[(*n_) - 1] == CMPLX (0.0, 0.0) && (*info_) == 0) {
19186 (*info_) = (*n_);
19187 goto _l0;
19188 }
19189
19190 _l0:;
19191 __calls[120].calls++;
19192 return 0;
19193 }
19194
19195 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_
19196 du_, complex_16 _p_ du2_, int_4 _p_ ipiv_, complex_16 _p_ b_, int_4 _p_ ldb_, int_4 _p_ info_)
19197 { // ** body not listed **
19302 }
19303
19304 int_4 _zlacgv (int_4 _p_ n_, complex_16 _p_ x_, int_4 _p_ incx_)
19305 { // ** body not listed **
19326 }
19327
19328 int_4 _zlacon (int_4 _p_ n_, complex_16 _p_ v_, complex_16 _p_ x_, real_8 _p_ est_, int_4 _p_ kase_)
19329 { // ** body not listed **
19425 }
19426
19427 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
19428 _p_ ldb_)
19429 { // ** body not listed **
19455 }
19456
19457 complex_16 _zladiv (complex_16 _p_ x_, complex_16 _p_ y_)
19458 { // ** body not listed **
19467 }
19468
19469 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_
19470 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_
19471 info_)
19472 { // ** body not listed **
19648 }
19649
19650 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_)
19651 { // ** body not listed **
19702 }
19703
19704 real_8 _zlanhs (char _p_ norm_, int_4 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
19705 { // ** body not listed **
19757 }
19758
19759 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_,
19760 complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_)
19761 { // ** body not listed **
19782 }
19783
19784 int_4 _zlarfg (int_4 _p_ n_, complex_16 _p_ alpha_, complex_16 _p_ x_, int_4 _p_ incx_, complex_16 _p_ tau_)
19785 { // ** body not listed **
19840 }
19841
19842 int_4 _zlarnv (int_4 _p_ idist_, int_4 _p_ iseed_, int_4 _p_ n_, complex_16 _p_ x_)
19843 { // ** body not listed **
19879 }
19880
19881 int_4 _zlartg (complex_16 _p_ f_, complex_16 _p_ g_, real_8 _p_ cs_, complex_16 _p_ sn_, complex_16 _p_ r_)
19882 { // ** body not listed **
19928 }
19929
19930 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
19931 _p_ n_, complex_16 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
19932 { // ** body not listed **
20083 }
20084
20085 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_
20086 a_, int_4 _p_ lda_)
20087 { // ** body not listed **
20122 }
20123
20124 int_4 _zlassq (int_4 _p_ n_, complex_16 _p_ x_, int_4 _p_ incx_, real_8 _p_ scale_, real_8 _p_ sumsq_)
20125 { // ** body not listed **
20156 }
20157
20158 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
20159 _p_ incx_)
20160 { // ** body not listed **
20201 }
20202
20203 int_4 _zlatrs (char _p_ uplo_, char _p_ trans_, char _p_ diag_, char _p_ normin_, int_4 _p_ n_, complex_16 _p_ a_,
20204 int_4 _p_ lda_, complex_16 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_)
20205 { // ** body not listed **
20678 }
20679
20680 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_,
20681 complex_16 _p_ s_)
20682 { // ** body not listed **
20718 }
20719
20720 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_
20721 , 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_
20722 work_, real_8 _p_ rwork_, int_4 _p_ info_)
20723 { // ** body not listed **
20850 }
20851
20852 if (leftv_) {
20853 is_ = 1;
20854 for (ki_ = 1; ki_ <= (*n_); (ki_)++) {
20855 if (somev_) {
20856 if (! (select_[ki_ - 1])) {
20857 goto _l130;
20858 }
20859 }
20860 smin_ = _max (ulp_ * (({complex_16 cdum__1_ = t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))]; _abs ((real_8) (creal (cdum__1_)))
20861 + _abs (cimag (cdum__1_));})), smlnum_);
20862 work_[(*n_) - 1] = CMPLX (1.0, 0.0);
20863 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20864 work_[k_ - 1] = -__dconjg (&t_[ki_ - 1 + ((*ldt_)) * ((k_ - 1))]);
20865 }
20866 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20867 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] - t_[ki_ - 1 + ((*ldt_)) * ((ki_ - 1))];
20868 if (({complex_16 cdum__1_ = t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))]; _abs ((real_8) (creal (cdum__1_))) + _abs (cimag
20869 (cdum__1_));}) < smin_) {
20870 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = smin_;
20871 }
20872 }
20873 if (ki_ < (*n_)) {
20874 (void) _zlatrs (_dc_11, _dc_14, _dc_12, _dc_35, ($6_ = (*n_) - ki_, &$6_), &t_[ki_ + 1 - 1 + ((*ldt_)) * ((ki_ + 1 -
20875 1))], ldt_, &work_[ki_ + 1 - 1], &scale_, rwork_, info_);
20876 work_[ki_ - 1] = scale_;
20877 }
20878 if (! (over_)) {
20879 (void) _zcopy (($7_ = (*n_) - ki_ + 1, &$7_), &work_[ki_ - 1], &_k1, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
20880 ii_ = _izamax (($8_ = (*n_) - ki_ + 1, &$8_), &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1) + ki_ - 1;
20881 remax_ = 1.0 / ({complex_16 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((is_ - 1))]; _abs ((real_8) (creal (cdum__1_))) +
20882 _abs (cimag (cdum__1_));});
20883 (void) _zdscal (($9_ = (*n_) - ki_ + 1, &$9_), &remax_, &vl_[ki_ - 1 + ((*ldvl_)) * ((is_ - 1))], &_k1);
20884 for (k_ = 1; k_ <= ki_ - 1; (k_)++) {
20885 vl_[k_ - 1 + ((*ldvl_)) * ((is_ - 1))] = CMPLX (0.0, 0.0);
20886 }
20887 }
20888 else {
20889 if (ki_ < (*n_)) {
20890 (void) _zgemv (_dc_7, n_, ($10_ = (*n_) - ki_, &$10_), ($11_ = CMPLX (1.0, 0.0), &$11_), &vl_[0 + ((*ldvl_)) * ((ki_ +
20891 1 - 1))], ldvl_, &work_[ki_ + 1 - 1], &_k1, ($12_ = (complex_16) (scale_), &$12_), &vl_[0 + ((*ldvl_)) * ((ki_ - 1))],
20892 &_k1);
20893 }
20894 ii_ = _izamax (n_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
20895 remax_ = 1.0 / ({complex_16 cdum__1_ = vl_[ii_ - 1 + ((*ldvl_)) * ((ki_ - 1))]; _abs ((real_8) (creal (cdum__1_))) +
20896 _abs (cimag (cdum__1_));});
20897 (void) _zdscal (n_, &remax_, &vl_[0 + ((*ldvl_)) * ((ki_ - 1))], &_k1);
20898 }
20899 for (k_ = ki_ + 1; k_ <= (*n_); (k_)++) {
20900 t_[k_ - 1 + ((*ldt_)) * ((k_ - 1))] = work_[k_ + (*n_) - 1];
20901 }
20902 is_ = is_ + 1;
20903 _l130:;
20904 }
20905
20906 }
20907
20908 _l0:;
20909 __calls[139].calls++;
20910 return 0;
20911 }
20912
20913 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_,
20914 int_4 _p_ ifst_, int_4 _p_ ilst_, int_4 _p_ info_)
20915 { // ** body not listed **
20983 }
20984
20985 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_,
20986 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_
20987 work_, int_4 _p_ lwork_, int_4 _p_ info_)
20988 { // ** body not listed **
21086 }
21087
21088 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
21089 _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_)
21090 { // ** body not listed **
21282 }
21283
21284 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_
21285 lda_, complex_16 _p_ tau_, complex_16 _p_ c_, int_4 _p_ ldc_, complex_16 _p_ work_, int_4 _p_ info_)
21286 { // ** body not listed **
21374 }
21375
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
