☰  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-fftpack4.c

 //! @file mathlib-fftpack4.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
 //
 //! 68 subprograms from FFTPACK4
 
 // Compiled from Fortran source code by VIF.
 // Selected subprograms are:
 //
 // CFFTB  COSQB1 DCOSQF DSINQB DZFFTI PASSB3 PASSF5 RADF3  RFFTF1 SINT
 // CFFTB1 COSQF  DCOSQI DSINQF EZFFT1 PASSB4 RADB2  RADF4  RFFTI  SINT1
 // CFFTF  COSQF1 DCOST  DSINQI EZFFTB PASSB5 RADB3  RADF5  RFFTI1 SINTI
 // CFFTF1 COSQI  DCOSTI DSINT  EZFFTF PASSF  RADB4  RADFG  SINQB  ZFFTB
 // CFFTI  COST   DFFTB  DSINTI EZFFTI PASSF2 RADB5  RFFTB  SINQF  ZFFTF
 // CFFTI1 COSTI  DFFTF  DZFFTB PASSB  PASSF3 RADBG  RFFTB1 SINQI  ZFFTI
 // COSQB  DCOSQB DFFTI  DZFFTF PASSB2 PASSF4 RADF2  RFFTF
 
 // Source: netlib.org/fftpack
 // 
 // DFFTPACK is a Double precision clone by Hugh C. Pumphrey of 
 // FFTPACK version 4 april 1985.
 // 
 // The original FFTPACK was public domain, so dfftpack is public domain
 // too. It is released in the hope it will be useful to someone. There is 
 // no warranty of any sort covering this software.
 
 // The license for FFTPACK4 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] = {
   {"cfftb1", 0}, // subroutine
   {"cfftb", 0}, // subroutine
   {"cfftf1", 0}, // subroutine
   {"cfftf", 0}, // subroutine
   {"cffti1", 0}, // subroutine
   {"cffti", 0}, // subroutine
   {"cosqb1", 0}, // subroutine
   {"cosqb", 0}, // subroutine
   {"cosqf1", 0}, // subroutine
   {"cosqf", 0}, // subroutine
   {"cosqi", 0}, // subroutine
   {"cost", 0}, // subroutine
   {"costi", 0}, // subroutine
   {"dcosqb", 0}, // subroutine
   {"dcosqf", 0}, // subroutine
   {"dcosqi", 0}, // subroutine
   {"dcost", 0}, // subroutine
   {"dcosti", 0}, // subroutine
   {"dfftb", 0}, // subroutine
   {"dfftf", 0}, // subroutine
   {"dffti", 0}, // subroutine
   {"dsinqb", 0}, // subroutine
   {"dsinqf", 0}, // subroutine
   {"dsinqi", 0}, // subroutine
   {"dsint", 0}, // subroutine
   {"dsinti", 0}, // subroutine
   {"dzfftb", 0}, // subroutine
   {"dzfftf", 0}, // subroutine
   {"dzffti", 0}, // subroutine
   {"ezfft1", 0}, // subroutine
   {"ezfftb", 0}, // subroutine
   {"ezfftf", 0}, // subroutine
   {"ezffti", 0}, // subroutine
   {"passb2", 0}, // subroutine
   {"passb3", 0}, // subroutine
   {"passb4", 0}, // subroutine
   {"passb5", 0}, // subroutine
   {"passb", 0}, // subroutine
   {"passf2", 0}, // subroutine
   {"passf3", 0}, // subroutine
   {"passf4", 0}, // subroutine
   {"passf5", 0}, // subroutine
   {"passf", 0}, // subroutine
   {"radb2", 0}, // subroutine
   {"radb3", 0}, // subroutine
   {"radb4", 0}, // subroutine
   {"radb5", 0}, // subroutine
   {"radbg", 0}, // subroutine
   {"radf2", 0}, // subroutine
   {"radf3", 0}, // subroutine
   {"radf4", 0}, // subroutine
   {"radf5", 0}, // subroutine
   {"radfg", 0}, // subroutine
   {"rfftb1", 0}, // subroutine
   {"rfftb", 0}, // subroutine
   {"rfftf1", 0}, // subroutine
   {"rfftf", 0}, // subroutine
   {"rffti1", 0}, // subroutine
   {"rffti", 0}, // subroutine
   {"sinqb", 0}, // subroutine
   {"sinqf", 0}, // subroutine
   {"sinqi", 0}, // subroutine
   {"sint1", 0}, // subroutine
   {"sint", 0}, // subroutine
   {"sinti", 0}, // subroutine
   {"zfftb", 0}, // subroutine
   {"zfftf", 0}, // subroutine
   {"zffti", 0}, // subroutine
   {NULL, 0}
 };
 
 wsave_);
 wsave_);
 wsave_);
 wsave_);
 wa2_);
 wa2_, real_8 _p_ wa3_);
 wa2_, real_8 _p_ wa3_, real_8 _p_ wa4_);
 real_8 _p_ c1_, real_8 _p_ c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_);
 wa2_);
 wa2_, real_8 _p_ wa3_);
 wa2_, real_8 _p_ wa3_, real_8 _p_ wa4_);
 real_8 _p_ c1_, real_8 _p_ c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_);
 );
 , real_8 _p_ wa3_);
 , real_8 _p_ wa3_, real_8 _p_ wa4_);
 real_8 _p_ c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_);
 );
 , real_8 _p_ wa3_);
 , real_8 _p_ wa3_, real_8 _p_ wa4_);
 real_8 _p_ c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_);
 ;
 int_4 _cfftb1 (int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ ch_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _cfftb (int_4 _p_ n_, complex_8 _p_ c_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _cfftf1 (int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ ch_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _cfftf (int_4 _p_ n_, complex_8 _p_ c_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _cffti1 (int_4 _p_ n_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _cffti (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _cosqb1 (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ w_, real_8 _p_ xh_)
 { // ** body not listed **
 }
 
 int_4 _cosqb (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _cosqf1 (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ w_, real_8 _p_ xh_)
 { // ** body not listed **
 }
 
 int_4 _cosqf (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _cosqi (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _cost (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _costi (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dcosqb (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dcosqf (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dcosqi (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dcost (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dcosti (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dfftb (int_4 _p_ n_, real_8 _p_ r_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dfftf (int_4 _p_ n_, real_8 _p_ r_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dffti (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dsinqb (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dsinqf (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dsinqi (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dsint (int_4 _p_ n_, real_8 _p_ x_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dsinti (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dzfftb (int_4 _p_ n_, real_8 _p_ r_, real_8 _p_ azero_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dzfftf (int_4 _p_ n_, real_8 _p_ r_, real_8 _p_ azero_, real_8 _p_ a_, real_8 _p_ b_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _dzffti (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _ezfft1 (int_4 _p_ n_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _ezfftb (int_4 _p_ n_, real_4 _p_ r_, real_4 _p_ azero_, real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _ezfftf (int_4 _p_ n_, real_4 _p_ r_, real_4 _p_ azero_, real_4 _p_ a_, real_4 _p_ b_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _ezffti (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _passb2 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_)
 { // ** body not listed **
 }
 
 int_4 _passb3 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_)
 { // ** body not listed **
 }
 
 int_4 _passb4 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_)
 { // ** body not listed **
 }
 
 int_4 _passb5 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_, real_8 _p_ wa4_)
 { // ** body not listed **
 }
 
 int_4 _passb (int_4 _p_ nac_, int_4 _p_ ido_, int_4 _p_ ip_, int_4 _p_ l1_, int_4 _p_ idl1_, real_8 _p_ cc_, real_8 _p_ 
 c1_, real_8 _p_ c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_)
 { // ** body not listed **
 }
 
 int_4 _passf2 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_)
 { // ** body not listed **
 }
 
 int_4 _passf3 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_)
 { // ** body not listed **
 }
 
 int_4 _passf4 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_)
 { // ** body not listed **
 }
 
 int_4 _passf5 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_, real_8 _p_ wa4_)
 { // ** body not listed **
 }
 
 int_4 _passf (int_4 _p_ nac_, int_4 _p_ ido_, int_4 _p_ ip_, int_4 _p_ l1_, int_4 _p_ idl1_, real_8 _p_ cc_, real_8 _p_ 
 c1_, real_8 _p_ c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_)
 { // ** body not listed **
 }
 
 int_4 _radb2 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_)
 { // ** body not listed **
 }
 int_4 _radb3 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_)
 { // ** body not listed **
 }
 
 int_4 _radb4 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_)
 { // ** body not listed **
 }
 
 int_4 _radb5 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_, real_8 _p_ wa4_)
 { // ** body not listed **
 }
 
 int_4 _radbg (int_4 _p_ ido_, int_4 _p_ ip_, int_4 _p_ l1_, int_4 _p_ idl1_, real_8 _p_ cc_, real_8 _p_ c1_, real_8 _p_ 
 c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_)
 { // ** body not listed **
 }
 
 int_4 _radf2 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_)
 { // ** body not listed **
 }
 int_4 _radf3 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_)
 { // ** body not listed **
 }
 
 int_4 _radf4 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_)
 { // ** body not listed **
 }
 
 int_4 _radf5 (int_4 _p_ ido_, int_4 _p_ l1_, real_8 _p_ cc_, real_8 _p_ ch_, real_8 _p_ wa1_, real_8 _p_ wa2_, real_8 
 _p_ wa3_, real_8 _p_ wa4_)
 { // ** body not listed **
 }
 
 int_4 _radfg (int_4 _p_ ido_, int_4 _p_ ip_, int_4 _p_ l1_, int_4 _p_ idl1_, real_8 _p_ cc_, real_8 _p_ c1_, real_8 _p_ 
 c2_, real_8 _p_ ch_, real_8 _p_ ch2_, real_8 _p_ wa_)
 { // ** body not listed **
 }
 
 int_4 _rfftb1 (int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ ch_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _rfftb (int_4 _p_ n_, real_4 _p_ r_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _rfftf1 (int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ ch_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _rfftf (int_4 _p_ n_, real_4 _p_ r_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _rffti1 (int_4 _p_ n_, real_8 _p_ wa_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _rffti (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _sinqb (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _sinqf (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _sinqi (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _sint1 (int_4 _p_ n_, real_8 _p_ war_, real_8 _p_ was_, real_8 _p_ xh_, real_8 _p_ x_, int_4 _p_ ifac_)
 { // ** body not listed **
 }
 
 int_4 _sint (int_4 _p_ n_, real_4 _p_ x_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _sinti (int_4 _p_ n_, real_4 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _zfftb (int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _zfftf (int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 
 int_4 _zffti (int_4 _p_ n_, real_8 _p_ wsave_)
 { // ** body not listed **
 }
 

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