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.10
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 };
178 _p_ ierr_);
180 ;
182 real_8 _p_ z_);
184 _p_ e2_, logical_4 _p_ matz_, real_8 _p_ z_);
186 real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_);
188 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
189 _p_ rv5_);
191 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_);
193 real_8 _p_ zr_, real_8 _p_ zi_);
195 real_8 _p_ scale_);
198 _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_);
200 int_4 _p_ matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_);
202 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_);
204 real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
206 real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_);
208 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_,
209 real_8 _p_ rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_);
211 int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
213 int_4 _p_ int_);
215 real_8 _p_ hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_);
217 real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
219 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_)
220 ;
222 real_8 _p_ wr_, real_8 _p_ wi_, int_4 _p_ ierr_);
224 , real_8 _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_);
226 real_8 _p_ ortr_, real_8 _p_ orti_);
230 real_8 _p_ z_);
233 real_8 _p_ z_);
236 _p_ ierr_);
238 _p_ ierr_);
240 real_8 _p_ wi_, real_8 _p_ z_, int_4 _p_ ierr_);
242 real_8 _p_ wi_, int_4 _p_ ierr_);
244 real_8 _p_ zi_);
246 real_8 _p_ zr_, real_8 _p_ zi_);
248 real_8 _p_ tau_);
250 real_8 _p_ e2_, real_8 _p_ tau_);
254 int_4 _p_ ierr_, real_8 _p_ rv1_);
256 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_
257 );
259 _p_ b_, int_4 _p_ ierr_, real_8 _p_ rv1_);
261 real_8 _p_ z_);
264 real_8 _p_ z_);
267 , real_8 _p_ z_, int_4 _p_ ierr_);
269 real_8 _p_ beta_, logical_4 _p_ matz_, real_8 _p_ z_);
271 real_8 _p_ beta_, real_8 _p_ z_);
273 real_8 _p_ w_, int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_);
279 real_8 _p_ z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_);
281 real_4 _p_ z_, int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_);
283 real_8 _p_ beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_);
285 real_8 _p_ z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_);
287 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
289 _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 real_8 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
295 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
297 _p_ fwork_, int_4 _p_ iwork_, int_4 _p_ ierr_);
299 _p_ z_, real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_);
301 _p_ ierr_, int_4 _p_ m_, logical_4 _p_ type_);
303 _p_ ierr_);
305 _p_ fv1_, int_4 _p_ ierr_);
307 int_4 _p_ job_, int_4 _p_ info_);
309 real_8 _p_ u_, logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_);
311 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_,
312 real_8 _p_ rv4_, real_8 _p_ rv6_);
322 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
323 _p_ rv5_);
325 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_
326 rv1_, real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_);
327 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_)
328 { // ** body not listed **
362 }
363 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_)
364 { // ** body not listed **
496 }
497 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_
498 z_)
499 { // ** body not listed **
537 }
538 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_,
539 logical_4 _p_ matz_, real_8 _p_ z_)
540 { // ** body not listed **
766 }
767 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_,
768 real_8 _p_ z_, int_4 _p_ ierr_, int_4 _p_ nv_, real_8 _p_ rv_, real_8 _p_ rv6_)
769 { // ** body not listed **
1020 }
1021 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_
1022 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_)
1023 { // ** body not listed **
1244 }
1245 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_,
1246 int_4 _p_ nv_, real_8 _p_ rv_)
1247 { // ** body not listed **
1445 }
1446 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_
1447 zr_, real_8 _p_ zi_)
1448 { // ** body not listed **
1490 }
1491 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_
1492 scale_)
1493 { // ** body not listed **
1633 }
1634 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_)
1635 { // ** body not listed **
1647 }
1648 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_
1649 , real_4 _p_ zr_, real_4 _p_ zi_, real_4 _p_ fv1_, real_4 _p_ fv2_, real_4 _p_ fv3_, int_4 _p_ ierr_)
1650 { // ** body not listed **
1674 }
1675 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_
1676 matz_, real_8 _p_ zr_, real_8 _p_ zi_, real_8 _p_ fv1_, int_4 _p_ iv1_, int_4 _p_ ierr_)
1677 { // ** body not listed **
1698 }
1699 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_
1700 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_)
1701 { // ** body not listed **
1725 }
1726 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_,
1727 real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
1728 { // ** body not listed **
1757 }
1758 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_,
1759 real_8 _p_ zi_, real_8 _p_ fv1_, real_8 _p_ fv2_, real_8 _p_ fm1_, int_4 _p_ ierr_)
1760 { // ** body not listed **
1789 }
1790 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
1791 _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_
1792 rm2_, real_8 _p_ rv1_, real_8 _p_ rv2_)
1793 { // ** body not listed **
1988 }
1989 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_
1990 m_, real_8 _p_ zr_, real_8 _p_ zi_)
1991 { // ** body not listed **
2036 }
2037 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_
2038 int_)
2039 { // ** body not listed **
2116 }
2117 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_
2118 hi_, real_8 _p_ wr_, real_8 _p_ wi_, real_8 _p_ zr_, real_8 _p_ zi_, int_4 _p_ ierr_)
2119 { // ** body not listed **
2428 }
2429 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_
2430 wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
2431 { // ** body not listed **
2600 }
2601 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
2602 _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_)
2603 { // ** body not listed **
2953 }
2954 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_
2955 wr_, real_8 _p_ wi_, int_4 _p_ ierr_)
2956 { // ** body not listed **
3137 }
3138 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
3139 _p_ orti_, int_4 _p_ m_, real_8 _p_ zr_, real_8 _p_ zi_)
3140 { // ** body not listed **
3182 }
3183 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_
3184 ortr_, real_8 _p_ orti_)
3185 { // ** body not listed **
3265 }
3266 int_4 _csroot (real_8 _p_ xr_, real_8 _p_ xi_, real_8 _p_ yr_, real_8 _p_ yi_)
3267 { // ** body not listed **
3290 }
3291 int_4 _eisdoc (void)
3292 { // ** body not listed **
3295 }
3296 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_
3297 z_)
3298 { // ** body not listed **
3336 }
3337 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_)
3338 { // ** body not listed **
3398 }
3399 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_
3400 z_)
3401 { // ** body not listed **
3433 }
3434 real_8 _epslon (real_8 _p_ x_)
3435 { // ** body not listed **
3449 }
3450 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_)
3451 { // ** body not listed **
3496 }
3497 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_)
3498 { // ** body not listed **
3531 }
3532 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_
3533 , real_8 _p_ z_, int_4 _p_ ierr_)
3534 { // ** body not listed **
3979 }
3980 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_,
3981 int_4 _p_ ierr_)
3982 { // ** body not listed **
4192 }
4193 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_
4194 zi_)
4195 { // ** body not listed **
4239 }
4240 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_
4241 zr_, real_8 _p_ zi_)
4242 { // ** body not listed **
4286 }
4287 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_
4288 tau_)
4289 { // ** body not listed **
4412 }
4413 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_
4414 e2_, real_8 _p_ tau_)
4415 { // ** body not listed **
4530 }
4531 int_4 _imtql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
4532 { // ** body not listed **
4621 }
4622 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_)
4623 { // ** body not listed **
4727 }
4728 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_
4729 ierr_, real_8 _p_ rv1_)
4730 { // ** body not listed **
4837 }
4838 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
4839 _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_)
4840 { // ** body not listed **
5224 }
5225 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_,
5226 int_4 _p_ ierr_, real_8 _p_ rv1_)
5227 { // ** body not listed **
5497 }
5498 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_
5499 z_)
5500 { // ** body not listed **
5535 }
5536 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_)
5537 { // ** body not listed **
5593 }
5594 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_
5595 z_)
5596 { // ** body not listed **
5633 }
5634 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_)
5635 { // ** body not listed **
5761 }
5762 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
5763 _p_ z_, int_4 _p_ ierr_)
5764 { // ** body not listed **
6101 }
6102 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_
6103 beta_, logical_4 _p_ matz_, real_8 _p_ z_)
6104 { // ** body not listed **
6355 }
6356 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_
6357 beta_, real_8 _p_ z_)
6358 { // ** body not listed **
6581 }
6582 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_
6583 , int_4 _p_ ind_, real_8 _p_ bd_, logical_4 _p_ type_, int_4 _p_ idef_, int_4 _p_ ierr_)
6584 { // ** body not listed **
6750 }
6751 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_)
6752 { // ** body not listed **
6776 }
6777 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_)
6778 { // ** body not listed **
6802 }
6803 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_)
6804 { // ** body not listed **
6881 }
6882 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_)
6883 { // ** body not listed **
6960 }
6961 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_
6962 z_, int_4 _p_ iv1_, real_8 _p_ fv1_, int_4 _p_ ierr_)
6963 { // ** body not listed **
6988 }
6989 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_,
6990 int_4 _p_ iv1_, real_4 _p_ fv1_, int_4 _p_ ierr_)
6991 { // ** body not listed **
7016 }
7017 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_
7018 beta_, int_4 _p_ matz_, real_8 _p_ z_, int_4 _p_ ierr_)
7019 { // ** body not listed **
7048 }
7049 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_
7050 z_, real_8 _p_ ort_, real_8 _p_ fv1_, int_4 _p_ ierr_)
7051 { // ** body not listed **
7076 }
7077 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_,
7078 real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
7079 { // ** body not listed **
7113 }
7114 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_,
7115 real_8 _p_ fv2_, int_4 _p_ ierr_)
7116 { // ** body not listed **
7135 }
7136 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_,
7137 real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
7138 { // ** body not listed **
7165 }
7166 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_,
7167 real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
7168 { // ** body not listed **
7195 }
7196 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_,
7197 real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
7198 { // ** body not listed **
7225 }
7226 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_,
7227 int_4 _p_ iwork_, int_4 _p_ ierr_)
7228 { // ** body not listed **
7257 }
7258 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_,
7259 real_8 _p_ fv1_, real_8 _p_ fv2_, int_4 _p_ ierr_)
7260 { // ** body not listed **
7295 }
7296 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_
7297 , int_4 _p_ m_, logical_4 _p_ type_)
7298 { // ** body not listed **
7522 }
7523 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_)
7524 { // ** body not listed **
7548 }
7549 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_,
7550 int_4 _p_ ierr_)
7551 { // ** body not listed **
7576 }
7577 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_
7578 job_, int_4 _p_ info_)
7579 { // ** body not listed **
7622 }
7623 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_
7624 , logical_4 _p_ matv_, real_8 _p_ v_, int_4 _p_ ierr_, real_8 _p_ rv1_)
7625 { // ** body not listed **
7934 }
7935 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_,
7936 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_,
7937 real_8 _p_ rv6_)
7938 { // ** body not listed **
8129 }
8130 int_4 _tql1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, int_4 _p_ ierr_)
8131 { // ** body not listed **
8234 }
8235 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_)
8236 { // ** body not listed **
8354 }
8355 int_4 _tqlrat (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e2_, int_4 _p_ ierr_)
8356 { // ** body not listed **
8472 }
8473 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_)
8474 { // ** body not listed **
8503 }
8504 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_)
8505 { // ** body not listed **
8541 }
8542 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_)
8543 { // ** body not listed **
8631 }
8632 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_)
8633 { // ** body not listed **
8750 }
8751 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_)
8752 { // ** body not listed **
8832 }
8833 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_
8834 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_)
8835 { // ** body not listed **
9108 }
9109 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_
9110 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_,
9111 real_8 _p_ rv2_, real_8 _p_ rv3_, real_8 _p_ rv4_, real_8 _p_ rv5_, real_8 _p_ rv6_)
9112 { // ** body not listed **
9449 }
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
