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