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  //! 24 subprograms from SLATEC-SDRIVE
    24  
    25  // Compiled from Fortran source code by VIF.
    26  // Selected subprograms are:
    27  //
    28  // CDCOR  CDRIV1 CDSTP  DDNTL  DDRIV2 DDSTP  SDCOR  SDPST  SDRIV2 SDSTP
    29  // CDNTL  CDRIV2 CDZRO  DDPST  DDRIV3 DDZRO  SDNTL  SDRIV1 SDRIV3 SDZRO
    30  // CDPST  CDRIV3 DDCOR  DDRIV1
    31  
    32  // SDRIVE routines from SLATEC.
    33  // 
    34  // Source: netlib.org/slatec/
    35  
    36  // The license for SLATEC Fortran source code is:
    37  //
    38  // The SLATEC Common Mathematical Library was developed at
    39  // US government research laboratories and is in the public domain.
    40  // 
    41  // Repository: http://www.netlib.org/slatec/
    42  // 
    43  // The SLATEC common mathematical library is issued by the following
    44  // 
    45  //         Air Force Weapons Laboratory, Albuquerque
    46  //         Lawrence Livermore National Laboratory, Livermore
    47  //         Los Alamos National Laboratory, Los Alamos
    48  //         National Institute of Standards and Technology, Washington
    49  //         National Energy Research Supercomputer Center, Livermore
    50  //         Oak Ridge National Laboratory, Oak Ridge
    51  //         Sandia National Laboratories, Albuquerque
    52  //         Sandia National Laboratories, Livermore
    53  // 
    54  // All questions concerning the distribution of the library should be
    55  // directed to the NATIONAL ENERGY SOFTWARE CENTER, 9700 Cass Ave.,
    56  // Argonne, Illinois  60439, and not to the authors of the subprograms.
    57  // 
    58  //                  * * * * * Notice * * * * *
    59  // 
    60  // This material was prepared as an account of work sponsored by the
    61  // United States Government.  Neither the United States, nor the
    62  // Department of Energy, nor the Department of Defense, nor any of
    63  // their employees, nor any of their contractors, subcontractors, or
    64  // their employees, makes any warranty, expressed or implied, or
    65  // assumes any legal liability or responsibility for the accuracy,
    66  // completeness, or usefulness of any information, apparatus, product,
    67  // or process disclosed, or represents that its use would not infringe
    68  // upon privately owned rights.
    69  
    70  /*
    71  Generated by VIF - experimental VIntage Fortran compiler.
    72  VIF release 1.1.7
    73  */
   152  static CALLS __calls[__ncalls] = {
   153    {"cdcor", 0}, // subroutine
   154    {"cdntl", 0}, // subroutine
   155    {"cdpst", 0}, // subroutine
   156    {"cdriv1", 0}, // subroutine
   157    {"cdriv2", 0}, // subroutine
   158    {"cdriv3", 0}, // subroutine
   159    {"cdstp", 0}, // subroutine
   160    {"cdzro", 0}, // subroutine
   161    {"ddcor", 0}, // subroutine
   162    {"ddntl", 0}, // subroutine
   163    {"ddpst", 0}, // subroutine
   164    {"ddriv1", 0}, // subroutine
   165    {"ddriv2", 0}, // subroutine
   166    {"ddriv3", 0}, // subroutine
   167    {"ddstp", 0}, // subroutine
   168    {"ddzro", 0}, // subroutine
   169    {"sdcor", 0}, // subroutine
   170    {"sdntl", 0}, // subroutine
   171    {"sdpst", 0}, // subroutine
   172    {"sdriv1", 0}, // subroutine
   173    {"sdriv2", 0}, // subroutine
   174    {"sdriv3", 0}, // subroutine
   175    {"sdstp", 0}, // subroutine
   176    {"sdzro", 0}, // subroutine
   177    {NULL, 0}
   178  };
   180  _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_ 
   181  nde_, int_4 _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, 
   182  logical_4 _p_ evalfa_, complex_8 _p_ save1_, complex_8 _p_ save2_, complex_8 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_);
   184  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 
   185  _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)(), 
   186  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_ 
   187  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_
   188  , 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_ 
   189  tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ jstate_);
   191  (*_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 
   192  _p_ nq_, complex_8 _p_ save2_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ 
   193  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_, 
   194  logical_4 _p_ ier_, int_4 _p_ ipvt_, complex_8 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_);
   196  mstate_, real_4 _p_ eps_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
   198  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_, 
   199  int_4 _p_ iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_);
   201  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 
   202  _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_, 
   203  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_ 
   204  mxstep_, real_4 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
   206  ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, 
   207  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_ 
   208  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_, 
   209  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_
   210  , complex_8 _p_ a_, logical_4 _p_ convrg_, complex_8 _p_ dfdy_, real_4 _p_ el_, complex_8 _p_ fac_, real_4 _p_ hold_, 
   211  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_ 
   212  rmax_, complex_8 _p_ save1_, complex_8 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ 
   213  mtrsv_, int_4 _p_ mxrdsv_);
   215  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_, 
   216  real_4 _p_ fc_, complex_8 _p_ y_);
   218  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_
   219  , 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_ 
   220  evalfa_, real_8 _p_ save1_, real_8 _p_ save2_, real_8 _p_ a_, real_8 _p_ d_, int_4 _p_ jstate_);
   222  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 
   223  _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 
   224  _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 
   225  _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_, 
   226  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_
   227  , int_4 _p_ iswflg_, int_4 _p_ jstate_);
   229  (*_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 
   230  _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 
   231  _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_, 
   232  int_4 _p_ ipvt_, real_8 _p_ save1_, int_4 _p_ iswflg_, real_8 _p_ bnd_, int_4 _p_ jstate_);
   234  mstate_, real_8 _p_ eps_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
   236  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 
   237  _p_ iwork_, int_4 _p_ leniw_, real_8 (*_g)(), int_4 _p_ ierflg_);
   239  , 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_ 
   240  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_ 
   241  lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, 
   242  real_8 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
   244  ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, 
   245  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_, 
   246  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_ 
   247  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_, 
   248  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_ 
   249  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_, 
   250  real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_);
   252  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 
   253  _p_ fc_, real_8 _p_ y_);
   255  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_
   256  , 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_ 
   257  evalfa_, real_4 _p_ save1_, real_4 _p_ save2_, real_4 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_);
   259  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 
   260  _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 
   261  _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 
   262  _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_, 
   263  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_
   264  , int_4 _p_ iswflg_, int_4 _p_ jstate_);
   266  (*_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 
   267  _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 
   268  _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_, 
   269  int_4 _p_ ipvt_, real_4 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_);
   271  mstate_, real_4 _p_ eps_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
   273  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 
   274  _p_ iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_);
   276  , 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_ 
   277  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_ 
   278  lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, 
   279  real_4 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
   281  ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, 
   282  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_, 
   283  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_ 
   284  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_, 
   285  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_ 
   286  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_, 
   287  real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_);
   289  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 
   290  _p_ fc_, real_4 _p_ y_);
   291  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_, 
   292  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 
   293  _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, logical_4 _p_ 
   294  evalfa_, complex_8 _p_ save1_, complex_8 _p_ save2_, complex_8 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_)
   295  { // ** body not listed **
   500  }
   501  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_, 
   502  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_
   503  , 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 
   504  _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_, 
   505  complex_8 _p_ yh_, complex_8 _p_ a_, logical_4 _p_ convrg_, real_4 _p_ el_, complex_8 _p_ fac_, logical_4 _p_ ier_, 
   506  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_ 
   507  tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ jstate_)
   508  { // ** body not listed **
   656  }
   657  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 
   658  _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_ 
   659  save2_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, real_4 _p_ uround_, 
   660  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 
   661  _p_ ipvt_, complex_8 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_)
   662  { // ** body not listed **
   965  }
   966  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_, 
   967  real_4 _p_ eps_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
   968  { // ** body not listed **
  1162  }
  1163  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_, 
  1164  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_ 
  1165  iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_)
  1166  { // ** body not listed **
  1309  }
  1310  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 
  1311  _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_, 
  1312  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_
  1313  , int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_4 
  1314  (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
  1315  { // ** body not listed **
  3145  }
  3146  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_, 
  3147  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_ 
  3148  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 
  3149  _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_, 
  3150  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_ 
  3151  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_, 
  3152  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_ 
  3153  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_
  3154  )
  3155  { // ** body not listed **
  3562  }
  3563  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_ 
  3564  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_
  3565  , complex_8 _p_ y_)
  3566  { // ** body not listed **
  3638  }
  3639  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_, 
  3640  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 
  3641  _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_, 
  3642  real_8 _p_ save1_, real_8 _p_ save2_, real_8 _p_ a_, real_8 _p_ d_, int_4 _p_ jstate_)
  3643  { // ** body not listed **
  3848  }
  3849  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_, 
  3850  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_
  3851  , 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_, 
  3852  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_, 
  3853  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_ 
  3854  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 
  3855  _p_ iswflg_, int_4 _p_ jstate_)
  3856  { // ** body not listed **
  4004  }
  4005  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 
  4006  _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_ 
  4007  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_ 
  4008  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_ 
  4009  save1_, int_4 _p_ iswflg_, real_8 _p_ bnd_, int_4 _p_ jstate_)
  4010  { // ** body not listed **
  4318  }
  4319  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 
  4320  _p_ eps_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
  4321  { // ** body not listed **
  4515  }
  4516  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 
  4517  _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_, 
  4518  int_4 _p_ leniw_, real_8 (*_g)(), int_4 _p_ ierflg_)
  4519  { // ** body not listed **
  4662  }
  4663  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 
  4664  _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_, 
  4665  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_, 
  4666  int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_8 
  4667  (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
  4668  { // ** body not listed **
  6437  }
  6438  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_, 
  6439  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_ 
  6440  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 
  6441  _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_, 
  6442  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_, 
  6443  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_ 
  6444  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_, 
  6445  real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_)
  6446  { // ** body not listed **
  6853  }
  6854  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_ 
  6855  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_, 
  6856  real_8 _p_ y_)
  6857  { // ** body not listed **
  6929  }
  6930  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_, 
  6931  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 
  6932  _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_, 
  6933  real_4 _p_ save1_, real_4 _p_ save2_, real_4 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_)
  6934  { // ** body not listed **
  7139  }
  7140  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_, 
  7141  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_
  7142  , 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_, 
  7143  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_, 
  7144  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_ 
  7145  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 
  7146  _p_ iswflg_, int_4 _p_ jstate_)
  7147  { // ** body not listed **
  7295  }
  7296  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 
  7297  _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_ 
  7298  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_ 
  7299  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_ 
  7300  save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_)
  7301  { // ** body not listed **
  7609  }
  7610  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 
  7611  _p_ eps_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
  7612  { // ** body not listed **
  7806  }
  7807  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 
  7808  _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_, 
  7809  int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_)
  7810  { // ** body not listed **
  7953  }
  7954  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 
  7955  _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_, 
  7956  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_, 
  7957  int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_4 
  7958  (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
  7959  { // ** body not listed **
  9728  }
  9729  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_, 
  9730  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_ 
  9731  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 
  9732  _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_, 
  9733  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_, 
  9734  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_ 
  9735  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_, 
  9736  real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_)
  9737  { // ** body not listed **
 10144  }
 10145  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_ 
 10146  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_, 
 10147  real_4 _p_ y_)
 10148  { // ** body not listed **
 10220  }


© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)