PLaSK library
Loading...
Searching...
No Matches
rfft1i.c
Go to the documentation of this file.
1/* rfft1i.f -- translated by f2c (version 20100827).
2 You must link the resulting object file with libf2c:
3 on Microsoft Windows system, link with libf2c.lib;
4 on Linux or Unix systems, link with .../path/to/libf2c.a -lm
5 or, if you install libf2c.a in a standard place, with -lf2c -lm
6 -- in that order, at the end of the command line, as in
7 cc *.o -lf2c -lm
8 Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
9
10 http://www.netlib.org/f2c/libf2c.zip
11*/
12
13#include "f2c.h"
14
15/* Table of constant values */
16
17static integer c__3 = 3;
18
19/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
20/* * * */
21/* * copyright (c) 2011 by UCAR * */
22/* * * */
23/* * University Corporation for Atmospheric Research * */
24/* * * */
25/* * all rights reserved * */
26/* * * */
27/* * FFTPACK version 5.1 * */
28/* * * */
29/* * A Fortran Package of Fast Fourier * */
30/* * * */
31/* * Subroutines and Example Programs * */
32/* * * */
33/* * by * */
34/* * * */
35/* * Paul Swarztrauber and Dick Valent * */
36/* * * */
37/* * of * */
38/* * * */
39/* * the National Center for Atmospheric Research * */
40/* * * */
41/* * Boulder, Colorado (80307) U.S.A. * */
42/* * * */
43/* * which is sponsored by * */
44/* * * */
45/* * the National Science Foundation * */
46/* * * */
47/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
48
49/* Subroutine */ int rfft1i_(integer *n, doublereal *wsave, integer *lensav,
50 integer *ier)
51{
52 /* Builtin functions */
53 double log(doublereal);
54
55 /* Local variables */
56 extern /* Subroutine */ int rffti1_(integer *, doublereal *, doublereal *)
57 , xerfft_(char *, integer *, ftnlen);
58
59
60 /* Parameter adjustments */
61 --wsave;
62
63 /* Function Body */
64 *ier = 0;
65
66 if (*lensav < *n + (integer) (log((doublereal) (*n)) / log(2.)) + 4) {
67 *ier = 2;
68 xerfft_("RFFT1I ", &c__3, (ftnlen)7);
69 }
70
71 if (*n == 1) {
72 return 0;
73 }
74
75 rffti1_(n, &wsave[1], &wsave[*n + 1]);
76 return 0;
77} /* rfft1i_ */
78