summaryrefslogtreecommitdiff
path: root/src/fftw3/dft/vrank-geq1.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/fftw3/dft/vrank-geq1.c')
-rw-r--r--src/fftw3/dft/vrank-geq1.c215
1 files changed, 0 insertions, 215 deletions
diff --git a/src/fftw3/dft/vrank-geq1.c b/src/fftw3/dft/vrank-geq1.c
deleted file mode 100644
index 92a5776..0000000
--- a/src/fftw3/dft/vrank-geq1.c
+++ /dev/null
@@ -1,215 +0,0 @@
-/*
- * Copyright (c) 2003 Matteo Frigo
- * Copyright (c) 2003 Massachusetts Institute of Technology
- *
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
- *
- */
-
-/* $Id: vrank-geq1.c,v 1.1 2008/10/17 06:11:08 scuri Exp $ */
-
-
-/* Plans for handling vector transform loops. These are *just* the
- loops, and rely on child plans for the actual DFTs.
-
- They form a wrapper around solvers that don't have apply functions
- for non-null vectors.
-
- vrank-geq1 plans also recursively handle the case of multi-dimensional
- vectors, obviating the need for most solvers to deal with this. We
- can also play games here, such as reordering the vector loops.
-
- Each vrank-geq1 plan reduces the vector rank by 1, picking out a
- dimension determined by the vecloop_dim field of the solver. */
-
-#include "dft.h"
-
-typedef struct {
- solver super;
- int vecloop_dim;
- const int *buddies;
- int nbuddies;
-} S;
-
-typedef struct {
- plan_dft super;
-
- plan *cld;
- int vl;
- int ivs, ovs;
- const S *solver;
-} P;
-
-static void apply(const plan *ego_, R *ri, R *ii, R *ro, R *io)
-{
- const P *ego = (const P *) ego_;
- int i, vl = ego->vl;
- int ivs = ego->ivs, ovs = ego->ovs;
- dftapply cldapply = ((plan_dft *) ego->cld)->apply;
-
- for (i = 0; i < vl; ++i) {
- cldapply(ego->cld,
- ri + i * ivs, ii + i * ivs, ro + i * ovs, io + i * ovs);
- }
-}
-
-static void awake(plan *ego_, int flg)
-{
- P *ego = (P *) ego_;
- AWAKE(ego->cld, flg);
-}
-
-static void destroy(plan *ego_)
-{
- P *ego = (P *) ego_;
- X(plan_destroy_internal)(ego->cld);
-}
-
-static void print(const plan *ego_, printer *p)
-{
- const P *ego = (const P *) ego_;
- const S *s = ego->solver;
- p->print(p, "(dft-vrank>=1-x%d/%d%(%p%))",
- ego->vl, s->vecloop_dim, ego->cld);
-}
-
-static int pickdim(const S *ego, const tensor *vecsz, int oop, int *dp)
-{
- return X(pickdim)(ego->vecloop_dim, ego->buddies, ego->nbuddies,
- vecsz, oop, dp);
-}
-
-static int applicable0(const solver *ego_, const problem *p_, int *dp)
-{
- if (DFTP(p_)) {
- const S *ego = (const S *) ego_;
- const problem_dft *p = (const problem_dft *) p_;
-
- return (1
- && FINITE_RNK(p->vecsz->rnk)
- && p->vecsz->rnk > 0
- && pickdim(ego, p->vecsz, p->ri != p->ro, dp)
- );
- }
-
- return 0;
-}
-
-static int applicable(const solver *ego_, const problem *p_,
- const planner *plnr, int *dp)
-{
- const S *ego = (const S *)ego_;
- const problem_dft *p;
-
- if (!applicable0(ego_, p_, dp)) return 0;
-
- /* fftw2 behavior */
- if (NO_VRANK_SPLITSP(plnr) && (ego->vecloop_dim != ego->buddies[0]))
- return 0;
-
- p = (const problem_dft *) p_;
-
- if (NO_UGLYP(plnr)) {
- /* Heuristic: if the transform is multi-dimensional, and the
- vector stride is less than the transform size, then we
- probably want to use a rank>=2 plan first in order to combine
- this vector with the transform-dimension vectors. */
- {
- iodim *d = p->vecsz->dims + *dp;
- if (1
- && p->sz->rnk > 1
- && X(imin)(X(iabs)(d->is), X(iabs)(d->os))
- < X(tensor_max_index)(p->sz)
- )
- return 0;
- }
-
- /* Heuristic: don't use a vrank-geq1 for rank-0 vrank-1
- transforms, since this case is better handled by rank-0
- solvers. */
- if (p->sz->rnk == 0 && p->vecsz->rnk == 1) return 0;
-
- if (NONTHREADED_ICKYP(plnr)) return 0; /* prefer threaded version */
- }
-
- return 1;
-}
-
-static plan *mkplan(const solver *ego_, const problem *p_, planner *plnr)
-{
- const S *ego = (const S *) ego_;
- const problem_dft *p;
- P *pln;
- plan *cld;
- int vdim;
- iodim *d;
-
- static const plan_adt padt = {
- X(dft_solve), awake, print, destroy
- };
-
- if (!applicable(ego_, p_, plnr, &vdim))
- return (plan *) 0;
- p = (const problem_dft *) p_;
-
- d = p->vecsz->dims + vdim;
-
- A(d->n > 1);
- cld = X(mkplan_d)(plnr,
- X(mkproblem_dft_d)(
- X(tensor_copy)(p->sz),
- X(tensor_copy_except)(p->vecsz, vdim),
- TAINT(p->ri, d->is), TAINT(p->ii, d->is),
- TAINT(p->ro, d->os), TAINT(p->io, d->os)));
- if (!cld) return (plan *) 0;
-
- pln = MKPLAN_DFT(P, &padt, apply);
-
- pln->cld = cld;
- pln->vl = d->n;
- pln->ivs = d->is;
- pln->ovs = d->os;
-
- pln->solver = ego;
- X(ops_zero)(&pln->super.super.ops);
- pln->super.super.ops.other = 3.14159; /* magic to prefer codelet loops */
- X(ops_madd2)(pln->vl, &cld->ops, &pln->super.super.ops);
- pln->super.super.pcost = pln->vl * cld->pcost;
-
- return &(pln->super.super);
-}
-
-static solver *mksolver(int vecloop_dim, const int *buddies, int nbuddies)
-{
- static const solver_adt sadt = { mkplan };
- S *slv = MKSOLVER(S, &sadt);
- slv->vecloop_dim = vecloop_dim;
- slv->buddies = buddies;
- slv->nbuddies = nbuddies;
- return &(slv->super);
-}
-
-void X(dft_vrank_geq1_register)(planner *p)
-{
- int i;
-
- /* FIXME: Should we try other vecloop_dim values? */
- static const int buddies[] = { 1, -1 };
-
- const int nbuddies = sizeof(buddies) / sizeof(buddies[0]);
-
- for (i = 0; i < nbuddies; ++i)
- REGISTER_SOLVER(p, mksolver(buddies[i], buddies, nbuddies));
-}
generated by cgit v1.2.3 (git 2.39.1) at 2024-11-22 03:31:53 +0000