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 "api.h" Chris@19: #include "rdft.h" Chris@19: Chris@19: #define N0(nembed)((nembed) ? (nembed) : n) Chris@19: Chris@19: X(plan) X(plan_many_r2r)(int rank, const int *n, Chris@19: int howmany, Chris@19: R *in, const int *inembed, Chris@19: int istride, int idist, Chris@19: R *out, const int *onembed, Chris@19: int ostride, int odist, Chris@19: const X(r2r_kind) * kind, unsigned flags) Chris@19: { Chris@19: X(plan) p; Chris@19: rdft_kind *k; Chris@19: Chris@19: if (!X(many_kosherp)(rank, n, howmany)) return 0; Chris@19: Chris@19: k = X(map_r2r_kind)(rank, kind); Chris@19: p = X(mkapiplan)( Chris@19: 0, flags, Chris@19: X(mkproblem_rdft_d)(X(mktensor_rowmajor)(rank, n, Chris@19: N0(inembed), N0(onembed), Chris@19: istride, ostride), Chris@19: X(mktensor_1d)(howmany, idist, odist), Chris@19: TAINT_UNALIGNED(in, flags), Chris@19: TAINT_UNALIGNED(out, flags), k)); Chris@19: X(ifree0)(k); Chris@19: return p; Chris@19: }