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.12
    74  */
   153  static CALLS __calls[__ncalls] = {
   154    {"cdcor", 0}, // subroutine
   155    {"cdcst", 0}, // subroutine
   156    {"cdntl", 0}, // subroutine
   157    {"cdntp", 0}, // subroutine
   158    {"cdpsc", 0}, // subroutine
   159    {"cdpst", 0}, // subroutine
   160    {"cdriv1", 0}, // subroutine
   161    {"cdriv2", 0}, // subroutine
   162    {"cdriv3", 0}, // subroutine
   163    {"cdscl", 0}, // subroutine
   164    {"cdstp", 0}, // subroutine
   165    {"cdzro", 0}, // subroutine
   166    {"ddcor", 0}, // subroutine
   167    {"ddcst", 0}, // subroutine
   168    {"ddntl", 0}, // subroutine
   169    {"ddntp", 0}, // subroutine
   170    {"ddpsc", 0}, // subroutine
   171    {"ddpst", 0}, // subroutine
   172    {"ddriv1", 0}, // subroutine
   173    {"ddriv2", 0}, // subroutine
   174    {"ddriv3", 0}, // subroutine
   175    {"ddscl", 0}, // subroutine
   176    {"ddstp", 0}, // subroutine
   177    {"ddzro", 0}, // subroutine
   178    {"sdcor", 0}, // subroutine
   179    {"sdcst", 0}, // subroutine
   180    {"sdntl", 0}, // subroutine
   181    {"sdntp", 0}, // subroutine
   182    {"sdpsc", 0}, // subroutine
   183    {"sdpst", 0}, // subroutine
   184    {"sdriv1", 0}, // subroutine
   185    {"sdriv2", 0}, // subroutine
   186    {"sdriv3", 0}, // subroutine
   187    {"sdscl", 0}, // subroutine
   188    {"sdstp", 0}, // subroutine
   189    {"sdzro", 0}, // subroutine
   190    {NULL, 0}
   191  };
   192  
   194  _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_ 
   195  nde_, int_4 _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, 
   196  logical_4 _p_ evalfa_, complex_8 _p_ save1_, complex_8 _p_ save2_, complex_8 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_);
   199  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 
   200  _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)(), 
   201  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_ 
   202  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_
   203  , 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_ 
   204  tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ jstate_);
   206  complex_8 _p_ yh_, complex_8 _p_ y_);
   209  (*_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 
   210  _p_ nq_, complex_8 _p_ save2_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ 
   211  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_, 
   212  logical_4 _p_ ier_, int_4 _p_ ipvt_, complex_8 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_);
   214  mstate_, real_4 _p_ eps_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
   216  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_, 
   217  int_4 _p_ iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_);
   219  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 
   220  _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_, 
   221  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_ 
   222  mxstep_, real_4 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
   224  real_4 _p_ rh_, complex_8 _p_ yh_);
   226  ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, 
   227  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_ 
   228  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_, 
   229  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_
   230  , complex_8 _p_ a_, logical_4 _p_ convrg_, complex_8 _p_ dfdy_, real_4 _p_ el_, complex_8 _p_ fac_, real_4 _p_ hold_, 
   231  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_ 
   232  rmax_, complex_8 _p_ save1_, complex_8 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ 
   233  mtrsv_, int_4 _p_ mxrdsv_);
   235  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_, 
   236  real_4 _p_ fc_, complex_8 _p_ y_);
   238  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_
   239  , 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_ 
   240  evalfa_, real_8 _p_ save1_, real_8 _p_ save2_, real_8 _p_ a_, real_8 _p_ d_, int_4 _p_ jstate_);
   243  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 
   244  _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 
   245  _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 
   246  _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_, 
   247  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_
   248  , int_4 _p_ iswflg_, int_4 _p_ jstate_);
   250  real_8 _p_ yh_, real_8 _p_ y_);
   253  (*_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 
   254  _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 
   255  _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_, 
   256  int_4 _p_ ipvt_, real_8 _p_ save1_, int_4 _p_ iswflg_, real_8 _p_ bnd_, int_4 _p_ jstate_);
   258  mstate_, real_8 _p_ eps_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
   260  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 
   261  _p_ iwork_, int_4 _p_ leniw_, real_8 (*_g)(), int_4 _p_ ierflg_);
   263  , 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_ 
   264  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_ 
   265  lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, 
   266  real_8 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
   268  real_8 _p_ rh_, real_8 _p_ yh_);
   270  ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, 
   271  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_, 
   272  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_ 
   273  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_, 
   274  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_ 
   275  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_, 
   276  real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_);
   278  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 
   279  _p_ fc_, real_8 _p_ y_);
   281  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_
   282  , 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_ 
   283  evalfa_, real_4 _p_ save1_, real_4 _p_ save2_, real_4 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_);
   286  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 
   287  _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 
   288  _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 
   289  _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_, 
   290  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_
   291  , int_4 _p_ iswflg_, int_4 _p_ jstate_);
   293  real_4 _p_ yh_, real_4 _p_ y_);
   296  (*_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 
   297  _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 
   298  _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_, 
   299  int_4 _p_ ipvt_, real_4 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_);
   301  mstate_, real_4 _p_ eps_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_);
   303  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 
   304  _p_ iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_);
   306  , 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_ 
   307  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_ 
   308  lenw_, int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, 
   309  real_4 (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_);
   311  real_4 _p_ rh_, real_4 _p_ yh_);
   313  ierror_, int_4 (*_jacobn)(), int_4 _p_ matdim_, int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ miter_, int_4 _p_ ml_, 
   314  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_, 
   315  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_ 
   316  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_, 
   317  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_ 
   318  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_, 
   319  real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_);
   321  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 
   322  _p_ fc_, real_4 _p_ y_);
   323  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_, 
   324  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 
   325  _p_ nq_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, logical_4 _p_ 
   326  evalfa_, complex_8 _p_ save1_, complex_8 _p_ save2_, complex_8 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_)
   327  { // ** body not listed **
   532  }
   533  
   534  int_4 _cdcst (int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ iswflg_, real_4 _p_ el_, real_4 _p_ tq_)
   535  { // ** body not listed **
   621  }
   622  
   623  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_, 
   624  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_
   625  , 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 
   626  _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_, 
   627  complex_8 _p_ yh_, complex_8 _p_ a_, logical_4 _p_ convrg_, real_4 _p_ el_, complex_8 _p_ fac_, logical_4 _p_ ier_, 
   628  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_ 
   629  tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ jstate_)
   630  { // ** body not listed **
   778  }
   779  
   780  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_ 
   781  yh_, complex_8 _p_ y_)
   782  { // ** body not listed **
   823  }
   824  
   825  int_4 _cdpsc (int_4 _p_ ksgn_, int_4 _p_ n_, int_4 _p_ nq_, complex_8 _p_ yh_)
   826  { // ** body not listed **
   850  }
   851  
   852  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 
   853  _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_ 
   854  save2_, real_4 _p_ t_, int_4 (*_users)(), complex_8 _p_ y_, complex_8 _p_ yh_, complex_8 _p_ ywt_, real_4 _p_ uround_, 
   855  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 
   856  _p_ ipvt_, complex_8 _p_ save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_)
   857  { // ** body not listed **
  1160  }
  1161  
  1162  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_, 
  1163  real_4 _p_ eps_, complex_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
  1164  { // ** body not listed **
  1358  }
  1359  
  1360  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_, 
  1361  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_ 
  1362  iwork_, int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_)
  1363  { // ** body not listed **
  1506  }
  1507  
  1508  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 
  1509  _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_, 
  1510  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_
  1511  , int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_4 
  1512  (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
  1513  { // ** body not listed **
  3343  }
  3344  
  3345  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 
  3346  _p_ rh_, complex_8 _p_ yh_)
  3347  { // ** body not listed **
  3368  }
  3369  
  3370  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_, 
  3371  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_ 
  3372  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 
  3373  _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_, 
  3374  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_ 
  3375  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_, 
  3376  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_ 
  3377  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_
  3378  )
  3379  { // ** body not listed **
  3786  }
  3787  
  3788  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_ 
  3789  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_
  3790  , complex_8 _p_ y_)
  3791  { // ** body not listed **
  3863  }
  3864  
  3865  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_, 
  3866  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 
  3867  _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_, 
  3868  real_8 _p_ save1_, real_8 _p_ save2_, real_8 _p_ a_, real_8 _p_ d_, int_4 _p_ jstate_)
  3869  { // ** body not listed **
  4074  }
  4075  
  4076  int_4 _ddcst (int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ iswflg_, real_8 _p_ el_, real_8 _p_ tq_)
  4077  { // ** body not listed **
  4163  }
  4164  
  4165  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_, 
  4166  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_
  4167  , 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_, 
  4168  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_, 
  4169  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_ 
  4170  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 
  4171  _p_ iswflg_, int_4 _p_ jstate_)
  4172  { // ** body not listed **
  4320  }
  4321  
  4322  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_
  4323  , real_8 _p_ y_)
  4324  { // ** body not listed **
  4365  }
  4366  
  4367  int_4 _ddpsc (int_4 _p_ ksgn_, int_4 _p_ n_, int_4 _p_ nq_, real_8 _p_ yh_)
  4368  { // ** body not listed **
  4392  }
  4393  
  4394  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 
  4395  _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_ 
  4396  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_ 
  4397  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_ 
  4398  save1_, int_4 _p_ iswflg_, real_8 _p_ bnd_, int_4 _p_ jstate_)
  4399  { // ** body not listed **
  4707  }
  4708  
  4709  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 
  4710  _p_ eps_, real_8 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
  4711  { // ** body not listed **
  4905  }
  4906  
  4907  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 
  4908  _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_, 
  4909  int_4 _p_ leniw_, real_8 (*_g)(), int_4 _p_ ierflg_)
  4910  { // ** body not listed **
  5053  }
  5054  
  5055  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 
  5056  _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_, 
  5057  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_, 
  5058  int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_8 
  5059  (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
  5060  { // ** body not listed **
  6829  }
  6830  
  6831  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 
  6832  _p_ rh_, real_8 _p_ yh_)
  6833  { // ** body not listed **
  6854  }
  6855  
  6856  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_, 
  6857  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_ 
  6858  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 
  6859  _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_, 
  6860  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_, 
  6861  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_ 
  6862  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_, 
  6863  real_8 _p_ save2_, real_8 _p_ tq_, real_8 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_)
  6864  { // ** body not listed **
  7271  }
  7272  
  7273  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_ 
  7274  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_, 
  7275  real_8 _p_ y_)
  7276  { // ** body not listed **
  7348  }
  7349  
  7350  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_, 
  7351  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 
  7352  _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_, 
  7353  real_4 _p_ save1_, real_4 _p_ save2_, real_4 _p_ a_, real_4 _p_ d_, int_4 _p_ jstate_)
  7354  { // ** body not listed **
  7559  }
  7560  
  7561  int_4 _sdcst (int_4 _p_ maxord_, int_4 _p_ mint_, int_4 _p_ iswflg_, real_4 _p_ el_, real_4 _p_ tq_)
  7562  { // ** body not listed **
  7648  }
  7649  
  7650  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_, 
  7651  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_
  7652  , 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_, 
  7653  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_, 
  7654  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_ 
  7655  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 
  7656  _p_ iswflg_, int_4 _p_ jstate_)
  7657  { // ** body not listed **
  7805  }
  7806  
  7807  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_
  7808  , real_4 _p_ y_)
  7809  { // ** body not listed **
  7850  }
  7851  
  7852  int_4 _sdpsc (int_4 _p_ ksgn_, int_4 _p_ n_, int_4 _p_ nq_, real_4 _p_ yh_)
  7853  { // ** body not listed **
  7877  }
  7878  
  7879  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 
  7880  _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_ 
  7881  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_ 
  7882  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_ 
  7883  save1_, int_4 _p_ iswflg_, real_4 _p_ bnd_, int_4 _p_ jstate_)
  7884  { // ** body not listed **
  8192  }
  8193  
  8194  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 
  8195  _p_ eps_, real_4 _p_ work_, int_4 _p_ lenw_, int_4 _p_ ierflg_)
  8196  { // ** body not listed **
  8390  }
  8391  
  8392  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 
  8393  _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_, 
  8394  int_4 _p_ leniw_, real_4 (*_g)(), int_4 _p_ ierflg_)
  8395  { // ** body not listed **
  8538  }
  8539  
  8540  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 
  8541  _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_, 
  8542  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_, 
  8543  int_4 _p_ iwork_, int_4 _p_ leniw_, int_4 (*_jacobn)(), int_4 (*_fa)(), int_4 _p_ nde_, int_4 _p_ mxstep_, real_4 
  8544  (*_g)(), int_4 (*_users)(), int_4 _p_ ierflg_)
  8545  { // ** body not listed **
 10314  }
 10315  
 10316  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 
 10317  _p_ rh_, real_4 _p_ yh_)
 10318  { // ** body not listed **
 10339  }
 10340  
 10341  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_, 
 10342  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_ 
 10343  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 
 10344  _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_, 
 10345  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_, 
 10346  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_ 
 10347  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_, 
 10348  real_4 _p_ save2_, real_4 _p_ tq_, real_4 _p_ trend_, int_4 _p_ iswflg_, int_4 _p_ mtrsv_, int_4 _p_ mxrdsv_)
 10349  { // ** body not listed **
 10756  }
 10757  
 10758  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_ 
 10759  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_, 
 10760  real_4 _p_ y_)
 10761  { // ** body not listed **
 10833  }
 10834  


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