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)
|