cannam@95: /* See bench.c.  We keep a few common subroutines in this file so
cannam@95:    that they can be re-used in the MPI test program. */
cannam@95: 
cannam@95: #include <math.h>
cannam@95: #include <stdio.h>
cannam@95: #include <string.h>
cannam@95: #include "fftw-bench.h"
cannam@95: 
cannam@95: #ifdef _OPENMP
cannam@95: #  include <omp.h>
cannam@95: #endif
cannam@95: 
cannam@95: #ifdef HAVE_SMP
cannam@95: int threads_ok = 1;
cannam@95: #endif
cannam@95: 
cannam@95: FFTW(plan) the_plan = 0;
cannam@95: 
cannam@95: static const char *wisdat = "wis.dat";
cannam@95: unsigned the_flags = 0;
cannam@95: int paranoid = 0;
cannam@95: int usewisdom = 0;
cannam@95: int havewisdom = 0;
cannam@95: int nthreads = 1;
cannam@95: int amnesia = 0;
cannam@95: 
cannam@95: extern void install_hook(void);  /* in hook.c */
cannam@95: extern void uninstall_hook(void);  /* in hook.c */
cannam@95: 
cannam@95: #ifdef FFTW_RANDOM_ESTIMATOR
cannam@95: extern unsigned FFTW(random_estimate_seed);
cannam@95: #endif
cannam@95: 
cannam@95: void useropt(const char *arg)
cannam@95: {
cannam@95:      int x;
cannam@95:      double y;
cannam@95: 
cannam@95:      if (!strcmp(arg, "patient")) the_flags |= FFTW_PATIENT;
cannam@95:      else if (!strcmp(arg, "estimate")) the_flags |= FFTW_ESTIMATE;
cannam@95:      else if (!strcmp(arg, "estimatepat")) the_flags |= FFTW_ESTIMATE_PATIENT;
cannam@95:      else if (!strcmp(arg, "exhaustive")) the_flags |= FFTW_EXHAUSTIVE;
cannam@95:      else if (!strcmp(arg, "unaligned")) the_flags |= FFTW_UNALIGNED;
cannam@95:      else if (!strcmp(arg, "nosimd")) the_flags |= FFTW_NO_SIMD;
cannam@95:      else if (!strcmp(arg, "noindirectop")) the_flags |= FFTW_NO_INDIRECT_OP;
cannam@95:      else if (!strcmp(arg, "wisdom-only")) the_flags |= FFTW_WISDOM_ONLY;
cannam@95:      else if (sscanf(arg, "flag=%d", &x) == 1) the_flags |= x;
cannam@95:      else if (sscanf(arg, "bflag=%d", &x) == 1) the_flags |= 1U << x;
cannam@95:      else if (!strcmp(arg, "paranoid")) paranoid = 1;
cannam@95:      else if (!strcmp(arg, "wisdom")) usewisdom = 1;
cannam@95:      else if (!strcmp(arg, "amnesia")) amnesia = 1;
cannam@95:      else if (sscanf(arg, "nthreads=%d", &x) == 1) nthreads = x;
cannam@95: #ifdef FFTW_RANDOM_ESTIMATOR
cannam@95:      else if (sscanf(arg, "eseed=%d", &x) == 1) FFTW(random_estimate_seed) = x;
cannam@95: #endif
cannam@95:      else if (sscanf(arg, "timelimit=%lg", &y) == 1) {
cannam@95: 	  FFTW(set_timelimit)(y);
cannam@95:      }
cannam@95: 
cannam@95:      else fprintf(stderr, "unknown user option: %s.  Ignoring.\n", arg);
cannam@95: }
cannam@95: 
cannam@95: void rdwisdom(void)
cannam@95: {
cannam@95:      FILE *f;
cannam@95:      double tim;
cannam@95:      int success = 0;
cannam@95: 
cannam@95:      if (havewisdom) return;
cannam@95: 
cannam@95: #ifdef HAVE_SMP
cannam@95:      if (threads_ok) {
cannam@95: 	  BENCH_ASSERT(FFTW(init_threads)());
cannam@95: 	  FFTW(plan_with_nthreads)(nthreads);
cannam@95: #ifdef _OPENMP
cannam@95: 	  omp_set_num_threads(nthreads);
cannam@95: #endif
cannam@95:      }
cannam@95:      else if (nthreads > 1 && verbose > 1) {
cannam@95: 	  fprintf(stderr, "bench: WARNING - nthreads = %d, but threads not supported\n", nthreads);
cannam@95: 	  nthreads = 1;
cannam@95:      }
cannam@95: #endif
cannam@95: 
cannam@95:      if (!usewisdom) return;
cannam@95: 
cannam@95:      timer_start(USER_TIMER);
cannam@95:      if ((f = fopen(wisdat, "r"))) {
cannam@95: 	  if (!import_wisdom(f))
cannam@95: 	       fprintf(stderr, "bench: ERROR reading wisdom\n");
cannam@95: 	  else
cannam@95: 	       success = 1;
cannam@95: 	  fclose(f);
cannam@95:      }
cannam@95:      tim = timer_stop(USER_TIMER);
cannam@95: 
cannam@95:      if (success) {
cannam@95: 	  if (verbose > 1) printf("READ WISDOM (%g seconds): ", tim);
cannam@95: 	  
cannam@95: 	  if (verbose > 3)
cannam@95: 	       export_wisdom(stdout);
cannam@95: 	  if (verbose > 1)
cannam@95: 	       printf("\n");
cannam@95:      }
cannam@95:      havewisdom = 1;
cannam@95: }
cannam@95: 
cannam@95: void wrwisdom(void)
cannam@95: {
cannam@95:      FILE *f;
cannam@95:      double tim;
cannam@95:      if (!havewisdom) return;
cannam@95: 
cannam@95:      timer_start(USER_TIMER);
cannam@95:      if ((f = fopen(wisdat, "w"))) {
cannam@95: 	  export_wisdom(f);
cannam@95: 	  fclose(f);
cannam@95:      }
cannam@95:      tim = timer_stop(USER_TIMER);
cannam@95:      if (verbose > 1) printf("write wisdom took %g seconds\n", tim);
cannam@95: }
cannam@95: 
cannam@95: static unsigned preserve_input_flags(bench_problem *p)
cannam@95: {
cannam@95:      /*
cannam@95:       * fftw3 cannot preserve input for multidimensional c2r transforms.
cannam@95:       * Enforce FFTW_DESTROY_INPUT
cannam@95:       */
cannam@95:      if (p->kind == PROBLEM_REAL && 
cannam@95: 	 p->sign > 0 && 
cannam@95: 	 !p->in_place && 
cannam@95: 	 p->sz->rnk > 1)
cannam@95: 	  p->destroy_input = 1;
cannam@95: 
cannam@95:      if (p->destroy_input)
cannam@95: 	  return FFTW_DESTROY_INPUT;
cannam@95:      else
cannam@95: 	  return FFTW_PRESERVE_INPUT;
cannam@95: }
cannam@95: 
cannam@95: int can_do(bench_problem *p)
cannam@95: {
cannam@95:      double tim;
cannam@95: 
cannam@95:      if (verbose > 2 && p->pstring)
cannam@95: 	  printf("Planning %s...\n", p->pstring);
cannam@95:      rdwisdom();
cannam@95: 
cannam@95:      timer_start(USER_TIMER);
cannam@95:      the_plan = mkplan(p, preserve_input_flags(p) | the_flags | FFTW_ESTIMATE);
cannam@95:      tim = timer_stop(USER_TIMER);
cannam@95:      if (verbose > 2) printf("estimate-planner time: %g s\n", tim);
cannam@95: 
cannam@95:      if (the_plan) {
cannam@95: 	  FFTW(destroy_plan)(the_plan);
cannam@95: 	  return 1;
cannam@95:      }
cannam@95:      return 0;
cannam@95: }
cannam@95: 
cannam@95: void setup(bench_problem *p)
cannam@95: {
cannam@95:      double tim;
cannam@95: 
cannam@95:      if (amnesia) {
cannam@95: 	  FFTW(forget_wisdom)();
cannam@95: 	  havewisdom = 0;
cannam@95:      }
cannam@95: 
cannam@95:      /* Regression test: check that fftw_malloc exists and links
cannam@95:       * properly */
cannam@95:      FFTW(free(FFTW(malloc(42))));
cannam@95: 
cannam@95:      rdwisdom();
cannam@95:      install_hook();
cannam@95: 
cannam@95: #ifdef HAVE_SMP
cannam@95:      if (verbose > 1 && nthreads > 1) printf("NTHREADS = %d\n", nthreads);
cannam@95: #endif
cannam@95: 
cannam@95:      timer_start(USER_TIMER);
cannam@95:      the_plan = mkplan(p, preserve_input_flags(p) | the_flags);
cannam@95:      tim = timer_stop(USER_TIMER);
cannam@95:      if (verbose > 1) printf("planner time: %g s\n", tim);
cannam@95: 
cannam@95:      BENCH_ASSERT(the_plan);
cannam@95:      
cannam@95:      {
cannam@95: 	  double add, mul, nfma, cost, pcost;
cannam@95: 	  FFTW(flops)(the_plan, &add, &mul, &nfma);
cannam@95: 	  cost = FFTW(estimate_cost)(the_plan);
cannam@95: 	  pcost = FFTW(cost)(the_plan);
cannam@95: 	  if (verbose > 1) {
cannam@95: 	       FFTW(print_plan)(the_plan);
cannam@95: 	       printf("\n");
cannam@95: 	       printf("flops: %0.0f add, %0.0f mul, %0.0f fma\n",
cannam@95: 		      add, mul, nfma);
cannam@95: 	       printf("estimated cost: %f, pcost = %f\n", cost, pcost);
cannam@95: 	  }
cannam@95:      }
cannam@95: }
cannam@95: 
cannam@95: 
cannam@95: void doit(int iter, bench_problem *p)
cannam@95: {
cannam@95:      int i;
cannam@95:      FFTW(plan) q = the_plan;
cannam@95: 
cannam@95:      UNUSED(p);
cannam@95:      for (i = 0; i < iter; ++i) 
cannam@95: 	  FFTW(execute)(q);
cannam@95: }
cannam@95: 
cannam@95: void done(bench_problem *p)
cannam@95: {
cannam@95:      UNUSED(p);
cannam@95: 
cannam@95:      FFTW(destroy_plan)(the_plan);
cannam@95:      uninstall_hook();
cannam@95: }
cannam@95: 
cannam@95: void cleanup(void)
cannam@95: {
cannam@95:      initial_cleanup();
cannam@95: 
cannam@95:      wrwisdom();
cannam@95: #ifdef HAVE_SMP
cannam@95:      FFTW(cleanup_threads)();
cannam@95: #else
cannam@95:      FFTW(cleanup)();
cannam@95: #endif
cannam@95: 
cannam@95: #    ifdef FFTW_DEBUG_MALLOC
cannam@95:      {
cannam@95: 	  /* undocumented memory checker */
cannam@95: 	  FFTW_EXTERN void FFTW(malloc_print_minfo)(int v);
cannam@95: 	  FFTW(malloc_print_minfo)(verbose);
cannam@95:      }
cannam@95: #    endif
cannam@95: 
cannam@95:      final_cleanup();
cannam@95: }