mathlib-slatec-depac.c

     1  //! @file mathlib-slatec-depac.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  //! 8 subprograms from SLATEC-DEPAC
    24  
    25  // Compiled from Fortran source code by VIF.
    26  // Selected subprograms are:
    27  //
    28  // DDEABM DDERKF DEABM  DERKF  DINTP  DSTEPS SINTRP STEPS
    29  
    30  // The license for SLATEC Fortran source code is:
    31  //
    32  // The SLATEC Common Mathematical Library was developed at
    33  // US government research laboratories and is in the public domain.
    34  // 
    35  // Repository: http://www.netlib.org/slatec/
    36  // 
    37  // The SLATEC common mathematical library is issued by the following
    38  // 
    39  //         Air Force Weapons Laboratory, Albuquerque
    40  //         Lawrence Livermore National Laboratory, Livermore
    41  //         Los Alamos National Laboratory, Los Alamos
    42  //         National Institute of Standards and Technology, Washington
    43  //         National Energy Research Supercomputer Center, Livermore
    44  //         Oak Ridge National Laboratory, Oak Ridge
    45  //         Sandia National Laboratories, Albuquerque
    46  //         Sandia National Laboratories, Livermore
    47  // 
    48  // All questions concerning the distribution of the library should be
    49  // directed to the NATIONAL ENERGY SOFTWARE CENTER, 9700 Cass Ave.,
    50  // Argonne, Illinois  60439, and not to the authors of the subprograms.
    51  // 
    52  //                  * * * * * Notice * * * * *
    53  // 
    54  // This material was prepared as an account of work sponsored by the
    55  // United States Government.  Neither the United States, nor the
    56  // Department of Energy, nor the Department of Defense, nor any of
    57  // their employees, nor any of their contractors, subcontractors, or
    58  // their employees, makes any warranty, expressed or implied, or
    59  // assumes any legal liability or responsibility for the accuracy,
    60  // completeness, or usefulness of any information, apparatus, product,
    61  // or process disclosed, or represents that its use would not infringe
    62  // upon privately owned rights.
    63  
    64  /*
    65  Generated by VIF - experimental VIntage Fortran compiler.
    66  VIF release 1.2.17
    67  */
    68  #if defined (__GNUC__)
    69  #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
    70  #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
    71  #if (__GNUC__ >= 14)
    72  #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
    73  #pragma GCC diagnostic ignored "-Wimplicit-int"
    74  #pragma GCC diagnostic ignored "-Wint-conversion"
    75  #pragma GCC diagnostic ignored "-Wreturn-mismatch"
    76  #endif
    77  #else
    78  #error VIF requires GCC
    79  #endif
   119  static CALLS __calls[__ncalls] = {
   120    {"ddeabm", 0}, // subroutine
   121    {"dderkf", 0}, // subroutine
   122    {"deabm", 0}, // subroutine
   123    {"derkf", 0}, // subroutine
   124    {"dintp", 0}, // subroutine
   125    {"dsteps", 0}, // subroutine
   126    {"sintrp", 0}, // subroutine
   127    {"steps", 0}, // subroutine
   128    {NULL, 0}
   129  };
   130  
   132  info_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ idid_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 
   133  _p_ liw_, real_8 _p_ rpar_, int_4 _p_ ipar_);
   135  info_, real_8 _p_ rtol_, real_8 _p_ atol_, int_4 _p_ idid_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 
   136  _p_ liw_, real_8 _p_ rpar_, int_4 _p_ ipar_);
   138  , real_4 _p_ rtol_, real_4 _p_ atol_, int_4 _p_ idid_, real_4 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ 
   139  liw_, real_4 _p_ rpar_, int_4 _p_ ipar_);
   141  , real_4 _p_ rtol_, real_4 _p_ atol_, int_4 _p_ idid_, real_4 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ 
   142  liw_, real_4 _p_ rpar_, int_4 _p_ ipar_);
   144  neqn_, int_4 _p_ kold_, real_8 _p_ phi_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_8 _p_ gi_, real_8 _p_ 
   145  alpha_, real_8 _p_ og_, real_8 _p_ ow_, real_8 _p_ ox_, real_8 _p_ oy_);
   147  real_8 _p_ wt_, logical_4 _p_ start_, real_8 _p_ hold_, int_4 _p_ k_, int_4 _p_ kold_, logical_4 _p_ crash_, real_8 _p_ 
   148  phi_, real_8 _p_ p_, real_8 _p_ yp_, real_8 _p_ psi_, real_8 _p_ alpha_, real_8 _p_ beta_, real_8 _p_ sig_, real_8 _p_ 
   149  v_, real_8 _p_ w_, real_8 _p_ g_, logical_4 _p_ phase1_, int_4 _p_ ns_, logical_4 _p_ nornd_, int_4 _p_ ksteps_, real_8 
   150  _p_ twou_, real_8 _p_ fouru_, real_8 _p_ xold_, int_4 _p_ kprev_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_8 
   151  _p_ gi_, real_8 _p_ rpar_, int_4 _p_ ipar_);
   153  neqn_, int_4 _p_ kold_, real_4 _p_ phi_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_4 _p_ gi_, real_4 _p_ 
   154  alpha_, real_4 _p_ og_, real_4 _p_ ow_, real_4 _p_ ox_, real_4 _p_ oy_);
   156  real_4 _p_ wt_, logical_4 _p_ start_, real_4 _p_ hold_, int_4 _p_ k_, int_4 _p_ kold_, logical_4 _p_ crash_, real_4 _p_ 
   157  phi_, real_4 _p_ p_, real_4 _p_ yp_, real_4 _p_ psi_, real_4 _p_ alpha_, real_4 _p_ beta_, real_4 _p_ sig_, real_4 _p_ 
   158  v_, real_4 _p_ w_, real_4 _p_ g_, logical_4 _p_ phase1_, int_4 _p_ ns_, logical_4 _p_ nornd_, int_4 _p_ ksteps_, real_4 
   159  _p_ twou_, real_4 _p_ fouru_, real_4 _p_ xold_, int_4 _p_ kprev_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_4 
   160  _p_ gi_, real_4 _p_ rpar_, int_4 _p_ ipar_);
   161  int_4 _ddeabm (real_4 (*_df)(), int_4 _p_ neq_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ tout_, int_4 _p_ info_, real_8 
   162  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ idid_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, 
   163  real_8 _p_ rpar_, int_4 _p_ ipar_)
   164  { // ** body not listed **
   376  }
   377  
   378  int_4 _dderkf (real_4 (*_df)(), int_4 _p_ neq_, real_8 _p_ t_, real_8 _p_ y_, real_8 _p_ tout_, int_4 _p_ info_, real_8 
   379  _p_ rtol_, real_8 _p_ atol_, int_4 _p_ idid_, real_8 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, 
   380  real_8 _p_ rpar_, int_4 _p_ ipar_)
   381  { // ** body not listed **
   565  }
   566  
   567  int_4 _deabm (real_4 (*_f)(), int_4 _p_ neq_, real_4 _p_ t_, real_4 _p_ y_, real_4 _p_ tout_, int_4 _p_ info_, real_4 
   568  _p_ rtol_, real_4 _p_ atol_, int_4 _p_ idid_, real_4 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, 
   569  real_4 _p_ rpar_, int_4 _p_ ipar_)
   570  { // ** body not listed **
   782  }
   783  
   784  int_4 _derkf (real_4 (*_f)(), int_4 _p_ neq_, real_4 _p_ t_, real_4 _p_ y_, real_4 _p_ tout_, int_4 _p_ info_, real_4 
   785  _p_ rtol_, real_4 _p_ atol_, int_4 _p_ idid_, real_4 _p_ rwork_, int_4 _p_ lrw_, int_4 _p_ iwork_, int_4 _p_ liw_, 
   786  real_4 _p_ rpar_, int_4 _p_ ipar_)
   787  { // ** body not listed **
   971  }
   972  
   973  int_4 _dintp (real_8 _p_ x_, real_8 _p_ y_, real_8 _p_ xout_, real_8 _p_ yout_, real_8 _p_ ypout_, int_4 _p_ neqn_, 
   974  int_4 _p_ kold_, real_8 _p_ phi_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_8 _p_ gi_, real_8 _p_ alpha_, 
   975  real_8 _p_ og_, real_8 _p_ ow_, real_8 _p_ ox_, real_8 _p_ oy_)
   976  { // ** body not listed **
  1059  }
  1060  
  1061  int_4 _dsteps (int_4 (*_df)(), int_4 _p_ neqn_, real_8 _p_ y_, real_8 _p_ x_, real_8 _p_ h_, real_8 _p_ eps_, real_8 
  1062  _p_ wt_, logical_4 _p_ start_, real_8 _p_ hold_, int_4 _p_ k_, int_4 _p_ kold_, logical_4 _p_ crash_, real_8 _p_ phi_, 
  1063  real_8 _p_ p_, real_8 _p_ yp_, real_8 _p_ psi_, real_8 _p_ alpha_, real_8 _p_ beta_, real_8 _p_ sig_, real_8 _p_ v_, 
  1064  real_8 _p_ w_, real_8 _p_ g_, logical_4 _p_ phase1_, int_4 _p_ ns_, logical_4 _p_ nornd_, int_4 _p_ ksteps_, real_8 _p_ 
  1065  twou_, real_8 _p_ fouru_, real_8 _p_ xold_, int_4 _p_ kprev_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_8 _p_ 
  1066  gi_, real_8 _p_ rpar_, int_4 _p_ ipar_)
  1067  { // ** body not listed **
  1631  }
  1632  
  1633  int_4 _sintrp (real_4 _p_ x_, real_4 _p_ y_, real_4 _p_ xout_, real_4 _p_ yout_, real_4 _p_ ypout_, int_4 _p_ neqn_, 
  1634  int_4 _p_ kold_, real_4 _p_ phi_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_4 _p_ gi_, real_4 _p_ alpha_, 
  1635  real_4 _p_ og_, real_4 _p_ ow_, real_4 _p_ ox_, real_4 _p_ oy_)
  1636  { // ** body not listed **
  1719  }
  1720  
  1721  int_4 _steps (int_4 (*_f)(), int_4 _p_ neqn_, real_4 _p_ y_, real_4 _p_ x_, real_4 _p_ h_, real_4 _p_ eps_, real_4 _p_ 
  1722  wt_, logical_4 _p_ start_, real_4 _p_ hold_, int_4 _p_ k_, int_4 _p_ kold_, logical_4 _p_ crash_, real_4 _p_ phi_, 
  1723  real_4 _p_ p_, real_4 _p_ yp_, real_4 _p_ psi_, real_4 _p_ alpha_, real_4 _p_ beta_, real_4 _p_ sig_, real_4 _p_ v_, 
  1724  real_4 _p_ w_, real_4 _p_ g_, logical_4 _p_ phase1_, int_4 _p_ ns_, logical_4 _p_ nornd_, int_4 _p_ ksteps_, real_4 _p_ 
  1725  twou_, real_4 _p_ fouru_, real_4 _p_ xold_, int_4 _p_ kprev_, int_4 _p_ ivc_, int_4 _p_ iv_, int_4 _p_ kgi_, real_4 _p_ 
  1726  gi_, real_4 _p_ rpar_, int_4 _p_ ipar_)
  1727  { // ** body not listed **
  2291  }
  2292  


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