mathlib-eispack.c

     1  //! @file mathlib-eispack.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  //! 81 subprograms from EISPACK
    24  
    25  // Compiled from Fortran source code by VIF.
    26  // Selected subprograms are:
    27  //
    28  // BAKVEC CDIV   COMHES ELMBAK HTRIB3 MINFIT RATQR  RGORT  RSPP   TQLRAT
    29  // BALANC CG     COMLR  ELMHES HTRIBK ORTBAK REBAK  RS     RST    TRBAK1
    30  // BALBAK CGLR   COMLR2 ELTRAN HTRID3 ORTHES REBAKB RSB    RT     TRBAK3
    31  // BANDR  CGQR   COMQR  EPSLON HTRIDI ORTRAN REDUC  RSG    SSPEV  TRED1
    32  // BANDV  CH     COMQR2 FIGI   IMTQL1 QZHES  REDUC2 RSGAB  SVD    TRED2
    33  // BISECT CH3    CORTB  FIGI2  IMTQL2 QZIT   RG     RSGBA  TINVIT TRED3
    34  // BQR    CINVIT CORTH  HQR    IMTQLV QZVAL  RGELM  RSM    TQL1   TRIDIB
    35  // CBABK2 COMBAK CSROOT HQR2   INVIT  QZVEC  RGG    RSP    TQL2   TSTURM
    36  // CBAL
    37  
    38  // EISPACK is a collection of double-precision Fortran subroutines that
    39  // compute the eigenvalues and eigenvectors of nine classes of matrices:
    40  // complex general, complex Hermitian, real general, real symmetric, real
    41  // symmetric banded, real symmetric tridiagonal, special real
    42  // tridiagonal, generalized real, and generalized real symmetric matices.
    43  // In addition, two routines are included that use singular value
    44  // decomposition to solve certain least-squares problems. 
    45  // 
    46  // Originally written around 1972–1973,[1] EISPACK, like LINPACK and MINPACK, 
    47  // originated from Argonne National Laboratory, has always been free, and 
    48  // aims to be portable, robust and reliable. The library drew heavily on 
    49  // algorithms developed by James Wilkinson, which were originally implemented 
    50  // in ALGOL. Brian Smith led a team at Argonne developing EISPACK, initially 
    51  // by translating these algorithms into FORTRAN. Jack Dongarra joined the 
    52  // team as an undergraduate intern at Argonne, and later went on to create 
    53  // LAPACK, which has largely superseded EISPACK and LINPACK. 
    54  // 
    55  // This version of EISPACK is adapted for VIF.
    56  // Adaptations are trivial, concerning choosing alternative names for
    57  // variable names that coincided with intrinsic functions.
    58  // 
    59  // Source: netlib.org/eispack/
    60  
    61  // The license for EISPACK Fortran source code is:
    62  //
    63  // EISPACK, like LINPACK and MINPACK, originated from Argonne National Laboratory 
    64  // and has always been free. EISPACK did not have a copyright notice or involved 
    65  // any kind of license. The term "open source" was not yet widely used at the time.
    66  // 
    67  // Source: netlib.org/eispack/
    68  // 
    69  // This particular version is John Burkardt's adaptation of original EISPACK code.
    70  // John Burkardt released his version under GNU LGPL.
    71  
    72  /*
    73  Generated by VIF - experimental VIntage Fortran compiler.
    74  VIF release 1.1.7
    75  */
    92  static CALLS __calls[__ncalls] = {
    93    {"bakvec", 0}, // subroutine
    94    {"balanc", 0}, // subroutine
    95    {"balbak", 0}, // subroutine
    96    {"bandr", 0}, // subroutine
    97    {"bandv", 0}, // subroutine
    98    {"bisect", 0}, // subroutine
    99    {"bqr", 0}, // subroutine
   100    {"cbabk2", 0}, // subroutine
   101    {"cbal", 0}, // subroutine
   102    {"cdiv", 0}, // subroutine
   103    {"cg", 0}, // subroutine
   104    {"cglr", 0}, // subroutine
   105    {"cgqr", 0}, // subroutine
   106    {"ch3", 0}, // subroutine
   107    {"ch", 0}, // subroutine
   108    {"cinvit", 0}, // subroutine
   109    {"combak", 0}, // subroutine
   110    {"comhes", 0}, // subroutine
   111    {"comlr2", 0}, // subroutine
   112    {"comlr", 0}, // subroutine
   113    {"comqr2", 0}, // subroutine
   114    {"comqr", 0}, // subroutine
   115    {"cortb", 0}, // subroutine
   116    {"corth", 0}, // subroutine
   117    {"csroot", 0}, // subroutine
   118    {"elmbak", 0}, // subroutine
   119    {"elmhes", 0}, // subroutine
   120    {"eltran", 0}, // subroutine
   121    {"epslon", 0}, // real*8 function
   122    {"figi2", 0}, // subroutine
   123    {"figi", 0}, // subroutine
   124    {"hqr2", 0}, // subroutine
   125    {"hqr", 0}, // subroutine
   126    {"htrib3", 0}, // subroutine
   127    {"htribk", 0}, // subroutine
   128    {"htrid3", 0}, // subroutine
   129    {"htridi", 0}, // subroutine
   130    {"imtql1", 0}, // subroutine
   131    {"imtql2", 0}, // subroutine
   132    {"imtqlv", 0}, // subroutine
   133    {"invit", 0}, // subroutine
   134    {"minfit", 0}, // subroutine
   135    {"ortbak", 0}, // subroutine
   136    {"orthes", 0}, // subroutine
   137    {"ortran", 0}, // subroutine
   138    {"qzhes", 0}, // subroutine
   139    {"qzit", 0}, // subroutine
   140    {"qzval", 0}, // subroutine
   141    {"qzvec", 0}, // subroutine
   142    {"ratqr", 0}, // subroutine
   143    {"rebakb", 0}, // subroutine
   144    {"rebak", 0}, // subroutine
   145    {"reduc2", 0}, // subroutine
   146    {"reduc", 0}, // subroutine
   147    {"rgelm", 0}, // subroutine
   148    {"rg", 0}, // subroutine
   149    {"rgg", 0}, // subroutine
   150    {"rgort", 0}, // subroutine
   151    {"rsb", 0}, // subroutine
   152    {"rs", 0}, // subroutine
   153    {"rsgab", 0}, // subroutine
   154    {"rsgba", 0}, // subroutine
   155    {"rsg", 0}, // subroutine
   156    {"rsm", 0}, // subroutine
   157    {"rsp", 0}, // subroutine
   158    {"rspp", 0}, // subroutine
   159    {"rst", 0}, // subroutine
   160    {"rt", 0}, // subroutine
   161    {"sspev", 0}, // subroutine
   162    {"svd", 0}, // subroutine
   163    {"tinvit", 0}, // subroutine
   164    {"tql1", 0}, // subroutine
   165    {"tql2", 0}, // subroutine
   166    {"tqlrat", 0}, // subroutine
   167    {"trbak1", 0}, // subroutine
   168    {"trbak3", 0}, // subroutine
   169    {"tred1", 0}, // subroutine
   170    {"tred2", 0}, // subroutine
   171    {"tred3", 0}, // subroutine
   172    {"tridib", 0}, // subroutine
   173    {"tsturm", 0}, // subroutine
   174    {NULL, 0}
   175  };
   177  _p_ ierr_);
   179  ;
   181  real_8 _p_ z_);
   183  _p_ e2_, logical_4 _p_ matz_, real_8 _p_ z_);
   185  real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_);
   187  real_8 _p_ ub_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ w_, int_4 _p_ ind_, int_4 _p_ ierr_, real_8 _p_ rv4_, real_8 
   188  _p_ rv5_);
   190  _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_);
   192  real_8 _p_ zr_, real_8 _p_ zi_);
   194  real_8 _p_ scale_);
   197  _p_ matz_, real_4 _p_ zr_, real_4 _p_ zi_, real_4 _p_ fv1_, real_4 _p_ fv2_, real_4 _p_ fv3_, int_4 _p_ ierr_);
   199  int_4 _p_ matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_);
   201  int_4 _p_ matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fv3_, int_4 _p_ ierr_);
   203  real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
   205  real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
   207  logical_4 _p_ select_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_, real_8 _p_ rm1_, 
   208  real_8 _p_ rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_);
   210  int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
   212  int_4 _p_ int_);
   214  real_8 _p_ hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_);
   216  real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
   218  orti_, real_8 _p_ hr_, real_8 _p_ hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_)
   219  ;
   221  real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
   223  , real_8 _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
   225  real_8 _p_ ortr_, real_8 _p_ orti_);
   228  real_8 _p_ z_);
   231  real_8 _p_ z_);
   234  _p_ ierr_);
   236  _p_ ierr_);
   238  real_8 _p_ wi_, real_8 _p_ z_, int_4 _p_ ierr_);
   240  real_8 _p_ wi_, int_4 _p_ ierr_);
   242  real_8 _p_ zi_);
   244  real_8 _p_ zr_, real_8 _p_ zi_);
   246  real_8 _p_ tau_);
   248  real_8 _p_ e2_, real_8 _p_ tau_);
   252  int_4 _p_ ierr_, real_8 _p_ rv1_);
   254  select_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ z_, int_4 _p_ ierr_, real_8 _p_ rm1_, real_8 _p_ rv1_, real_8 _p_ rv2_
   255  );
   257  _p_ b_, int_4 _p_ ierr_, real_8 _p_ rv1_);
   259  real_8 _p_ z_);
   262  real_8 _p_ z_);
   265  , real_8 _p_ z_, int_4 _p_ ierr_);
   267  real_8 _p_ beta_, logical_4 _p_ matz_, real_8 _p_ z_);
   269  real_8 _p_ beta_, real_8 _p_ z_);
   271  real_8 _p_ w_, int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_);
   277  real_8 _p_ z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_);
   279  real_4 _p_ z_, int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_);
   281  real_8 _p_ beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_);
   283  real_8 _p_ z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_);
   285  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   287  _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   289  real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   291  real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   293  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   295  _p_ fwork_, int_4 _p_ iwork_, int_4 _p_ ierr_);
   297  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   299  _p_ ierr_, int_4 _p_ m_, logical_4 _p_ type_);
   301  _p_ ierr_);
   303  _p_ fv1_, int_4 _p_ ierr_);
   305  int_4 _p_ job_, int_4 _p_ info_);
   307  real_8 _p_ u_, logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_);
   309  real_8 _p_ w_, int_4 _p_ ind_, real_8 _p_ z_, int_4 _p_ ierr_, real_8 _p_ rv1_, real_8 _p_ rv2_, real_8 _p_ rv3_, 
   310  real_8 _p_ rv4_, real_8 _p_ rv6_);
   320  real_8 _p_ ub_, int_4 _p_ m11_, int_4 _p_ m_, real_8 _p_ w_, int_4 _p_ ind_, int_4 _p_ ierr_, real_8 _p_ rv4_, real_8 
   321  _p_ rv5_);
   323  real_8 _p_ lb_, real_8 _p_ ub_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ ierr_, real_8 _p_ 
   324  rv1_, real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_);
   325  int_4 _bakvec (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ t_, real_8 _p_ e_, int_4 _p_ m_, real_8 _p_ z_, int_4 _p_ ierr_)
   326  { // ** body not listed **
   360  }
   361  int_4 _balanc (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ scale_)
   362  { // ** body not listed **
   494  }
   495  int_4 _balbak (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ scale_, int_4 _p_ m_, real_8 _p_ 
   496  z_)
   497  { // ** body not listed **
   535  }
   536  int_4 _bandr (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ mb_, real_8 _p_ a_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, 
   537  logical_4 _p_ matz_, real_8 _p_ z_)
   538  { // ** body not listed **
   764  }
   765  int_4 _bandv (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ mbw_, real_8 _p_ a_, real_8 _p_ e21_, int_4 _p_ m_, real_8 _p_ w_, 
   766  real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_)
   767  { // ** body not listed **
  1018  }
  1019  int_4 _bisect (int_4 _p_ n_, real_8 _p_ eps1_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, real_8 _p_ lb_, real_8 _p_ 
  1020  ub_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ w_, int_4 _p_ ind_, int_4 _p_ ierr_, real_8 _p_ rv4_, real_8 _p_ rv5_)
  1021  { // ** body not listed **
  1242  }
  1243  int_4 _bqr (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ mb_, real_8 _p_ a_, real_8 _p_ t_, real_8 _p_ r_, int_4 _p_ ierr_, 
  1244  int_4 _p_ nv_, real_8 _p_ rv_)
  1245  { // ** body not listed **
  1443  }
  1444  int_4 _cbabk2 (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ scale_, int_4 _p_ m_, real_8 _p_ 
  1445  zr_, real_8 _p_ zi_)
  1446  { // ** body not listed **
  1488  }
  1489  int_4 _cbal (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ 
  1490  scale_)
  1491  { // ** body not listed **
  1631  }
  1632  int_4 _cdiv (real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ br_, real_8 _p_ bi_, real_8 _p_ cr_, real_8 _p_ ci_)
  1633  { // ** body not listed **
  1645  }
  1646  int_4 _cg (int_4 _p_ nm_, int_4 _p_ n_, real_4 _p_ ar_, real_4 _p_ ai_, real_4 _p_ wr_, real_4 _p_ wi_, int_4 _p_ matz_
  1647  , real_4 _p_ zr_, real_4 _p_ zi_, real_4 _p_ fv1_, real_4 _p_ fv2_, real_4 _p_ fv3_, int_4 _p_ ierr_)
  1648  { // ** body not listed **
  1672  }
  1673  int_4 _cglr (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ 
  1674  matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_)
  1675  { // ** body not listed **
  1696  }
  1697  int_4 _cgqr (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ 
  1698  matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fv3_, int_4 _p_ ierr_)
  1699  { // ** body not listed **
  1723  }
  1724  int_4 _ch3 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ zr_, real_8 _p_ zi_, 
  1725  real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
  1726  { // ** body not listed **
  1755  }
  1756  int_4 _ch (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ zr_, 
  1757  real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
  1758  { // ** body not listed **
  1787  }
  1788  int_4 _cinvit (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ wr_, real_8 _p_ wi_, logical_4 
  1789  _p_ select_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_, real_8 _p_ rm1_, real_8 _p_ 
  1790  rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_)
  1791  { // ** body not listed **
  1986  }
  1987  int_4 _combak (int_4 _p_ nm_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ ar_, real_8 _p_ ai_, int_4 _p_ int_, int_4 _p_ 
  1988  m_, real_8 _p_ zr_, real_8 _p_ zi_)
  1989  { // ** body not listed **
  2034  }
  2035  int_4 _comhes (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ ar_, real_8 _p_ ai_, int_4 _p_ 
  2036  int_)
  2037  { // ** body not listed **
  2114  }
  2115  int_4 _comlr2 (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, int_4 _p_ int_, real_8 _p_ hr_, real_8 _p_ 
  2116  hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_)
  2117  { // ** body not listed **
  2426  }
  2427  int_4 _comlr (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ hr_, real_8 _p_ hi_, real_8 _p_ 
  2428  wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
  2429  { // ** body not listed **
  2598  }
  2599  int_4 _comqr2 (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ ortr_, real_8 _p_ orti_, real_8 
  2600  _p_ hr_, real_8 _p_ hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_)
  2601  { // ** body not listed **
  2951  }
  2952  int_4 _comqr (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ hr_, real_8 _p_ hi_, real_8 _p_ 
  2953  wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
  2954  { // ** body not listed **
  3135  }
  3136  int_4 _cortb (int_4 _p_ nm_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ ortr_, real_8 
  3137  _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_)
  3138  { // ** body not listed **
  3180  }
  3181  int_4 _corth (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ 
  3182  ortr_, real_8 _p_ orti_)
  3183  { // ** body not listed **
  3263  }
  3264  int_4 _csroot (real_8 _p_ xr_, real_8 _p_ xi_, real_8 _p_ yr_, real_8 _p_ yi_)
  3265  { // ** body not listed **
  3288  }
  3289  int_4 _elmbak (int_4 _p_ nm_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ a_, int_4 _p_ int_, int_4 _p_ m_, real_8 _p_ 
  3290  z_)
  3291  { // ** body not listed **
  3329  }
  3330  int_4 _elmhes (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ a_, int_4 _p_ int_)
  3331  { // ** body not listed **
  3391  }
  3392  int_4 _eltran (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ a_, int_4 _p_ int_, real_8 _p_ 
  3393  z_)
  3394  { // ** body not listed **
  3426  }
  3427  real_8 _epslon (real_8 _p_ x_)
  3428  { // ** body not listed **
  3442  }
  3443  int_4 _figi2 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ t_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ z_, int_4 _p_ ierr_)
  3444  { // ** body not listed **
  3489  }
  3490  int_4 _figi (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ t_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, int_4 _p_ ierr_)
  3491  { // ** body not listed **
  3524  }
  3525  int_4 _hqr2 (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ h_, real_8 _p_ wr_, real_8 _p_ wi_
  3526  , real_8 _p_ z_, int_4 _p_ ierr_)
  3527  { // ** body not listed **
  3972  }
  3973  int_4 _hqr (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ h_, real_8 _p_ wr_, real_8 _p_ wi_, 
  3974  int_4 _p_ ierr_)
  3975  { // ** body not listed **
  4185  }
  4186  int_4 _htrib3 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ tau_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ 
  4187  zi_)
  4188  { // ** body not listed **
  4232  }
  4233  int_4 _htribk (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ tau_, int_4 _p_ m_, real_8 _p_ 
  4234  zr_, real_8 _p_ zi_)
  4235  { // ** body not listed **
  4279  }
  4280  int_4 _htrid3 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, real_8 _p_ 
  4281  tau_)
  4282  { // ** body not listed **
  4405  }
  4406  int_4 _htridi (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ ar_, real_8 _p_ ai_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ 
  4407  e2_, real_8 _p_ tau_)
  4408  { // ** body not listed **
  4523  }
  4524  int_4 _imtql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
  4525  { // ** body not listed **
  4614  }
  4615  int_4 _imtql2 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ z_, int_4 _p_ ierr_)
  4616  { // ** body not listed **
  4720  }
  4721  int_4 _imtqlv (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, real_8 _p_ w_, int_4 _p_ ind_, int_4 _p_ 
  4722  ierr_, real_8 _p_ rv1_)
  4723  { // ** body not listed **
  4830  }
  4831  int_4 _invit (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ wr_, real_8 _p_ wi_, logical_4 _p_ select_, int_4 
  4832  _p_ mm_, int_4 _p_ m_, real_8 _p_ z_, int_4 _p_ ierr_, real_8 _p_ rm1_, real_8 _p_ rv1_, real_8 _p_ rv2_)
  4833  { // ** body not listed **
  5217  }
  5218  int_4 _minfit (int_4 _p_ nm_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ ip_, real_8 _p_ b_, 
  5219  int_4 _p_ ierr_, real_8 _p_ rv1_)
  5220  { // ** body not listed **
  5490  }
  5491  int_4 _ortbak (int_4 _p_ nm_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ a_, real_8 _p_ ort_, int_4 _p_ m_, real_8 _p_ 
  5492  z_)
  5493  { // ** body not listed **
  5528  }
  5529  int_4 _orthes (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ a_, real_8 _p_ ort_)
  5530  { // ** body not listed **
  5586  }
  5587  int_4 _ortran (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ low_, int_4 _p_ igh_, real_8 _p_ a_, real_8 _p_ ort_, real_8 _p_ 
  5588  z_)
  5589  { // ** body not listed **
  5626  }
  5627  int_4 _qzhes (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, logical_4 _p_ matz_, real_8 _p_ z_)
  5628  { // ** body not listed **
  5754  }
  5755  int_4 _qzit (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ eps1_, logical_4 _p_ matz_, real_8 
  5756  _p_ z_, int_4 _p_ ierr_)
  5757  { // ** body not listed **
  6094  }
  6095  int_4 _qzval (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ alfr_, real_8 _p_ alfi_, real_8 _p_ 
  6096  beta_, logical_4 _p_ matz_, real_8 _p_ z_)
  6097  { // ** body not listed **
  6348  }
  6349  int_4 _qzvec (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ alfr_, real_8 _p_ alfi_, real_8 _p_ 
  6350  beta_, real_8 _p_ z_)
  6351  { // ** body not listed **
  6574  }
  6575  int_4 _ratqr (int_4 _p_ n_, real_8 _p_ eps1_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, int_4 _p_ m_, real_8 _p_ w_
  6576  , int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_)
  6577  { // ** body not listed **
  6743  }
  6744  int_4 _rebakb (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ b_, real_8 _p_ dl_, int_4 _p_ m_, real_8 _p_ z_)
  6745  { // ** body not listed **
  6769  }
  6770  int_4 _rebak (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ b_, real_8 _p_ dl_, int_4 _p_ m_, real_8 _p_ z_)
  6771  { // ** body not listed **
  6795  }
  6796  int_4 _reduc2 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ dl_, int_4 _p_ ierr_)
  6797  { // ** body not listed **
  6874  }
  6875  int_4 _reduc (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ dl_, int_4 _p_ ierr_)
  6876  { // ** body not listed **
  6953  }
  6954  int_4 _rgelm (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ matz_, real_8 _p_ 
  6955  z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_)
  6956  { // ** body not listed **
  6981  }
  6982  int_4 _rg (int_4 _p_ nm_, int_4 _p_ n_, real_4 _p_ a_, real_4 _p_ wr_, real_4 _p_ wi_, int_4 _p_ matz_, real_4 _p_ z_, 
  6983  int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_)
  6984  { // ** body not listed **
  7009  }
  7010  int_4 _rgg (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ alfr_, real_8 _p_ alfi_, real_8 _p_ 
  7011  beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_)
  7012  { // ** body not listed **
  7041  }
  7042  int_4 _rgort (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ matz_, real_8 _p_ 
  7043  z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_)
  7044  { // ** body not listed **
  7069  }
  7070  int_4 _rsb (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ mb_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, 
  7071  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7072  { // ** body not listed **
  7106  }
  7107  int_4 _rs (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, real_8 _p_ fv1_, 
  7108  real_8 _p_ fv2_, int_4 _p_ ierr_)
  7109  { // ** body not listed **
  7128  }
  7129  int_4 _rsgab (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, 
  7130  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7131  { // ** body not listed **
  7158  }
  7159  int_4 _rsgba (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, 
  7160  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7161  { // ** body not listed **
  7188  }
  7189  int_4 _rsg (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, 
  7190  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7191  { // ** body not listed **
  7218  }
  7219  int_4 _rsm (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ m_, real_8 _p_ z_, real_8 _p_ fwork_, 
  7220  int_4 _p_ iwork_, int_4 _p_ ierr_)
  7221  { // ** body not listed **
  7250  }
  7251  int_4 _rsp (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ nv_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, 
  7252  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7253  { // ** body not listed **
  7288  }
  7289  int_4 _rspp (int_4 _p_ n_, int_4 _p_ nv_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_
  7290  , int_4 _p_ m_, logical_4 _p_ type_)
  7291  { // ** body not listed **
  7515  }
  7516  int_4 _rst (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ w_, real_8 _p_ e_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_)
  7517  { // ** body not listed **
  7541  }
  7542  int_4 _rt (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_, int_4 _p_ matz_, real_8 _p_ z_, real_8 _p_ fv1_, 
  7543  int_4 _p_ ierr_)
  7544  { // ** body not listed **
  7569  }
  7570  int_4 _sspev (real_4 _p_ a_, int_4 _p_ n_, real_4 _p_ e_, real_4 _p_ v_, int_4 _p_ ldv_, real_4 _p_ work_, int_4 _p_ 
  7571  job_, int_4 _p_ info_)
  7572  { // ** body not listed **
  7615  }
  7616  int_4 _svd (int_4 _p_ nm_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_, logical_4 _p_ matu_, real_8 _p_ u_
  7617  , logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_)
  7618  { // ** body not listed **
  7927  }
  7928  int_4 _tinvit (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, int_4 _p_ m_, real_8 _p_ w_, 
  7929  int_4 _p_ ind_, real_8 _p_ z_, int_4 _p_ ierr_, real_8 _p_ rv1_, real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, 
  7930  real_8 _p_ rv6_)
  7931  { // ** body not listed **
  8122  }
  8123  int_4 _tql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
  8124  { // ** body not listed **
  8227  }
  8228  int_4 _tql2 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ z_, int_4 _p_ ierr_)
  8229  { // ** body not listed **
  8347  }
  8348  int_4 _tqlrat (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e2_, int_4 _p_ ierr_)
  8349  { // ** body not listed **
  8465  }
  8466  int_4 _trbak1 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ e_, int_4 _p_ m_, real_8 _p_ z_)
  8467  { // ** body not listed **
  8496  }
  8497  int_4 _trbak3 (int_4 _p_ nm_, int_4 _p_ n_, int_4 _p_ nv_, real_8 _p_ a_, int_4 _p_ m_, real_8 _p_ z_)
  8498  { // ** body not listed **
  8534  }
  8535  int_4 _tred1 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_)
  8536  { // ** body not listed **
  8624  }
  8625  int_4 _tred2 (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ z_)
  8626  { // ** body not listed **
  8743  }
  8744  int_4 _tred3 (int_4 _p_ n_, int_4 _p_ nv_, real_8 _p_ a_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_)
  8745  { // ** body not listed **
  8825  }
  8826  int_4 _tridib (int_4 _p_ n_, real_8 _p_ eps1_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, real_8 _p_ lb_, real_8 _p_ 
  8827  ub_, int_4 _p_ m11_, int_4 _p_ m_, real_8 _p_ w_, int_4 _p_ ind_, int_4 _p_ ierr_, real_8 _p_ rv4_, real_8 _p_ rv5_)
  8828  { // ** body not listed **
  9101  }
  9102  int_4 _tsturm (int_4 _p_ nm_, int_4 _p_ n_, real_8 _p_ eps1_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ e2_, real_8 _p_ 
  9103  lb_, real_8 _p_ ub_, int_4 _p_ mm_, int_4 _p_ m_, real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ ierr_, real_8 _p_ rv1_, 
  9104  real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_)
  9105  { // ** body not listed **
  9442  }


© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)