Mercurial > hg > sv-dependency-builds

annotate src/fftw-3.3.8/libbench2/allocate.c @ 83:ae30d91d2ffe

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Replace these with versions built using an older toolset (so as to avoid ABI compatibilities when linking on Ubuntu 14.04 for packaging purposes)
author Chris Cannam
date Fri, 07 Feb 2020 11:51:13 +0000
parents d0c2a83c1364
children
rev   line source
Chris@82 1 /* not worth copyrighting */
Chris@82 2
Chris@82 3
Chris@82 4 #include "libbench2/bench.h"
Chris@82 5
Chris@82 6 static void bounds(bench_problem *p, int *ilb, int *iub, int *olb, int *oub)
Chris@82 7 {
Chris@82 8 bench_tensor *t = tensor_append(p->sz, p->vecsz);
Chris@82 9 tensor_ibounds(t, ilb, iub);
Chris@82 10 tensor_obounds(t, olb, oub);
Chris@82 11 tensor_destroy(t);
Chris@82 12 }
Chris@82 13
Chris@82 14 /*
Chris@82 15 * Allocate I/O arrays for a problem.
Chris@82 16 *
Chris@82 17 * This is the default routine that can be overridden by the user in
Chris@82 18 * complicated cases.
Chris@82 19 */
Chris@82 20 void problem_alloc(bench_problem *p)
Chris@82 21 {
Chris@82 22 int ilb, iub, olb, oub;
Chris@82 23 int isz, osz;
Chris@82 24
Chris@82 25 bounds(p, &ilb, &iub, &olb, &oub);
Chris@82 26 isz = iub - ilb;
Chris@82 27 osz = oub - olb;
Chris@82 28
Chris@82 29 if (p->kind == PROBLEM_COMPLEX) {
Chris@82 30 bench_complex *in, *out;
Chris@82 31
Chris@82 32 p->iphyssz = isz;
Chris@82 33 p->inphys = in = (bench_complex *) bench_malloc(isz * sizeof(bench_complex));
Chris@82 34 p->in = in - ilb;
Chris@82 35
Chris@82 36 if (p->in_place) {
Chris@82 37 p->out = p->in;
Chris@82 38 p->outphys = p->inphys;
Chris@82 39 p->ophyssz = p->iphyssz;
Chris@82 40 } else {
Chris@82 41 p->ophyssz = osz;
Chris@82 42 p->outphys = out = (bench_complex *) bench_malloc(osz * sizeof(bench_complex));
Chris@82 43 p->out = out - olb;
Chris@82 44 }
Chris@82 45 } else if (p->kind == PROBLEM_R2R) {
Chris@82 46 bench_real *in, *out;
Chris@82 47
Chris@82 48 p->iphyssz = isz;
Chris@82 49 p->inphys = in = (bench_real *) bench_malloc(isz * sizeof(bench_real));
Chris@82 50 p->in = in - ilb;
Chris@82 51
Chris@82 52 if (p->in_place) {
Chris@82 53 p->out = p->in;
Chris@82 54 p->outphys = p->inphys;
Chris@82 55 p->ophyssz = p->iphyssz;
Chris@82 56 } else {
Chris@82 57 p->ophyssz = osz;
Chris@82 58 p->outphys = out = (bench_real *) bench_malloc(osz * sizeof(bench_real));
Chris@82 59 p->out = out - olb;
Chris@82 60 }
Chris@82 61 } else if (p->kind == PROBLEM_REAL && p->sign < 0) { /* R2HC */
Chris@82 62 bench_real *in;
Chris@82 63 bench_complex *out;
Chris@82 64
Chris@82 65 isz = isz > osz*2 ? isz : osz*2;
Chris@82 66 p->iphyssz = isz;
Chris@82 67 p->inphys = in = (bench_real *) bench_malloc(p->iphyssz * sizeof(bench_real));
Chris@82 68 p->in = in - ilb;
Chris@82 69
Chris@82 70 if (p->in_place) {
Chris@82 71 p->out = p->in;
Chris@82 72 p->outphys = p->inphys;
Chris@82 73 p->ophyssz = p->iphyssz / 2;
Chris@82 74 } else {
Chris@82 75 p->ophyssz = osz;
Chris@82 76 p->outphys = out = (bench_complex *) bench_malloc(osz * sizeof(bench_complex));
Chris@82 77 p->out = out - olb;
Chris@82 78 }
Chris@82 79 } else if (p->kind == PROBLEM_REAL && p->sign > 0) { /* HC2R */
Chris@82 80 bench_real *out;
Chris@82 81 bench_complex *in;
Chris@82 82
Chris@82 83 osz = osz > isz*2 ? osz : isz*2;
Chris@82 84 p->ophyssz = osz;
Chris@82 85 p->outphys = out = (bench_real *) bench_malloc(p->ophyssz * sizeof(bench_real));
Chris@82 86 p->out = out - olb;
Chris@82 87
Chris@82 88 if (p->in_place) {
Chris@82 89 p->in = p->out;
Chris@82 90 p->inphys = p->outphys;
Chris@82 91 p->iphyssz = p->ophyssz / 2;
Chris@82 92 } else {
Chris@82 93 p->iphyssz = isz;
Chris@82 94 p->inphys = in = (bench_complex *) bench_malloc(isz * sizeof(bench_complex));
Chris@82 95 p->in = in - ilb;
Chris@82 96 }
Chris@82 97 } else {
Chris@82 98 BENCH_ASSERT(0); /* TODO */
Chris@82 99 }
Chris@82 100 }
Chris@82 101
Chris@82 102 void problem_free(bench_problem *p)
Chris@82 103 {
Chris@82 104 if (p->outphys && p->outphys != p->inphys)
Chris@82 105 bench_free(p->outphys);
Chris@82 106 if (p->inphys)
Chris@82 107 bench_free(p->inphys);
Chris@82 108 tensor_destroy(p->sz);
Chris@82 109 tensor_destroy(p->vecsz);
Chris@82 110 }