mathlib-slatec-sdrive.c
1 //! @file mathlib-slatec-sdrive.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 //! 36 subprograms from SLATEC-SDRIVE
24
25 // Compiled from Fortran source code by VIF.
26 // Selected subprograms are:
27 //
28 // CDCOR CDPSC CDRIV3 DDCOR DDPSC DDRIV3 SDCOR SDNTP SDRIV1 SDSCL
29 // CDCST CDPST CDSCL DDCST DDPST DDSCL SDCST SDPSC SDRIV2 SDSTP
30 // CDNTL CDRIV1 CDSTP DDNTL DDRIV1 DDSTP SDNTL SDPST SDRIV3 SDZRO
31 // CDNTP CDRIV2 CDZRO DDNTP DDRIV2 DDZRO
32
33 // SDRIVE routines from SLATEC.
34 //
35 // Source: netlib.org/slatec/
36
37 // The license for SLATEC Fortran source code is:
38 //
39 // The SLATEC Common Mathematical Library was developed at
40 // US government research laboratories and is in the public domain.
41 //
42 // Repository: http://www.netlib.org/slatec/
43 //
44 // The SLATEC common mathematical library is issued by the following
45 //
46 // Air Force Weapons Laboratory, Albuquerque
47 // Lawrence Livermore National Laboratory, Livermore
48 // Los Alamos National Laboratory, Los Alamos
49 // National Institute of Standards and Technology, Washington
50 // National Energy Research Supercomputer Center, Livermore
51 // Oak Ridge National Laboratory, Oak Ridge
52 // Sandia National Laboratories, Albuquerque
53 // Sandia National Laboratories, Livermore
54 //
55 // All questions concerning the distribution of the library should be
56 // directed to the NATIONAL ENERGY SOFTWARE CENTER, 9700 Cass Ave.,
57 // Argonne, Illinois 60439, and not to the authors of the subprograms.
58 //
59 // * * * * * Notice * * * * *
60 //
61 // This material was prepared as an account of work sponsored by the
62 // United States Government. Neither the United States, nor the
63 // Department of Energy, nor the Department of Defense, nor any of
64 // their employees, nor any of their contractors, subcontractors, or
65 // their employees, makes any warranty, expressed or implied, or
66 // assumes any legal liability or responsibility for the accuracy,
67 // completeness, or usefulness of any information, apparatus, product,
68 // or process disclosed, or represents that its use would not infringe
69 // upon privately owned rights.
70
71 /*
72 Generated by VIF - experimental VIntage Fortran compiler.
73 VIF release 1.2.17
74 */
75 #if defined (__GNUC__)
76 #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
77 #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
78 #if (__GNUC__ >= 14)
79 #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
80 #pragma GCC diagnostic ignored "-Wimplicit-int"
81 #pragma GCC diagnostic ignored "-Wint-conversion"
82 #pragma GCC diagnostic ignored "-Wreturn-mismatch"
83 #endif
84 #else
85 #error VIF requires GCC
86 #endif
165 static CALLS __calls[__ncalls] = {
166 {"cdcor", 0}, // subroutine
167 {"cdcst", 0}, // subroutine
168 {"cdntl", 0}, // subroutine
169 {"cdntp", 0}, // subroutine
170 {"cdpsc", 0}, // subroutine
171 {"cdpst", 0}, // subroutine
172 {"cdriv1", 0}, // subroutine
173 {"cdriv2", 0}, // subroutine
174 {"cdriv3", 0}, // subroutine
175 {"cdscl", 0}, // subroutine
176 {"cdstp", 0}, // subroutine
177 {"cdzro", 0}, // subroutine
178 {"ddcor", 0}, // subroutine
179 {"ddcst", 0}, // subroutine
180 {"ddntl", 0}, // subroutine
181 {"ddntp", 0}, // subroutine
182 {"ddpsc", 0}, // subroutine
183 {"ddpst", 0}, // subroutine
184 {"ddriv1", 0}, // subroutine
185 {"ddriv2", 0}, // subroutine
186 {"ddriv3", 0}, // subroutine
187 {"ddscl", 0}, // subroutine
188 {"ddstp", 0}, // subroutine
189 {"ddzro", 0}, // subroutine
190 {"sdcor", 0}, // subroutine
191 {"sdcst", 0}, // subroutine
192 {"sdntl", 0}, // subroutine
193 {"sdntp", 0}, // subroutine
194 {"sdpsc", 0}, // subroutine
195 {"sdpst", 0}, // subroutine
196 {"sdriv1", 0}, // subroutine
197 {"sdriv2", 0}, // subroutine
198 {"sdriv3", 0}, // subroutine
199 {"sdscl", 0}, // subroutine
200 {"sdstp", 0}, // subroutine
201 {"sdzro", 0}, // subroutine
202 {NULL, 0}
203 };
204
206 _p_ impl_, int_4 _p_ ipvt_, int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_
207 nde_, int_4 _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_,
208 logical_4 _p_ evalfa_, complex_8 _p_ save1_, complex_8 _p_ save2_, complex_8 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_);
211 impl_, int_4 _p_ jtask_, int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4
212 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, complex_8 _p_ save1_, real_4 _p_ t_, real_4 _p_ uround_, int_4 (*_users)(),
213 complex_8 _p_ y_, complex_8 _p_ ywt_, real_4 _p_ h_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_, real_4 _p_
214 rc_, complex_8 _p_ yh_, complex_8 _p_ a_, logical_4 _p_ convrg_, real_4 _p_ el_, complex_8 _p_ fac_, logical_4 _p_ ier_
215 , int_4 _p_ ipvt_, int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rh_, real_4 _p_ rmax_, complex_8 _p_ save2_, real_4 _p_
216 tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ jstate_);
218 complex_8 _p_ yh_, complex_8 _p_ y_);
221 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4
222 _p_ nq_, complex_8 _p_ save2_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_
223 ywt_, real_4 _p_ uround_, int_4 _p_ nfe_, int_4 _p_ nje_, complex_8 _p_ a_, complex_8 _p_ dfdy_, complex_8 _p_ fac_,
224 logical_4 _p_ ier_, int_4 _p_ ipvt_, complex_8 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_);
226 mstate_, real_4 _p_ eps_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
228 mstate_, int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ mint_, complex_8 _p_ work_, int_4 _p_ lenw_,
229 int_4 _p_ iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_);
231 tout_, int_4 _p_ ntask_, int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ ierror_, int_4 _p_ mint_, int_4
232 _p_ miter_, int_4 _p_ impl_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ mxord_, real_4 _p_ hmax_, complex_8 _p_ work_,
233 int_4 _p_ lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_
234 mxstep_, real_4 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
236 real_4 _p_ rh_, complex_8 _p_ yh_);
238 ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_,
239 int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, complex_8 _p_ ywt_, real_4 _p_ uround_, real_4 (*_users)(), real_4 _p_
240 avgh_, real_4 _p_ avgord_, real_4 _p_ h_, real_4 _p_ hused_, int_4 _p_ jtask_, int_4 _p_ mntold_, int_4 _p_ mtrold_,
241 int_4 _p_ nfe_, int_4 _p_ nje_, int_4 _p_ nqused_, int_4 _p_ nstep_, real_4 _p_ t_, complex_8 _p_ y_, complex_8 _p_ yh_
242 , complex_8 _p_ a_, logical_4 _p_ convrg_, complex_8 _p_ dfdy_, real_4 _p_ el_, complex_8 _p_ fac_, real_4 _p_ hold_,
243 int_4 _p_ ipvt_, int_4 _p_ jstate_, int_4 _p_ jstepl_, int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rc_, real_4 _p_
244 rmax_, complex_8 _p_ save1_, complex_8 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_
245 mtrsv_, int_4 _p_ mxrdsv_);
247 real_4 _p_ re_, real_4 _p_ t_, complex_8 _p_ yh_, real_4 _p_ uround_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ fb_,
248 real_4 _p_ fc_, complex_8 _p_ y_);
250 impl_, int_4 _p_ ipvt_, int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_
251 , int_4 _p_ nq_, real_8 _p_ t_, int_4 (*_users)(), real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ ywt_, logical_4 _p_
252 evalfa_, real_8 _p_ save1_, real_8 _p_ save2_, real_8 _p_ a_, real_8 _p_ d_, int_4 _p_ jstate_);
255 impl_, int_4 _p_ jtask_, int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4
256 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, real_8 _p_ save1_, real_8 _p_ t_, real_8 _p_ uround_, int_4 (*_users)(), real_8
257 _p_ y_, real_8 _p_ ywt_, real_8 _p_ h_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_, real_8 _p_ rc_, real_8
258 _p_ yh_, real_8 _p_ a_, logical_4 _p_ convrg_, real_8 _p_ el_, real_8 _p_ fac_, logical_4 _p_ ier_, int_4 _p_ ipvt_,
259 int_4 _p_ nq_, int_4 _p_ nwait_, real_8 _p_ rh_, real_8 _p_ rmax_, real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_
260 , int_4 _p_ iswflg_, int_4 _p_ jstate_);
262 real_8 _p_ yh_, real_8 _p_ y_);
265 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4
266 _p_ nq_, real_8 _p_ save2_, real_8 _p_ t_, int_4 (*_users)(), real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ ywt_, real_8
267 _p_ uround_, int_4 _p_ nfe_, int_4 _p_ nje_, real_8 _p_ a_, real_8 _p_ dfdy_, real_8 _p_ fac_, logical_4 _p_ ier_,
268 int_4 _p_ ipvt_, real_8 _p_ save1_, int_4 _p_ iswflg_, real_8 _p_ bnd_, int_4 _p_ jstate_);
270 mstate_, real_8 _p_ eps_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
272 mstate_, int_4 _p_ nroot_, real_8 _p_ eps_, real_8 _p_ ewt_, int_4 _p_ mint_, real_8 _p_ work_, int_4 _p_ lenw_, int_4
273 _p_ iwork_, int_4 _p_ leniw_, real_8 (*_g)(), int_4 _p_ ierflg_);
275 , int_4 _p_ ntask_, int_4 _p_ nroot_, real_8 _p_ eps_, real_8 _p_ ewt_, int_4 _p_ ierror_, int_4 _p_ mint_, int_4 _p_
276 miter_, int_4 _p_ impl_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ mxord_, real_8 _p_ hmax_, real_8 _p_ work_, int_4 _p_
277 lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_,
278 real_8 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
280 real_8 _p_ rh_, real_8 _p_ yh_);
282 ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_,
283 int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, real_8 _p_ ywt_, real_8 _p_ uround_, real_4 (*_users)(), real_8 _p_ avgh_,
284 real_8 _p_ avgord_, real_8 _p_ h_, real_8 _p_ hused_, int_4 _p_ jtask_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_
285 nfe_, int_4 _p_ nje_, int_4 _p_ nqused_, int_4 _p_ nstep_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ a_,
286 logical_4 _p_ convrg_, real_8 _p_ dfdy_, real_8 _p_ el_, real_8 _p_ fac_, real_8 _p_ hold_, int_4 _p_ ipvt_, int_4 _p_
287 jstate_, int_4 _p_ jstepl_, int_4 _p_ nq_, int_4 _p_ nwait_, real_8 _p_ rc_, real_8 _p_ rmax_, real_8 _p_ save1_,
288 real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_);
290 real_8 _p_ re_, real_8 _p_ t_, real_8 _p_ yh_, real_8 _p_ uround_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ fb_, real_8
291 _p_ fc_, real_8 _p_ y_);
293 impl_, int_4 _p_ ipvt_, int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_
294 , int_4 _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), real_4 _p_ y_, real_4 _p_ yh_, real_4 _p_ ywt_, logical_4 _p_
295 evalfa_, real_4 _p_ save1_, real_4 _p_ save2_, real_4 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_);
298 impl_, int_4 _p_ jtask_, int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4
299 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, real_4 _p_ save1_, real_4 _p_ t_, real_4 _p_ uround_, int_4 (*_users)(), real_4
300 _p_ y_, real_4 _p_ ywt_, real_4 _p_ h_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_, real_4 _p_ rc_, real_4
301 _p_ yh_, real_4 _p_ a_, logical_4 _p_ convrg_, real_4 _p_ el_, real_4 _p_ fac_, logical_4 _p_ ier_, int_4 _p_ ipvt_,
302 int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rh_, real_4 _p_ rmax_, real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_
303 , int_4 _p_ iswflg_, int_4 _p_ jstate_);
305 real_4 _p_ yh_, real_4 _p_ y_);
308 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4
309 _p_ nq_, real_4 _p_ save2_, real_4 _p_ t_, int_4 (*_users)(), real_4 _p_ y_, real_4 _p_ yh_, real_4 _p_ ywt_, real_4
310 _p_ uround_, int_4 _p_ nfe_, int_4 _p_ nje_, real_4 _p_ a_, real_4 _p_ dfdy_, real_4 _p_ fac_, logical_4 _p_ ier_,
311 int_4 _p_ ipvt_, real_4 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_);
313 mstate_, real_4 _p_ eps_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
315 mstate_, int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ mint_, real_4 _p_ work_, int_4 _p_ lenw_, int_4
316 _p_ iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_);
318 , int_4 _p_ ntask_, int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ ierror_, int_4 _p_ mint_, int_4 _p_
319 miter_, int_4 _p_ impl_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ mxord_, real_4 _p_ hmax_, real_4 _p_ work_, int_4 _p_
320 lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_,
321 real_4 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
323 real_4 _p_ rh_, real_4 _p_ yh_);
325 ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_,
326 int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, real_4 _p_ ywt_, real_4 _p_ uround_, real_4 (*_users)(), real_4 _p_ avgh_,
327 real_4 _p_ avgord_, real_4 _p_ h_, real_4 _p_ hused_, int_4 _p_ jtask_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_
328 nfe_, int_4 _p_ nje_, int_4 _p_ nqused_, int_4 _p_ nstep_, real_4 _p_ t_, real_4 _p_ y_, real_4 _p_ yh_, real_4 _p_ a_,
329 logical_4 _p_ convrg_, real_4 _p_ dfdy_, real_4 _p_ el_, real_4 _p_ fac_, real_4 _p_ hold_, int_4 _p_ ipvt_, int_4 _p_
330 jstate_, int_4 _p_ jstepl_, int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rc_, real_4 _p_ rmax_, real_4 _p_ save1_,
331 real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_);
333 real_4 _p_ re_, real_4 _p_ t_, real_4 _p_ yh_, real_4 _p_ uround_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ fb_, real_4
334 _p_ fc_, real_4 _p_ y_);
335 int_4 _cdcor (complex_8 _p_ dfdy_, real_4 _p_ el_, int_4 (*_fa)(), real_4 _p_ h_, int_4 _p_ ierror_, int_4 _p_ impl_,
336 int_4 _p_ ipvt_, int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4
337 _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, logical_4 _p_
338 evalfa_, complex_8 _p_ save1_, complex_8 _p_ save2_, complex_8 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_)
339 { // ** body not listed **
544 }
545
546 int_4 _cdcst (int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ iswflg_, real_4 _p_ el_, real_4 _p_ tq_)
547 { // ** body not listed **
633 }
634
635 int_4 _cdntl (real_4 _p_ eps_, int_4 (*_f)(), int_4 (*_fa)(), real_4 _p_ hmax_, real_4 _p_ hold_, int_4 _p_ impl_,
636 int_4 _p_ jtask_, int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_
637 , int_4 _p_ n_, int_4 _p_ nde_, complex_8 _p_ save1_, real_4 _p_ t_, real_4 _p_ uround_, int_4 (*_users)(), complex_8
638 _p_ y_, complex_8 _p_ ywt_, real_4 _p_ h_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_, real_4 _p_ rc_,
639 complex_8 _p_ yh_, complex_8 _p_ a_, logical_4 _p_ convrg_, real_4 _p_ el_, complex_8 _p_ fac_, logical_4 _p_ ier_,
640 int_4 _p_ ipvt_, int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rh_, real_4 _p_ rmax_, complex_8 _p_ save2_, real_4 _p_
641 tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ jstate_)
642 { // ** body not listed **
790 }
791
792 int_4 _cdntp (real_4 _p_ h_, int_4 _p_ k_, int_4 _p_ n_, int_4 _p_ nq_, real_4 _p_ t_, real_4 _p_ tout_, complex_8 _p_
793 yh_, complex_8 _p_ y_)
794 { // ** body not listed **
835 }
836
837 int_4 _cdpsc (int_4 _p_ ksgn_, int_4 _p_ n_, int_4 _p_ nq_, complex_8 _p_ yh_)
838 { // ** body not listed **
862 }
863
864 int_4 _cdpst (real_4 _p_ el_, int_4 (*_f)(), int_4 (*_fa)(), real_4 _p_ h_, int_4 _p_ impl_, int_4 (*_jacobn)(), int_4
865 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4 _p_ nq_, complex_8 _p_
866 save2_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, real_4 _p_ uround_,
867 int_4 _p_ nfe_, int_4 _p_ nje_, complex_8 _p_ a_, complex_8 _p_ dfdy_, complex_8 _p_ fac_, logical_4 _p_ ier_, int_4
868 _p_ ipvt_, complex_8 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_)
869 { // ** body not listed **
1172 }
1173
1174 int_4 _cdriv1 (int_4 _p_ n_, real_4 _p_ t_, complex_8 _p_ y_, real_4 (*_f)(), real_4 _p_ tout_, int_4 _p_ mstate_,
1175 real_4 _p_ eps_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
1176 { // ** body not listed **
1370 }
1371
1372 int_4 _cdriv2 (int_4 _p_ n_, real_4 _p_ t_, complex_8 _p_ y_, real_4 (*_f)(), real_4 _p_ tout_, int_4 _p_ mstate_,
1373 int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ mint_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_
1374 iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_)
1375 { // ** body not listed **
1518 }
1519
1520 int_4 _cdriv3 (int_4 _p_ n_, real_4 _p_ t_, complex_8 _p_ y_, int_4 (*_f)(), int_4 _p_ nstate_, real_4 _p_ tout_, int_4
1521 _p_ ntask_, int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ ierror_, int_4 _p_ mint_, int_4 _p_ miter_,
1522 int_4 _p_ impl_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ mxord_, real_4 _p_ hmax_, complex_8 _p_ work_, int_4 _p_ lenw_
1523 , int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_4
1524 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
1525 { // ** body not listed **
3355 }
3356
3357 int_4 _cdscl (real_4 _p_ hmax_, int_4 _p_ n_, int_4 _p_ nq_, real_4 _p_ rmax_, real_4 _p_ h_, real_4 _p_ rc_, real_4
3358 _p_ rh_, complex_8 _p_ yh_)
3359 { // ** body not listed **
3380 }
3381
3382 int_4 _cdstp (real_4 _p_ eps_, int_4 (*_f)(), real_4 (*_fa)(), real_4 _p_ hmax_, int_4 _p_ impl_, int_4 _p_ ierror_,
3383 int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_
3384 mu_, int_4 _p_ n_, int_4 _p_ nde_, complex_8 _p_ ywt_, real_4 _p_ uround_, real_4 (*_users)(), real_4 _p_ avgh_, real_4
3385 _p_ avgord_, real_4 _p_ h_, real_4 _p_ hused_, int_4 _p_ jtask_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_,
3386 int_4 _p_ nje_, int_4 _p_ nqused_, int_4 _p_ nstep_, real_4 _p_ t_, complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_
3387 a_, logical_4 _p_ convrg_, complex_8 _p_ dfdy_, real_4 _p_ el_, complex_8 _p_ fac_, real_4 _p_ hold_, int_4 _p_ ipvt_,
3388 int_4 _p_ jstate_, int_4 _p_ jstepl_, int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rc_, real_4 _p_ rmax_, complex_8 _p_
3389 save1_, complex_8 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_
3390 )
3391 { // ** body not listed **
3798 }
3799
3800 int_4 _cdzro (real_4 _p_ ae_, real_4 (*_f)(), real_4 _p_ h_, int_4 _p_ n_, int_4 _p_ nq_, int_4 _p_ iroot_, real_4 _p_
3801 re_, real_4 _p_ t_, complex_8 _p_ yh_, real_4 _p_ uround_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ fb_, real_4 _p_ fc_
3802 , complex_8 _p_ y_)
3803 { // ** body not listed **
3875 }
3876
3877 int_4 _ddcor (real_8 _p_ dfdy_, real_8 _p_ el_, int_4 (*_fa)(), real_8 _p_ h_, int_4 _p_ ierror_, int_4 _p_ impl_,
3878 int_4 _p_ ipvt_, int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4
3879 _p_ nq_, real_8 _p_ t_, int_4 (*_users)(), real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ ywt_, logical_4 _p_ evalfa_,
3880 real_8 _p_ save1_, real_8 _p_ save2_, real_8 _p_ a_, real_8 _p_ d_, int_4 _p_ jstate_)
3881 { // ** body not listed **
4086 }
4087
4088 int_4 _ddcst (int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ iswflg_, real_8 _p_ el_, real_8 _p_ tq_)
4089 { // ** body not listed **
4175 }
4176
4177 int_4 _ddntl (real_8 _p_ eps_, int_4 (*_f)(), int_4 (*_fa)(), real_8 _p_ hmax_, real_8 _p_ hold_, int_4 _p_ impl_,
4178 int_4 _p_ jtask_, int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_
4179 , int_4 _p_ n_, int_4 _p_ nde_, real_8 _p_ save1_, real_8 _p_ t_, real_8 _p_ uround_, int_4 (*_users)(), real_8 _p_ y_,
4180 real_8 _p_ ywt_, real_8 _p_ h_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_, real_8 _p_ rc_, real_8 _p_ yh_,
4181 real_8 _p_ a_, logical_4 _p_ convrg_, real_8 _p_ el_, real_8 _p_ fac_, logical_4 _p_ ier_, int_4 _p_ ipvt_, int_4 _p_
4182 nq_, int_4 _p_ nwait_, real_8 _p_ rh_, real_8 _p_ rmax_, real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4
4183 _p_ iswflg_, int_4 _p_ jstate_)
4184 { // ** body not listed **
4332 }
4333
4334 int_4 _ddntp (real_8 _p_ h_, int_4 _p_ k_, int_4 _p_ n_, int_4 _p_ nq_, real_8 _p_ t_, real_8 _p_ tout_, real_8 _p_ yh_
4335 , real_8 _p_ y_)
4336 { // ** body not listed **
4377 }
4378
4379 int_4 _ddpsc (int_4 _p_ ksgn_, int_4 _p_ n_, int_4 _p_ nq_, real_8 _p_ yh_)
4380 { // ** body not listed **
4404 }
4405
4406 int_4 _ddpst (real_8 _p_ el_, int_4 (*_f)(), int_4 (*_fa)(), real_8 _p_ h_, int_4 _p_ impl_, int_4 (*_jacobn)(), int_4
4407 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4 _p_ nq_, real_8 _p_
4408 save2_, real_8 _p_ t_, int_4 (*_users)(), real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ ywt_, real_8 _p_ uround_, int_4 _p_
4409 nfe_, int_4 _p_ nje_, real_8 _p_ a_, real_8 _p_ dfdy_, real_8 _p_ fac_, logical_4 _p_ ier_, int_4 _p_ ipvt_, real_8 _p_
4410 save1_, int_4 _p_ iswflg_, real_8 _p_ bnd_, int_4 _p_ jstate_)
4411 { // ** body not listed **
4719 }
4720
4721 int_4 _ddriv1 (int_4 _p_ n_, real_8 _p_ t_, real_8 _p_ y_, real_4 (*_f)(), real_8 _p_ tout_, int_4 _p_ mstate_, real_8
4722 _p_ eps_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
4723 { // ** body not listed **
4917 }
4918
4919 int_4 _ddriv2 (int_4 _p_ n_, real_8 _p_ t_, real_8 _p_ y_, real_4 (*_f)(), real_8 _p_ tout_, int_4 _p_ mstate_, int_4
4920 _p_ nroot_, real_8 _p_ eps_, real_8 _p_ ewt_, int_4 _p_ mint_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ iwork_,
4921 int_4 _p_ leniw_, real_8 (*_g)(), int_4 _p_ ierflg_)
4922 { // ** body not listed **
5065 }
5066
5067 int_4 _ddriv3 (int_4 _p_ n_, real_8 _p_ t_, real_8 _p_ y_, int_4 (*_f)(), int_4 _p_ nstate_, real_8 _p_ tout_, int_4
5068 _p_ ntask_, int_4 _p_ nroot_, real_8 _p_ eps_, real_8 _p_ ewt_, int_4 _p_ ierror_, int_4 _p_ mint_, int_4 _p_ miter_,
5069 int_4 _p_ impl_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ mxord_, real_8 _p_ hmax_, real_8 _p_ work_, int_4 _p_ lenw_,
5070 int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_8
5071 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
5072 { // ** body not listed **
6841 }
6842
6843 int_4 _ddscl (real_8 _p_ hmax_, int_4 _p_ n_, int_4 _p_ nq_, real_8 _p_ rmax_, real_8 _p_ h_, real_8 _p_ rc_, real_8
6844 _p_ rh_, real_8 _p_ yh_)
6845 { // ** body not listed **
6866 }
6867
6868 int_4 _ddstp (real_8 _p_ eps_, int_4 (*_f)(), real_4 (*_fa)(), real_8 _p_ hmax_, int_4 _p_ impl_, int_4 _p_ ierror_,
6869 int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_
6870 mu_, int_4 _p_ n_, int_4 _p_ nde_, real_8 _p_ ywt_, real_8 _p_ uround_, real_4 (*_users)(), real_8 _p_ avgh_, real_8
6871 _p_ avgord_, real_8 _p_ h_, real_8 _p_ hused_, int_4 _p_ jtask_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_,
6872 int_4 _p_ nje_, int_4 _p_ nqused_, int_4 _p_ nstep_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ a_,
6873 logical_4 _p_ convrg_, real_8 _p_ dfdy_, real_8 _p_ el_, real_8 _p_ fac_, real_8 _p_ hold_, int_4 _p_ ipvt_, int_4 _p_
6874 jstate_, int_4 _p_ jstepl_, int_4 _p_ nq_, int_4 _p_ nwait_, real_8 _p_ rc_, real_8 _p_ rmax_, real_8 _p_ save1_,
6875 real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_)
6876 { // ** body not listed **
7283 }
7284
7285 int_4 _ddzro (real_8 _p_ ae_, real_8 (*_f)(), real_8 _p_ h_, int_4 _p_ n_, int_4 _p_ nq_, int_4 _p_ iroot_, real_8 _p_
7286 re_, real_8 _p_ t_, real_8 _p_ yh_, real_8 _p_ uround_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ fb_, real_8 _p_ fc_,
7287 real_8 _p_ y_)
7288 { // ** body not listed **
7360 }
7361
7362 int_4 _sdcor (real_4 _p_ dfdy_, real_4 _p_ el_, int_4 (*_fa)(), real_4 _p_ h_, int_4 _p_ ierror_, int_4 _p_ impl_,
7363 int_4 _p_ ipvt_, int_4 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4
7364 _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), real_4 _p_ y_, real_4 _p_ yh_, real_4 _p_ ywt_, logical_4 _p_ evalfa_,
7365 real_4 _p_ save1_, real_4 _p_ save2_, real_4 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_)
7366 { // ** body not listed **
7571 }
7572
7573 int_4 _sdcst (int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ iswflg_, real_4 _p_ el_, real_4 _p_ tq_)
7574 { // ** body not listed **
7660 }
7661
7662 int_4 _sdntl (real_4 _p_ eps_, int_4 (*_f)(), int_4 (*_fa)(), real_4 _p_ hmax_, real_4 _p_ hold_, int_4 _p_ impl_,
7663 int_4 _p_ jtask_, int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_
7664 , int_4 _p_ n_, int_4 _p_ nde_, real_4 _p_ save1_, real_4 _p_ t_, real_4 _p_ uround_, int_4 (*_users)(), real_4 _p_ y_,
7665 real_4 _p_ ywt_, real_4 _p_ h_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_, real_4 _p_ rc_, real_4 _p_ yh_,
7666 real_4 _p_ a_, logical_4 _p_ convrg_, real_4 _p_ el_, real_4 _p_ fac_, logical_4 _p_ ier_, int_4 _p_ ipvt_, int_4 _p_
7667 nq_, int_4 _p_ nwait_, real_4 _p_ rh_, real_4 _p_ rmax_, real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4
7668 _p_ iswflg_, int_4 _p_ jstate_)
7669 { // ** body not listed **
7817 }
7818
7819 int_4 _sdntp (real_4 _p_ h_, int_4 _p_ k_, int_4 _p_ n_, int_4 _p_ nq_, real_4 _p_ t_, real_4 _p_ tout_, real_4 _p_ yh_
7820 , real_4 _p_ y_)
7821 { // ** body not listed **
7862 }
7863
7864 int_4 _sdpsc (int_4 _p_ ksgn_, int_4 _p_ n_, int_4 _p_ nq_, real_4 _p_ yh_)
7865 { // ** body not listed **
7889 }
7890
7891 int_4 _sdpst (real_4 _p_ el_, int_4 (*_f)(), int_4 (*_fa)(), real_4 _p_ h_, int_4 _p_ impl_, int_4 (*_jacobn)(), int_4
7892 _p_ matdim_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ n_, int_4 _p_ nde_, int_4 _p_ nq_, real_4 _p_
7893 save2_, real_4 _p_ t_, int_4 (*_users)(), real_4 _p_ y_, real_4 _p_ yh_, real_4 _p_ ywt_, real_4 _p_ uround_, int_4 _p_
7894 nfe_, int_4 _p_ nje_, real_4 _p_ a_, real_4 _p_ dfdy_, real_4 _p_ fac_, logical_4 _p_ ier_, int_4 _p_ ipvt_, real_4 _p_
7895 save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_)
7896 { // ** body not listed **
8204 }
8205
8206 int_4 _sdriv1 (int_4 _p_ n_, real_4 _p_ t_, real_4 _p_ y_, real_4 (*_f)(), real_4 _p_ tout_, int_4 _p_ mstate_, real_4
8207 _p_ eps_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
8208 { // ** body not listed **
8402 }
8403
8404 int_4 _sdriv2 (int_4 _p_ n_, real_4 _p_ t_, real_4 _p_ y_, real_4 (*_f)(), real_4 _p_ tout_, int_4 _p_ mstate_, int_4
8405 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ mint_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ iwork_,
8406 int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_)
8407 { // ** body not listed **
8550 }
8551
8552 int_4 _sdriv3 (int_4 _p_ n_, real_4 _p_ t_, real_4 _p_ y_, int_4 (*_f)(), int_4 _p_ nstate_, real_4 _p_ tout_, int_4
8553 _p_ ntask_, int_4 _p_ nroot_, real_4 _p_ eps_, real_4 _p_ ewt_, int_4 _p_ ierror_, int_4 _p_ mint_, int_4 _p_ miter_,
8554 int_4 _p_ impl_, int_4 _p_ ml_, int_4 _p_ mu_, int_4 _p_ mxord_, real_4 _p_ hmax_, real_4 _p_ work_, int_4 _p_ lenw_,
8555 int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_4
8556 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
8557 { // ** body not listed **
10326 }
10327
10328 int_4 _sdscl (real_4 _p_ hmax_, int_4 _p_ n_, int_4 _p_ nq_, real_4 _p_ rmax_, real_4 _p_ h_, real_4 _p_ rc_, real_4
10329 _p_ rh_, real_4 _p_ yh_)
10330 { // ** body not listed **
10351 }
10352
10353 int_4 _sdstp (real_4 _p_ eps_, int_4 (*_f)(), real_4 (*_fa)(), real_4 _p_ hmax_, int_4 _p_ impl_, int_4 _p_ ierror_,
10354 int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_
10355 mu_, int_4 _p_ n_, int_4 _p_ nde_, real_4 _p_ ywt_, real_4 _p_ uround_, real_4 (*_users)(), real_4 _p_ avgh_, real_4
10356 _p_ avgord_, real_4 _p_ h_, real_4 _p_ hused_, int_4 _p_ jtask_, int_4 _p_ mntold_, int_4 _p_ mtrold_, int_4 _p_ nfe_,
10357 int_4 _p_ nje_, int_4 _p_ nqused_, int_4 _p_ nstep_, real_4 _p_ t_, real_4 _p_ y_, real_4 _p_ yh_, real_4 _p_ a_,
10358 logical_4 _p_ convrg_, real_4 _p_ dfdy_, real_4 _p_ el_, real_4 _p_ fac_, real_4 _p_ hold_, int_4 _p_ ipvt_, int_4 _p_
10359 jstate_, int_4 _p_ jstepl_, int_4 _p_ nq_, int_4 _p_ nwait_, real_4 _p_ rc_, real_4 _p_ rmax_, real_4 _p_ save1_,
10360 real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_)
10361 { // ** body not listed **
10768 }
10769
10770 int_4 _sdzro (real_4 _p_ ae_, real_4 (*_f)(), real_4 _p_ h_, int_4 _p_ n_, int_4 _p_ nq_, int_4 _p_ iroot_, real_4 _p_
10771 re_, real_4 _p_ t_, real_4 _p_ yh_, real_4 _p_ uround_, real_4 _p_ b_, real_4 _p_ c_, real_4 _p_ fb_, real_4 _p_ fc_,
10772 real_4 _p_ y_)
10773 { // ** body not listed **
10845 }
10846
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|