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  //! 82 subprograms from EISPACK
    24  
    25  // Compiled from Fortran source code by VIF.
    26  // Selected subprograms are:
    27  //
    28  // BAKVEC CDIV   COMLR  ELMBAK HTRIB3 MINFIT RATQR  RGORT  RSPP   TQLRAT
    29  // BALANC CG     COMLR2 ELMHES HTRIBK ORTBAK REBAK  RS     RST    TRBAK1
    30  // BALBAK CGLR   COMQR  ELTRAN HTRID3 ORTHES REBAKB RSB    RT     TRBAK3
    31  // BANDR  CGQR   COMQR2 EPSLON HTRIDI ORTRAN REDUC  RSG    SSPEV  TRED1
    32  // BANDV  CH     CORTB  FIGI   IMTQL1 QZHES  REDUC2 RSGAB  SVD    TRED2
    33  // BISECT CH3    CORTH  FIGI2  IMTQL2 QZIT   RG     RSGBA  TINVIT TRED3
    34  // BQR    CINVIT CSROOT HQR    IMTQLV QZVAL  RGELM  RSM    TQL1   TRIDIB
    35  // CBABK2 COMBAK EISDOC HQR2   INVIT  QZVEC  RGG    RSP    TQL2   TSTURM
    36  // CBAL   COMHES
    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.2.12
    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    {"eisdoc", 0}, // subroutine
   119    {"elmbak", 0}, // subroutine
   120    {"elmhes", 0}, // subroutine
   121    {"eltran", 0}, // subroutine
   122    {"epslon", 0}, // real*8 function
   123    {"figi2", 0}, // subroutine
   124    {"figi", 0}, // subroutine
   125    {"hqr2", 0}, // subroutine
   126    {"hqr", 0}, // subroutine
   127    {"htrib3", 0}, // subroutine
   128    {"htribk", 0}, // subroutine
   129    {"htrid3", 0}, // subroutine
   130    {"htridi", 0}, // subroutine
   131    {"imtql1", 0}, // subroutine
   132    {"imtql2", 0}, // subroutine
   133    {"imtqlv", 0}, // subroutine
   134    {"invit", 0}, // subroutine
   135    {"minfit", 0}, // subroutine
   136    {"ortbak", 0}, // subroutine
   137    {"orthes", 0}, // subroutine
   138    {"ortran", 0}, // subroutine
   139    {"qzhes", 0}, // subroutine
   140    {"qzit", 0}, // subroutine
   141    {"qzval", 0}, // subroutine
   142    {"qzvec", 0}, // subroutine
   143    {"ratqr", 0}, // subroutine
   144    {"rebakb", 0}, // subroutine
   145    {"rebak", 0}, // subroutine
   146    {"reduc2", 0}, // subroutine
   147    {"reduc", 0}, // subroutine
   148    {"rgelm", 0}, // subroutine
   149    {"rg", 0}, // subroutine
   150    {"rgg", 0}, // subroutine
   151    {"rgort", 0}, // subroutine
   152    {"rsb", 0}, // subroutine
   153    {"rs", 0}, // subroutine
   154    {"rsgab", 0}, // subroutine
   155    {"rsgba", 0}, // subroutine
   156    {"rsg", 0}, // subroutine
   157    {"rsm", 0}, // subroutine
   158    {"rsp", 0}, // subroutine
   159    {"rspp", 0}, // subroutine
   160    {"rst", 0}, // subroutine
   161    {"rt", 0}, // subroutine
   162    {"sspev", 0}, // subroutine
   163    {"svd", 0}, // subroutine
   164    {"tinvit", 0}, // subroutine
   165    {"tql1", 0}, // subroutine
   166    {"tql2", 0}, // subroutine
   167    {"tqlrat", 0}, // subroutine
   168    {"trbak1", 0}, // subroutine
   169    {"trbak3", 0}, // subroutine
   170    {"tred1", 0}, // subroutine
   171    {"tred2", 0}, // subroutine
   172    {"tred3", 0}, // subroutine
   173    {"tridib", 0}, // subroutine
   174    {"tsturm", 0}, // subroutine
   175    {NULL, 0}
   176  };
   177  
   179  _p_ ierr_);
   181  ;
   183  real_8 _p_ z_);
   185  _p_ e2_, logical_4 _p_ matz_, real_8 _p_ z_);
   187  real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_);
   189  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 
   190  _p_ rv5_);
   192  _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_);
   194  real_8 _p_ zr_, real_8 _p_ zi_);
   196  real_8 _p_ scale_);
   199  _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_);
   201  int_4 _p_ matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_);
   203  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_);
   205  real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
   207  real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
   209  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_, 
   210  real_8 _p_ rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_);
   212  int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
   214  int_4 _p_ int_);
   216  real_8 _p_ hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_);
   218  real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
   220  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_)
   221  ;
   223  real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
   225  , real_8 _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
   227  real_8 _p_ ortr_, real_8 _p_ orti_);
   231  real_8 _p_ z_);
   234  real_8 _p_ z_);
   237  _p_ ierr_);
   239  _p_ ierr_);
   241  real_8 _p_ wi_, real_8 _p_ z_, int_4 _p_ ierr_);
   243  real_8 _p_ wi_, int_4 _p_ ierr_);
   245  real_8 _p_ zi_);
   247  real_8 _p_ zr_, real_8 _p_ zi_);
   249  real_8 _p_ tau_);
   251  real_8 _p_ e2_, real_8 _p_ tau_);
   255  int_4 _p_ ierr_, real_8 _p_ rv1_);
   257  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_
   258  );
   260  _p_ b_, int_4 _p_ ierr_, real_8 _p_ rv1_);
   262  real_8 _p_ z_);
   265  real_8 _p_ z_);
   268  , real_8 _p_ z_, int_4 _p_ ierr_);
   270  real_8 _p_ beta_, logical_4 _p_ matz_, real_8 _p_ z_);
   272  real_8 _p_ beta_, real_8 _p_ z_);
   274  real_8 _p_ w_, int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_);
   280  real_8 _p_ z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_);
   282  real_4 _p_ z_, int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_);
   284  real_8 _p_ beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_);
   286  real_8 _p_ z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_);
   288  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   290  _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   292  real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   294  real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   296  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   298  _p_ fwork_, int_4 _p_ iwork_, int_4 _p_ ierr_);
   300  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   302  _p_ ierr_, int_4 _p_ m_, logical_4 _p_ type_);
   304  _p_ ierr_);
   306  _p_ fv1_, int_4 _p_ ierr_);
   308  int_4 _p_ job_, int_4 _p_ info_);
   310  real_8 _p_ u_, logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_);
   312  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_, 
   313  real_8 _p_ rv4_, real_8 _p_ rv6_);
   323  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 
   324  _p_ rv5_);
   326  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_ 
   327  rv1_, real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_);
   328  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_)
   329  { // ** body not listed **
   363  }
   364  
   365  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_)
   366  { // ** body not listed **
   498  }
   499  
   500  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_ 
   501  z_)
   502  { // ** body not listed **
   540  }
   541  
   542  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_, 
   543  logical_4 _p_ matz_, real_8 _p_ z_)
   544  { // ** body not listed **
   770  }
   771  
   772  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_, 
   773  real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_)
   774  { // ** body not listed **
  1025  }
  1026  
  1027  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_ 
  1028  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_)
  1029  { // ** body not listed **
  1250  }
  1251  
  1252  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_, 
  1253  int_4 _p_ nv_, real_8 _p_ rv_)
  1254  { // ** body not listed **
  1452  }
  1453  
  1454  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_ 
  1455  zr_, real_8 _p_ zi_)
  1456  { // ** body not listed **
  1498  }
  1499  
  1500  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_ 
  1501  scale_)
  1502  { // ** body not listed **
  1642  }
  1643  
  1644  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_)
  1645  { // ** body not listed **
  1657  }
  1658  
  1659  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_
  1660  , real_4 _p_ zr_, real_4 _p_ zi_, real_4 _p_ fv1_, real_4 _p_ fv2_, real_4 _p_ fv3_, int_4 _p_ ierr_)
  1661  { // ** body not listed **
  1685  }
  1686  
  1687  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_ 
  1688  matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_)
  1689  { // ** body not listed **
  1710  }
  1711  
  1712  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_ 
  1713  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_)
  1714  { // ** body not listed **
  1738  }
  1739  
  1740  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_, 
  1741  real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
  1742  { // ** body not listed **
  1771  }
  1772  
  1773  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_, 
  1774  real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
  1775  { // ** body not listed **
  1804  }
  1805  
  1806  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 
  1807  _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_ 
  1808  rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_)
  1809  { // ** body not listed **
  2004  }
  2005  
  2006  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_ 
  2007  m_, real_8 _p_ zr_, real_8 _p_ zi_)
  2008  { // ** body not listed **
  2053  }
  2054  
  2055  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_ 
  2056  int_)
  2057  { // ** body not listed **
  2134  }
  2135  
  2136  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_ 
  2137  hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_)
  2138  { // ** body not listed **
  2447  }
  2448  
  2449  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_ 
  2450  wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
  2451  { // ** body not listed **
  2620  }
  2621  
  2622  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 
  2623  _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_)
  2624  { // ** body not listed **
  2974  }
  2975  
  2976  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_ 
  2977  wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
  2978  { // ** body not listed **
  3159  }
  3160  
  3161  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 
  3162  _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_)
  3163  { // ** body not listed **
  3205  }
  3206  
  3207  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_ 
  3208  ortr_, real_8 _p_ orti_)
  3209  { // ** body not listed **
  3289  }
  3290  
  3291  int_4 _csroot (real_8 _p_ xr_, real_8 _p_ xi_, real_8 _p_ yr_, real_8 _p_ yi_)
  3292  { // ** body not listed **
  3315  }
  3316  
  3317  int_4 _eisdoc (void)
  3318  { // ** body not listed **
  3321  }
  3322  
  3323  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_ 
  3324  z_)
  3325  { // ** body not listed **
  3363  }
  3364  
  3365  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_)
  3366  { // ** body not listed **
  3426  }
  3427  
  3428  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_ 
  3429  z_)
  3430  { // ** body not listed **
  3462  }
  3463  
  3464  real_8 _epslon (real_8 _p_ x_)
  3465  { // ** body not listed **
  3479  }
  3480  
  3481  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_)
  3482  { // ** body not listed **
  3527  }
  3528  
  3529  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_)
  3530  { // ** body not listed **
  3563  }
  3564  
  3565  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_
  3566  , real_8 _p_ z_, int_4 _p_ ierr_)
  3567  { // ** body not listed **
  4012  }
  4013  
  4014  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_, 
  4015  int_4 _p_ ierr_)
  4016  { // ** body not listed **
  4226  }
  4227  
  4228  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_ 
  4229  zi_)
  4230  { // ** body not listed **
  4274  }
  4275  
  4276  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_ 
  4277  zr_, real_8 _p_ zi_)
  4278  { // ** body not listed **
  4322  }
  4323  
  4324  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_ 
  4325  tau_)
  4326  { // ** body not listed **
  4449  }
  4450  
  4451  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_ 
  4452  e2_, real_8 _p_ tau_)
  4453  { // ** body not listed **
  4568  }
  4569  
  4570  int_4 _imtql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
  4571  { // ** body not listed **
  4660  }
  4661  
  4662  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_)
  4663  { // ** body not listed **
  4767  }
  4768  
  4769  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_ 
  4770  ierr_, real_8 _p_ rv1_)
  4771  { // ** body not listed **
  4878  }
  4879  
  4880  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 
  4881  _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_)
  4882  { // ** body not listed **
  5266  }
  5267  
  5268  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_, 
  5269  int_4 _p_ ierr_, real_8 _p_ rv1_)
  5270  { // ** body not listed **
  5540  }
  5541  
  5542  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_ 
  5543  z_)
  5544  { // ** body not listed **
  5579  }
  5580  
  5581  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_)
  5582  { // ** body not listed **
  5638  }
  5639  
  5640  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_ 
  5641  z_)
  5642  { // ** body not listed **
  5679  }
  5680  
  5681  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_)
  5682  { // ** body not listed **
  5808  }
  5809  
  5810  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 
  5811  _p_ z_, int_4 _p_ ierr_)
  5812  { // ** body not listed **
  6149  }
  6150  
  6151  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_ 
  6152  beta_, logical_4 _p_ matz_, real_8 _p_ z_)
  6153  { // ** body not listed **
  6404  }
  6405  
  6406  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_ 
  6407  beta_, real_8 _p_ z_)
  6408  { // ** body not listed **
  6631  }
  6632  
  6633  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_
  6634  , int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_)
  6635  { // ** body not listed **
  6801  }
  6802  
  6803  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_)
  6804  { // ** body not listed **
  6828  }
  6829  
  6830  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_)
  6831  { // ** body not listed **
  6855  }
  6856  
  6857  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_)
  6858  { // ** body not listed **
  6935  }
  6936  
  6937  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_)
  6938  { // ** body not listed **
  7015  }
  7016  
  7017  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_ 
  7018  z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_)
  7019  { // ** body not listed **
  7044  }
  7045  
  7046  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_, 
  7047  int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_)
  7048  { // ** body not listed **
  7073  }
  7074  
  7075  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_ 
  7076  beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_)
  7077  { // ** body not listed **
  7106  }
  7107  
  7108  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_ 
  7109  z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_)
  7110  { // ** body not listed **
  7135  }
  7136  
  7137  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_, 
  7138  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7139  { // ** body not listed **
  7173  }
  7174  
  7175  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_, 
  7176  real_8 _p_ fv2_, int_4 _p_ ierr_)
  7177  { // ** body not listed **
  7196  }
  7197  
  7198  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_, 
  7199  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7200  { // ** body not listed **
  7227  }
  7228  
  7229  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_, 
  7230  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7231  { // ** body not listed **
  7258  }
  7259  
  7260  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_, 
  7261  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7262  { // ** body not listed **
  7289  }
  7290  
  7291  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_, 
  7292  int_4 _p_ iwork_, int_4 _p_ ierr_)
  7293  { // ** body not listed **
  7322  }
  7323  
  7324  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_, 
  7325  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7326  { // ** body not listed **
  7361  }
  7362  
  7363  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_
  7364  , int_4 _p_ m_, logical_4 _p_ type_)
  7365  { // ** body not listed **
  7589  }
  7590  
  7591  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_)
  7592  { // ** body not listed **
  7616  }
  7617  
  7618  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_, 
  7619  int_4 _p_ ierr_)
  7620  { // ** body not listed **
  7645  }
  7646  
  7647  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_ 
  7648  job_, int_4 _p_ info_)
  7649  { // ** body not listed **
  7692  }
  7693  
  7694  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_
  7695  , logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_)
  7696  { // ** body not listed **
  8005  }
  8006  
  8007  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_, 
  8008  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_, 
  8009  real_8 _p_ rv6_)
  8010  { // ** body not listed **
  8201  }
  8202  
  8203  int_4 _tql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
  8204  { // ** body not listed **
  8307  }
  8308  
  8309  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_)
  8310  { // ** body not listed **
  8428  }
  8429  
  8430  int_4 _tqlrat (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e2_, int_4 _p_ ierr_)
  8431  { // ** body not listed **
  8547  }
  8548  
  8549  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_)
  8550  { // ** body not listed **
  8579  }
  8580  
  8581  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_)
  8582  { // ** body not listed **
  8618  }
  8619  
  8620  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_)
  8621  { // ** body not listed **
  8709  }
  8710  
  8711  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_)
  8712  { // ** body not listed **
  8829  }
  8830  
  8831  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_)
  8832  { // ** body not listed **
  8912  }
  8913  
  8914  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_ 
  8915  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_)
  8916  { // ** body not listed **
  9189  }
  9190  
  9191  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_ 
  9192  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_, 
  9193  real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_)
  9194  { // ** body not listed **
  9531  }
  9532  


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