annotate fft/fftw/fftw-3.3.4/kernel/plan.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
Chris@19 22 #include "ifftw.h"
Chris@19 23
Chris@19 24 /* "Plan: To bother about the best method of accomplishing an
Chris@19 25 accidental result." (Ambrose Bierce, The Enlarged Devil's
Chris@19 26 Dictionary). */
Chris@19 27
Chris@19 28 plan *X(mkplan)(size_t size, const plan_adt *adt)
Chris@19 29 {
Chris@19 30 plan *p = (plan *)MALLOC(size, PLANS);
Chris@19 31
Chris@19 32 A(adt->destroy);
Chris@19 33 p->adt = adt;
Chris@19 34 X(ops_zero)(&p->ops);
Chris@19 35 p->pcost = 0.0;
Chris@19 36 p->wakefulness = SLEEPY;
Chris@19 37 p->could_prune_now_p = 0;
Chris@19 38
Chris@19 39 return p;
Chris@19 40 }
Chris@19 41
Chris@19 42 /*
Chris@19 43 * destroy a plan
Chris@19 44 */
Chris@19 45 void X(plan_destroy_internal)(plan *ego)
Chris@19 46 {
Chris@19 47 if (ego) {
Chris@19 48 A(ego->wakefulness == SLEEPY);
Chris@19 49 ego->adt->destroy(ego);
Chris@19 50 X(ifree)(ego);
Chris@19 51 }
Chris@19 52 }
Chris@19 53
Chris@19 54 /* dummy destroy routine for plans with no local state */
Chris@19 55 void X(plan_null_destroy)(plan *ego)
Chris@19 56 {
Chris@19 57 UNUSED(ego);
Chris@19 58 /* nothing */
Chris@19 59 }
Chris@19 60
Chris@19 61 void X(plan_awake)(plan *ego, enum wakefulness wakefulness)
Chris@19 62 {
Chris@19 63 if (ego) {
Chris@19 64 A(((wakefulness == SLEEPY) ^ (ego->wakefulness == SLEEPY)));
Chris@19 65
Chris@19 66 ego->adt->awake(ego, wakefulness);
Chris@19 67 ego->wakefulness = wakefulness;
Chris@19 68 }
Chris@19 69 }
Chris@19 70