cannam@127: /* not worth copyrighting */
cannam@127: 
cannam@127: #include "bench.h"
cannam@127: #include <math.h>
cannam@127: 
cannam@127: double mflops(const bench_problem *p, double t)
cannam@127: {
cannam@127:      int size = tensor_sz(p->sz);
cannam@127:      int vsize = tensor_sz(p->vecsz);
cannam@127: 
cannam@127:      if (size <= 1) /* a copy: just return reals copied / time */
cannam@127: 	  switch (p->kind) {
cannam@127: 	      case PROBLEM_COMPLEX:
cannam@127: 		   return (2.0 * size * vsize / (t * 1.0e6));
cannam@127: 	      case PROBLEM_REAL:
cannam@127: 	      case PROBLEM_R2R:
cannam@127: 		   return (1.0 * size * vsize / (t * 1.0e6));
cannam@127: 	  }
cannam@127: 
cannam@127:      switch (p->kind) {
cannam@127: 	 case PROBLEM_COMPLEX:
cannam@127: 	      return (5.0 * size * vsize * log((double)size) / 
cannam@127: 		      (log(2.0) * t * 1.0e6));
cannam@127: 	 case PROBLEM_REAL:
cannam@127: 	 case PROBLEM_R2R:
cannam@127: 	      return (2.5 * vsize * size * log((double) size) / 
cannam@127: 		      (log(2.0) * t * 1.0e6));
cannam@127:      }
cannam@127:      BENCH_ASSERT(0 /* can't happen */);
cannam@127:      return 0.0;
cannam@127: }
cannam@127: