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)