mathlib-toeplitz.c
1 //! @file mathlib-toeplitz.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 //! Subprograms from TOEPLITZ
24
25 // Compiled from Fortran source code by VIF.
26
27 // TOEPLITZ solves a variety of Toeplitz and circulant linear systems. The package can also handle
28 // circulant Toeplitz systems, and some other more complicated but related forms. The package was
29 // written in the early 1980's by a joint working group of American and Soviet mathematicians.
30 //
31 // The package lives on for instance as a Python package.
32 //
33 // Reference:
34 // Oleg Arushanian, MK Samarin, Valentin Voevodin, Evgeny Tyrtyshnikov,
35 // Burton Garbow, James Boyle, Wayne Cowell, Kenneth Dritz.
36 // The TOEPLITZ Package User's Guide.
37 // Argonne National Laboratory.
38 // ANL-83-16, 1983.
39
40 // The license for TOEPLITZ Fortran source code is:
41 //
42 // TOEPLITZ can be obtained from:
43 // https://people.sc.fsu.edu/~jburkardt/f77_src/toeplitz/toeplitz.html
44 // distributed under MIT license.
45 //
46 // The original code can be obtained from:
47 // netlib.org/toeplitz
48 //
49 //
50 // MIT License
51 //
52 // Permission is hereby granted, free of charge, to any person obtaining a copy
53 // of this software and associated documentation files (the "Software"), to deal
54 // in the Software without restriction, including without limitation the rights
55 // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
56 // copies of the Software, and to permit persons to whom the Software is
57 // furnished to do so, subject to the following conditions:
58 //
59 // The above copyright notice and this permission notice shall be included in all
60 // copies or substantial portions of the Software.
61 //
62 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
63 // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
64 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
65 // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
66 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
67 // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
68 // SOFTWARE.
69 //
70 //
71 // Reference:
72 //
73 // Oleg Arushanian, MK Samarin, Valentin Voevodin, Evgeny Tyrtyshnikov,
74 // Burton Garbow, James Boyle, Wayne Cowell, Kenneth Dritz.
75 // The TOEPLITZ Package User's Guide.
76 // Argonne National Laboratory.
77 // ANL-83-16, 1983.
78
79 // VIF THU 20 NOV 2025 22:22:03 ** MARCEL ** TOEPLITZ PAGE 00001
80
81
82
83 // M M A RRRRRR CCCCC EEEEEEE L
84 // MM MM A A R R C C E L
85 // M M M M A A R R C E L
86 // M M M A A RRRRRR C EEEEE L
87 // M M AAAAAAA R R C E L
88 // M M A A R R C C E L
89 // M M A A R R CCCCC EEEEEEE LLLLLLL
90
91 // TTTTTTT OOOOO EEEEEEE PPPPPP L III TTTTTTT ZZZZZZZ
92 // T O O E P P L I T Z
93 // T O O E P P L I T Z
94 // T O O EEEEE PPPPPP L I T Z
95 // T O O E P L I T Z
96 // T O O E P L I T *
97 // T OOOOO EEEEEEE P LLLLLLL III T ZZZZZZZ
98
99
100
101
102 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** DIAGNOSTICS PAGE 00002
103
104 // ** linker ** object size 101184 bytes
105 // ** statistics ** 19 subprograms, no errors, no warnings
106
107
108 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** JOB CONTROL PAGE 00003
109
110
111 // Line JCL *...*....1....|....2....|....3....|....4....|....5....|....6....|....7..*.|....8
112
113
114
115 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00004
116
117 /*
118 Generated by VIF - experimental VIntage Fortran compiler.
119 VIF release 1.4.1
120 */
121
122 #if defined (__GNUC__)
123 #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
124 #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
125 #if (__GNUC__ >= 14)
126 #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
127 #pragma GCC diagnostic ignored "-Wimplicit-int"
128 #pragma GCC diagnostic ignored "-Wint-conversion"
129 #pragma GCC diagnostic ignored "-Wreturn-mismatch"
130 #endif
131 #else
132 #error VIF requires GCC
133 #endif
134
135 #include <vif.h>
136
137 static int_4 _km1 = -1, _k0 = 0, _k1 = 1;
138
139
140
141
142
143 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** TYPEDEFS PAGE 00005
144
145 static FORMAT *__fmt_a = NULL;
146
147
148 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** FREQUENCY TABLE PAGE 00006
149
150 #define __ncalls 20
151 static CALLS __calls[__ncalls] = {
152 { // ** body not listed **"cccslz", 0}, // subroutine
154 { // ** body not listed **"ccslz", 0}, // subroutine
156 { // ** body not listed **"cgslz", 0}, // subroutine
158 { // ** body not listed **"cqrz", 0}, // subroutine
160 { // ** body not listed **"ctgslz", 0}, // subroutine
162 { // ** body not listed **"salwz", 0}, // subroutine
164 { // ** body not listed **"tgsld", 0}, // subroutine
166 { // ** body not listed **"tgslz", 0}, // subroutine
168 { // ** body not listed **"tsld", 0}, // subroutine
170 { // ** body not listed **"tslz", 0}, // subroutine
172 };
173
174
175
176 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** COMMON BLOCKS PAGE 00007
177
178
179
180 // VIF THU 20 NOV 2025 22:22:03 ** GLOBAL SCOPE ** PROTOTYPES PAGE 00008
181
182
183 prototype int_4 _cccslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
184 k_, int_4 _p_ lda_);
185 prototype int_4 _ccgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
186 k_, int_4 _p_ lda_);
187 prototype int_4 _ccslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
188 lda_);
189 prototype int_4 _cctslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
190 k_, int_4 _p_ lda_);
191 prototype int_4 _cgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
192 lda_);
193 prototype int_4 _cqrd (real_8 _p_ a_, real_8 _p_ q_, real_8 _p_ s_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ ldq_, int_4
194 _p_ lds_);
195 prototype int_4 _cqrz (complex_16 _p_ a_, complex_16 _p_ q_, complex_16 _p_ s_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
196 ldq_, int_4 _p_ lds_);
197 prototype int_4 _cslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_);
198 prototype int_4 _ctgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
199 k_, int_4 _p_ lda_);
200 prototype int_4 _ctslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
201 lda_);
202 prototype int_4 _salwz (complex_16 _p_ a_, complex_16 _p_ r1_, complex_16 _p_ r2_, int_4 _p_ m_, int_4 _p_ l_, int_4
203 _p_ lda_, int_4 _p_ job_);
204 prototype int_4 _tgsld1 (real_8 _p_ a1_, real_8 _p_ a2_, real_8 _p_ b_, real_8 _p_ x_, real_8 _p_ c1_, real_8 _p_ c2_,
205 real_8 _p_ r1_, real_8 _p_ r2_, real_8 _p_ r3_, real_8 _p_ r5_, real_8 _p_ r6_, real_8 _p_ r_, int_4 _p_ m_, int_4 _p_
206 l_, int_4 _p_ lda_);
207 prototype int_4 _tgsld (real_8 _p_ a_, real_8 _p_ x_, real_8 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_);
208 prototype int_4 _tgslz1 (complex_16 _p_ a1_, complex_16 _p_ a2_, complex_16 _p_ b_, complex_16 _p_ x_, complex_16 _p_
209 c1_, complex_16 _p_ c2_, complex_16 _p_ r1_, complex_16 _p_ r2_, complex_16 _p_ r3_, complex_16 _p_ r5_, complex_16 _p_
210 r6_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_);
211 prototype int_4 _tgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_
212 lda_);
213 prototype int_4 _tsld1 (real_8 _p_ a1_, real_8 _p_ a2_, real_8 _p_ b_, real_8 _p_ x_, real_8 _p_ c1_, real_8 _p_ c2_,
214 int_4 _p_ m_);
215 prototype int_4 _tsld (real_8 _p_ a_, real_8 _p_ x_, real_8 _p_ r_, int_4 _p_ m_);
216 prototype int_4 _tslz1 (complex_16 _p_ a1_, complex_16 _p_ a2_, complex_16 _p_ b_, complex_16 _p_ x_, complex_16 _p_
217 c1_, complex_16 _p_ c2_, int_4 _p_ m_);
218 prototype int_4 _tslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_);
219
220
221
222
223
224
225
226
227
228
229
230
231
232 // VIF THU 20 NOV 2025 22:22:03 ** CCCSLZ ** CCCSLZ PAGE 00009
233
234 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
235 // S S U U B B R R O O U U T I NN N E
236 // S U U B B R R O O U U T I N N N E
237 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
238 // S U U B B R R O O U U T I N N N E
239 // S S U U B B R R O O U U T I N NN E
240 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
241
242 // CCCCC CCCCC CCCCC SSSSS L ZZZZZZZ
243 // C C C C C C S S L Z
244 // C C C S L Z
245 // C C C SSSSS L Z
246 // C C C S L Z
247 // C C C C C C S S L *
248 // CCCCC CCCCC CCCCC SSSSS LLLLLLL ZZZZZZZ
249
250
251 // VIF THU 20 NOV 2025 22:22:03 ** CCCSLZ ** DIAGNOSTICS PAGE 00010
252
253
254 // VIF THU 20 NOV 2025 22:22:03 ** CCCSLZ ** GENERATED CODE PAGE 00011
255
256 int_4 _cccslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ k_, int_4
257 _p_ lda_)
258 { // ** body not listed **
279 }
280
281 // VIF THU 20 NOV 2025 22:22:03 ** CCGSLZ ** CCGSLZ PAGE 00012
282
283 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
284 // S S U U B B R R O O U U T I NN N E
285 // S U U B B R R O O U U T I N N N E
286 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
287 // S U U B B R R O O U U T I N N N E
288 // S S U U B B R R O O U U T I N NN E
289 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
290
291 // CCCCC CCCCC GGGGG SSSSS L ZZZZZZZ
292 // C C C C G G S S L Z
293 // C C G S L Z
294 // C C G GGGG SSSSS L Z
295 // C C G G S L Z
296 // C C C C G G S S L *
297 // CCCCC CCCCC GGGGG SSSSS LLLLLLL ZZZZZZZ
298
299
300 // VIF THU 20 NOV 2025 22:22:03 ** CCGSLZ ** DIAGNOSTICS PAGE 00013
301
302
303 // VIF THU 20 NOV 2025 22:22:03 ** CCGSLZ ** GENERATED CODE PAGE 00014
304
305 int_4 _ccgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ k_, int_4
306 _p_ lda_)
307 { // ** body not listed **
330 }
331
332 // VIF THU 20 NOV 2025 22:22:03 ** CCSLZ ** CCSLZ PAGE 00015
333
334 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
335 // S S U U B B R R O O U U T I NN N E
336 // S U U B B R R O O U U T I N N N E
337 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
338 // S U U B B R R O O U U T I N N N E
339 // S S U U B B R R O O U U T I N NN E
340 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
341
342 // CCCCC CCCCC SSSSS L ZZZZZZZ
343 // C C C C S S L Z
344 // C C S L Z
345 // C C SSSSS L Z
346 // C C S L Z
347 // C C C C S S L *
348 // CCCCC CCCCC SSSSS LLLLLLL ZZZZZZZ
349
350
351 // VIF THU 20 NOV 2025 22:22:03 ** CCSLZ ** DIAGNOSTICS PAGE 00016
352
353
354 // VIF THU 20 NOV 2025 22:22:03 ** CCSLZ ** GENERATED CODE PAGE 00017
355
356 int_4 _ccslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_)
357 { // ** body not listed **
374 }
375
376 // VIF THU 20 NOV 2025 22:22:03 ** CCTSLZ ** CCTSLZ PAGE 00018
377
378 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
379 // S S U U B B R R O O U U T I NN N E
380 // S U U B B R R O O U U T I N N N E
381 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
382 // S U U B B R R O O U U T I N N N E
383 // S S U U B B R R O O U U T I N NN E
384 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
385
386 // CCCCC CCCCC TTTTTTT SSSSS L ZZZZZZZ
387 // C C C C T S S L Z
388 // C C T S L Z
389 // C C T SSSSS L Z
390 // C C T S L Z
391 // C C C C T S S L *
392 // CCCCC CCCCC T SSSSS LLLLLLL ZZZZZZZ
393
394
395 // VIF THU 20 NOV 2025 22:22:03 ** CCTSLZ ** DIAGNOSTICS PAGE 00019
396
397
398 // VIF THU 20 NOV 2025 22:22:03 ** CCTSLZ ** GENERATED CODE PAGE 00020
399
400 int_4 _cctslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ k_, int_4
401 _p_ lda_)
402 { // ** body not listed **
425 }
426
427 // VIF THU 20 NOV 2025 22:22:03 ** CGSLZ ** CGSLZ PAGE 00021
428
429 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
430 // S S U U B B R R O O U U T I NN N E
431 // S U U B B R R O O U U T I N N N E
432 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
433 // S U U B B R R O O U U T I N N N E
434 // S S U U B B R R O O U U T I N NN E
435 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
436
437 // CCCCC GGGGG SSSSS L ZZZZZZZ
438 // C C G G S S L Z
439 // C G S L Z
440 // C G GGGG SSSSS L Z
441 // C G G S L Z
442 // C C G G S S L *
443 // CCCCC GGGGG SSSSS LLLLLLL ZZZZZZZ
444
445
446 // VIF THU 20 NOV 2025 22:22:03 ** CGSLZ ** DIAGNOSTICS PAGE 00022
447
448
449 // VIF THU 20 NOV 2025 22:22:03 ** CGSLZ ** GENERATED CODE PAGE 00023
450
451 int_4 _cgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_)
452 { // ** body not listed **
471 }
472
473 // VIF THU 20 NOV 2025 22:22:03 ** CQRD ** CQRD PAGE 00024
474
475 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
476 // S S U U B B R R O O U U T I NN N E
477 // S U U B B R R O O U U T I N N N E
478 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
479 // S U U B B R R O O U U T I N N N E
480 // S S U U B B R R O O U U T I N NN E
481 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
482
483 // CCCCC QQQQQ RRRRRR DDDDDD
484 // C C Q Q R R D D
485 // C Q Q R R D D
486 // C Q Q RRRRRR D D
487 // C Q Q Q R R D D
488 // C C Q Q R R D D
489 // CCCCC QQQQ Q R R DDDDDD
490
491
492 // VIF THU 20 NOV 2025 22:22:03 ** CQRD ** DIAGNOSTICS PAGE 00025
493
494
495 // VIF THU 20 NOV 2025 22:22:03 ** CQRD ** GENERATED CODE PAGE 00026
496
497 int_4 _cqrd (real_8 _p_ a_, real_8 _p_ q_, real_8 _p_ s_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ ldq_, int_4 _p_ lds_)
498 { // ** body not listed **
546 }
547 // VIF THU 20 NOV 2025 22:22:03 ** CQRD ** GENERATED CODE PAGE 00027
548
549
550 // VIF THU 20 NOV 2025 22:22:03 ** CQRZ ** CQRZ PAGE 00028
551
552 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
553 // S S U U B B R R O O U U T I NN N E
554 // S U U B B R R O O U U T I N N N E
555 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
556 // S U U B B R R O O U U T I N N N E
557 // S S U U B B R R O O U U T I N NN E
558 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
559
560 // CCCCC QQQQQ RRRRRR ZZZZZZZ
561 // C C Q Q R R Z
562 // C Q Q R R Z
563 // C Q Q RRRRRR Z
564 // C Q Q Q R R Z
565 // C C Q Q R R *
566 // CCCCC QQQQ Q R R ZZZZZZZ
567
568
569 // VIF THU 20 NOV 2025 22:22:03 ** CQRZ ** DIAGNOSTICS PAGE 00029
570
571
572 // VIF THU 20 NOV 2025 22:22:03 ** CQRZ ** GENERATED CODE PAGE 00030
573
574 int_4 _cqrz (complex_16 _p_ a_, complex_16 _p_ q_, complex_16 _p_ s_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ ldq_, int_4
575 _p_ lds_)
576 { // ** body not listed **
627 }
628
629 // VIF THU 20 NOV 2025 22:22:03 ** CSLZ ** CSLZ PAGE 00032
630
631 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
632 // S S U U B B R R O O U U T I NN N E
633 // S U U B B R R O O U U T I N N N E
634 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
635 // S U U B B R R O O U U T I N N N E
636 // S S U U B B R R O O U U T I N NN E
637 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
638
639 // CCCCC SSSSS L ZZZZZZZ
640 // C C S S L Z
641 // C S L Z
642 // C SSSSS L Z
643 // C S L Z
644 // C C S S L *
645 // CCCCC SSSSS LLLLLLL ZZZZZZZ
646
647
648 // VIF THU 20 NOV 2025 22:22:03 ** CSLZ ** DIAGNOSTICS PAGE 00033
649
650
651 // VIF THU 20 NOV 2025 22:22:03 ** CSLZ ** GENERATED CODE PAGE 00034
652
653 int_4 _cslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_)
654 { // ** body not listed **
699 }
700
701 // VIF THU 20 NOV 2025 22:22:03 ** CTGSLZ ** CTGSLZ PAGE 00035
702
703 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
704 // S S U U B B R R O O U U T I NN N E
705 // S U U B B R R O O U U T I N N N E
706 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
707 // S U U B B R R O O U U T I N N N E
708 // S S U U B B R R O O U U T I N NN E
709 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
710
711 // CCCCC TTTTTTT GGGGG SSSSS L ZZZZZZZ
712 // C C T G G S S L Z
713 // C T G S L Z
714 // C T G GGGG SSSSS L Z
715 // C T G G S L Z
716 // C C T G G S S L *
717 // CCCCC T GGGGG SSSSS LLLLLLL ZZZZZZZ
718
719
720 // VIF THU 20 NOV 2025 22:22:03 ** CTGSLZ ** DIAGNOSTICS PAGE 00036
721
722
723 // VIF THU 20 NOV 2025 22:22:03 ** CTGSLZ ** GENERATED CODE PAGE 00037
724
725 int_4 _ctgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ k_, int_4
726 _p_ lda_)
727 { // ** body not listed **
750 }
751
752 // VIF THU 20 NOV 2025 22:22:03 ** CTSLZ ** CTSLZ PAGE 00038
753
754 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
755 // S S U U B B R R O O U U T I NN N E
756 // S U U B B R R O O U U T I N N N E
757 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
758 // S U U B B R R O O U U T I N N N E
759 // S S U U B B R R O O U U T I N NN E
760 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
761
762 // CCCCC TTTTTTT SSSSS L ZZZZZZZ
763 // C C T S S L Z
764 // C T S L Z
765 // C T SSSSS L Z
766 // C T S L Z
767 // C C T S S L *
768 // CCCCC T SSSSS LLLLLLL ZZZZZZZ
769
770
771 // VIF THU 20 NOV 2025 22:22:03 ** CTSLZ ** DIAGNOSTICS PAGE 00039
772
773
774 // VIF THU 20 NOV 2025 22:22:03 ** CTSLZ ** GENERATED CODE PAGE 00040
775
776 int_4 _ctslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_)
777 { // ** body not listed **
795 }
796
797 // VIF THU 20 NOV 2025 22:22:03 ** SALWZ ** SALWZ PAGE 00041
798
799 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
800 // S S U U B B R R O O U U T I NN N E
801 // S U U B B R R O O U U T I N N N E
802 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
803 // S U U B B R R O O U U T I N N N E
804 // S S U U B B R R O O U U T I N NN E
805 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
806
807 // SSSSS A L W W ZZZZZZZ
808 // S S A A L W W Z
809 // S A A L W W Z
810 // SSSSS A A L W W W Z
811 // S AAAAAAA L W W W Z
812 // S S A A L W W W *
813 // SSSSS A A LLLLLLL WW WW ZZZZZZZ
814
815
816 // VIF THU 20 NOV 2025 22:22:03 ** SALWZ ** DIAGNOSTICS PAGE 00042
817
818
819 // VIF THU 20 NOV 2025 22:22:03 ** SALWZ ** GENERATED CODE PAGE 00043
820
821 int_4 _salwz (complex_16 _p_ a_, complex_16 _p_ r1_, complex_16 _p_ r2_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_,
822 int_4 _p_ job_)
823 { // ** body not listed **
859 }
860
861 // VIF THU 20 NOV 2025 22:22:03 ** TGSLD1 ** TGSLD1 PAGE 00044
862
863 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
864 // S S U U B B R R O O U U T I NN N E
865 // S U U B B R R O O U U T I N N N E
866 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
867 // S U U B B R R O O U U T I N N N E
868 // S S U U B B R R O O U U T I N NN E
869 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
870
871 // TTTTTTT GGGGG SSSSS L DDDDDD 1
872 // T G G S S L D D 11
873 // T G S L D D 1 1
874 // T G GGGG SSSSS L D D 1
875 // T G G S L D D 1
876 // T G G S S L D D 1
877 // T GGGGG SSSSS LLLLLLL DDDDDD 11111
878
879
880 // VIF THU 20 NOV 2025 22:22:03 ** TGSLD1 ** DIAGNOSTICS PAGE 00045
881
882
883 // VIF THU 20 NOV 2025 22:22:03 ** TGSLD1 ** GENERATED CODE PAGE 00046
884
885 int_4 _tgsld1 (real_8 _p_ a1_, real_8 _p_ a2_, real_8 _p_ b_, real_8 _p_ x_, real_8 _p_ c1_, real_8 _p_ c2_, real_8 _p_
886 r1_, real_8 _p_ r2_, real_8 _p_ r3_, real_8 _p_ r5_, real_8 _p_ r6_, real_8 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4
887 _p_ lda_)
888 { // ** body not listed **
1043 }
1044
1045 // VIF THU 20 NOV 2025 22:22:03 ** TGSLD ** TGSLD PAGE 00050
1046
1047 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1048 // S S U U B B R R O O U U T I NN N E
1049 // S U U B B R R O O U U T I N N N E
1050 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1051 // S U U B B R R O O U U T I N N N E
1052 // S S U U B B R R O O U U T I N NN E
1053 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1054
1055 // TTTTTTT GGGGG SSSSS L DDDDDD
1056 // T G G S S L D D
1057 // T G S L D D
1058 // T G GGGG SSSSS L D D
1059 // T G G S L D D
1060 // T G G S S L D D
1061 // T GGGGG SSSSS LLLLLLL DDDDDD
1062
1063
1064 // VIF THU 20 NOV 2025 22:22:03 ** TGSLD ** DIAGNOSTICS PAGE 00051
1065
1066
1067 // VIF THU 20 NOV 2025 22:22:03 ** TGSLD ** GENERATED CODE PAGE 00052
1068
1069 int_4 _tgsld (real_8 _p_ a_, real_8 _p_ x_, real_8 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_)
1070 { // ** body not listed **
1084 }
1085
1086 // VIF THU 20 NOV 2025 22:22:03 ** TGSLZ1 ** TGSLZ1 PAGE 00053
1087
1088 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1089 // S S U U B B R R O O U U T I NN N E
1090 // S U U B B R R O O U U T I N N N E
1091 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1092 // S U U B B R R O O U U T I N N N E
1093 // S S U U B B R R O O U U T I N NN E
1094 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1095
1096 // TTTTTTT GGGGG SSSSS L ZZZZZZZ 1
1097 // T G G S S L Z 11
1098 // T G S L Z 1 1
1099 // T G GGGG SSSSS L Z 1
1100 // T G G S L Z 1
1101 // T G G S S L * 1
1102 // T GGGGG SSSSS LLLLLLL ZZZZZZZ 11111
1103
1104
1105 // VIF THU 20 NOV 2025 22:22:03 ** TGSLZ1 ** DIAGNOSTICS PAGE 00054
1106
1107
1108 // VIF THU 20 NOV 2025 22:22:03 ** TGSLZ1 ** GENERATED CODE PAGE 00055
1109
1110 int_4 _tgslz1 (complex_16 _p_ a1_, complex_16 _p_ a2_, complex_16 _p_ b_, complex_16 _p_ x_, complex_16 _p_ c1_,
1111 complex_16 _p_ c2_, complex_16 _p_ r1_, complex_16 _p_ r2_, complex_16 _p_ r3_, complex_16 _p_ r5_, complex_16 _p_ r6_,
1112 complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_)
1113 { // ** body not listed **
1268 }
1269
1270 // VIF THU 20 NOV 2025 22:22:03 ** TGSLZ ** TGSLZ PAGE 00059
1271
1272 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1273 // S S U U B B R R O O U U T I NN N E
1274 // S U U B B R R O O U U T I N N N E
1275 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1276 // S U U B B R R O O U U T I N N N E
1277 // S S U U B B R R O O U U T I N NN E
1278 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1279
1280 // TTTTTTT GGGGG SSSSS L ZZZZZZZ
1281 // T G G S S L Z
1282 // T G S L Z
1283 // T G GGGG SSSSS L Z
1284 // T G G S L Z
1285 // T G G S S L *
1286 // T GGGGG SSSSS LLLLLLL ZZZZZZZ
1287
1288
1289 // VIF THU 20 NOV 2025 22:22:03 ** TGSLZ ** DIAGNOSTICS PAGE 00060
1290
1291
1292 // VIF THU 20 NOV 2025 22:22:03 ** TGSLZ ** GENERATED CODE PAGE 00061
1293
1294 int_4 _tgslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_, int_4 _p_ l_, int_4 _p_ lda_)
1295 { // ** body not listed **
1309 }
1310
1311 // VIF THU 20 NOV 2025 22:22:03 ** TSLD1 ** TSLD1 PAGE 00062
1312
1313 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1314 // S S U U B B R R O O U U T I NN N E
1315 // S U U B B R R O O U U T I N N N E
1316 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1317 // S U U B B R R O O U U T I N N N E
1318 // S S U U B B R R O O U U T I N NN E
1319 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1320
1321 // TTTTTTT SSSSS L DDDDDD 1
1322 // T S S L D D 11
1323 // T S L D D 1 1
1324 // T SSSSS L D D 1
1325 // T S L D D 1
1326 // T S S L D D 1
1327 // T SSSSS LLLLLLL DDDDDD 11111
1328
1329
1330 // VIF THU 20 NOV 2025 22:22:03 ** TSLD1 ** DIAGNOSTICS PAGE 00063
1331
1332
1333 // VIF THU 20 NOV 2025 22:22:03 ** TSLD1 ** GENERATED CODE PAGE 00064
1334
1335 int_4 _tsld1 (real_8 _p_ a1_, real_8 _p_ a2_, real_8 _p_ b_, real_8 _p_ x_, real_8 _p_ c1_, real_8 _p_ c2_, int_4 _p_
1336 m_)
1337 { // ** body not listed **
1392 }
1393
1394 // VIF THU 20 NOV 2025 22:22:03 ** TSLD ** TSLD PAGE 00066
1395
1396 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1397 // S S U U B B R R O O U U T I NN N E
1398 // S U U B B R R O O U U T I N N N E
1399 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1400 // S U U B B R R O O U U T I N N N E
1401 // S S U U B B R R O O U U T I N NN E
1402 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1403
1404 // TTTTTTT SSSSS L DDDDDD
1405 // T S S L D D
1406 // T S L D D
1407 // T SSSSS L D D
1408 // T S L D D
1409 // T S S L D D
1410 // T SSSSS LLLLLLL DDDDDD
1411
1412
1413 // VIF THU 20 NOV 2025 22:22:03 ** TSLD ** DIAGNOSTICS PAGE 00067
1414
1415
1416 // VIF THU 20 NOV 2025 22:22:03 ** TSLD ** GENERATED CODE PAGE 00068
1417
1418 int_4 _tsld (real_8 _p_ a_, real_8 _p_ x_, real_8 _p_ r_, int_4 _p_ m_)
1419 { // ** body not listed **
1423 }
1424
1425 // VIF THU 20 NOV 2025 22:22:03 ** TSLZ1 ** TSLZ1 PAGE 00069
1426
1427 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1428 // S S U U B B R R O O U U T I NN N E
1429 // S U U B B R R O O U U T I N N N E
1430 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1431 // S U U B B R R O O U U T I N N N E
1432 // S S U U B B R R O O U U T I N NN E
1433 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1434
1435 // TTTTTTT SSSSS L ZZZZZZZ 1
1436 // T S S L Z 11
1437 // T S L Z 1 1
1438 // T SSSSS L Z 1
1439 // T S L Z 1
1440 // T S S L * 1
1441 // T SSSSS LLLLLLL ZZZZZZZ 11111
1442
1443
1444 // VIF THU 20 NOV 2025 22:22:03 ** TSLZ1 ** DIAGNOSTICS PAGE 00070
1445
1446
1447 // VIF THU 20 NOV 2025 22:22:03 ** TSLZ1 ** GENERATED CODE PAGE 00071
1448
1449 int_4 _tslz1 (complex_16 _p_ a1_, complex_16 _p_ a2_, complex_16 _p_ b_, complex_16 _p_ x_, complex_16 _p_ c1_,
1450 complex_16 _p_ c2_, int_4 _p_ m_)
1451 { // ** body not listed **
1506 }
1507
1508 // VIF THU 20 NOV 2025 22:22:03 ** TSLZ ** TSLZ PAGE 00073
1509
1510 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
1511 // S S U U B B R R O O U U T I NN N E
1512 // S U U B B R R O O U U T I N N N E
1513 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
1514 // S U U B B R R O O U U T I N N N E
1515 // S S U U B B R R O O U U T I N NN E
1516 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
1517
1518 // TTTTTTT SSSSS L ZZZZZZZ
1519 // T S S L Z
1520 // T S L Z
1521 // T SSSSS L Z
1522 // T S L Z
1523 // T S S L *
1524 // T SSSSS LLLLLLL ZZZZZZZ
1525
1526
1527 // VIF THU 20 NOV 2025 22:22:03 ** TSLZ ** DIAGNOSTICS PAGE 00074
1528
1529
1530 // VIF THU 20 NOV 2025 22:22:03 ** TSLZ ** GENERATED CODE PAGE 00075
1531
1532 int_4 _tslz (complex_16 _p_ a_, complex_16 _p_ x_, complex_16 _p_ r_, int_4 _p_ m_)
1533 { // ** body not listed **
1537 }
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|