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.12
    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  
   594  static struct {
   595    int_4 iplost_, iesp_, istatc_, iys_, iba_, ibian_, ibjan_, ibjgp_, ipian_, ipjan_, ipjgp_, ipigp_, ipr_, ipc_, ipic_, 
   596    ipisp_, iprsp_, ipa_, lenyh_, lenyhm_, lenwk_, lreq_, lrat_, lrest_, lwmin_, moss_, msbj_, nslj_, ngp_, nlu_, nnz_, 
   597    nsp_, nzl_, nzu_, ilss_[34];
   598    real_8 rlss_[6], con0_, conmin_, ccmxj_, psmall_, rbig_, seth_;
   599  } dlss01_;
   600  static struct {
   601    int_4 iownd_[6], iowns_[6], icf_, ierpj_, iersl_, jcur_, jstart_, kflag_, l_, lyh_, lewt_, lacor_, lsavf_, lwm_, liwm_, 
   602    meth_, miter_, maxord_, maxcor_, msbp_, mxncf_, n_, nq_, nst_, nfe_, nje_, nqu_, init_, mxstep_, mxhnil_, nhnil_, 
   603    nslast_, nyh_, ils_[37], ialth_, ipup_, lmax_, meo_, nqnyh_, nslp_;
   604    real_8 rowns_[209], ccmax_, el0_, h_, hmin_, hmxi_, hu_, rc_, tn_, uround_, rls_[218], conit_, crate_, el_[13], 
   605    elco_[156], hold_, rmax_, tesco_[36];
   606  } dls001_;
   607  static struct {
   608    int_4 insufr_, insufi_, ixpr_, iowns2_[2], jtyp_, mused_, mxordn_, mxords_, iownd2_[3], ilsa_[9], icount_, irflag_;
   609    real_8 tsw_, rowns2_[20], pdnorm_, rownd2_, rlsa_[22], cm1_[12], cm2_[5], pdest_, pdlast_, ratio_;
   610  } dlsa01_;
   611  static struct {
   612    int_4 lg0_, lg1_, lgx_, iownr3_[2], irfnd_, itaskc_, ngc_, nge_, iownd3_[3], imax_, last_, idum3_[4], ilsr_[9];
   613    real_8 rownr3_[2], t0_, tlast_, toutc_, alpha_, x2_, rdum3_[3], rlsr_[5];
   614  } dlsr01_;
   615  static struct {
   616    int_4 newt_, nsfi_, nslj_, njev_, ils2_[4];
   617    real_8 stifr_, rls2_;
   618  } dls002_;
   619  static struct {
   620    int_4 jpre_, jacflg_, locwp_, lociwp_, lsavx_, kmp_, maxl_, mnewt_, nni_, nli_, nps_, ncfn_, ncfl_, ilsp_[13];
   621    real_8 delt_, epcon_, sqrtn_, rsqrtn_, rlsp_[4];
   622  } dlpk01_;
   625  _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_, 
   626  int_4 _p_ flag_);
   628  _p_ ar_, int_4 _p_ p_, int_4 _p_ flag_);
   630  _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_, 
   631  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 
   632  _p_ flag_);
   634  _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_, 
   635  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 
   636  _p_ tmp_, int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ flag_);
   638  _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_, 
   639  real_8 _p_ tmp_);
   641  _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_, 
   642  real_8 _p_ tmp_);
   645  ydot_, int_4 _p_ mb_, int_4 _p_ nb_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_);
   647  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_);
   649  , real_8 _p_ ydot_, int_4 _p_ ier_, int_4 (*_res)(), int_4 (*_adda)());
   651  hl0_, real_8 _p_ wk_, int_4 (*_f)(), real_8 _p_ w_);
   653  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_ 
   654  hl0_, int_4 _p_ jpre_, int_4 _p_ ier_, int_4 _p_ npsl_);
   658  _p_ ier_);
   660  _p_ ewt_);
   667  );
   669  ipflag_, real_4 (*_f)(), int_4 (*_jac)());
   671  , int_4 _p_ ic_, int_4 _p_ jc_, int_4 _p_ ipflag_, real_4 (*_res)(), int_4 (*_jac)(), real_4 (*_adda)());
   673  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_ 
   674  temp_, real_8 _p_ h0_, int_4 _p_ niter_, int_4 _p_ ier_);
   676  , 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 
   677  _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_);
   679  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_, 
   680  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 
   681  _p_ jroot_);
   683  , 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 
   684  _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_);
   686  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_, 
   687  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_);
   689  _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_, 
   690  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_)
   691  ;
   693  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 
   694  _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_);
   696  _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_, 
   697  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_)
   698  ;
   700  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_, 
   701  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)(), 
   702  int_4 _p_ ng_, int_4 _p_ jroot_);
   704  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_, 
   705  int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_);
   708  ldhes_, int_4 _p_ kmp_, real_8 _p_ snormw_);
   710  , 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)(), 
   711  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 
   712  _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   714  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 
   715  (*_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_
   716  , real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   718  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 
   719  (*_adda)());
   721  savf_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)());
   723  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)());
   725  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_, 
   726  int_4 _p_ ipper_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)());
   728  ftem_, real_8 _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)());
   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  ftem_, real_8 _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)());
   735  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 
   736  (*_adda)());
   738  ftem_, real_8 _p_ savf_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_f)(), int_4 (*_jac)());
   740  real_8 _p_ g0_, real_8 _p_ g1_, real_8 _p_ gx_, int_4 _p_ jroot_, int_4 _p_ irt_);
   742  g0_, real_8 _p_ g1_, real_8 _p_ gx_, real_8 _p_ x_, int_4 _p_ jroot_);
   744  savf_, int_4 _p_ jok_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)());
   747  _p_ ip_);
   749  wm_, int_4 _p_ iwm_, real_4 (*_f)(), real_4 (*_psol)());
   753  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_ 
   754  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_ 
   755  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_)
   756  ;
   758  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_ 
   759  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_ 
   760  ipvt_, int_4 _p_ liom_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   768  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 
   769  (*_pjac)(), int_4 (*_slvs)());
   771  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 
   772  (*_pjac)(), int_4 (*_slvs)());
   774  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 
   775  (*_adda)(), int_4 (*_jac)(), int_4 (*_pjac)(), int_4 (*_slvs)());
   777  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 
   778  (*_jac)(), real_4 (*_psol)());
   780  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 
   781  (*_jac)(), real_4 (*_psol)());
   785  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 
   786  _p_ x_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_);
   791  int_4 _p_ jgp_, int_4 _p_ incl_, int_4 _p_ jdone_, int_4 _p_ ier_);
   793  _p_ nsp_, int_4 _p_ isp_, int_4 _p_ path_, int_4 _p_ flag_);
   795  head_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ flag_);
   797  head_, int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ tag_, int_4 _p_ flag_);
   799  int_4 _p_ mark_);
   801  _p_ last_, int_4 _p_ next_, int_4 _p_ mark_);
   803  int_4 _p_ next_, int_4 _p_ mark_);
   805  r_, logical_4 _p_ dflag_);
   807  i1_, int_4 _p_ i2_, int_4 _p_ nr_, real_8 _p_ r1_, real_8 _p_ r2_);
   808  int_4 _adjlr (int_4 _p_ n_, int_4 _p_ isp_, int_4 _p_ ldif_)
   809  { // ** body not listed **
   820  }
   821  
   822  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_, 
   823  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 
   824  _p_ flag_)
   825  { // ** body not listed **
   943  }
   944  
   945  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_, 
   946  int_4 _p_ p_, int_4 _p_ flag_)
   947  { // ** body not listed **
   989  }
   990  
   991  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_, 
   992  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_ 
   993  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_)
   994  { // ** body not listed **
  1383  }
  1384  
  1385  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_, 
  1386  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_ 
  1387  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_, 
  1388  int_4 _p_ irl_, int_4 _p_ jrl_, int_4 _p_ flag_)
  1389  { // ** body not listed **
  1542  }
  1543  
  1544  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_, 
  1545  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_ 
  1546  tmp_)
  1547  { // ** body not listed **
  1586  }
  1587  
  1588  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_, 
  1589  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_ 
  1590  tmp_)
  1591  { // ** body not listed **
  1629  }
  1630  
  1631  int_4 _cntnzu (int_4 _p_ n_, int_4 _p_ ia_, int_4 _p_ ja_, int_4 _p_ nzsut_)
  1632  { // ** body not listed **
  1672  }
  1673  
  1674  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 
  1675  _p_ mb_, int_4 _p_ nb_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_)
  1676  { // ** body not listed **
  1702  }
  1703  
  1704  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 
  1705  _p_ miter_, int_4 _p_ ml_, int_4 _p_ mu_, real_8 _p_ pw_, int_4 _p_ ipvt_, int_4 _p_ ier_)
  1706  { // ** body not listed **
  1754  }
  1755  
  1756  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 
  1757  _p_ ydot_, int_4 _p_ ier_, int_4 (*_res)(), int_4 (*_adda)())
  1758  { // ** body not listed **
  1813  }
  1814  
  1815  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 
  1816  _p_ wk_, int_4 (*_f)(), real_8 _p_ w_)
  1817  { // ** body not listed **
  1835  }
  1836  
  1837  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 
  1838  (*_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_ 
  1839  jpre_, int_4 _p_ ier_, int_4 _p_ npsl_)
  1840  { // ** body not listed **
  1899  }
  1900  
  1901  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_)
  1902  { // ** body not listed **
  1920  }
  1921  
  1922  int_4 _dcfode (int_4 _p_ meth_, real_8 _p_ elco_, real_8 _p_ tesco_)
  1923  { // ** body not listed **
  1998  }
  1999  
  2000  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_)
  2001  { // ** body not listed **
  2075  }
  2076  
  2077  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_)
  2078  { // ** body not listed **
  2109  }
  2110  
  2111  real_8 _dfnorm (int_4 _p_ n_, real_8 _p_ a_, real_8 _p_ w_)
  2112  { // ** body not listed **
  2127  }
  2128  
  2129  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_)
  2130  { // ** body not listed **
  2235  }
  2236  
  2237  int_4 _dhels (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, real_8 _p_ q_, real_8 _p_ b_)
  2238  { // ** body not listed **
  2260  }
  2261  
  2262  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_)
  2263  { // ** body not listed **
  2359  }
  2360  
  2361  int_4 _dhesl (real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ n_, int_4 _p_ ipvt_, real_8 _p_ b_)
  2362  { // ** body not listed **
  2390  }
  2391  
  2392  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_)
  2393  { // ** body not listed **
  2468  }
  2469  
  2470  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_, 
  2471  real_4 (*_f)(), int_4 (*_jac)())
  2472  { // ** body not listed **
  2512  }
  2513  
  2514  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 
  2515  _p_ ic_, int_4 _p_ jc_, int_4 _p_ ipflag_, real_4 (*_res)(), int_4 (*_jac)(), real_4 (*_adda)())
  2516  { // ** body not listed **
  2557  }
  2558  
  2559  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_ 
  2560  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 
  2561  _p_ h0_, int_4 _p_ niter_, int_4 _p_ ier_)
  2562  { // ** body not listed **
  2679  }
  2680  
  2681  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 
  2682  _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_, 
  2683  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ jt_)
  2684  { // ** body not listed **
  3566  }
  3567  
  3568  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 
  3569  _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_, 
  3570  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_)
  3571  { // ** body not listed **
  4560  }
  4561  
  4562  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 
  4563  _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_, 
  4564  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_)
  4565  { // ** body not listed **
  5314  }
  5315  
  5316  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 
  5317  _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_, 
  5318  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), int_4 _p_ mf_)
  5319  { // ** body not listed **
  6297  }
  6298  
  6299  int_4 _dlsoibt (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  6300  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_ 
  6301  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_)
  6302  { // ** body not listed **
  7110  }
  7111  
  7112  int_4 _dlsodi (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  7113  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_ 
  7114  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_)
  7115  { // ** body not listed **
  7934  }
  7935  
  7936  int_4 _dlsodis (int_4 (*_res)(), real_4 (*_adda)(), int_4 (*_jac)(), int_4 _p_ neq_, real_8 _p_ y_, real_8 _p_ ydoti_, 
  7937  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_ 
  7938  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_)
  7939  { // ** body not listed **
  9017  }
  9018  
  9019  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 
  9020  _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_, 
  9021  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_, 
  9022  int_4 _p_ jroot_)
  9023  { // ** body not listed **
  9968  }
  9969  
  9970  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 
  9971  _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_, 
  9972  int_4 _p_ iwork_, int_4 _p_ liw_, int_4 (*_jac)(), real_4 (*_psol)(), int_4 _p_ mf_)
  9973  { // ** body not listed **
 10830  }
 10831  
 10832  real_8 _dmnorm (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ w_)
 10833  { // ** body not listed **
 10844  }
 10845  
 10846  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 
 10847  _p_ kmp_, real_8 _p_ snormw_)
 10848  { // ** body not listed **
 10883  }
 10884  
 10885  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 
 10886  _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 
 10887  (*_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_ 
 10888  wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 10889  { // ** body not listed **
 10977  }
 10978  
 10979  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 
 10980  _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 
 10981  (*_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_ 
 10982  wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 10983  { // ** body not listed **
 11076  }
 11077  
 11078  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 
 11079  _p_ savr_, real_8 _p_ s_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 11080  { // ** body not listed **
 11210  }
 11211  
 11212  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_, 
 11213  real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)())
 11214  { // ** body not listed **
 11230  }
 11231  
 11232  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 
 11233  _p_ ia_, int_4 _p_ ja_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 _p_ ipper_, int_4 (*_f)(), int_4 (*_jac)())
 11234  { // ** body not listed **
 11493  }
 11494  
 11495  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 
 11496  _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_
 11497  , int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 11498  { // ** body not listed **
 11902  }
 11903  
 11904  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 
 11905  _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)())
 11906  { // ** body not listed **
 12053  }
 12054  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_, 
 12055  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)())
 12056  { // ** body not listed **
 12207  }
 12208  
 12209  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 
 12210  _p_ savf_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_f)(), int_4 (*_jac)())
 12211  { // ** body not listed **
 12337  }
 12338  
 12339  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 
 12340  _p_ savr_, real_8 _p_ s_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_res)(), int_4 (*_jac)(), int_4 (*_adda)())
 12341  { // ** body not listed **
 12451  }
 12452  
 12453  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 
 12454  _p_ savf_, real_8 _p_ wk_, int_4 _p_ iwk_, int_4 (*_f)(), int_4 (*_jac)())
 12455  { // ** body not listed **
 12624  }
 12625  
 12626  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_ 
 12627  g0_, real_8 _p_ g1_, real_8 _p_ gx_, int_4 _p_ jroot_, int_4 _p_ irt_)
 12628  { // ** body not listed **
 12770  }
 12771  
 12772  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_, 
 12773  real_8 _p_ g1_, real_8 _p_ gx_, real_8 _p_ x_, int_4 _p_ jroot_)
 12774  { // ** body not listed **
 12992  }
 12993  
 12994  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_, 
 12995  int_4 _p_ jok_, real_8 _p_ wm_, int_4 _p_ iwm_, real_4 (*_f)(), int_4 (*_jac)())
 12996  { // ** body not listed **
 13015  }
 13016  
 13017  int_4 _dslsbt (real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ x_, real_8 _p_ tem_)
 13018  { // ** body not listed **
 13028  }
 13029  
 13030  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_)
 13031  { // ** body not listed **
 13068  }
 13069  
 13070  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 
 13071  _p_ iwm_, real_4 (*_f)(), real_4 (*_psol)())
 13072  { // ** body not listed **
 13207  }
 13208  
 13209  int_4 _dsolss (real_8 _p_ wk_, int_4 _p_ iwk_, real_8 _p_ x_, real_8 _p_ tem_)
 13210  { // ** body not listed **
 13254  }
 13255  
 13256  int_4 _dsolsy (real_8 _p_ wm_, int_4 _p_ iwm_, real_8 _p_ x_, real_8 _p_ tem_)
 13257  { // ** body not listed **
 13303  }
 13304  
 13305  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 
 13306  _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 
 13307  _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 
 13308  _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_)
 13309  { // ** body not listed **
 13453  }
 13454  
 13455  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 
 13456  _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 
 13457  (*_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_ 
 13458  liom_, real_8 _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 13459  { // ** body not listed **
 13592  }
 13593  
 13594  int_4 _dsrcar (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13595  { // ** body not listed **
 13704  }
 13705  
 13706  int_4 _dsrckr (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13707  { // ** body not listed **
 13814  }
 13815  
 13816  int_4 _dsrcma (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13817  { // ** body not listed **
 13890  }
 13891  
 13892  int_4 _dsrcms (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13893  { // ** body not listed **
 13966  }
 13967  
 13968  int_4 _dsrcom (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 13969  { // ** body not listed **
 14012  }
 14013  
 14014  int_4 _dsrcpk (real_8 _p_ rsav_, int_4 _p_ isav_, int_4 _p_ job_)
 14015  { // ** body not listed **
 14080  }
 14081  
 14082  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 
 14083  _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 
 14084  (*_slvs)())
 14085  { // ** body not listed **
 14653  }
 14654  
 14655  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 
 14656  _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 
 14657  (*_slvs)())
 14658  { // ** body not listed **
 15067  }
 15068  
 15069  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 
 15070  _p_ savf_, real_8 _p_ savr_, real_8 _p_ acor_, real_8 _p_ wm_, int_4 _p_ iwm_, int_4 (*_res)(), real_4 (*_adda)(), 
 15071  int_4 (*_jac)(), int_4 (*_pjac)(), int_4 (*_slvs)())
 15072  { // ** body not listed **
 15482  }
 15483  
 15484  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 
 15485  _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 
 15486  (*_psol)())
 15487  { // ** body not listed **
 15909  }
 15910  
 15911  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 
 15912  _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 
 15913  (*_psol)())
 15914  { // ** body not listed **
 16394  }
 16395  
 16396  real_8 _dumach (void)
 16397  { // ** body not listed **
 16411  }
 16412  
 16413  int_4 _dumsum (real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ c_)
 16414  { // ** body not listed **
 16418  }
 16419  
 16420  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 
 16421  _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 
 16422  _p_ wp_, int_4 _p_ iwp_, real_8 _p_ wk_, int_4 _p_ iflag_)
 16423  { // ** body not listed **
 16461  }
 16462  
 16463  real_8 _dvnorm (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ w_)
 16464  { // ** body not listed **
 16476  }
 16477  
 16478  int_4 _iumach (void)
 16479  { // ** body not listed **
 16484  }
 16485  
 16486  int_4 _ixsav (int_4 _p_ ipar_, int_4 _p_ ivalue_, logical_4 _p_ iset_)
 16487  { // ** body not listed **
 16528  }
 16529  
 16530  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_ 
 16531  jgp_, int_4 _p_ incl_, int_4 _p_ jdone_, int_4 _p_ ier_)
 16532  { // ** body not listed **
 16581  }
 16582  
 16583  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_, 
 16584  int_4 _p_ isp_, int_4 _p_ path_, int_4 _p_ flag_)
 16585  { // ** body not listed **
 16630  }
 16631  
 16632  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_, 
 16633  int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ flag_)
 16634  { // ** body not listed **
 16677  }
 16678  
 16679  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_, 
 16680  int_4 _p_ last_, int_4 _p_ next_, int_4 _p_ mark_, int_4 _p_ tag_, int_4 _p_ flag_)
 16681  { // ** body not listed **
 16753  }
 16754  
 16755  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_ 
 16756  mark_)
 16757  { // ** body not listed **
 16799  }
 16800  
 16801  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_, 
 16802  int_4 _p_ next_, int_4 _p_ mark_)
 16803  { // ** body not listed **
 16891  }
 16892  
 16893  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_ 
 16894  next_, int_4 _p_ mark_)
 16895  { // ** body not listed **
 16982  }
 16983  
 16984  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_, 
 16985  logical_4 _p_ dflag_)
 16986  { // ** body not listed **
 17052  }
 17053  
 17054  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 
 17055  _p_ i2_, int_4 _p_ nr_, real_8 _p_ r1_, real_8 _p_ r2_)
 17056  { // ** body not listed **
 17306  }
 17307  


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