Chris@19: /* Chris@19: * Copyright (c) 2003, 2007-14 Matteo Frigo Chris@19: * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology Chris@19: * Chris@19: * This program is free software; you can redistribute it and/or modify Chris@19: * it under the terms of the GNU General Public License as published by Chris@19: * the Free Software Foundation; either version 2 of the License, or Chris@19: * (at your option) any later version. Chris@19: * Chris@19: * This program is distributed in the hope that it will be useful, Chris@19: * but WITHOUT ANY WARRANTY; without even the implied warranty of Chris@19: * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Chris@19: * GNU General Public License for more details. Chris@19: * Chris@19: * You should have received a copy of the GNU General Public License Chris@19: * along with this program; if not, write to the Free Software Chris@19: * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA Chris@19: * Chris@19: */ Chris@19: Chris@19: #include "dft.h" Chris@19: Chris@19: typedef void (*dftwapply)(const plan *ego, R *rio, R *iio); Chris@19: typedef struct ct_solver_s ct_solver; Chris@19: typedef plan *(*ct_mkinferior)(const ct_solver *ego, Chris@19: INT r, INT irs, INT ors, Chris@19: INT m, INT ms, Chris@19: INT v, INT ivs, INT ovs, Chris@19: INT mstart, INT mcount, Chris@19: R *rio, R *iio, planner *plnr); Chris@19: typedef int (*ct_force_vrecursion)(const ct_solver *ego, Chris@19: const problem_dft *p); Chris@19: Chris@19: typedef struct { Chris@19: plan super; Chris@19: dftwapply apply; Chris@19: } plan_dftw; Chris@19: Chris@19: extern plan *X(mkplan_dftw)(size_t size, const plan_adt *adt, dftwapply apply); Chris@19: Chris@19: #define MKPLAN_DFTW(type, adt, apply) \ Chris@19: (type *)X(mkplan_dftw)(sizeof(type), adt, apply) Chris@19: Chris@19: struct ct_solver_s { Chris@19: solver super; Chris@19: INT r; Chris@19: int dec; Chris@19: # define DECDIF 0 Chris@19: # define DECDIT 1 Chris@19: # define TRANSPOSE 2 Chris@19: ct_mkinferior mkcldw; Chris@19: ct_force_vrecursion force_vrecursionp; Chris@19: }; Chris@19: Chris@19: int X(ct_applicable)(const ct_solver *, const problem *, planner *); Chris@19: ct_solver *X(mksolver_ct)(size_t size, INT r, int dec, Chris@19: ct_mkinferior mkcldw, Chris@19: ct_force_vrecursion force_vrecursionp); Chris@19: extern ct_solver *(*X(mksolver_ct_hook))(size_t, INT, int, Chris@19: ct_mkinferior, ct_force_vrecursion); Chris@19: Chris@19: void X(regsolver_ct_directw)(planner *plnr, Chris@19: kdftw codelet, const ct_desc *desc, int dec); Chris@19: void X(regsolver_ct_directwbuf)(planner *plnr, Chris@19: kdftw codelet, const ct_desc *desc, int dec); Chris@19: solver *X(mksolver_ctsq)(kdftwsq codelet, const ct_desc *desc, int dec); Chris@19: void X(regsolver_ct_directwsq)(planner *plnr, kdftwsq codelet, Chris@19: const ct_desc *desc, int dec);