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