annotate fft/fftw/fftw-3.3.4/api/execute-dft-c2r.c @ 40:223f770b5341 kissfft-double tip

Try a double-precision kissfft
author Chris Cannam
date Wed, 07 Sep 2016 10:40:32 +0100
parents 26056e866c29
children
rev   line source
Chris@19 1 /*
Chris@19 2 * Copyright (c) 2003, 2007-14 Matteo Frigo
Chris@19 3 * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
Chris@19 4 *
Chris@19 5 * This program is free software; you can redistribute it and/or modify
Chris@19 6 * it under the terms of the GNU General Public License as published by
Chris@19 7 * the Free Software Foundation; either version 2 of the License, or
Chris@19 8 * (at your option) any later version.
Chris@19 9 *
Chris@19 10 * This program is distributed in the hope that it will be useful,
Chris@19 11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
Chris@19 12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
Chris@19 13 * GNU General Public License for more details.
Chris@19 14 *
Chris@19 15 * You should have received a copy of the GNU General Public License
Chris@19 16 * along with this program; if not, write to the Free Software
Chris@19 17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
Chris@19 18 *
Chris@19 19 */
Chris@19 20
Chris@19 21 #include "api.h"
Chris@19 22 #include "rdft.h"
Chris@19 23
Chris@19 24 /* guru interface: requires care in alignment, r - i, etcetera. */
Chris@19 25 void X(execute_dft_c2r)(const X(plan) p, C *in, R *out)
Chris@19 26 {
Chris@19 27 plan_rdft2 *pln = (plan_rdft2 *) p->pln;
Chris@19 28 problem_rdft2 *prb = (problem_rdft2 *) p->prb;
Chris@19 29 pln->apply((plan *) pln, out, out + (prb->r1 - prb->r0), in[0], in[0]+1);
Chris@19 30 }