PLaSK library
Loading...
Searching...
No Matches
cfftmb.c
Go to the documentation of this file.
1/* cfftmb.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__6 = 6;
18static integer c__8 = 8;
19static integer c__10 = 10;
20static integer c_n1 = -1;
21
22/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
23/* * * */
24/* * copyright (c) 2011 by UCAR * */
25/* * * */
26/* * University Corporation for Atmospheric Research * */
27/* * * */
28/* * all rights reserved * */
29/* * * */
30/* * FFTPACK version 5.1 * */
31/* * * */
32/* * A Fortran Package of Fast Fourier * */
33/* * * */
34/* * Subroutines and Example Programs * */
35/* * * */
36/* * by * */
37/* * * */
38/* * Paul Swarztrauber and Dick Valent * */
39/* * * */
40/* * of * */
41/* * * */
42/* * the National Center for Atmospheric Research * */
43/* * * */
44/* * Boulder, Colorado (80307) U.S.A. * */
45/* * * */
46/* * which is sponsored by * */
47/* * * */
48/* * the National Science Foundation * */
49/* * * */
50/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
51
52/* Subroutine */ int cfftmb_(integer *lot, integer *jump, integer *n, integer
53 *inc, doublecomplex *c__, integer *lenc, doublereal *wsave, integer *
54 lensav, doublereal *work, integer *lenwrk, integer *ier)
55{
56 /* Builtin functions */
57 double log(doublereal);
58
59 /* Local variables */
60 integer iw1;
61 extern /* Subroutine */ int cmfm1b_(integer *, integer *, integer *,
62 integer *, doublecomplex *, doublereal *, doublereal *,
63 doublereal *, doublereal *);
64 extern logical xercon_(integer *, integer *, integer *, integer *);
65 extern /* Subroutine */ int xerfft_(char *, integer *, ftnlen);
66
67
68 /* Parameter adjustments */
69 --c__;
70 --wsave;
71 --work;
72
73 /* Function Body */
74 *ier = 0;
75
76 if (*lenc < (*lot - 1) * *jump + *inc * (*n - 1) + 1) {
77 *ier = 1;
78 xerfft_("CFFTMB ", &c__6, (ftnlen)7);
79 } else if (*lensav < (*n << 1) + (integer) (log((doublereal) (*n)) / log(
80 2.)) + 4) {
81 *ier = 2;
82 xerfft_("CFFTMB ", &c__8, (ftnlen)7);
83 } else if (*lenwrk < (*lot << 1) * *n) {
84 *ier = 3;
85 xerfft_("CFFTMB ", &c__10, (ftnlen)7);
86 } else if (! xercon_(inc, jump, n, lot)) {
87 *ier = 4;
88 xerfft_("CFFTMB ", &c_n1, (ftnlen)7);
89 }
90
91 if (*n == 1) {
92 return 0;
93 }
94
95 iw1 = *n + *n + 1;
96 cmfm1b_(lot, jump, n, inc, &c__[1], &work[1], &wsave[1], &wsave[iw1], &
97 wsave[iw1 + 1]);
98 return 0;
99} /* cfftmb_ */
100