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.10
    84  */
   506  static CALLS __calls[__ncalls] = {
   507    {"adjlr", 0}, // subroutine
   508    {"cdrv", 0}, // subroutine
   509    {"nroc", 0}, // subroutine
   510    {"nsfc", 0}, // subroutine
   511    {"nnfc", 0}, // subroutine
   512    {"nnsc", 0}, // subroutine
   513    {"nntc", 0}, // subroutine
   514    {"cntnzu", 0}, // subroutine
   515    {"daigbt", 0}, // subroutine
   516    {"dainvg", 0}, // subroutine
   517    {"dainvgs", 0}, // subroutine
   518    {"datp", 0}, // subroutine
   519    {"datv", 0}, // subroutine
   520    {"dbnorm", 0}, // real*8 function
   521    {"dcfode", 0}, // subroutine
   522    {"ddecbt", 0}, // subroutine
   523    {"dewset", 0}, // subroutine
   524    {"dfnorm", 0}, // real*8 function
   525    {"dhefa", 0}, // subroutine
   526    {"dhels", 0}, // subroutine
   527    {"dheqr", 0}, // subroutine
   528    {"dhesl", 0}, // subroutine
   529    {"dintdy", 0}, // subroutine
   530    {"diprep", 0}, // subroutine
   531    {"diprepi", 0}, // subroutine
   532    {"dlhin", 0}, // subroutine
   533    {"dlsoda", 0}, // subroutine
   534    {"dlsodar", 0}, // subroutine
   535    {"dlsode", 0}, // subroutine
   536    {"dlsodes", 0}, // subroutine
   537    {"dlsoibt", 0}, // subroutine
   538    {"dlsodi", 0}, // subroutine
   539    {"dlsodis", 0}, // subroutine
   540    {"dlsodkr", 0}, // subroutine
   541    {"dlsodpk", 0}, // subroutine
   542    {"dmnorm", 0}, // real*8 function
   543    {"dorthog", 0}, // subroutine
   544    {"dpcg", 0}, // subroutine
   545    {"dpcgs", 0}, // subroutine
   546    {"dpjibt", 0}, // subroutine
   547    {"dpkset", 0}, // subroutine
   548    {"dprep", 0}, // subroutine
   549    {"dprepi", 0}, // subroutine
   550    {"dprepj", 0}, // subroutine
   551    {"dprepji", 0}, // subroutine
   552    {"dprja", 0}, // subroutine
   553    {"dprjis", 0}, // subroutine
   554    {"dprjs", 0}, // subroutine
   555    {"drchek", 0}, // subroutine
   556    {"droots", 0}, // subroutine
   557    {"dsetpk", 0}, // subroutine
   558    {"dslsbt", 0}, // subroutine
   559    {"dsolbt", 0}, // subroutine
   560    {"dsolpk", 0}, // subroutine
   561    {"dsolss", 0}, // subroutine
   562    {"dsolsy", 0}, // subroutine
   563    {"dspigmr", 0}, // subroutine
   564    {"dspiom", 0}, // subroutine
   565    {"dsrcar", 0}, // subroutine
   566    {"dsrckr", 0}, // subroutine
   567    {"dsrcma", 0}, // subroutine
   568    {"dsrcms", 0}, // subroutine
   569    {"dsrcom", 0}, // subroutine
   570    {"dsrcpk", 0}, // subroutine
   571    {"dstoda", 0}, // subroutine
   572    {"dstode", 0}, // subroutine
   573    {"dstodi", 0}, // subroutine
   574    {"dstodpk", 0}, // subroutine
   575    {"dstoka", 0}, // subroutine
   576    {"dumach", 0}, // real*8 function
   577    {"dumsum", 0}, // subroutine
   578    {"dusol", 0}, // subroutine
   579    {"dvnorm", 0}, // real*8 function
   580    {"iumach", 0}, // integer*4 function
   581    {"ixsav", 0}, // integer*4 function
   582    {"jgroup", 0}, // subroutine
   583    {"odrv", 0}, // subroutine
   584    {"md", 0}, // subroutine
   585    {"mdi", 0}, // subroutine
   586    {"mdm", 0}, // subroutine
   587    {"mdp", 0}, // subroutine
   588    {"mdu", 0}, // subroutine
   589    {"sro", 0}, // subroutine
   590    {"xerrwd", 0}, // subroutine
   591    {NULL, 0}
   592  };
   593  static struct {
   594    int_4 iplost_, iesp_, istatc_, iys_, iba_, ibian_, ibjan_, ibjgp_, ipian_, ipjan_, ipjgp_, ipigp_, ipr_, ipc_, ipic_, 
   595    ipisp_, iprsp_, ipa_, lenyh_, lenyhm_, lenwk_, lreq_, lrat_, lrest_, lwmin_, moss_, msbj_, nslj_, ngp_, nlu_, nnz_, 
   596    nsp_, nzl_, nzu_, ilss_[34];
   597    real_8 rlss_[6], con0_, conmin_, ccmxj_, psmall_, rbig_, seth_;
   598  } dlss01_;
   599  static struct {
   600    int_4 iownd_[6], iowns_[6], icf_, ierpj_, iersl_, jcur_, jstart_, kflag_, l_, lyh_, lewt_, lacor_, lsavf_, lwm_, liwm_, 
   601    meth_, miter_, maxord_, maxcor_, msbp_, mxncf_, n_, nq_, nst_, nfe_, nje_, nqu_, init_, mxstep_, mxhnil_, nhnil_, 
   602    nslast_, nyh_, ils_[37], ialth_, ipup_, lmax_, meo_, nqnyh_, nslp_;
   603    real_8 rowns_[209], ccmax_, el0_, h_, hmin_, hmxi_, hu_, rc_, tn_, uround_, rls_[218], conit_, crate_, el_[13], 
   604    elco_[156], hold_, rmax_, tesco_[36];
   605  } dls001_;
   606  static struct {
   607    int_4 insufr_, insufi_, ixpr_, iowns2_[2], jtyp_, mused_, mxordn_, mxords_, iownd2_[3], ilsa_[9], icount_, irflag_;
   608    real_8 tsw_, rowns2_[20], pdnorm_, rownd2_, rlsa_[22], cm1_[12], cm2_[5], pdest_, pdlast_, ratio_;
   609  } dlsa01_;
   610  static struct {
   611    int_4 lg0_, lg1_, lgx_, iownr3_[2], irfnd_, itaskc_, ngc_, nge_, iownd3_[3], imax_, last_, idum3_[4], ilsr_[9];
   612    real_8 rownr3_[2], t0_, tlast_, toutc_, alpha_, x2_, rdum3_[3], rlsr_[5];
   613  } dlsr01_;
   614  static struct {
   615    int_4 newt_, nsfi_, nslj_, njev_, ils2_[4];
   616    real_8 stifr_, rls2_;
   617  } dls002_;
   618  static struct {
   619    int_4 jpre_, jacflg_, locwp_, lociwp_, lsavx_, kmp_, maxl_, mnewt_, nni_, nli_, nps_, ncfn_, ncfl_, ilsp_[13];
   620    real_8 delt_, epcon_, sqrtn_, rsqrtn_, rlsp_[4];
   621  } dlpk01_;
   624  _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_, 
   625  int_4 _p_ flag_);
   627  _p_ ar_, int_4 _p_ p_, int_4 _p_ flag_);
   629  _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_, 
   630  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 
   631  _p_ flag_);
   633  _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_, 
   634  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 
   635  _p_ tmp_, int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ flag_);
   637  _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_, 
   638  real_8 _p_ tmp_);
   640  _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_, 
   641  real_8 _p_ tmp_);
   644  ydot_, int_4 _p_ mb_, int_4 _p_ nb_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_);
   646  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_);
   648  , real_8 _p_ ydot_, int_4 _p_ ier_, int_4 (*_res)(), int_4 (*_adda)());
   650  hl0_, real_8 _p_ wk_, int_4 (*_f)(), real_8 _p_ w_);
   652  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_ 
   653  hl0_, int_4 _p_ jpre_, int_4 _p_ ier_, int_4 _p_ npsl_);
   657  _p_ ier_);
   659  _p_ ewt_);
   666  );
   668  ipflag_, real_4 (*_f)(), int_4 (*_jac)());
   670  , int_4 _p_ ic_, int_4 _p_ jc_, int_4 _p_ ipflag_, real_4 (*_res)(), int_4 (*_jac)(), real_4 (*_adda)());
   672  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_ 
   673  temp_, real_8 _p_ h0_, int_4 _p_ niter_, int_4 _p_ ier_);
   675  , 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 
   676  _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_);
   678  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_, 
   679  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 
   680  _p_ jroot_);
   682  , 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 
   683  _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_);
   685  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_, 
   686  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_);
   688  _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_, 
   689  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_)
   690  ;
   692  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 
   693  _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_);
   695  _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_, 
   696  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_)
   697  ;
   699  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_, 
   700  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)(), 
   701  int_4 _p_ ng_, int_4 _p_ jroot_);
   703  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_, 
   704  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_);
   707  ldhes_, int_4 _p_ kmp_, real_8 _p_ snormw_);
   709  , 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)(), 
   710  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 
   711  _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   713  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 
   714  (*_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_
   715  , real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   717  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 
   718  (*_adda)());
   720  savf_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)());
   722  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)());
   724  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_, 
   725  int_4 _p_ ipper_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)());
   727  ftem_, real_8 _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)());
   729  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 
   730  (*_adda)());
   732  ftem_, real_8 _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)());
   734  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 
   735  (*_adda)());
   737  ftem_, real_8 _p_ savf_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_f)(), int_4 (*_jac)());
   739  real_8 _p_ g0_, real_8 _p_ g1_, real_8 _p_ gx_, int_4 _p_ jroot_, int_4 _p_ irt_);
   741  g0_, real_8 _p_ g1_, real_8 _p_ gx_, real_8 _p_ x_, int_4 _p_ jroot_);
   743  savf_, int_4 _p_ jok_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)());
   746  _p_ ip_);
   748  wm_, int_4 _p_ iwm_, real_4 (*_f)(), real_4 (*_psol)());
   752  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_ 
   753  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_ 
   754  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_)
   755  ;
   757  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_ 
   758  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_ 
   759  ipvt_, int_4 _p_ liom_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   767  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 
   768  (*_pjac)(), int_4 (*_slvs)());
   770  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 
   771  (*_pjac)(), int_4 (*_slvs)());
   773  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 
   774  (*_adda)(), int_4 (*_jac)(), int_4 (*_pjac)(), int_4 (*_slvs)());
   776  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 
   777  (*_jac)(), real_4 (*_psol)());
   779  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 
   780  (*_jac)(), real_4 (*_psol)());
   784  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 
   785  _p_ x_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   790  int_4 _p_ jgp_, int_4 _p_ incl_, int_4 _p_ jdone_, int_4 _p_ ier_);
   792  _p_ nsp_, int_4 _p_ isp_, int_4 _p_ path_, int_4 _p_ flag_);
   794  head_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ flag_);
   796  head_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ tag_, int_4 _p_ flag_);
   798  int_4 _p_ mark_);
   800  _p_ last_, int_4 _p_ next_, int_4 _p_ mark_);
   802  int_4 _p_ next_, int_4 _p_ mark_);
   804  r_, logical_4 _p_ dflag_);
   806  i1_, int_4 _p_ i2_, int_4 _p_ nr_, real_8 _p_ r1_, real_8 _p_ r2_);
   807  int_4 _adjlr (int_4 _p_ n_, int_4 _p_ isp_, int_4 _p_ ldif_)
   808  { // ** body not listed **
   819  }
   820  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_, 
   821  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 
   822  _p_ flag_)
   823  { // ** body not listed **
   941  }
   942  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_, 
   943  int_4 _p_ p_, int_4 _p_ flag_)
   944  { // ** body not listed **
   986  }
   987  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_, 
   988  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_ 
   989  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_)
   990  { // ** body not listed **
  1379  }
  1380  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_, 
  1381  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_ 
  1382  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_, 
  1383  int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ flag_)
  1384  { // ** body not listed **
  1537  }
  1538  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_, 
  1539  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_ 
  1540  tmp_)
  1541  { // ** body not listed **
  1580  }
  1581  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_, 
  1582  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_ 
  1583  tmp_)
  1584  { // ** body not listed **
  1622  }
  1623  int_4 _cntnzu (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ nzsut_)
  1624  { // ** body not listed **
  1664  }
  1665  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 
  1666  _p_ mb_, int_4 _p_ nb_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_)
  1667  { // ** body not listed **
  1693  }
  1694  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 
  1695  _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_)
  1696  { // ** body not listed **
  1744  }
  1745  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 
  1746  _p_ ydot_, int_4 _p_ ier_, int_4 (*_res)(), int_4 (*_adda)())
  1747  { // ** body not listed **
  1802  }
  1803  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 
  1804  _p_ wk_, int_4 (*_f)(), real_8 _p_ w_)
  1805  { // ** body not listed **
  1823  }
  1824  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 
  1825  (*_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_ 
  1826  jpre_, int_4 _p_ ier_, int_4 _p_ npsl_)
  1827  { // ** body not listed **
  1886  }
  1887  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_)
  1888  { // ** body not listed **
  1906  }
  1907  int_4 _dcfode (int_4 _p_ meth_, real_8 _p_ elco_, real_8 _p_ tesco_)
  1908  { // ** body not listed **
  1983  }
  1984  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_)
  1985  { // ** body not listed **
  2059  }
  2060  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_)
  2061  { // ** body not listed **
  2092  }
  2093  real_8 _dfnorm (int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_)
  2094  { // ** body not listed **
  2109  }
  2110  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_)
  2111  { // ** body not listed **
  2216  }
  2217  int_4 _dhels (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, real_8 _p_ q_, real_8 _p_ b_)
  2218  { // ** body not listed **
  2240  }
  2241  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_)
  2242  { // ** body not listed **
  2338  }
  2339  int_4 _dhesl (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, int_4 _p_ ipvt_, real_8 _p_ b_)
  2340  { // ** body not listed **
  2368  }
  2369  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_)
  2370  { // ** body not listed **
  2445  }
  2446  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_, 
  2447  real_4 (*_f)(), int_4 (*_jac)())
  2448  { // ** body not listed **
  2488  }
  2489  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 
  2490  _p_ ic_, int_4 _p_ jc_, int_4 _p_ ipflag_, real_4 (*_res)(), int_4 (*_jac)(), real_4 (*_adda)())
  2491  { // ** body not listed **
  2532  }
  2533  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_ 
  2534  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 
  2535  _p_ h0_, int_4 _p_ niter_, int_4 _p_ ier_)
  2536  { // ** body not listed **
  2653  }
  2654  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 
  2655  _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_, 
  2656  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_)
  2657  { // ** body not listed **
  3539  }
  3540  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 
  3541  _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_, 
  3542  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_)
  3543  { // ** body not listed **
  4532  }
  4533  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 
  4534  _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_, 
  4535  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_)
  4536  { // ** body not listed **
  5285  }
  5286  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 
  5287  _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_, 
  5288  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_)
  5289  { // ** body not listed **
  6267  }
  6268  int_4 _dlsoibt (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  6269  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_ 
  6270  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_)
  6271  { // ** body not listed **
  7079  }
  7080  int_4 _dlsodi (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  7081  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_ 
  7082  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_)
  7083  { // ** body not listed **
  7902  }
  7903  int_4 _dlsodis (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  7904  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_ 
  7905  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_)
  7906  { // ** body not listed **
  8984  }
  8985  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 
  8986  _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_, 
  8987  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_, 
  8988  int_4 _p_ jroot_)
  8989  { // ** body not listed **
  9934  }
  9935  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 
  9936  _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_, 
  9937  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_)
  9938  { // ** body not listed **
 10795  }
 10796  real_8 _dmnorm (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ w_)
 10797  { // ** body not listed **
 10808  }
 10809  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 
 10810  _p_ kmp_, real_8 _p_ snormw_)
 10811  { // ** body not listed **
 10846  }
 10847  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 
 10848  _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 
 10849  (*_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_ 
 10850  wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 10851  { // ** body not listed **
 10939  }
 10940  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 
 10941  _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 
 10942  (*_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_ 
 10943  wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 10944  { // ** body not listed **
 11037  }
 11038  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 
 11039  _p_ savr_, real_8 _p_ s_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 11040  { // ** body not listed **
 11170  }
 11171  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_, 
 11172  real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)())
 11173  { // ** body not listed **
 11189  }
 11190  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 
 11191  _p_ ia_, int_4 _p_ ja_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 _p_ ipper_, int_4 (*_f)(), int_4 (*_jac)())
 11192  { // ** body not listed **
 11451  }
 11452  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 
 11453  _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_
 11454  , int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 11455  { // ** body not listed **
 11859  }
 11860  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 
 11861  _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)())
 11862  { // ** body not listed **
 12009  }
 12010  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_, 
 12011  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)())
 12012  { // ** body not listed **
 12163  }
 12164  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 
 12165  _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)())
 12166  { // ** body not listed **
 12292  }
 12293  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 
 12294  _p_ savr_, real_8 _p_ s_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 12295  { // ** body not listed **
 12405  }
 12406  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 
 12407  _p_ savf_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_f)(), int_4 (*_jac)())
 12408  { // ** body not listed **
 12577  }
 12578  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_ 
 12579  g0_, real_8 _p_ g1_, real_8 _p_ gx_, int_4 _p_ jroot_, int_4 _p_ irt_)
 12580  { // ** body not listed **
 12722  }
 12723  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_, 
 12724  real_8 _p_ g1_, real_8 _p_ gx_, real_8 _p_ x_, int_4 _p_ jroot_)
 12725  { // ** body not listed **
 12943  }
 12944  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_, 
 12945  int_4 _p_ jok_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)())
 12946  { // ** body not listed **
 12965  }
 12966  int_4 _dslsbt (real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ x_, real_8 _p_ tem_)
 12967  { // ** body not listed **
 12977  }
 12978  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_)
 12979  { // ** body not listed **
 13016  }
 13017  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 
 13018  _p_ iwm_, real_4 (*_f)(), real_4 (*_psol)())
 13019  { // ** body not listed **
 13154  }
 13155  int_4 _dsolss (real_8 _p_ wk_, int_4 _p_ iwk_, real_8 _p_ x_, real_8 _p_ tem_)
 13156  { // ** body not listed **
 13200  }
 13201  int_4 _dsolsy (real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ x_, real_8 _p_ tem_)
 13202  { // ** body not listed **
 13248  }
 13249  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 
 13250  _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 
 13251  _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 
 13252  _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_)
 13253  { // ** body not listed **
 13397  }
 13398  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 
 13399  _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 
 13400  (*_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_ 
 13401  liom_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 13402  { // ** body not listed **
 13535  }
 13536  int_4 _dsrcar (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13537  { // ** body not listed **
 13646  }
 13647  int_4 _dsrckr (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13648  { // ** body not listed **
 13755  }
 13756  int_4 _dsrcma (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13757  { // ** body not listed **
 13830  }
 13831  int_4 _dsrcms (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13832  { // ** body not listed **
 13905  }
 13906  int_4 _dsrcom (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13907  { // ** body not listed **
 13950  }
 13951  int_4 _dsrcpk (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13952  { // ** body not listed **
 14017  }
 14018  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 
 14019  _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 
 14020  (*_slvs)())
 14021  { // ** body not listed **
 14589  }
 14590  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 
 14591  _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 
 14592  (*_slvs)())
 14593  { // ** body not listed **
 15002  }
 15003  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 
 15004  _p_ savf_, real_8 _p_ savr_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), real_4 (*_adda)(), 
 15005  int_4 (*_jac)(), int_4 (*_pjac)(), int_4 (*_slvs)())
 15006  { // ** body not listed **
 15416  }
 15417  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 
 15418  _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 
 15419  (*_psol)())
 15420  { // ** body not listed **
 15842  }
 15843  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 
 15844  _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 
 15845  (*_psol)())
 15846  { // ** body not listed **
 16326  }
 16327  real_8 _dumach (void)
 16328  { // ** body not listed **
 16342  }
 16343  int_4 _dumsum (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_)
 16344  { // ** body not listed **
 16348  }
 16349  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 
 16350  _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 
 16351  _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 16352  { // ** body not listed **
 16390  }
 16391  real_8 _dvnorm (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ w_)
 16392  { // ** body not listed **
 16404  }
 16405  int_4 _iumach (void)
 16406  { // ** body not listed **
 16411  }
 16412  int_4 _ixsav (int_4 _p_ ipar_, int_4 _p_ ivalue_, logical_4 _p_ iset_)
 16413  { // ** body not listed **
 16454  }
 16455  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_ 
 16456  jgp_, int_4 _p_ incl_, int_4 _p_ jdone_, int_4 _p_ ier_)
 16457  { // ** body not listed **
 16506  }
 16507  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_, 
 16508  int_4 _p_ isp_, int_4 _p_ path_, int_4 _p_ flag_)
 16509  { // ** body not listed **
 16554  }
 16555  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_, 
 16556  int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ flag_)
 16557  { // ** body not listed **
 16600  }
 16601  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_, 
 16602  int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ tag_, int_4 _p_ flag_)
 16603  { // ** body not listed **
 16675  }
 16676  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_ 
 16677  mark_)
 16678  { // ** body not listed **
 16720  }
 16721  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_, 
 16722  int_4 _p_ next_, int_4 _p_ mark_)
 16723  { // ** body not listed **
 16811  }
 16812  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_ 
 16813  next_, int_4 _p_ mark_)
 16814  { // ** body not listed **
 16901  }
 16902  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_, 
 16903  logical_4 _p_ dflag_)
 16904  { // ** body not listed **
 16970  }
 16971  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 
 16972  _p_ i2_, int_4 _p_ nr_, real_8 _p_ r1_, real_8 _p_ r2_)
 16973  { // ** body not listed **
 17223  }


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