PLaSK library
Loading...
Searching...
No Matches
mcfti1.c
Go to the documentation of this file.
1/* mcfti1.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/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
16/* * * */
17/* * copyright (c) 2011 by UCAR * */
18/* * * */
19/* * University Corporation for Atmospheric Research * */
20/* * * */
21/* * all rights reserved * */
22/* * * */
23/* * FFTPACK version 5.1 * */
24/* * * */
25/* * A Fortran Package of Fast Fourier * */
26/* * * */
27/* * Subroutines and Example Programs * */
28/* * * */
29/* * by * */
30/* * * */
31/* * Paul Swarztrauber and Dick Valent * */
32/* * * */
33/* * of * */
34/* * * */
35/* * the National Center for Atmospheric Research * */
36/* * * */
37/* * Boulder, Colorado (80307) U.S.A. * */
38/* * * */
39/* * which is sponsored by * */
40/* * * */
41/* * the National Science Foundation * */
42/* * * */
43/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
44
45/* Subroutine */ int mcfti1_(integer *n, doublereal *wa, doublereal *fnf,
46 doublereal *fac)
47{
48 /* System generated locals */
49 integer i__1;
50
51 /* Local variables */
52 integer k1, l1, l2, nf, ip, iw, ido;
53 extern /* Subroutine */ int tables_(integer *, integer *, doublereal *),
54 factor_(integer *, integer *, doublereal *);
55
56
57 /* Parameter adjustments */
58 --fac;
59 --wa;
60
61 /* Function Body */
62 factor_(n, &nf, &fac[1]);
63 *fnf = (doublereal) nf;
64 iw = 1;
65 l1 = 1;
66 i__1 = nf;
67 for (k1 = 1; k1 <= i__1; ++k1) {
68 ip = (integer) fac[k1];
69 l2 = l1 * ip;
70 ido = *n / l2;
71 tables_(&ido, &ip, &wa[iw]);
72 iw += (ip - 1) * (ido + ido);
73 l1 = l2;
74/* L110: */
75 }
76 return 0;
77} /* mcfti1_ */
78