diff options
Diffstat (limited to 'src/fftw3/rdft/rdft2-tensor-max-index.c')
-rw-r--r-- | src/fftw3/rdft/rdft2-tensor-max-index.c | 44 |
1 files changed, 0 insertions, 44 deletions
diff --git a/src/fftw3/rdft/rdft2-tensor-max-index.c b/src/fftw3/rdft/rdft2-tensor-max-index.c deleted file mode 100644 index 0d29af3..0000000 --- a/src/fftw3/rdft/rdft2-tensor-max-index.c +++ /dev/null @@ -1,44 +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: rdft2-tensor-max-index.c,v 1.1 2008/10/17 06:11:29 scuri Exp $ */ - -#include "rdft.h" - -/* like X(tensor_max_index), but takes into account the special n/2+1 - final dimension for the complex output/input of an R2HC/HC2R transform. */ -int X(rdft2_tensor_max_index)(const tensor *sz, rdft_kind k) -{ - int i; - int n = 0; - - A(FINITE_RNK(sz->rnk)); - for (i = 0; i + 1 < sz->rnk; ++i) { - const iodim *p = sz->dims + i; - n += (p->n - 1) * X(imax)(X(iabs)(p->is), X(iabs)(p->os)); - } - if (i < sz->rnk) { - const iodim *p = sz->dims + i; - int is, os; - X(rdft2_strides)(k, p, &is, &os); - n += X(imax)((p->n - 1) * X(iabs)(is), (p->n/2) * X(iabs)(os)); - } - return n; -} |