☰  Menu Navigation ‣ Home ‣ About Marcel ‣ Privacy and cookies Algol 68 Genie ‣ The Algol 68 Genie project ‣ Downloads and tutorials Vintage Fortran ‣ Vintage Fortran project Blog ‣ Eighty years of freedom 1945-2025 ‣ Thermodynamics ex machina ‣ Wireless headsets, HiFi and Linux ‣ Snoopy calendar 1969-2025 ‣ Wonderful WiFi World ‣ …

mathlib-pchip.c

 //! @file mathlib-pchip.c
 //
 //! @section copyright
 //
 // This file is part of VIF - vintage fortran compiler.
 // Copyright 2020-2025 J. Marcel van der Veer <algol68g@xs4all.nl>.
 //
 //! @section license
 //
 // This program is free software; you can redistribute it and/or modify it 
 // under the terms of the gnu general public license as published by the 
 // free software foundation; either version 3 of the license, or 
 // (at your option) any later version.
 //
 // This program is distributed in the hope that it will be useful, but 
 // without any warranty; without even the implied warranty of merchantability 
 // or fitness for a particular purpose. See the GNU general public license for 
 // more details. you should have received a copy of the GNU general public 
 // license along with this program. If not, see <http://www.gnu.org/licenses/>.
 //
 //! @Synopsis
 //
 //! 41 subprograms from PCHIP
 
 // Compiled from Fortran source code by VIF.
 // Selected subprograms are:
 //
 // CHFCM  DCHFDV DPCHCE DPCHDF DPCHIC DPCHSP PCHCE  PCHDF  PCHIA  PCHKT
 // CHFDV  DCHFEV DPCHCI DPCHFD DPCHID DPCHST PCHCI  PCHDOC PCHIC  PCHSP
 // CHFEV  DCHFIE DPCHCM DPCHFE DPCHIM DPCHSW PCHCM  PCHFD  PCHID  PCHST
 // CHFIE  DPCHBS DPCHCS DPCHIA DPCHKT PCHBS  PCHCS  PCHFE  PCHIM  PCHSW
 // DCHFCM
 
 // PCHIP routines from SLATEC.
 // 
 // Source: netlib.org/slatec/
 
 // The license for PCHIP Fortran source code is:
 //
 // The SLATEC Common Mathematical Library was developed at
 // US government research laboratories and is in the public domain.
 // 
 // Repository: http://www.netlib.org/slatec/
 // 
 // The SLATEC common mathematical library is issued by the following
 // 
 //         Air Force Weapons Laboratory, Albuquerque
 //         Lawrence Livermore National Laboratory, Livermore
 //         Los Alamos National Laboratory, Los Alamos
 //         National Institute of Standards and Technology, Washington
 //         National Energy Research Supercomputer Center, Livermore
 //         Oak Ridge National Laboratory, Oak Ridge
 //         Sandia National Laboratories, Albuquerque
 //         Sandia National Laboratories, Livermore
 // 
 // All questions concerning the distribution of the library should be
 // directed to the NATIONAL ENERGY SOFTWARE CENTER, 9700 Cass Ave.,
 // Argonne, Illinois  60439, and not to the authors of the subprograms.
 // 
 //                  * * * * * Notice * * * * *
 // 
 // This material was prepared as an account of work sponsored by the
 // United States Government.  Neither the United States, nor the
 // Department of Energy, nor the Department of Defense, nor any of
 // their employees, nor any of their contractors, subcontractors, or
 // their employees, makes any warranty, expressed or implied, or
 // assumes any legal liability or responsibility for the accuracy,
 // completeness, or usefulness of any information, apparatus, product,
 // or process disclosed, or represents that its use would not infringe
 // upon privately owned rights.
 
 /*
 Generated by VIF - experimental VIntage Fortran compiler.
 VIF release 1.2.17
 */
 #if defined (__GNUC__)
 #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
 #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
 #if (__GNUC__ >= 14)
 #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
 #pragma GCC diagnostic ignored "-Wimplicit-int"
 #pragma GCC diagnostic ignored "-Wint-conversion"
 #pragma GCC diagnostic ignored "-Wreturn-mismatch"
 #endif
 #else
 #error VIF requires GCC
 #endif
 static CALLS __calls[__ncalls] = {
   {"chfcm", 0}, // integer*4 function
   {"chfdv", 0}, // subroutine
   {"chfev", 0}, // subroutine
   {"chfie", 0}, // real*4 function
   {"dchfcm", 0}, // integer*4 function
   {"dchfdv", 0}, // subroutine
   {"dchfev", 0}, // subroutine
   {"dchfie", 0}, // real*8 function
   {"dpchbs", 0}, // subroutine
   {"dpchce", 0}, // subroutine
   {"dpchci", 0}, // subroutine
   {"dpchcm", 0}, // subroutine
   {"dpchcs", 0}, // subroutine
   {"dpchdf", 0}, // real*8 function
   {"dpchfd", 0}, // subroutine
   {"dpchfe", 0}, // subroutine
   {"dpchia", 0}, // real*8 function
   {"dpchic", 0}, // subroutine
   {"dpchid", 0}, // real*8 function
   {"dpchim", 0}, // subroutine
   {"dpchkt", 0}, // subroutine
   {"dpchsp", 0}, // subroutine
   {"dpchst", 0}, // real*8 function
   {"dpchsw", 0}, // subroutine
   {"pchbs", 0}, // subroutine
   {"pchce", 0}, // subroutine
   {"pchci", 0}, // subroutine
   {"pchcm", 0}, // subroutine
   {"pchcs", 0}, // subroutine
   {"pchdf", 0}, // real*4 function
   {"pchdoc", 0}, // subroutine
   {"pchfd", 0}, // subroutine
   {"pchfe", 0}, // subroutine
   {"pchia", 0}, // real*4 function
   {"pchic", 0}, // subroutine
   {"pchid", 0}, // real*4 function
   {"pchim", 0}, // subroutine
   {"pchkt", 0}, // subroutine
   {"pchsp", 0}, // subroutine
   {"pchst", 0}, // real*4 function
   {"pchsw", 0}, // subroutine
   {NULL, 0}
 };
 
 int_4 _p_ ne_, real_4 _p_ xe_, real_4 _p_ fe_, real_4 _p_ de_, int_4 _p_ next_, int_4 _p_ ierr_);
 int_4 _p_ ne_, real_4 _p_ xe_, real_4 _p_ fe_, int_4 _p_ next_, int_4 _p_ ierr_);
 , real_4 _p_ a_, real_4 _p_ b_);
 , int_4 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, real_8 _p_ de_, int_4 _p_ next_, int_4 _p_ ierr_);
 , int_4 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, int_4 _p_ next_, int_4 _p_ ierr_);
 d2_, real_8 _p_ a_, real_8 _p_ b_);
 , int_4 _p_ nknots_, real_8 _p_ t_, real_8 _p_ bcoef_, int_4 _p_ ndim_, int_4 _p_ kord_, int_4 _p_ ierr_);
 real_8 _p_ d_, int_4 _p_ incfd_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ismon_, int_4 _p_ ierr_);
 incfd_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, real_8 _p_ de_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, int_4 _p_ ierr_);
 skip_, real_8 _p_ a_, real_8 _p_ b_, int_4 _p_ ierr_);
 real_8 _p_ d_, int_4 _p_ incfd_, real_8 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ia_, int_4 _p_ ib_, int_4 _p_ ierr_);
 int_4 _p_ incfd_, real_8 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_);
 _p_ slope_, int_4 _p_ ierr_);
 int_4 _p_ nknots_, real_4 _p_ t_, real_4 _p_ bcoef_, int_4 _p_ ndim_, int_4 _p_ kord_, int_4 _p_ ierr_);
 real_4 _p_ d_, int_4 _p_ incfd_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ismon_, int_4 _p_ ierr_);
 incfd_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ne_, real_4 _p_ xe_, real_4 _p_ fe_, real_4 _p_ de_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ne_, real_4 _p_ xe_, real_4 _p_ fe_, int_4 _p_ ierr_);
 skip_, real_4 _p_ a_, real_4 _p_ b_, int_4 _p_ ierr_);
 real_4 _p_ d_, int_4 _p_ incfd_, real_4 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_);
 skip_, int_4 _p_ ia_, int_4 _p_ ib_, int_4 _p_ ierr_);
 _p_ incfd_, real_4 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_);
 slope_, int_4 _p_ ierr_);
 int_4 _chfcm (real_4 _p_ d1_, real_4 _p_ d2_, real_4 _p_ delta_)
 { // ** body not listed **
 }
 
 int_4 _chfdv (real_4 _p_ x1_, real_4 _p_ x2_, real_4 _p_ f1_, real_4 _p_ f2_, real_4 _p_ d1_, real_4 _p_ d2_, int_4 _p_ 
 ne_, real_4 _p_ xe_, real_4 _p_ fe_, real_4 _p_ de_, int_4 _p_ next_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _chfev (real_4 _p_ x1_, real_4 _p_ x2_, real_4 _p_ f1_, real_4 _p_ f2_, real_4 _p_ d1_, real_4 _p_ d2_, int_4 _p_ 
 ne_, real_4 _p_ xe_, real_4 _p_ fe_, int_4 _p_ next_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_4 _chfie (real_4 _p_ x1_, real_4 _p_ x2_, real_4 _p_ f1_, real_4 _p_ f2_, real_4 _p_ d1_, real_4 _p_ d2_, real_4 
 _p_ a_, real_4 _p_ b_)
 { // ** body not listed **
 }
 
 int_4 _dchfcm (real_8 _p_ d1_, real_8 _p_ d2_, real_8 _p_ delta_)
 { // ** body not listed **
 }
 
 int_4 _dchfdv (real_8 _p_ x1_, real_8 _p_ x2_, real_8 _p_ f1_, real_8 _p_ f2_, real_8 _p_ d1_, real_8 _p_ d2_, int_4 
 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, real_8 _p_ de_, int_4 _p_ next_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dchfev (real_8 _p_ x1_, real_8 _p_ x2_, real_8 _p_ f1_, real_8 _p_ f2_, real_8 _p_ d1_, real_8 _p_ d2_, int_4 
 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, int_4 _p_ next_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_8 _dchfie (real_8 _p_ x1_, real_8 _p_ x2_, real_8 _p_ f1_, real_8 _p_ f2_, real_8 _p_ d1_, real_8 _p_ d2_, real_8 
 _p_ a_, real_8 _p_ b_)
 { // ** body not listed **
 }
 
 int_4 _dpchbs (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, int_4 _p_ knotyp_, int_4 
 _p_ nknots_, real_8 _p_ t_, real_8 _p_ bcoef_, int_4 _p_ ndim_, int_4 _p_ kord_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchce (int_4 _p_ ic_, real_8 _p_ vc_, int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ h_, real_8 _p_ slope_, real_8 _p_ 
 d_, int_4 _p_ incfd_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchci (int_4 _p_ n_, real_8 _p_ h_, real_8 _p_ slope_, real_8 _p_ d_, int_4 _p_ incfd_)
 { // ** body not listed **
 }
 
 int_4 _dpchcm (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ismon_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchcs (real_8 _p_ switch_, int_4 _p_ n_, real_8 _p_ h_, real_8 _p_ slope_, real_8 _p_ d_, int_4 _p_ incfd_, 
 int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_8 _dpchdf (int_4 _p_ k_, real_8 _p_ x_, real_8 _p_ s_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchfd (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, real_8 _p_ de_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchfe (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ne_, real_8 _p_ xe_, real_8 _p_ fe_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_8 _dpchia (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, 
 real_8 _p_ a_, real_8 _p_ b_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchic (int_4 _p_ ic_, real_8 _p_ vc_, real_8 _p_ switch_, int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 
 _p_ d_, int_4 _p_ incfd_, real_8 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_8 _dpchid (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ia_, int_4 _p_ ib_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 int_4 _dpchim (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ incfd_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _dpchkt (int_4 _p_ n_, real_8 _p_ x_, int_4 _p_ knotyp_, real_8 _p_ t_)
 { // ** body not listed **
 }
 
 int_4 _dpchsp (int_4 _p_ ic_, real_8 _p_ vc_, int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ f_, real_8 _p_ d_, int_4 _p_ 
 incfd_, real_8 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_8 _dpchst (real_8 _p_ arg1_, real_8 _p_ arg2_)
 { // ** body not listed **
 }
 
 int_4 _dpchsw (real_8 _p_ dfmax_, int_4 _p_ iextrm_, real_8 _p_ d1_, real_8 _p_ d2_, real_8 _p_ h_, real_8 _p_ slope_, 
 int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchbs (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, int_4 _p_ knotyp_, int_4 _p_ 
 nknots_, real_4 _p_ t_, real_4 _p_ bcoef_, int_4 _p_ ndim_, int_4 _p_ kord_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchce (int_4 _p_ ic_, real_4 _p_ vc_, int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ h_, real_4 _p_ slope_, real_4 _p_ 
 d_, int_4 _p_ incfd_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchci (int_4 _p_ n_, real_4 _p_ h_, real_4 _p_ slope_, real_4 _p_ d_, int_4 _p_ incfd_)
 { // ** body not listed **
 }
 
 int_4 _pchcm (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ismon_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchcs (real_4 _p_ switch_, int_4 _p_ n_, real_4 _p_ h_, real_4 _p_ slope_, real_4 _p_ d_, int_4 _p_ incfd_, 
 int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_4 _pchdf (int_4 _p_ k_, real_4 _p_ x_, real_4 _p_ s_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchdoc (void)
 { // ** body not listed **
 }
 
 int_4 _pchfd (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ne_, real_4 _p_ xe_, real_4 _p_ fe_, real_4 _p_ de_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchfe (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ne_, real_4 _p_ xe_, real_4 _p_ fe_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_4 _pchia (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, real_4 
 _p_ a_, real_4 _p_ b_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchic (int_4 _p_ ic_, real_4 _p_ vc_, real_4 _p_ switch_, int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ 
 d_, int_4 _p_ incfd_, real_4 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_4 _pchid (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, logical_4 _p_ skip_, int_4 
 _p_ ia_, int_4 _p_ ib_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 int_4 _pchim (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ incfd_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 int_4 _pchkt (int_4 _p_ n_, real_4 _p_ x_, int_4 _p_ knotyp_, real_4 _p_ t_)
 { // ** body not listed **
 }
 
 int_4 _pchsp (int_4 _p_ ic_, real_4 _p_ vc_, int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ f_, real_4 _p_ d_, int_4 _p_ 
 incfd_, real_4 _p_ wk_, int_4 _p_ nwk_, int_4 _p_ ierr_)
 { // ** body not listed **
 }
 
 real_4 _pchst (real_4 _p_ arg1_, real_4 _p_ arg2_)
 { // ** body not listed **
 }
 
 int_4 _pchsw (real_4 _p_ dfmax_, int_4 _p_ iextrm_, real_4 _p_ d1_, real_4 _p_ d2_, real_4 _p_ h_, real_4 _p_ slope_, 
 int_4 _p_ ierr_)
 { // ** body not listed **
 }
 

Deze website is gearchiveerd door de KB, nationale bibliotheek. This website is archived by the National Library of the Netherlands.