Chris@42: /* Chris@42: * Copyright (c) 2001 Matteo Frigo Chris@42: * Copyright (c) 2001 Massachusetts Institute of Technology Chris@42: * Chris@42: * This program is free software; you can redistribute it and/or modify Chris@42: * it under the terms of the GNU General Public License as published by Chris@42: * the Free Software Foundation; either version 2 of the License, or Chris@42: * (at your option) any later version. Chris@42: * Chris@42: * This program is distributed in the hope that it will be useful, Chris@42: * but WITHOUT ANY WARRANTY; without even the implied warranty of Chris@42: * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the Chris@42: * GNU General Public License for more details. Chris@42: * Chris@42: * You should have received a copy of the GNU General Public License Chris@42: * along with this program; if not, write to the Free Software Chris@42: * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA Chris@42: * Chris@42: */ Chris@42: Chris@42: Chris@42: #include "bench.h" Chris@42: Chris@42: /* set I/O arrays to zero. Default routine */ Chris@42: void problem_zero(bench_problem *p) Chris@42: { Chris@42: bench_complex czero = {0, 0}; Chris@42: if (p->kind == PROBLEM_COMPLEX) { Chris@42: caset((bench_complex *) p->inphys, p->iphyssz, czero); Chris@42: caset((bench_complex *) p->outphys, p->ophyssz, czero); Chris@42: } else if (p->kind == PROBLEM_R2R) { Chris@42: aset((bench_real *) p->inphys, p->iphyssz, 0.0); Chris@42: aset((bench_real *) p->outphys, p->ophyssz, 0.0); Chris@42: } else if (p->kind == PROBLEM_REAL && p->sign < 0) { Chris@42: aset((bench_real *) p->inphys, p->iphyssz, 0.0); Chris@42: caset((bench_complex *) p->outphys, p->ophyssz, czero); Chris@42: } else if (p->kind == PROBLEM_REAL && p->sign > 0) { Chris@42: caset((bench_complex *) p->inphys, p->iphyssz, czero); Chris@42: aset((bench_real *) p->outphys, p->ophyssz, 0.0); Chris@42: } else { Chris@42: BENCH_ASSERT(0); /* TODO */ Chris@42: } Chris@42: }