mathlib-slatec-dassl.c
1 //! @file mathlib-slatec-dassl.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 SLATEC-DASSL
24
25 // Compiled from Fortran source code by VIF.
26
27 // The license for SLATEC Fortran source code is:
28 //
29 // The SLATEC Common Mathematical Library was developed at
30 // US government research laboratories and is in the public domain.
31 //
32 // Repository: http://www.netlib.org/slatec/
33 //
34 // The SLATEC common mathematical library is issued by the following
35 //
36 // Air Force Weapons Laboratory, Albuquerque
37 // Lawrence Livermore National Laboratory, Livermore
38 // Los Alamos National Laboratory, Los Alamos
39 // National Institute of Standards and Technology, Washington
40 // National Energy Research Supercomputer Center, Livermore
41 // Oak Ridge National Laboratory, Oak Ridge
42 // Sandia National Laboratories, Albuquerque
43 // Sandia National Laboratories, Livermore
44 //
45 // All questions concerning the distribution of the library should be
46 // directed to the NATIONAL ENERGY SOFTWARE CENTER, 9700 Cass Ave.,
47 // Argonne, Illinois 60439, and not to the authors of the subprograms.
48 //
49 // * * * * * Notice * * * * *
50 //
51 // This material was prepared as an account of work sponsored by the
52 // United States Government. Neither the United States, nor the
53 // Department of Energy, nor the Department of Defense, nor any of
54 // their employees, nor any of their contractors, subcontractors, or
55 // their employees, makes any warranty, expressed or implied, or
56 // assumes any legal liability or responsibility for the accuracy,
57 // completeness, or usefulness of any information, apparatus, product,
58 // or process disclosed, or represents that its use would not infringe
59 // upon privately owned rights.
60
61 // VIF THU 20 NOV 2025 22:19:14 ** MARCEL ** SLATEC DASSL PAGE 00001
62
63
64
65 // M M A RRRRRR CCCCC EEEEEEE L
66 // MM MM A A R R C C E L
67 // M M M M A A R R C E L
68 // M M M A A RRRRRR C EEEEE L
69 // M M AAAAAAA R R C E L
70 // M M A A R R C C E L
71 // M M A A R R CCCCC EEEEEEE LLLLLLL
72
73 // SSSSS L A TTTTTTT EEEEEEE CCCCC DDDDDD A SSSSS SSSSS L
74 // S S L A A T E C C D D A A S S S S L
75 // S L A A T E C D D A A S S L
76 // SSSSS L A A T EEEEE C ***** D D A A SSSSS SSSSS L
77 // S L AAAAAAA T E C D D AAAAAAA S S L
78 // S S L A A T E C C D D A A S S S S L
79 // SSSSS LLLLLLL A A T EEEEEEE CCCCC DDDDDD A A SSSSS SSSSS LLLLLLL
80
81
82
83
84 // VIF THU 20 NOV 2025 22:19:14 ** GLOBAL SCOPE ** DIAGNOSTICS PAGE 00002
85
86 // ** linker ** object size 426736 bytes
87 // ** statistics ** 16 subprograms, no errors, no warnings
88
89
90 // VIF THU 20 NOV 2025 22:19:14 ** GLOBAL SCOPE ** JOB CONTROL PAGE 00003
91
92
93 // Line JCL *...*....1....|....2....|....3....|....4....|....5....|....6....|....7..*.|....8
94
95
96
97 // VIF THU 20 NOV 2025 22:19:14 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00004
98
99 /*
100 Generated by VIF - experimental VIntage Fortran compiler.
101 VIF release 1.4.1
102 */
103
104 #if defined (__GNUC__)
105 #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
106 #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
107 #if (__GNUC__ >= 14)
108 #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
109 #pragma GCC diagnostic ignored "-Wimplicit-int"
110 #pragma GCC diagnostic ignored "-Wint-conversion"
111 #pragma GCC diagnostic ignored "-Wreturn-mismatch"
112 #endif
113 #else
114 #error VIF requires GCC
115 #endif
116
117 #include <vif.h>
118
119 static int_4 _km1 = -1, _k0 = 0, _k1 = 1;
120
121
122 #define _dc_0 "SLATEC"
123 #define _dc_1 "DDASSL"
124 #define _dc_2 "THE LAST STEP TERMINATED WITH A NEGATIVE VALUE OF IDID = "
125 #define _dc_3 " AND NO APPROPRIATE ACTION WAS TAKEN. "
126 #define _dc_4 "RUN TERMINATED"
127 #define _dc_5 "AT CURRENT T = "
128 #define _dc_6 " 500 STEPS TAKEN ON THIS "
129 #define _dc_7 "CALL BEFORE REACHING TOUT"
130 #define _dc_8 "AT T = "
131 #define _dc_9 " TOO MUCH ACCURACY REQUESTED FOR "
132 #define _dc_10 "PRECISION OF MACHINE. RTOL AND ATOL WERE INCREASED TO "
133 #define _dc_11 "APPROPRIATE VALUES"
134 #define _dc_12 " SOME ELEMENT OF WT HAS BECOME .LE. "
135 #define _dc_13 "0.0"
136 #define _dc_14 " AND STEPSIZE H = "
137 #define _dc_15 " THE ERROR TEST FAILED REPEATEDLY OR WITH ABS(H)=HMIN"
138 #define _dc_16 " THE CORRECTOR FAILED TO CONVERGE REPEATEDLY OR WITH "
139 #define _dc_17 "ABS(H)=HMIN"
140 #define _dc_18 " THE ITERATION MATRIX IS SINGULAR"
141 #define _dc_19 " THE CORRECTOR COULD NOT CONVERGE. ALSO, THE ERROR TEST "
142 #define _dc_20 "FAILED REPEATEDLY."
143 #define _dc_21 " THE CORRECTOR COULD NOT CONVERGE BECAUSE IRES WAS EQUAL "
144 #define _dc_22 "TO MINUS ONE"
145 #define _dc_23 " IRES WAS EQUAL TO MINUS TWO"
146 #define _dc_24 " THE INITIAL YPRIME COULD NOT BE COMPUTED"
147 #define _dc_25 "SOME ELEMENT OF INFO VECTOR IS NOT ZERO OR ONE"
148 #define _dc_26 "NEQ = "
149 #define _dc_27 " .LE. 0"
150 #define _dc_28 "MAXORD = "
151 #define _dc_29 " NOT IN RANGE"
152 // VIF THU 20 NOV 2025 22:19:14 ** GLOBAL SCOPE ** DEFINITIONS PAGE 00005
153
154 #define _dc_30 "RWORK LENGTH NEEDED, LENRW = "
155 #define _dc_31 ", EXCEEDS LRW = "
156 #define _dc_32 "IWORK LENGTH NEEDED, LENIW = "
157 #define _dc_33 ", EXCEEDS LIW = "
158 #define _dc_34 "SOME ELEMENT OF RTOL IS .LT. 0"
159 #define _dc_35 "SOME ELEMENT OF ATOL IS .LT. 0"
160 #define _dc_36 "ALL ELEMENTS OF RTOL AND ATOL ARE ZERO"
161 #define _dc_37 "INFO(4) = 1 AND TSTOP = "
162 #define _dc_38 " BEHIND TOUT = "
163 #define _dc_39 "HMAX = "
164 #define _dc_40 " .LT. 0.0"
165 #define _dc_41 "TOUT = "
166 #define _dc_42 " BEHIND T = "
167 #define _dc_43 "INFO(8)=1 AND H0=0.0"
168 #define _dc_44 "SOME ELEMENT OF WT IS .LE. 0.0"
169 #define _dc_45 " TOO CLOSE TO T = "
170 #define _dc_46 " TO START INTEGRATION"
171 #define _dc_47 "INFO(4)=1 AND TSTOP = "
172 #define _dc_48 "ML = "
173 #define _dc_49 " ILLEGAL. EITHER .LT. 0 OR .GT. NEQ"
174 #define _dc_50 "MU = "
175 #define _dc_51 "TOUT = T = "
176 #define _dc_52 "REPEATED OCCURRENCES OF ILLEGAL INPUT$$"
177 #define _dc_53 "RUN TERMINATED. APPARENT INFINITE LOOP"
178 #define _dc_54 "SDASSL"
179
180
181
182 // VIF THU 20 NOV 2025 22:19:14 ** GLOBAL SCOPE ** TYPEDEFS PAGE 00006
183
184 static FORMAT *__fmt_a = NULL;
185 typedef char char_3[4];
186 typedef char char_7[8];
187 typedef char char_15[16];
188 typedef char char_31[32];
189 typedef char char_63[64];
190 typedef char char_127[128];
191 typedef char char_255[256];
192 typedef char char_511[512];
193 typedef char char_1023[1024];
194
195
196 // VIF THU 20 NOV 2025 22:19:14 ** GLOBAL SCOPE ** FREQUENCY TABLE PAGE 00007
197
198 #define __ncalls 17
199 static CALLS __calls[__ncalls] = {
200 { // ** body not listed **"ddaini", 0}, // subroutine
202 { // ** body not listed **"ddanrm", 0}, // real*8 function
204 { // ** body not listed **"ddassl", 0}, // subroutine
206 { // ** body not listed **"ddatrp", 0}, // subroutine
208 { // ** body not listed **"sdaini", 0}, // subroutine
210 { // ** body not listed **"sdanrm", 0}, // real*4 function
212 { // ** body not listed **"sdassl", 0}, // subroutine
214 { // ** body not listed **"sdatrp", 0}, // subroutine
216 { // ** body not listed **NULL, 0}
217 };
315 { // ** body not listed **
514 }
515
516 // VIF THU 20 NOV 2025 22:19:14 ** DDAJAC ** DDAJAC PAGE 00017
517
518 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
519 // S S U U B B R R O O U U T I NN N E
520 // S U U B B R R O O U U T I N N N E
521 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
522 // S U U B B R R O O U U T I N N N E
523 // S S U U B B R R O O U U T I N NN E
524 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
525
526 // DDDDDD DDDDDD A J A CCCCC
527 // D D D D A A J A A C C
528 // D D D D A A J A A C
529 // D D D D A A J A A C
530 // D D D D AAAAAAA J J AAAAAAA C
531 // D D D D A A J J A A C C
532 // DDDDDD DDDDDD A A JJJJJ A A CCCCC
533
534
535 // VIF THU 20 NOV 2025 22:19:14 ** DDAJAC ** DIAGNOSTICS PAGE 00018
536
537
538 // VIF THU 20 NOV 2025 22:19:14 ** DDAJAC ** GENERATED CODE PAGE 00019
539
540 int_4 _ddajac (int_4 _p_ neq_, real_8 _p_ x_, real_8 _p_ y_, real_8 _p_ yprime_, real_8 _p_ delta_, real_8 _p_ cj_,
541 real_8 _p_ h_, int_4 _p_ ier_, real_8 _p_ wt_, real_8 _p_ e_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4
542 _p_ ires_, real_8 _p_ uround_, int_4 (*_jac)(), real_8 _p_ rpar_, int_4 _p_ ipar_, int_4 _p_ ntemp_)
543 { // ** body not listed **
655 }
656
657 // VIF THU 20 NOV 2025 22:19:14 ** DDANRM ** DDANRM PAGE 00022
658
659 // RRRRRR EEEEEEE A L 88888
660 // R R E A A L * * 8 8
661 // R R E A A L * * 8 8
662 // RRRRRR EEEEE A A L ******* 88888
663 // R R E AAAAAAA L * * 8 8
664 // R R E A A L * * 8 8
665 // R R EEEEEEE A A LLLLLLL 88888
666
667 // FFFFFFF U U N N CCCCC TTTTTTT III OOOOO N N
668 // F U U NN N C C T I O O NN N
669 // F U U N N N C T I O O N N N
670 // FFFFF U U N N N C T I O O N N N
671 // F U U N N N C T I O O N N N
672 // F U U N NN C C T I O O N NN
673 // F UUUUU N N CCCCC T III OOOOO N N
674
675 // DDDDDD DDDDDD A N N RRRRRR M M
676 // D D D D A A NN N R R MM MM
677 // D D D D A A N N N R R M M M M
678 // D D D D A A N N N RRRRRR M M M
679 // D D D D AAAAAAA N N N R R M M
680 // D D D D A A N NN R R M M
681 // DDDDDD DDDDDD A A N N R R M M
682
683
684 // VIF THU 20 NOV 2025 22:19:14 ** DDANRM ** DIAGNOSTICS PAGE 00023
685
686
687 // VIF THU 20 NOV 2025 22:19:14 ** DDANRM ** GENERATED CODE PAGE 00024
688
689 real_8 _ddanrm (int_4 _p_ neq_, real_8 _p_ v_, real_8 _p_ wt_, real_8 _p_ rpar_, int_4 _p_ ipar_)
690 { // ** body not listed **
713 }
714
715 // VIF THU 20 NOV 2025 22:19:14 ** DDASLV ** DDASLV PAGE 00025
716
717 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
718 // S S U U B B R R O O U U T I NN N E
719 // S U U B B R R O O U U T I N N N E
720 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
721 // S U U B B R R O O U U T I N N N E
722 // S S U U B B R R O O U U T I N NN E
723 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
724
725 // DDDDDD DDDDDD A SSSSS L V V
726 // D D D D A A S S L V V
727 // D D D D A A S L V V
728 // D D D D A A SSSSS L V V
729 // D D D D AAAAAAA S L V V
730 // D D D D A A S S L V V
731 // DDDDDD DDDDDD A A SSSSS LLLLLLL V
732
733
734 // VIF THU 20 NOV 2025 22:19:14 ** DDASLV ** DIAGNOSTICS PAGE 00026
735
736
737 // VIF THU 20 NOV 2025 22:19:14 ** DDASLV ** GENERATED CODE PAGE 00027
738
739 int_4 _ddaslv (int_4 _p_ neq_, real_8 _p_ delta_, real_8 _p_ wm_, int_4 _p_ iwm_)
740 { // ** body not listed **
764 }
765
766 // VIF THU 20 NOV 2025 22:19:14 ** DDASSL ** DDASSL PAGE 00028
767
768 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
769 // S S U U B B R R O O U U T I NN N E
770 // S U U B B R R O O U U T I N N N E
771 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
772 // S U U B B R R O O U U T I N N N E
773 // S S U U B B R R O O U U T I N NN E
774 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
775
776 // DDDDDD DDDDDD A SSSSS SSSSS L
777 // D D D D A A S S S S L
778 // D D D D A A S S L
779 // D D D D A A SSSSS SSSSS L
780 // D D D D AAAAAAA S S L
781 // D D D D A A S S S S L
782 // DDDDDD DDDDDD A A SSSSS SSSSS LLLLLLL
783
784
785 // VIF THU 20 NOV 2025 22:19:14 ** DDASSL ** DIAGNOSTICS PAGE 00029
786
787
788 // VIF THU 20 NOV 2025 22:19:14 ** DDASSL ** GENERATED CODE PAGE 00030
789
790 int_4 _ddassl (real_4 (*_res)(), int_4 _p_ neq_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ yprime_, real_8 _p_ tout_,
791 int_4 _p_ info_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ idid_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_
792 iwork_, int_4 _p_ liw_, real_8 _p_ rpar_, int_4 _p_ ipar_, int_4 (*_jac)())
793 { // ** body not listed **
2864 }
2865
2866 // VIF THU 20 NOV 2025 22:19:14 ** DDASTP ** DDASTP PAGE 00070
2867
2868 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
2869 // S S U U B B R R O O U U T I NN N E
2870 // S U U B B R R O O U U T I N N N E
2871 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
2872 // S U U B B R R O O U U T I N N N E
2873 // S S U U B B R R O O U U T I N NN E
2874 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
2875
2876 // DDDDDD DDDDDD A SSSSS TTTTTTT PPPPPP
2877 // D D D D A A S S T P P
2878 // D D D D A A S T P P
2879 // D D D D A A SSSSS T PPPPPP
2880 // D D D D AAAAAAA S T P
2881 // D D D D A A S S T P
2882 // DDDDDD DDDDDD A A SSSSS T P
2883
2884
2885 // VIF THU 20 NOV 2025 22:19:14 ** DDASTP ** DIAGNOSTICS PAGE 00071
2886
2887
2888 // VIF THU 20 NOV 2025 22:19:14 ** DDASTP ** GENERATED CODE PAGE 00072
2889
2890 int_4 _ddastp (real_8 _p_ x_, real_8 _p_ y_, real_8 _p_ yprime_, int_4 _p_ neq_, int_4 (*_res)(), int_4 (*_jac)(),
2891 real_8 _p_ h_, real_8 _p_ wt_, int_4 _p_ jstart_, int_4 _p_ idid_, real_8 _p_ rpar_, int_4 _p_ ipar_, real_8 _p_ phi_,
2892 real_8 _p_ delta_, real_8 _p_ e_, real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ alpha_, real_8 _p_ beta_, real_8 _p_
2893 gamma_, real_8 _p_ psi_, real_8 _p_ sigma_, real_8 _p_ cj_, real_8 _p_ cjold_, real_8 _p_ hold_, real_8 _p_ s_, real_8
2894 _p_ hmin_, real_8 _p_ uround_, int_4 _p_ iphase_, int_4 _p_ jcalc_, int_4 _p_ k_, int_4 _p_ kold_, int_4 _p_ ns_, int_4
2895 _p_ nonneg_, int_4 _p_ ntemp_)
2896 { // ** body not listed **
3351 }
3352
3353 // VIF THU 20 NOV 2025 22:19:14 ** DDATRP ** DDATRP PAGE 00081
3354
3355 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3356 // S S U U B B R R O O U U T I NN N E
3357 // S U U B B R R O O U U T I N N N E
3358 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3359 // S U U B B R R O O U U T I N N N E
3360 // S S U U B B R R O O U U T I N NN E
3361 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3362
3363 // DDDDDD DDDDDD A TTTTTTT RRRRRR PPPPPP
3364 // D D D D A A T R R P P
3365 // D D D D A A T R R P P
3366 // D D D D A A T RRRRRR PPPPPP
3367 // D D D D AAAAAAA T R R P
3368 // D D D D A A T R R P
3369 // DDDDDD DDDDDD A A T R R P
3370
3371
3372 // VIF THU 20 NOV 2025 22:19:14 ** DDATRP ** DIAGNOSTICS PAGE 00082
3373
3374
3375 // VIF THU 20 NOV 2025 22:19:14 ** DDATRP ** GENERATED CODE PAGE 00083
3376
3377 int_4 _ddatrp (real_8 _p_ x_, real_8 _p_ xout_, real_8 _p_ yout_, real_8 _p_ ypout_, int_4 _p_ neq_, int_4 _p_ kold_,
3378 real_8 _p_ phi_, real_8 _p_ psi_)
3379 { // ** body not listed **
3402 }
3403
3404 // VIF THU 20 NOV 2025 22:19:14 ** DDAWTS ** DDAWTS PAGE 00084
3405
3406 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3407 // S S U U B B R R O O U U T I NN N E
3408 // S U U B B R R O O U U T I N N N E
3409 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3410 // S U U B B R R O O U U T I N N N E
3411 // S S U U B B R R O O U U T I N NN E
3412 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3413
3414 // DDDDDD DDDDDD A W W TTTTTTT SSSSS
3415 // D D D D A A W W T S S
3416 // D D D D A A W W T S
3417 // D D D D A A W W W T SSSSS
3418 // D D D D AAAAAAA W W W T S
3419 // D D D D A A W W W T S S
3420 // DDDDDD DDDDDD A A WW WW T SSSSS
3421
3422
3423 // VIF THU 20 NOV 2025 22:19:14 ** DDAWTS ** DIAGNOSTICS PAGE 00085
3424
3425
3426 // VIF THU 20 NOV 2025 22:19:14 ** DDAWTS ** GENERATED CODE PAGE 00086
3427
3428 int_4 _ddawts (int_4 _p_ neq_, int_4 _p_ iwt_, real_8 _p_ rtol_, real_8 _p_ atol_, real_8 _p_ y_, real_8 _p_ wt_,
3429 real_8 _p_ rpar_, int_4 _p_ ipar_)
3430 { // ** body not listed **
3446 }
3447
3448 // VIF THU 20 NOV 2025 22:19:14 ** SDAINI ** SDAINI PAGE 00087
3449
3450 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3451 // S S U U B B R R O O U U T I NN N E
3452 // S U U B B R R O O U U T I N N N E
3453 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3454 // S U U B B R R O O U U T I N N N E
3455 // S S U U B B R R O O U U T I N NN E
3456 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3457
3458 // SSSSS DDDDDD A III N N III
3459 // S S D D A A I NN N I
3460 // S D D A A I N N N I
3461 // SSSSS D D A A I N N N I
3462 // S D D AAAAAAA I N N N I
3463 // S S D D A A I N NN I
3464 // SSSSS DDDDDD A A III N N III
3465
3466
3467 // VIF THU 20 NOV 2025 22:19:14 ** SDAINI ** DIAGNOSTICS PAGE 00088
3468
3469
3470 // VIF THU 20 NOV 2025 22:19:14 ** SDAINI ** GENERATED CODE PAGE 00089
3471
3472 int_4 _sdaini (real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ yprime_, int_4 _p_ neq_, int_4 (*_res)(), int_4 (*_jac)(),
3473 real_4 _p_ h_, real_4 _p_ wt_, int_4 _p_ idid_, real_4 _p_ rpar_, int_4 _p_ ipar_, real_4 _p_ phi_, real_4 _p_ delta_,
3474 real_4 _p_ e_, real_4 _p_ wm_, int_4 _p_ iwm_, real_4 _p_ hmin_, real_4 _p_ uround_, int_4 _p_ nonneg_, int_4 _p_
3475 ntemp_)
3476 { // ** body not listed **
3675 }
3676
3677 // VIF THU 20 NOV 2025 22:19:14 ** SDAJAC ** SDAJAC PAGE 00093
3678
3679 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3680 // S S U U B B R R O O U U T I NN N E
3681 // S U U B B R R O O U U T I N N N E
3682 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3683 // S U U B B R R O O U U T I N N N E
3684 // S S U U B B R R O O U U T I N NN E
3685 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3686
3687 // SSSSS DDDDDD A J A CCCCC
3688 // S S D D A A J A A C C
3689 // S D D A A J A A C
3690 // SSSSS D D A A J A A C
3691 // S D D AAAAAAA J J AAAAAAA C
3692 // S S D D A A J J A A C C
3693 // SSSSS DDDDDD A A JJJJJ A A CCCCC
3694
3695
3696 // VIF THU 20 NOV 2025 22:19:14 ** SDAJAC ** DIAGNOSTICS PAGE 00094
3697
3698
3699 // VIF THU 20 NOV 2025 22:19:14 ** SDAJAC ** GENERATED CODE PAGE 00095
3700
3701 int_4 _sdajac (int_4 _p_ neq_, real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ yprime_, real_4 _p_ delta_, real_4 _p_ cj_,
3702 real_4 _p_ h_, int_4 _p_ ier_, real_4 _p_ wt_, real_4 _p_ e_, real_4 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4
3703 _p_ ires_, real_4 _p_ uround_, int_4 (*_jac)(), real_4 _p_ rpar_, int_4 _p_ ipar_, int_4 _p_ ntemp_)
3704 { // ** body not listed **
3816 }
3817
3818 // VIF THU 20 NOV 2025 22:19:14 ** SDANRM ** SDANRM PAGE 00098
3819
3820 // RRRRRR EEEEEEE A L 4
3821 // R R E A A L * * 4 4
3822 // R R E A A L * * 4 4
3823 // RRRRRR EEEEE A A L ******* 4444444
3824 // R R E AAAAAAA L * * 4
3825 // R R E A A L * * 4
3826 // R R EEEEEEE A A LLLLLLL 4
3827
3828 // FFFFFFF U U N N CCCCC TTTTTTT III OOOOO N N
3829 // F U U NN N C C T I O O NN N
3830 // F U U N N N C T I O O N N N
3831 // FFFFF U U N N N C T I O O N N N
3832 // F U U N N N C T I O O N N N
3833 // F U U N NN C C T I O O N NN
3834 // F UUUUU N N CCCCC T III OOOOO N N
3835
3836 // SSSSS DDDDDD A N N RRRRRR M M
3837 // S S D D A A NN N R R MM MM
3838 // S D D A A N N N R R M M M M
3839 // SSSSS D D A A N N N RRRRRR M M M
3840 // S D D AAAAAAA N N N R R M M
3841 // S S D D A A N NN R R M M
3842 // SSSSS DDDDDD A A N N R R M M
3843
3844
3845 // VIF THU 20 NOV 2025 22:19:14 ** SDANRM ** DIAGNOSTICS PAGE 00099
3846
3847
3848 // VIF THU 20 NOV 2025 22:19:14 ** SDANRM ** GENERATED CODE PAGE 00100
3849
3850 real_4 _sdanrm (int_4 _p_ neq_, real_4 _p_ v_, real_4 _p_ wt_, real_4 _p_ rpar_, int_4 _p_ ipar_)
3851 { // ** body not listed **
3874 }
3875
3876 // VIF THU 20 NOV 2025 22:19:14 ** SDASLV ** SDASLV PAGE 00101
3877
3878 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3879 // S S U U B B R R O O U U T I NN N E
3880 // S U U B B R R O O U U T I N N N E
3881 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3882 // S U U B B R R O O U U T I N N N E
3883 // S S U U B B R R O O U U T I N NN E
3884 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3885
3886 // SSSSS DDDDDD A SSSSS L V V
3887 // S S D D A A S S L V V
3888 // S D D A A S L V V
3889 // SSSSS D D A A SSSSS L V V
3890 // S D D AAAAAAA S L V V
3891 // S S D D A A S S L V V
3892 // SSSSS DDDDDD A A SSSSS LLLLLLL V
3893
3894
3895 // VIF THU 20 NOV 2025 22:19:14 ** SDASLV ** DIAGNOSTICS PAGE 00102
3896
3897
3898 // VIF THU 20 NOV 2025 22:19:14 ** SDASLV ** GENERATED CODE PAGE 00103
3899
3900 int_4 _sdaslv (int_4 _p_ neq_, real_4 _p_ delta_, real_4 _p_ wm_, int_4 _p_ iwm_)
3901 { // ** body not listed **
3925 }
3926
3927 // VIF THU 20 NOV 2025 22:19:14 ** SDASSL ** SDASSL PAGE 00104
3928
3929 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
3930 // S S U U B B R R O O U U T I NN N E
3931 // S U U B B R R O O U U T I N N N E
3932 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
3933 // S U U B B R R O O U U T I N N N E
3934 // S S U U B B R R O O U U T I N NN E
3935 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
3936
3937 // SSSSS DDDDDD A SSSSS SSSSS L
3938 // S S D D A A S S S S L
3939 // S D D A A S S L
3940 // SSSSS D D A A SSSSS SSSSS L
3941 // S D D AAAAAAA S S L
3942 // S S D D A A S S S S L
3943 // SSSSS DDDDDD A A SSSSS SSSSS LLLLLLL
3944
3945
3946 // VIF THU 20 NOV 2025 22:19:14 ** SDASSL ** DIAGNOSTICS PAGE 00105
3947
3948
3949 // VIF THU 20 NOV 2025 22:19:14 ** SDASSL ** GENERATED CODE PAGE 00106
3950
3951 int_4 _sdassl (real_4 (*_res)(), int_4 _p_ neq_, real_4 _p_ t_, real_4 _p_ y_, real_4 _p_ yprime_, real_4 _p_ tout_,
3952 int_4 _p_ info_, real_4 _p_ rtol_, real_4 _p_ atol_, int_4 _p_ idid_, real_4 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_
3953 iwork_, int_4 _p_ liw_, real_4 _p_ rpar_, int_4 _p_ ipar_, int_4 (*_jac)())
3954 { // ** body not listed **
6025 }
6026
6027 // VIF THU 20 NOV 2025 22:19:14 ** SDASTP ** SDASTP PAGE 00146
6028
6029 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
6030 // S S U U B B R R O O U U T I NN N E
6031 // S U U B B R R O O U U T I N N N E
6032 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
6033 // S U U B B R R O O U U T I N N N E
6034 // S S U U B B R R O O U U T I N NN E
6035 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
6036
6037 // SSSSS DDDDDD A SSSSS TTTTTTT PPPPPP
6038 // S S D D A A S S T P P
6039 // S D D A A S T P P
6040 // SSSSS D D A A SSSSS T PPPPPP
6041 // S D D AAAAAAA S T P
6042 // S S D D A A S S T P
6043 // SSSSS DDDDDD A A SSSSS T P
6044
6045
6046 // VIF THU 20 NOV 2025 22:19:14 ** SDASTP ** DIAGNOSTICS PAGE 00147
6047
6048
6049 // VIF THU 20 NOV 2025 22:19:14 ** SDASTP ** GENERATED CODE PAGE 00148
6050
6051 int_4 _sdastp (real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ yprime_, int_4 _p_ neq_, int_4 (*_res)(), int_4 (*_jac)(),
6052 real_4 _p_ h_, real_4 _p_ wt_, int_4 _p_ jstart_, int_4 _p_ idid_, real_4 _p_ rpar_, int_4 _p_ ipar_, real_4 _p_ phi_,
6053 real_4 _p_ delta_, real_4 _p_ e_, real_4 _p_ wm_, int_4 _p_ iwm_, real_4 _p_ alpha_, real_4 _p_ beta_, real_4 _p_
6054 gamma_, real_4 _p_ psi_, real_4 _p_ sigma_, real_4 _p_ cj_, real_4 _p_ cjold_, real_4 _p_ hold_, real_4 _p_ s_, real_4
6055 _p_ hmin_, real_4 _p_ uround_, int_4 _p_ iphase_, int_4 _p_ jcalc_, int_4 _p_ k_, int_4 _p_ kold_, int_4 _p_ ns_, int_4
6056 _p_ nonneg_, int_4 _p_ ntemp_)
6057 { // ** body not listed **
6512 }
6513
6514 // VIF THU 20 NOV 2025 22:19:14 ** SDATRP ** SDATRP PAGE 00157
6515
6516 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
6517 // S S U U B B R R O O U U T I NN N E
6518 // S U U B B R R O O U U T I N N N E
6519 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
6520 // S U U B B R R O O U U T I N N N E
6521 // S S U U B B R R O O U U T I N NN E
6522 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
6523
6524 // SSSSS DDDDDD A TTTTTTT RRRRRR PPPPPP
6525 // S S D D A A T R R P P
6526 // S D D A A T R R P P
6527 // SSSSS D D A A T RRRRRR PPPPPP
6528 // S D D AAAAAAA T R R P
6529 // S S D D A A T R R P
6530 // SSSSS DDDDDD A A T R R P
6531
6532
6533 // VIF THU 20 NOV 2025 22:19:14 ** SDATRP ** DIAGNOSTICS PAGE 00158
6534
6535
6536 // VIF THU 20 NOV 2025 22:19:14 ** SDATRP ** GENERATED CODE PAGE 00159
6537
6538 int_4 _sdatrp (real_4 _p_ x_, real_4 _p_ xout_, real_4 _p_ yout_, real_4 _p_ ypout_, int_4 _p_ neq_, int_4 _p_ kold_,
6539 real_4 _p_ phi_, real_4 _p_ psi_)
6540 { // ** body not listed **
6563 }
6564
6565 // VIF THU 20 NOV 2025 22:19:14 ** SDAWTS ** SDAWTS PAGE 00160
6566
6567 // SSSSS U U BBBBBB RRRRRR OOOOO U U TTTTTTT III N N EEEEEEE
6568 // S S U U B B R R O O U U T I NN N E
6569 // S U U B B R R O O U U T I N N N E
6570 // SSSSS U U BBBBBB RRRRRR O O U U T I N N N EEEEE
6571 // S U U B B R R O O U U T I N N N E
6572 // S S U U B B R R O O U U T I N NN E
6573 // SSSSS UUUUU BBBBBB R R OOOOO UUUUU T III N N EEEEEEE
6574
6575 // SSSSS DDDDDD A W W TTTTTTT SSSSS
6576 // S S D D A A W W T S S
6577 // S D D A A W W T S
6578 // SSSSS D D A A W W W T SSSSS
6579 // S D D AAAAAAA W W W T S
6580 // S S D D A A W W W T S S
6581 // SSSSS DDDDDD A A WW WW T SSSSS
6582
6583
6584 // VIF THU 20 NOV 2025 22:19:14 ** SDAWTS ** DIAGNOSTICS PAGE 00161
6585
6586
6587 // VIF THU 20 NOV 2025 22:19:14 ** SDAWTS ** GENERATED CODE PAGE 00162
6588
6589 int_4 _sdawts (int_4 _p_ neq_, int_4 _p_ iwt_, real_4 _p_ rtol_, real_4 _p_ atol_, real_4 _p_ y_, real_4 _p_ wt_,
6590 real_4 _p_ rpar_, int_4 _p_ ipar_)
6591 { // ** body not listed **
6607 }
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
