mathlib-odepack.c

     1  //! @file mathlib-odepack.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  //! 84 subprograms from ODEPACK
    24  
    25  // Compiled from Fortran source code by VIF.
    26  // Selected subprograms are:
    27  //
    28  // ADJLR  DCFODE DIPREP DLSODK DPREP  DROOTS DSPIOM DSTODE IUMACH NNFC
    29  // CDRV   DDECBT DIPREP DLSODP DPREPI DSETPK DSRCAR DSTODI IXSAV  NNSC
    30  // CNTNZU DEWSET DLHIN  DLSOIB DPREPJ DSLSBT DSRCKR DSTODP JGROUP NNTC
    31  // DAIGBT DFNORM DLSODA DMNORM DPREPJ DSOLBT DSRCMA DSTOKA MD     NROC
    32  // DAINVG DHEFA  DLSODA DORTHO DPRJA  DSOLPK DSRCMS DUMACH MDI    NSFC
    33  // DAINVG DHELS  DLSODE DPCG   DPRJIS DSOLSS DSRCOM DUMSUM MDM    ODRV
    34  // DATP   DHEQR  DLSODE DPCGS  DPRJS  DSOLSY DSRCPK DUSOL  MDP    SRO
    35  // DATV   DHESL  DLSODI DPJIBT DRCHEK DSPIGM DSTODA DVNORM MDU    XERRWD
    36  // DBNORM DINTDY DLSODI DPKSET
    37  
    38  // ODEPACK - A Systematized Collection of ODE Solvers
    39  // 
    40  // Source: netlib.org/odepack
    41  // 
    42  // Original author is Alan C. Hindmarsh.
    43  // 
    44  // Work performed under the auspices of the U.S. Department of Energy
    45  // by the Lawrence Livermore National Laboratory under contract
    46  // No. W-7405-Eng-48, and supported (formerly) by the DOE Office of 
    47  // Energy Research, Applied Mathematical Sciences Research Program.
    48  // 
    49  
    50  // The license for ODEPACK Fortran source code is:
    51  //
    52  // NOTICE
    53  // 
    54  // This work was produced at the University of California, Lawrence Livermore
    55  // National Laboratory (UC LLNL) under contract no. W-7405-ENG-48 (Contract
    56  // 48) between the U.S. Department of Energy (DOE) and The Regents of the
    57  // University of California (University) for the operation of UC LLNL. The
    58  // rights of the Federal Government are reserved under Contract 48 subject
    59  // to the restrictions agreed upon by the DOE and University as allowed
    60  // under DOE Acquisition Letter 97-1.
    61  // 
    62  // The ODEPACK package has been declared to be in the Public Domain.
    63  // 
    64  // DISCLAIMER
    65  // 
    66  // This work was prepared as an account of work sponsored by an agency of
    67  // the United States Government. Neither the United States Government nor the
    68  // University of California nor any of their employees, makes any warranty,
    69  // express or implied, or assumes any liability or responsibility for the
    70  // accuracy, completeness, or usefulness of any information, apparatus,
    71  // product, or process disclosed, or represents that its use would not
    72  // infringe privately-owned rights. Reference herein to any specific
    73  // commercial products, process, or service by trade name, trademark,
    74  // manufacturer or otherwise does not necessarily constitute or imply
    75  // its endorsement, recommendation, or favoring by the United States
    76  // Government or the University of California. The views and opinions of
    77  // authors expressed herein do not necessarily state or reflect those of
    78  // the United States Government or the University of California, and shall
    79  // not be used for advertising or product endorsement purposes.
    80  
    81  /*
    82  Generated by VIF - experimental VIntage Fortran compiler.
    83  VIF release 1.2.17
    84  */
    85  #if defined (__GNUC__)
    86  #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
    87  #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
    88  #if (__GNUC__ >= 14)
    89  #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
    90  #pragma GCC diagnostic ignored "-Wimplicit-int"
    91  #pragma GCC diagnostic ignored "-Wint-conversion"
    92  #pragma GCC diagnostic ignored "-Wreturn-mismatch"
    93  #endif
    94  #else
    95  #error VIF requires GCC
    96  #endif
   518  static CALLS __calls[__ncalls] = {
   519    {"adjlr", 0}, // subroutine
   520    {"cdrv", 0}, // subroutine
   521    {"nroc", 0}, // subroutine
   522    {"nsfc", 0}, // subroutine
   523    {"nnfc", 0}, // subroutine
   524    {"nnsc", 0}, // subroutine
   525    {"nntc", 0}, // subroutine
   526    {"cntnzu", 0}, // subroutine
   527    {"daigbt", 0}, // subroutine
   528    {"dainvg", 0}, // subroutine
   529    {"dainvgs", 0}, // subroutine
   530    {"datp", 0}, // subroutine
   531    {"datv", 0}, // subroutine
   532    {"dbnorm", 0}, // real*8 function
   533    {"dcfode", 0}, // subroutine
   534    {"ddecbt", 0}, // subroutine
   535    {"dewset", 0}, // subroutine
   536    {"dfnorm", 0}, // real*8 function
   537    {"dhefa", 0}, // subroutine
   538    {"dhels", 0}, // subroutine
   539    {"dheqr", 0}, // subroutine
   540    {"dhesl", 0}, // subroutine
   541    {"dintdy", 0}, // subroutine
   542    {"diprep", 0}, // subroutine
   543    {"diprepi", 0}, // subroutine
   544    {"dlhin", 0}, // subroutine
   545    {"dlsoda", 0}, // subroutine
   546    {"dlsodar", 0}, // subroutine
   547    {"dlsode", 0}, // subroutine
   548    {"dlsodes", 0}, // subroutine
   549    {"dlsoibt", 0}, // subroutine
   550    {"dlsodi", 0}, // subroutine
   551    {"dlsodis", 0}, // subroutine
   552    {"dlsodkr", 0}, // subroutine
   553    {"dlsodpk", 0}, // subroutine
   554    {"dmnorm", 0}, // real*8 function
   555    {"dorthog", 0}, // subroutine
   556    {"dpcg", 0}, // subroutine
   557    {"dpcgs", 0}, // subroutine
   558    {"dpjibt", 0}, // subroutine
   559    {"dpkset", 0}, // subroutine
   560    {"dprep", 0}, // subroutine
   561    {"dprepi", 0}, // subroutine
   562    {"dprepj", 0}, // subroutine
   563    {"dprepji", 0}, // subroutine
   564    {"dprja", 0}, // subroutine
   565    {"dprjis", 0}, // subroutine
   566    {"dprjs", 0}, // subroutine
   567    {"drchek", 0}, // subroutine
   568    {"droots", 0}, // subroutine
   569    {"dsetpk", 0}, // subroutine
   570    {"dslsbt", 0}, // subroutine
   571    {"dsolbt", 0}, // subroutine
   572    {"dsolpk", 0}, // subroutine
   573    {"dsolss", 0}, // subroutine
   574    {"dsolsy", 0}, // subroutine
   575    {"dspigmr", 0}, // subroutine
   576    {"dspiom", 0}, // subroutine
   577    {"dsrcar", 0}, // subroutine
   578    {"dsrckr", 0}, // subroutine
   579    {"dsrcma", 0}, // subroutine
   580    {"dsrcms", 0}, // subroutine
   581    {"dsrcom", 0}, // subroutine
   582    {"dsrcpk", 0}, // subroutine
   583    {"dstoda", 0}, // subroutine
   584    {"dstode", 0}, // subroutine
   585    {"dstodi", 0}, // subroutine
   586    {"dstodpk", 0}, // subroutine
   587    {"dstoka", 0}, // subroutine
   588    {"dumach", 0}, // real*8 function
   589    {"dumsum", 0}, // subroutine
   590    {"dusol", 0}, // subroutine
   591    {"dvnorm", 0}, // real*8 function
   592    {"iumach", 0}, // integer*4 function
   593    {"ixsav", 0}, // integer*4 function
   594    {"jgroup", 0}, // subroutine
   595    {"odrv", 0}, // subroutine
   596    {"md", 0}, // subroutine
   597    {"mdi", 0}, // subroutine
   598    {"mdm", 0}, // subroutine
   599    {"mdp", 0}, // subroutine
   600    {"mdu", 0}, // subroutine
   601    {"sro", 0}, // subroutine
   602    {"xerrwd", 0}, // subroutine
   603    {NULL, 0}
   604  };
   605  
   606  static struct {
   607    int_4 iplost_, iesp_, istatc_, iys_, iba_, ibian_, ibjan_, ibjgp_, ipian_, ipjan_, ipjgp_, ipigp_, ipr_, ipc_, ipic_, 
   608    ipisp_, iprsp_, ipa_, lenyh_, lenyhm_, lenwk_, lreq_, lrat_, lrest_, lwmin_, moss_, msbj_, nslj_, ngp_, nlu_, nnz_, 
   609    nsp_, nzl_, nzu_, ilss_[34];
   610    real_8 rlss_[6], con0_, conmin_, ccmxj_, psmall_, rbig_, seth_;
   611  } dlss01_;
   612  static struct {
   613    int_4 iownd_[6], iowns_[6], icf_, ierpj_, iersl_, jcur_, jstart_, kflag_, l_, lyh_, lewt_, lacor_, lsavf_, lwm_, liwm_, 
   614    meth_, miter_, maxord_, maxcor_, msbp_, mxncf_, n_, nq_, nst_, nfe_, nje_, nqu_, init_, mxstep_, mxhnil_, nhnil_, 
   615    nslast_, nyh_, ils_[37], ialth_, ipup_, lmax_, meo_, nqnyh_, nslp_;
   616    real_8 rowns_[209], ccmax_, el0_, h_, hmin_, hmxi_, hu_, rc_, tn_, uround_, rls_[218], conit_, crate_, el_[13], 
   617    elco_[156], hold_, rmax_, tesco_[36];
   618  } dls001_;
   619  static struct {
   620    int_4 insufr_, insufi_, ixpr_, iowns2_[2], jtyp_, mused_, mxordn_, mxords_, iownd2_[3], ilsa_[9], icount_, irflag_;
   621    real_8 tsw_, rowns2_[20], pdnorm_, rownd2_, rlsa_[22], cm1_[12], cm2_[5], pdest_, pdlast_, ratio_;
   622  } dlsa01_;
   623  static struct {
   624    int_4 lg0_, lg1_, lgx_, iownr3_[2], irfnd_, itaskc_, ngc_, nge_, iownd3_[3], imax_, last_, idum3_[4], ilsr_[9];
   625    real_8 rownr3_[2], t0_, tlast_, toutc_, alpha_, x2_, rdum3_[3], rlsr_[5];
   626  } dlsr01_;
   627  static struct {
   628    int_4 newt_, nsfi_, nslj_, njev_, ils2_[4];
   629    real_8 stifr_, rls2_;
   630  } dls002_;
   631  static struct {
   632    int_4 jpre_, jacflg_, locwp_, lociwp_, lsavx_, kmp_, maxl_, mnewt_, nni_, nli_, nps_, ncfn_, ncfl_, ilsp_[13];
   633    real_8 delt_, epcon_, sqrtn_, rsqrtn_, rlsp_[4];
   634  } dlpk01_;
   637  _p_ a_, real_8 _p_ b_, real_8 _p_ z_, int_4 _p_ nsp_, int_4 _p_ isp_, real_8 _p_ rsp_, int_4 _p_ esp_, int_4 _p_ path_, 
   638  int_4 _p_ flag_);
   640  _p_ ar_, int_4 _p_ p_, int_4 _p_ flag_);
   642  _p_ il_, int_4 _p_ jl_, int_4 _p_ ijl_, int_4 _p_ jumax_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, int_4 _p_ q_, 
   643  int_4 _p_ ira_, int_4 _p_ jra_, int_4 _p_ irac_, int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ iru_, int_4 _p_ jru_, int_4 
   644  _p_ flag_);
   646  _p_ a_, real_8 _p_ z_, real_8 _p_ b_, int_4 _p_ lmax_, int_4 _p_ il_, int_4 _p_ jl_, int_4 _p_ ijl_, real_8 _p_ l_, 
   647  real_8 _p_ d_, int_4 _p_ umax_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, real_8 _p_ u_, real_8 _p_ row_, real_8 
   648  _p_ tmp_, int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ flag_);
   650  _p_ l_, real_8 _p_ d_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, real_8 _p_ u_, real_8 _p_ z_, real_8 _p_ b_, 
   651  real_8 _p_ tmp_);
   653  _p_ l_, real_8 _p_ d_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, real_8 _p_ u_, real_8 _p_ z_, real_8 _p_ b_, 
   654  real_8 _p_ tmp_);
   657  ydot_, int_4 _p_ mb_, int_4 _p_ nb_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_);
   659  ydot_, int_4 _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_);
   661  , real_8 _p_ ydot_, int_4 _p_ ier_, int_4 (*_res)(), int_4 (*_adda)());
   663  hl0_, real_8 _p_ wk_, int_4 (*_f)(), real_8 _p_ w_);
   665  ftem_, int_4 (*_f)(), int_4 (*_psol)(), real_8 _p_ z_, real_8 _p_ vtem_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ 
   666  hl0_, int_4 _p_ jpre_, int_4 _p_ ier_, int_4 _p_ npsl_);
   670  _p_ ier_);
   672  _p_ ewt_);
   679  );
   681  ipflag_, real_4 (*_f)(), int_4 (*_jac)());
   683  , int_4 _p_ ic_, int_4 _p_ jc_, int_4 _p_ ipflag_, real_4 (*_res)(), int_4 (*_jac)(), real_4 (*_adda)());
   685  real_8 _p_ tout_, real_8 _p_ uround_, real_8 _p_ ewt_, int_4 _p_ itol_, real_8 _p_ atol_, real_8 _p_ y_, real_8 _p_ 
   686  temp_, real_8 _p_ h0_, int_4 _p_ niter_, int_4 _p_ ier_);
   688  , real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 
   689  _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_);
   691  itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, 
   692  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_, real_4 (*_g)(), int_4 _p_ ng_, int_4 
   693  _p_ jroot_);
   695  , real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 
   696  _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_);
   698  itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, 
   699  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_);
   701  _p_ ydoti_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, 
   702  int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 _p_ mf_)
   703  ;
   705  ydoti_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 
   706  _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 _p_ mf_);
   708  _p_ ydoti_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, 
   709  int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 _p_ mf_)
   710  ;
   712  itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, 
   713  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_, real_4 (*_g)(), 
   714  int_4 _p_ ng_, int_4 _p_ jroot_);
   716  itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, 
   717  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_);
   720  ldhes_, int_4 _p_ kmp_, real_8 _p_ snormw_);
   722  , int_4 _p_ n_, int_4 _p_ maxl_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 _p_ mnewt_, real_4 (*_f)(), 
   723  int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ p_, real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ lpcg_, real_8 
   724  _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   726  wght_, int_4 _p_ n_, int_4 _p_ maxl_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 _p_ mnewt_, real_4 
   727  (*_f)(), int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ p_, real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ lpcg_
   728  , real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   730  rtem_, real_8 _p_ savr_, real_8 _p_ s_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4 (*_jac)(), int_4 
   731  (*_adda)());
   733  savf_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)());
   735  ftem_, int_4 _p_ ia_, int_4 _p_ ja_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 _p_ ipper_, int_4 (*_f)(), int_4 (*_jac)());
   737  ewt_, real_8 _p_ rtem_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ ic_, int_4 _p_ jc_, real_8 _p_ wk_, int_4 _p_ iwk_, 
   738  int_4 _p_ ipper_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)());
   740  ftem_, real_8 _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)());
   742  rtem_, real_8 _p_ savr_, real_8 _p_ s_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4 (*_jac)(), int_4 
   743  (*_adda)());
   745  ftem_, real_8 _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)());
   747  rtem_, real_8 _p_ savr_, real_8 _p_ s_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_res)(), int_4 (*_jac)(), int_4 
   748  (*_adda)());
   750  ftem_, real_8 _p_ savf_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_f)(), int_4 (*_jac)());
   752  real_8 _p_ g0_, real_8 _p_ g1_, real_8 _p_ gx_, int_4 _p_ jroot_, int_4 _p_ irt_);
   754  g0_, real_8 _p_ g1_, real_8 _p_ gx_, real_8 _p_ x_, int_4 _p_ jroot_);
   756  savf_, int_4 _p_ jok_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)());
   759  _p_ ip_);
   761  wm_, int_4 _p_ iwm_, real_4 (*_f)(), real_4 (*_psol)());
   765  wght_, int_4 _p_ n_, int_4 _p_ maxl_, int_4 _p_ maxlp1_, int_4 _p_ kmp_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ 
   766  jpre_, int_4 _p_ mnewt_, real_4 (*_f)(), int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ v_, real_8 _p_ 
   767  hes_, real_8 _p_ q_, int_4 _p_ lgmr_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, real_8 _p_ dl_, int_4 _p_ iflag_)
   768  ;
   770  wght_, int_4 _p_ n_, int_4 _p_ maxl_, int_4 _p_ kmp_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 _p_ 
   771  mnewt_, real_4 (*_f)(), int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ v_, real_8 _p_ hes_, int_4 _p_ 
   772  ipvt_, int_4 _p_ liom_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   780  ewt_, real_8 _p_ savf_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)(), int_4 
   781  (*_pjac)(), int_4 (*_slvs)());
   783  ewt_, real_8 _p_ savf_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)(), int_4 
   784  (*_pjac)(), int_4 (*_slvs)());
   786  ewt_, real_8 _p_ savf_, real_8 _p_ savr_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), real_4 
   787  (*_adda)(), int_4 (*_jac)(), int_4 (*_pjac)(), int_4 (*_slvs)());
   789  ewt_, real_8 _p_ savf_, real_8 _p_ savx_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 
   790  (*_jac)(), real_4 (*_psol)());
   792  ewt_, real_8 _p_ savf_, real_8 _p_ savx_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 
   793  (*_jac)(), real_4 (*_psol)());
   797  wght_, int_4 _p_ n_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ mnewt_, int_4 (*_psol)(), int_4 _p_ npsl_, real_8 
   798  _p_ x_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   803  int_4 _p_ jgp_, int_4 _p_ incl_, int_4 _p_ jdone_, int_4 _p_ ier_);
   805  _p_ nsp_, int_4 _p_ isp_, int_4 _p_ path_, int_4 _p_ flag_);
   807  head_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ flag_);
   809  head_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ tag_, int_4 _p_ flag_);
   811  int_4 _p_ mark_);
   813  _p_ last_, int_4 _p_ next_, int_4 _p_ mark_);
   815  int_4 _p_ next_, int_4 _p_ mark_);
   817  r_, logical_4 _p_ dflag_);
   819  i1_, int_4 _p_ i2_, int_4 _p_ nr_, real_8 _p_ r1_, real_8 _p_ r2_);
   820  int_4 _adjlr (int_4 _p_ n_, int_4 _p_ isp_, int_4 _p_ ldif_)
   821  { // ** body not listed **
   832  }
   833  
   834  int_4 _cdrv (int_4 _p_ n_, int_4 _p_ r_, int_4 _p_ c_, int_4 _p_ ic_, int_4 _p_ ia_, int_4 _p_ ja_, real_8 _p_ a_, 
   835  real_8 _p_ b_, real_8 _p_ z_, int_4 _p_ nsp_, int_4 _p_ isp_, real_8 _p_ rsp_, int_4 _p_ esp_, int_4 _p_ path_, int_4 
   836  _p_ flag_)
   837  { // ** body not listed **
   955  }
   956  
   957  int_4 _nroc (int_4 _p_ n_, int_4 _p_ ic_, int_4 _p_ ia_, int_4 _p_ ja_, real_8 _p_ a_, int_4 _p_ jar_, real_8 _p_ ar_, 
   958  int_4 _p_ p_, int_4 _p_ flag_)
   959  { // ** body not listed **
  1001  }
  1002  
  1003  int_4 _nsfc (int_4 _p_ n_, int_4 _p_ r_, int_4 _p_ ic_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ jlmax_, int_4 _p_ il_, 
  1004  int_4 _p_ jl_, int_4 _p_ ijl_, int_4 _p_ jumax_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, int_4 _p_ q_, int_4 _p_ 
  1005  ira_, int_4 _p_ jra_, int_4 _p_ irac_, int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ iru_, int_4 _p_ jru_, int_4 _p_ flag_)
  1006  { // ** body not listed **
  1395  }
  1396  
  1397  int_4 _nnfc (int_4 _p_ n_, int_4 _p_ r_, int_4 _p_ c_, int_4 _p_ ic_, int_4 _p_ ia_, int_4 _p_ ja_, real_8 _p_ a_, 
  1398  real_8 _p_ z_, real_8 _p_ b_, int_4 _p_ lmax_, int_4 _p_ il_, int_4 _p_ jl_, int_4 _p_ ijl_, real_8 _p_ l_, real_8 _p_ 
  1399  d_, int_4 _p_ umax_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, real_8 _p_ u_, real_8 _p_ row_, real_8 _p_ tmp_, 
  1400  int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ flag_)
  1401  { // ** body not listed **
  1554  }
  1555  
  1556  int_4 _nnsc (int_4 _p_ n_, int_4 _p_ r_, int_4 _p_ c_, int_4 _p_ il_, int_4 _p_ jl_, int_4 _p_ ijl_, real_8 _p_ l_, 
  1557  real_8 _p_ d_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, real_8 _p_ u_, real_8 _p_ z_, real_8 _p_ b_, real_8 _p_ 
  1558  tmp_)
  1559  { // ** body not listed **
  1598  }
  1599  
  1600  int_4 _nntc (int_4 _p_ n_, int_4 _p_ r_, int_4 _p_ c_, int_4 _p_ il_, int_4 _p_ jl_, int_4 _p_ ijl_, real_8 _p_ l_, 
  1601  real_8 _p_ d_, int_4 _p_ iu_, int_4 _p_ ju_, int_4 _p_ iju_, real_8 _p_ u_, real_8 _p_ z_, real_8 _p_ b_, real_8 _p_ 
  1602  tmp_)
  1603  { // ** body not listed **
  1641  }
  1642  
  1643  int_4 _cntnzu (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ nzsut_)
  1644  { // ** body not listed **
  1684  }
  1685  
  1686  int_4 _daigbt (int_4 (*_res)(), int_4 (*_adda)(), int_4 _p_ neq_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ ydot_, int_4 
  1687  _p_ mb_, int_4 _p_ nb_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_)
  1688  { // ** body not listed **
  1714  }
  1715  
  1716  int_4 _dainvg (int_4 (*_res)(), int_4 (*_adda)(), int_4 _p_ neq_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ ydot_, int_4 
  1717  _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_)
  1718  { // ** body not listed **
  1766  }
  1767  
  1768  int_4 _dainvgs (int_4 _p_ neq_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ wk_, int_4 _p_ iwk_, real_8 _p_ tem_, real_8 
  1769  _p_ ydot_, int_4 _p_ ier_, int_4 (*_res)(), int_4 (*_adda)())
  1770  { // ** body not listed **
  1825  }
  1826  
  1827  int_4 _datp (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ p_, real_8 _p_ wght_, real_8 _p_ hl0_, real_8 
  1828  _p_ wk_, int_4 (*_f)(), real_8 _p_ w_)
  1829  { // ** body not listed **
  1847  }
  1848  
  1849  int_4 _datv (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ v_, real_8 _p_ wght_, real_8 _p_ ftem_, int_4 
  1850  (*_f)(), int_4 (*_psol)(), real_8 _p_ z_, real_8 _p_ vtem_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ hl0_, int_4 _p_ 
  1851  jpre_, int_4 _p_ ier_, int_4 _p_ npsl_)
  1852  { // ** body not listed **
  1911  }
  1912  
  1913  real_8 _dbnorm (int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ nra_, int_4 _p_ ml_, int_4 _p_ mu_, real_8 _p_ w_)
  1914  { // ** body not listed **
  1932  }
  1933  
  1934  int_4 _dcfode (int_4 _p_ meth_, real_8 _p_ elco_, real_8 _p_ tesco_)
  1935  { // ** body not listed **
  2010  }
  2011  
  2012  int_4 _ddecbt (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, int_4 _p_ ip_, int_4 _p_ ier_)
  2013  { // ** body not listed **
  2087  }
  2088  
  2089  int_4 _dewset (int_4 _p_ n_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, real_8 _p_ ycur_, real_8 _p_ ewt_)
  2090  { // ** body not listed **
  2121  }
  2122  
  2123  real_8 _dfnorm (int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_)
  2124  { // ** body not listed **
  2139  }
  2140  
  2141  int_4 _dhefa (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, int_4 _p_ ipvt_, int_4 _p_ info_, int_4 _p_ job_)
  2142  { // ** body not listed **
  2247  }
  2248  
  2249  int_4 _dhels (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, real_8 _p_ q_, real_8 _p_ b_)
  2250  { // ** body not listed **
  2272  }
  2273  
  2274  int_4 _dheqr (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, real_8 _p_ q_, int_4 _p_ info_, int_4 _p_ ijob_)
  2275  { // ** body not listed **
  2371  }
  2372  
  2373  int_4 _dhesl (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, int_4 _p_ ipvt_, real_8 _p_ b_)
  2374  { // ** body not listed **
  2402  }
  2403  
  2404  int_4 _dintdy (real_8 _p_ t_, int_4 _p_ k_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ dky_, int_4 _p_ iflag_)
  2405  { // ** body not listed **
  2480  }
  2481  
  2482  int_4 _diprep (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ rwork_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ ipflag_, 
  2483  real_4 (*_f)(), int_4 (*_jac)())
  2484  { // ** body not listed **
  2524  }
  2525  
  2526  int_4 _diprepi (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ s_, real_8 _p_ rwork_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 
  2527  _p_ ic_, int_4 _p_ jc_, int_4 _p_ ipflag_, real_4 (*_res)(), int_4 (*_jac)(), real_4 (*_adda)())
  2528  { // ** body not listed **
  2569  }
  2570  
  2571  int_4 _dlhin (int_4 _p_ neq_, int_4 _p_ n_, real_8 _p_ t0_, real_8 _p_ y0_, real_8 _p_ ydot_, int_4 (*_f)(), real_8 _p_ 
  2572  tout_, real_8 _p_ uround_, real_8 _p_ ewt_, int_4 _p_ itol_, real_8 _p_ atol_, real_8 _p_ y_, real_8 _p_ temp_, real_8 
  2573  _p_ h0_, int_4 _p_ niter_, int_4 _p_ ier_)
  2574  { // ** body not listed **
  2691  }
  2692  
  2693  int_4 _dlsoda (int_4 (*_f)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 
  2694  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, 
  2695  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_)
  2696  { // ** body not listed **
  3578  }
  3579  
  3580  int_4 _dlsodar (int_4 (*_f)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 
  3581  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, 
  3582  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_, real_4 (*_g)(), int_4 _p_ ng_, int_4 _p_ jroot_)
  3583  { // ** body not listed **
  4572  }
  4573  
  4574  int_4 _dlsode (int_4 (*_f)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 
  4575  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, 
  4576  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_)
  4577  { // ** body not listed **
  5326  }
  5327  
  5328  int_4 _dlsodes (int_4 (*_f)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 
  5329  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, 
  5330  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_)
  5331  { // ** body not listed **
  6309  }
  6310  
  6311  int_4 _dlsoibt (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  6312  real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ 
  6313  istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 _p_ mf_)
  6314  { // ** body not listed **
  7122  }
  7123  
  7124  int_4 _dlsodi (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  7125  real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ 
  7126  istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 _p_ mf_)
  7127  { // ** body not listed **
  7946  }
  7947  
  7948  int_4 _dlsodis (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  7949  real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ 
  7950  istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 _p_ mf_)
  7951  { // ** body not listed **
  9029  }
  9030  
  9031  int_4 _dlsodkr (int_4 (*_f)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 
  9032  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, 
  9033  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_, real_4 (*_g)(), int_4 _p_ ng_, 
  9034  int_4 _p_ jroot_)
  9035  { // ** body not listed **
  9980  }
  9981  
  9982  int_4 _dlsodpk (int_4 (*_f)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ t_, real_8 _p_ tout_, int_4 _p_ itol_, real_8 
  9983  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ itask_, int_4 _p_ istate_, int_4 _p_ iopt_, real_8 _p_ rwork_, int_4 _p_ lrw_, 
  9984  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_)
  9985  { // ** body not listed **
 10842  }
 10843  
 10844  real_8 _dmnorm (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ w_)
 10845  { // ** body not listed **
 10856  }
 10857  
 10858  int_4 _dorthog (real_8 _p_ vnew_, real_8 _p_ v_, real_8 _p_ hes_, int_4 _p_ n_, int_4 _p_ ll_, int_4 _p_ ldhes_, int_4 
 10859  _p_ kmp_, real_8 _p_ snormw_)
 10860  { // ** body not listed **
 10895  }
 10896  
 10897  int_4 _dpcg (int_4 _p_ neq_, real_8 _p_ tn_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ r_, real_8 _p_ wght_, int_4 
 10898  _p_ n_, int_4 _p_ maxl_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 _p_ mnewt_, real_4 (*_f)(), int_4 
 10899  (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ p_, real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ lpcg_, real_8 _p_ 
 10900  wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 10901  { // ** body not listed **
 10989  }
 10990  
 10991  int_4 _dpcgs (int_4 _p_ neq_, real_8 _p_ tn_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ r_, real_8 _p_ wght_, int_4 
 10992  _p_ n_, int_4 _p_ maxl_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 _p_ mnewt_, real_4 (*_f)(), int_4 
 10993  (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ p_, real_8 _p_ w_, real_8 _p_ z_, int_4 _p_ lpcg_, real_8 _p_ 
 10994  wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 10995  { // ** body not listed **
 11088  }
 11089  
 11090  int_4 _dpjibt (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ ewt_, real_8 _p_ rtem_, real_8 
 11091  _p_ savr_, real_8 _p_ s_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 11092  { // ** body not listed **
 11222  }
 11223  
 11224  int_4 _dpkset (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ysv_, real_8 _p_ ewt_, real_8 _p_ ftem_, real_8 _p_ savf_, 
 11225  real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)())
 11226  { // ** body not listed **
 11242  }
 11243  
 11244  int_4 _dprep (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, real_8 _p_ savf_, real_8 _p_ ewt_, real_8 _p_ ftem_, int_4 
 11245  _p_ ia_, int_4 _p_ ja_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 _p_ ipper_, int_4 (*_f)(), int_4 (*_jac)())
 11246  { // ** body not listed **
 11505  }
 11506  
 11507  int_4 _dprepi (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ s_, real_8 _p_ yh_, real_8 _p_ savr_, real_8 _p_ ewt_, real_8 
 11508  _p_ rtem_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ ic_, int_4 _p_ jc_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 _p_ ipper_
 11509  , int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 11510  { // ** body not listed **
 11914  }
 11915  
 11916  int_4 _dprepj (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ ewt_, real_8 _p_ ftem_, real_8 
 11917  _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)())
 11918  { // ** body not listed **
 12065  }
 12066  int_4 _dprepji (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ ewt_, real_8 _p_ rtem_, 
 12067  real_8 _p_ savr_, real_8 _p_ s_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 12068  { // ** body not listed **
 12219  }
 12220  
 12221  int_4 _dprja (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ ewt_, real_8 _p_ ftem_, real_8 
 12222  _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)())
 12223  { // ** body not listed **
 12349  }
 12350  
 12351  int_4 _dprjis (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ ewt_, real_8 _p_ rtem_, real_8 
 12352  _p_ savr_, real_8 _p_ s_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 12353  { // ** body not listed **
 12463  }
 12464  
 12465  int_4 _dprjs (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ ewt_, real_8 _p_ ftem_, real_8 
 12466  _p_ savf_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_f)(), int_4 (*_jac)())
 12467  { // ** body not listed **
 12636  }
 12637  
 12638  int_4 _drchek (int_4 _p_ job_, int_4 (*_g)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ 
 12639  g0_, real_8 _p_ g1_, real_8 _p_ gx_, int_4 _p_ jroot_, int_4 _p_ irt_)
 12640  { // ** body not listed **
 12782  }
 12783  
 12784  int_4 _droots (int_4 _p_ ng_, real_8 _p_ hmin_, int_4 _p_ jflag_, real_8 _p_ x0_, real_8 _p_ x1_, real_8 _p_ g0_, 
 12785  real_8 _p_ g1_, real_8 _p_ gx_, real_8 _p_ x_, int_4 _p_ jroot_)
 12786  { // ** body not listed **
 13004  }
 13005  
 13006  int_4 _dsetpk (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ysv_, real_8 _p_ ewt_, real_8 _p_ ftem_, real_8 _p_ savf_, 
 13007  int_4 _p_ jok_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)())
 13008  { // ** body not listed **
 13027  }
 13028  
 13029  int_4 _dslsbt (real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ x_, real_8 _p_ tem_)
 13030  { // ** body not listed **
 13040  }
 13041  
 13042  int_4 _dsolbt (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_, real_8 _p_ y_, int_4 _p_ ip_)
 13043  { // ** body not listed **
 13080  }
 13081  
 13082  int_4 _dsolpk (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ x_, real_8 _p_ ewt_, real_8 _p_ wm_, int_4 
 13083  _p_ iwm_, real_4 (*_f)(), real_4 (*_psol)())
 13084  { // ** body not listed **
 13219  }
 13220  
 13221  int_4 _dsolss (real_8 _p_ wk_, int_4 _p_ iwk_, real_8 _p_ x_, real_8 _p_ tem_)
 13222  { // ** body not listed **
 13266  }
 13267  
 13268  int_4 _dsolsy (real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ x_, real_8 _p_ tem_)
 13269  { // ** body not listed **
 13315  }
 13316  
 13317  int_4 _dspigmr (int_4 _p_ neq_, real_8 _p_ tn_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ b_, real_8 _p_ wght_, int_4 
 13318  _p_ n_, int_4 _p_ maxl_, int_4 _p_ maxlp1_, int_4 _p_ kmp_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 
 13319  _p_ mnewt_, real_4 (*_f)(), int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ v_, real_8 _p_ hes_, real_8 
 13320  _p_ q_, int_4 _p_ lgmr_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, real_8 _p_ dl_, int_4 _p_ iflag_)
 13321  { // ** body not listed **
 13465  }
 13466  
 13467  int_4 _dspiom (int_4 _p_ neq_, real_8 _p_ tn_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ b_, real_8 _p_ wght_, int_4 
 13468  _p_ n_, int_4 _p_ maxl_, int_4 _p_ kmp_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ jpre_, int_4 _p_ mnewt_, real_4 
 13469  (*_f)(), int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 _p_ v_, real_8 _p_ hes_, int_4 _p_ ipvt_, int_4 _p_ 
 13470  liom_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 13471  { // ** body not listed **
 13604  }
 13605  
 13606  int_4 _dsrcar (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13607  { // ** body not listed **
 13716  }
 13717  
 13718  int_4 _dsrckr (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13719  { // ** body not listed **
 13826  }
 13827  
 13828  int_4 _dsrcma (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13829  { // ** body not listed **
 13902  }
 13903  
 13904  int_4 _dsrcms (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13905  { // ** body not listed **
 13978  }
 13979  
 13980  int_4 _dsrcom (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13981  { // ** body not listed **
 14024  }
 14025  
 14026  int_4 _dsrcpk (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 14027  { // ** body not listed **
 14092  }
 14093  
 14094  int_4 _dstoda (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ yh1_, real_8 _p_ ewt_, real_8 
 14095  _p_ savf_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)(), int_4 (*_pjac)(), int_4 
 14096  (*_slvs)())
 14097  { // ** body not listed **
 14665  }
 14666  
 14667  int_4 _dstode (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ yh1_, real_8 _p_ ewt_, real_8 
 14668  _p_ savf_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)(), int_4 (*_pjac)(), int_4 
 14669  (*_slvs)())
 14670  { // ** body not listed **
 15079  }
 15080  
 15081  int_4 _dstodi (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ yh1_, real_8 _p_ ewt_, real_8 
 15082  _p_ savf_, real_8 _p_ savr_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), real_4 (*_adda)(), 
 15083  int_4 (*_jac)(), int_4 (*_pjac)(), int_4 (*_slvs)())
 15084  { // ** body not listed **
 15494  }
 15495  
 15496  int_4 _dstodpk (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ yh1_, real_8 _p_ ewt_, real_8 
 15497  _p_ savf_, real_8 _p_ savx_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)(), real_4 
 15498  (*_psol)())
 15499  { // ** body not listed **
 15921  }
 15922  
 15923  int_4 _dstoka (int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ yh_, int_4 _p_ nyh_, real_8 _p_ yh1_, real_8 _p_ ewt_, real_8 
 15924  _p_ savf_, real_8 _p_ savx_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)(), real_4 
 15925  (*_psol)())
 15926  { // ** body not listed **
 16406  }
 16407  
 16408  real_8 _dumach (void)
 16409  { // ** body not listed **
 16423  }
 16424  
 16425  int_4 _dumsum (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_)
 16426  { // ** body not listed **
 16430  }
 16431  
 16432  int_4 _dusol (int_4 _p_ neq_, real_8 _p_ tn_, real_8 _p_ y_, real_8 _p_ savf_, real_8 _p_ b_, real_8 _p_ wght_, int_4 
 16433  _p_ n_, real_8 _p_ delta_, real_8 _p_ hl0_, int_4 _p_ mnewt_, int_4 (*_psol)(), int_4 _p_ npsl_, real_8 _p_ x_, real_8 
 16434  _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 16435  { // ** body not listed **
 16473  }
 16474  
 16475  real_8 _dvnorm (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ w_)
 16476  { // ** body not listed **
 16488  }
 16489  
 16490  int_4 _iumach (void)
 16491  { // ** body not listed **
 16496  }
 16497  
 16498  int_4 _ixsav (int_4 _p_ ipar_, int_4 _p_ ivalue_, logical_4 _p_ iset_)
 16499  { // ** body not listed **
 16540  }
 16541  
 16542  int_4 _jgroup (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ maxg_, int_4 _p_ ngrp_, int_4 _p_ igp_, int_4 _p_ 
 16543  jgp_, int_4 _p_ incl_, int_4 _p_ jdone_, int_4 _p_ ier_)
 16544  { // ** body not listed **
 16593  }
 16594  
 16595  int_4 _odrv (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, real_8 _p_ a_, int_4 _p_ p_, int_4 _p_ ip_, int_4 _p_ nsp_, 
 16596  int_4 _p_ isp_, int_4 _p_ path_, int_4 _p_ flag_)
 16597  { // ** body not listed **
 16642  }
 16643  
 16644  int_4 _md (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ max_, int_4 _p_ v_, int_4 _p_ l_, int_4 _p_ head_, 
 16645  int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ flag_)
 16646  { // ** body not listed **
 16689  }
 16690  
 16691  int_4 _mdi (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ max_, int_4 _p_ v_, int_4 _p_ l_, int_4 _p_ head_, 
 16692  int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ tag_, int_4 _p_ flag_)
 16693  { // ** body not listed **
 16765  }
 16766  
 16767  int_4 _mdm (int_4 _p_ vk_, int_4 _p_ tail_, int_4 _p_ v_, int_4 _p_ l_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ 
 16768  mark_)
 16769  { // ** body not listed **
 16811  }
 16812  
 16813  int_4 _mdp (int_4 _p_ k_, int_4 _p_ ek_, int_4 _p_ tail_, int_4 _p_ v_, int_4 _p_ l_, int_4 _p_ head_, int_4 _p_ last_, 
 16814  int_4 _p_ next_, int_4 _p_ mark_)
 16815  { // ** body not listed **
 16903  }
 16904  
 16905  int_4 _mdu (int_4 _p_ ek_, int_4 _p_ dmin_, int_4 _p_ v_, int_4 _p_ l_, int_4 _p_ head_, int_4 _p_ last_, int_4 _p_ 
 16906  next_, int_4 _p_ mark_)
 16907  { // ** body not listed **
 16994  }
 16995  
 16996  int_4 _sro (int_4 _p_ n_, int_4 _p_ ip_, int_4 _p_ ia_, int_4 _p_ ja_, real_8 _p_ a_, int_4 _p_ q_, int_4 _p_ r_, 
 16997  logical_4 _p_ dflag_)
 16998  { // ** body not listed **
 17064  }
 17065  
 17066  int_4 _xerrwd (char _p_ msg_, int_4 _p_ nmes_, int_4 _p_ nerr_, int_4 _p_ level_, int_4 _p_ ni_, int_4 _p_ i1_, int_4 
 17067  _p_ i2_, int_4 _p_ nr_, real_8 _p_ r1_, real_8 _p_ r2_)
 17068  { // ** body not listed **
 17318  }
 17319  


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