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: int no_speed_allocation = 0; /* 1 to not allocate array data in speed() */
Chris@42: 
Chris@42: void speed(const char *param, int setup_only)
Chris@42: {
Chris@42:      double *t;
Chris@42:      int iter = 0, k;
Chris@42:      bench_problem *p;
Chris@42:      double tmin, y;
Chris@42: 
Chris@42:      t = (double *) bench_malloc(time_repeat * sizeof(double));
Chris@42: 
Chris@42:      for (k = 0; k < time_repeat; ++k) 
Chris@42: 	  t[k] = 0;
Chris@42: 
Chris@42:      p = problem_parse(param);
Chris@42:      BENCH_ASSERT(can_do(p));
Chris@42:      if (!no_speed_allocation) {
Chris@42: 	  problem_alloc(p);
Chris@42: 	  problem_zero(p);
Chris@42:      }
Chris@42: 
Chris@42:      timer_start(LIBBENCH_TIMER);
Chris@42:      setup(p);
Chris@42:      p->setup_time = bench_cost_postprocess(timer_stop(LIBBENCH_TIMER));
Chris@42: 
Chris@42:      /* reset the input to zero again, because the planner in paranoid
Chris@42: 	mode sets it to random values, thus making the benchmark
Chris@42: 	diverge. */
Chris@42:      if (!no_speed_allocation) 
Chris@42: 	  problem_zero(p);
Chris@42:      
Chris@42:      if (setup_only)
Chris@42: 	  goto done;
Chris@42: 
Chris@42:  start_over:
Chris@42:      for (iter = 1; iter < (1<<30); iter *= 2) {
Chris@42: 	  tmin = 1.0e20;
Chris@42: 	  for (k = 0; k < time_repeat; ++k) {
Chris@42: 	       timer_start(LIBBENCH_TIMER);
Chris@42: 	       doit(iter, p);
Chris@42: 	       y = bench_cost_postprocess(timer_stop(LIBBENCH_TIMER));
Chris@42: 	       if (y < 0) /* yes, it happens */
Chris@42: 		    goto start_over;
Chris@42: 	       t[k] = y;
Chris@42: 	       if (y < tmin)
Chris@42: 		    tmin = y;
Chris@42: 	  }
Chris@42: 	  
Chris@42: 	  if (tmin >= time_min)
Chris@42: 	       goto done;
Chris@42:      }
Chris@42: 
Chris@42:      goto start_over; /* this also happens */
Chris@42: 
Chris@42:  done:
Chris@42:      done(p);
Chris@42: 
Chris@42:      if (iter) 
Chris@42: 	  for (k = 0; k < time_repeat; ++k) 
Chris@42: 	       t[k] /= iter;
Chris@42:      else
Chris@42: 	  for (k = 0; k < time_repeat; ++k) 
Chris@42: 	       t[k] = 0;
Chris@42: 
Chris@42:      report(p, t, time_repeat);
Chris@42: 
Chris@42:      if (!no_speed_allocation)
Chris@42: 	  problem_destroy(p);
Chris@42:      bench_free(t);
Chris@42:      return;
Chris@42: }