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.17
    75  */
    76  #if defined (__GNUC__)
    77  #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
    78  #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
    79  #if (__GNUC__ >= 14)
    80  #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
    81  #pragma GCC diagnostic ignored "-Wimplicit-int"
    82  #pragma GCC diagnostic ignored "-Wint-conversion"
    83  #pragma GCC diagnostic ignored "-Wreturn-mismatch"
    84  #endif
    85  #else
    86  #error VIF requires GCC
    87  #endif
   104  static CALLS __calls[__ncalls] = {
   105    {"bakvec", 0}, // subroutine
   106    {"balanc", 0}, // subroutine
   107    {"balbak", 0}, // subroutine
   108    {"bandr", 0}, // subroutine
   109    {"bandv", 0}, // subroutine
   110    {"bisect", 0}, // subroutine
   111    {"bqr", 0}, // subroutine
   112    {"cbabk2", 0}, // subroutine
   113    {"cbal", 0}, // subroutine
   114    {"cdiv", 0}, // subroutine
   115    {"cg", 0}, // subroutine
   116    {"cglr", 0}, // subroutine
   117    {"cgqr", 0}, // subroutine
   118    {"ch3", 0}, // subroutine
   119    {"ch", 0}, // subroutine
   120    {"cinvit", 0}, // subroutine
   121    {"combak", 0}, // subroutine
   122    {"comhes", 0}, // subroutine
   123    {"comlr2", 0}, // subroutine
   124    {"comlr", 0}, // subroutine
   125    {"comqr2", 0}, // subroutine
   126    {"comqr", 0}, // subroutine
   127    {"cortb", 0}, // subroutine
   128    {"corth", 0}, // subroutine
   129    {"csroot", 0}, // subroutine
   130    {"eisdoc", 0}, // subroutine
   131    {"elmbak", 0}, // subroutine
   132    {"elmhes", 0}, // subroutine
   133    {"eltran", 0}, // subroutine
   134    {"epslon", 0}, // real*8 function
   135    {"figi2", 0}, // subroutine
   136    {"figi", 0}, // subroutine
   137    {"hqr2", 0}, // subroutine
   138    {"hqr", 0}, // subroutine
   139    {"htrib3", 0}, // subroutine
   140    {"htribk", 0}, // subroutine
   141    {"htrid3", 0}, // subroutine
   142    {"htridi", 0}, // subroutine
   143    {"imtql1", 0}, // subroutine
   144    {"imtql2", 0}, // subroutine
   145    {"imtqlv", 0}, // subroutine
   146    {"invit", 0}, // subroutine
   147    {"minfit", 0}, // subroutine
   148    {"ortbak", 0}, // subroutine
   149    {"orthes", 0}, // subroutine
   150    {"ortran", 0}, // subroutine
   151    {"qzhes", 0}, // subroutine
   152    {"qzit", 0}, // subroutine
   153    {"qzval", 0}, // subroutine
   154    {"qzvec", 0}, // subroutine
   155    {"ratqr", 0}, // subroutine
   156    {"rebakb", 0}, // subroutine
   157    {"rebak", 0}, // subroutine
   158    {"reduc2", 0}, // subroutine
   159    {"reduc", 0}, // subroutine
   160    {"rgelm", 0}, // subroutine
   161    {"rg", 0}, // subroutine
   162    {"rgg", 0}, // subroutine
   163    {"rgort", 0}, // subroutine
   164    {"rsb", 0}, // subroutine
   165    {"rs", 0}, // subroutine
   166    {"rsgab", 0}, // subroutine
   167    {"rsgba", 0}, // subroutine
   168    {"rsg", 0}, // subroutine
   169    {"rsm", 0}, // subroutine
   170    {"rsp", 0}, // subroutine
   171    {"rspp", 0}, // subroutine
   172    {"rst", 0}, // subroutine
   173    {"rt", 0}, // subroutine
   174    {"sspev", 0}, // subroutine
   175    {"svd", 0}, // subroutine
   176    {"tinvit", 0}, // subroutine
   177    {"tql1", 0}, // subroutine
   178    {"tql2", 0}, // subroutine
   179    {"tqlrat", 0}, // subroutine
   180    {"trbak1", 0}, // subroutine
   181    {"trbak3", 0}, // subroutine
   182    {"tred1", 0}, // subroutine
   183    {"tred2", 0}, // subroutine
   184    {"tred3", 0}, // subroutine
   185    {"tridib", 0}, // subroutine
   186    {"tsturm", 0}, // subroutine
   187    {NULL, 0}
   188  };
   189  
   191  _p_ ierr_);
   193  ;
   195  real_8 _p_ z_);
   197  _p_ e2_, logical_4 _p_ matz_, real_8 _p_ z_);
   199  real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_);
   201  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 
   202  _p_ rv5_);
   204  _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_);
   206  real_8 _p_ zr_, real_8 _p_ zi_);
   208  real_8 _p_ scale_);
   211  _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_);
   213  int_4 _p_ matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_);
   215  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_);
   217  real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
   219  real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
   221  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_, 
   222  real_8 _p_ rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_);
   224  int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
   226  int_4 _p_ int_);
   228  real_8 _p_ hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_);
   230  real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
   232  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_)
   233  ;
   235  real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
   237  , real_8 _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
   239  real_8 _p_ ortr_, real_8 _p_ orti_);
   243  real_8 _p_ z_);
   246  real_8 _p_ z_);
   249  _p_ ierr_);
   251  _p_ ierr_);
   253  real_8 _p_ wi_, real_8 _p_ z_, int_4 _p_ ierr_);
   255  real_8 _p_ wi_, int_4 _p_ ierr_);
   257  real_8 _p_ zi_);
   259  real_8 _p_ zr_, real_8 _p_ zi_);
   261  real_8 _p_ tau_);
   263  real_8 _p_ e2_, real_8 _p_ tau_);
   267  int_4 _p_ ierr_, real_8 _p_ rv1_);
   269  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_
   270  );
   272  _p_ b_, int_4 _p_ ierr_, real_8 _p_ rv1_);
   274  real_8 _p_ z_);
   277  real_8 _p_ z_);
   280  , real_8 _p_ z_, int_4 _p_ ierr_);
   282  real_8 _p_ beta_, logical_4 _p_ matz_, real_8 _p_ z_);
   284  real_8 _p_ beta_, real_8 _p_ z_);
   286  real_8 _p_ w_, int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_);
   292  real_8 _p_ z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_);
   294  real_4 _p_ z_, int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_);
   296  real_8 _p_ beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_);
   298  real_8 _p_ z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_);
   300  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   302  _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   304  real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   306  real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   308  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   310  _p_ fwork_, int_4 _p_ iwork_, int_4 _p_ ierr_);
   312  _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
   314  _p_ ierr_, int_4 _p_ m_, logical_4 _p_ type_);
   316  _p_ ierr_);
   318  _p_ fv1_, int_4 _p_ ierr_);
   320  int_4 _p_ job_, int_4 _p_ info_);
   322  real_8 _p_ u_, logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_);
   324  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_, 
   325  real_8 _p_ rv4_, real_8 _p_ rv6_);
   335  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 
   336  _p_ rv5_);
   338  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_ 
   339  rv1_, real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_);
   340  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_)
   341  { // ** body not listed **
   375  }
   376  
   377  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_)
   378  { // ** body not listed **
   510  }
   511  
   512  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_ 
   513  z_)
   514  { // ** body not listed **
   552  }
   553  
   554  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_, 
   555  logical_4 _p_ matz_, real_8 _p_ z_)
   556  { // ** body not listed **
   782  }
   783  
   784  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_, 
   785  real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_)
   786  { // ** body not listed **
  1037  }
  1038  
  1039  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_ 
  1040  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_)
  1041  { // ** body not listed **
  1262  }
  1263  
  1264  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_, 
  1265  int_4 _p_ nv_, real_8 _p_ rv_)
  1266  { // ** body not listed **
  1464  }
  1465  
  1466  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_ 
  1467  zr_, real_8 _p_ zi_)
  1468  { // ** body not listed **
  1510  }
  1511  
  1512  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_ 
  1513  scale_)
  1514  { // ** body not listed **
  1654  }
  1655  
  1656  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_)
  1657  { // ** body not listed **
  1669  }
  1670  
  1671  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_
  1672  , real_4 _p_ zr_, real_4 _p_ zi_, real_4 _p_ fv1_, real_4 _p_ fv2_, real_4 _p_ fv3_, int_4 _p_ ierr_)
  1673  { // ** body not listed **
  1697  }
  1698  
  1699  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_ 
  1700  matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_)
  1701  { // ** body not listed **
  1722  }
  1723  
  1724  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_ 
  1725  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_)
  1726  { // ** body not listed **
  1750  }
  1751  
  1752  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_, 
  1753  real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
  1754  { // ** body not listed **
  1783  }
  1784  
  1785  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_, 
  1786  real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
  1787  { // ** body not listed **
  1816  }
  1817  
  1818  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 
  1819  _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_ 
  1820  rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_)
  1821  { // ** body not listed **
  2016  }
  2017  
  2018  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_ 
  2019  m_, real_8 _p_ zr_, real_8 _p_ zi_)
  2020  { // ** body not listed **
  2065  }
  2066  
  2067  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_ 
  2068  int_)
  2069  { // ** body not listed **
  2146  }
  2147  
  2148  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_ 
  2149  hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_)
  2150  { // ** body not listed **
  2459  }
  2460  
  2461  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_ 
  2462  wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
  2463  { // ** body not listed **
  2632  }
  2633  
  2634  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 
  2635  _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_)
  2636  { // ** body not listed **
  2986  }
  2987  
  2988  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_ 
  2989  wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
  2990  { // ** body not listed **
  3171  }
  3172  
  3173  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 
  3174  _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_)
  3175  { // ** body not listed **
  3217  }
  3218  
  3219  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_ 
  3220  ortr_, real_8 _p_ orti_)
  3221  { // ** body not listed **
  3301  }
  3302  
  3303  int_4 _csroot (real_8 _p_ xr_, real_8 _p_ xi_, real_8 _p_ yr_, real_8 _p_ yi_)
  3304  { // ** body not listed **
  3327  }
  3328  
  3329  int_4 _eisdoc (void)
  3330  { // ** body not listed **
  3333  }
  3334  
  3335  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_ 
  3336  z_)
  3337  { // ** body not listed **
  3375  }
  3376  
  3377  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_)
  3378  { // ** body not listed **
  3438  }
  3439  
  3440  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_ 
  3441  z_)
  3442  { // ** body not listed **
  3474  }
  3475  
  3476  real_8 _epslon (real_8 _p_ x_)
  3477  { // ** body not listed **
  3491  }
  3492  
  3493  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_)
  3494  { // ** body not listed **
  3539  }
  3540  
  3541  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_)
  3542  { // ** body not listed **
  3575  }
  3576  
  3577  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_
  3578  , real_8 _p_ z_, int_4 _p_ ierr_)
  3579  { // ** body not listed **
  4024  }
  4025  
  4026  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_, 
  4027  int_4 _p_ ierr_)
  4028  { // ** body not listed **
  4238  }
  4239  
  4240  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_ 
  4241  zi_)
  4242  { // ** body not listed **
  4286  }
  4287  
  4288  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_ 
  4289  zr_, real_8 _p_ zi_)
  4290  { // ** body not listed **
  4334  }
  4335  
  4336  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_ 
  4337  tau_)
  4338  { // ** body not listed **
  4461  }
  4462  
  4463  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_ 
  4464  e2_, real_8 _p_ tau_)
  4465  { // ** body not listed **
  4580  }
  4581  
  4582  int_4 _imtql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
  4583  { // ** body not listed **
  4672  }
  4673  
  4674  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_)
  4675  { // ** body not listed **
  4779  }
  4780  
  4781  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_ 
  4782  ierr_, real_8 _p_ rv1_)
  4783  { // ** body not listed **
  4890  }
  4891  
  4892  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 
  4893  _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_)
  4894  { // ** body not listed **
  5278  }
  5279  
  5280  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_, 
  5281  int_4 _p_ ierr_, real_8 _p_ rv1_)
  5282  { // ** body not listed **
  5552  }
  5553  
  5554  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_ 
  5555  z_)
  5556  { // ** body not listed **
  5591  }
  5592  
  5593  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_)
  5594  { // ** body not listed **
  5650  }
  5651  
  5652  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_ 
  5653  z_)
  5654  { // ** body not listed **
  5691  }
  5692  
  5693  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_)
  5694  { // ** body not listed **
  5820  }
  5821  
  5822  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 
  5823  _p_ z_, int_4 _p_ ierr_)
  5824  { // ** body not listed **
  6161  }
  6162  
  6163  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_ 
  6164  beta_, logical_4 _p_ matz_, real_8 _p_ z_)
  6165  { // ** body not listed **
  6416  }
  6417  
  6418  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_ 
  6419  beta_, real_8 _p_ z_)
  6420  { // ** body not listed **
  6643  }
  6644  
  6645  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_
  6646  , int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_)
  6647  { // ** body not listed **
  6813  }
  6814  
  6815  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_)
  6816  { // ** body not listed **
  6840  }
  6841  
  6842  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_)
  6843  { // ** body not listed **
  6867  }
  6868  
  6869  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_)
  6870  { // ** body not listed **
  6947  }
  6948  
  6949  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_)
  6950  { // ** body not listed **
  7027  }
  7028  
  7029  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_ 
  7030  z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_)
  7031  { // ** body not listed **
  7056  }
  7057  
  7058  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_, 
  7059  int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_)
  7060  { // ** body not listed **
  7085  }
  7086  
  7087  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_ 
  7088  beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_)
  7089  { // ** body not listed **
  7118  }
  7119  
  7120  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_ 
  7121  z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_)
  7122  { // ** body not listed **
  7147  }
  7148  
  7149  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_, 
  7150  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7151  { // ** body not listed **
  7185  }
  7186  
  7187  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_, 
  7188  real_8 _p_ fv2_, int_4 _p_ ierr_)
  7189  { // ** body not listed **
  7208  }
  7209  
  7210  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_, 
  7211  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7212  { // ** body not listed **
  7239  }
  7240  
  7241  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_, 
  7242  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7243  { // ** body not listed **
  7270  }
  7271  
  7272  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_, 
  7273  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7274  { // ** body not listed **
  7301  }
  7302  
  7303  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_, 
  7304  int_4 _p_ iwork_, int_4 _p_ ierr_)
  7305  { // ** body not listed **
  7334  }
  7335  
  7336  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_, 
  7337  real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
  7338  { // ** body not listed **
  7373  }
  7374  
  7375  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_
  7376  , int_4 _p_ m_, logical_4 _p_ type_)
  7377  { // ** body not listed **
  7601  }
  7602  
  7603  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_)
  7604  { // ** body not listed **
  7628  }
  7629  
  7630  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_, 
  7631  int_4 _p_ ierr_)
  7632  { // ** body not listed **
  7657  }
  7658  
  7659  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_ 
  7660  job_, int_4 _p_ info_)
  7661  { // ** body not listed **
  7704  }
  7705  
  7706  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_
  7707  , logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_)
  7708  { // ** body not listed **
  8017  }
  8018  
  8019  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_, 
  8020  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_, 
  8021  real_8 _p_ rv6_)
  8022  { // ** body not listed **
  8213  }
  8214  
  8215  int_4 _tql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
  8216  { // ** body not listed **
  8319  }
  8320  
  8321  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_)
  8322  { // ** body not listed **
  8440  }
  8441  
  8442  int_4 _tqlrat (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e2_, int_4 _p_ ierr_)
  8443  { // ** body not listed **
  8559  }
  8560  
  8561  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_)
  8562  { // ** body not listed **
  8591  }
  8592  
  8593  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_)
  8594  { // ** body not listed **
  8630  }
  8631  
  8632  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_)
  8633  { // ** body not listed **
  8721  }
  8722  
  8723  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_)
  8724  { // ** body not listed **
  8841  }
  8842  
  8843  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_)
  8844  { // ** body not listed **
  8924  }
  8925  
  8926  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_ 
  8927  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_)
  8928  { // ** body not listed **
  9201  }
  9202  
  9203  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_ 
  9204  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_, 
  9205  real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_)
  9206  { // ** body not listed **
  9543  }
  9544  


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