diff options
Diffstat (limited to 'Microbenchmarks/Arithmetic-Throughput')
| -rw-r--r-- | Microbenchmarks/Arithmetic-Throughput/Makefile | 46 | ||||
| -rw-r--r-- | Microbenchmarks/Arithmetic-Throughput/dpu/task.c | 100 | ||||
| -rw-r--r-- | Microbenchmarks/Arithmetic-Throughput/host/app.c | 224 | ||||
| -rwxr-xr-x | Microbenchmarks/Arithmetic-Throughput/run.sh | 17 | ||||
| -rwxr-xr-x | Microbenchmarks/Arithmetic-Throughput/support/common.h | 48 | ||||
| -rw-r--r-- | Microbenchmarks/Arithmetic-Throughput/support/cyclecount.h | 20 | ||||
| -rw-r--r-- | Microbenchmarks/Arithmetic-Throughput/support/params.h | 56 | ||||
| -rwxr-xr-x | Microbenchmarks/Arithmetic-Throughput/support/timer.h | 59 |
8 files changed, 570 insertions, 0 deletions
diff --git a/Microbenchmarks/Arithmetic-Throughput/Makefile b/Microbenchmarks/Arithmetic-Throughput/Makefile new file mode 100644 index 0000000..137e7b5 --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/Makefile @@ -0,0 +1,46 @@ +DPU_DIR := dpu +HOST_DIR := host +BUILDDIR ?= bin +NR_TASKLETS ?= 16 +BL ?= 10 +NR_DPUS ?= 1 +OP ?= ADD +TYPE ?= INT32 + +define conf_filename + ${BUILDDIR}/.NR_DPUS_$(1)_NR_TASKLETS_$(2)_BL_$(3)_$(4)_$(5).conf +endef +CONF := $(call conf_filename,${NR_DPUS},${NR_TASKLETS},${BL},${OP},${TYPE}) + +HOST_TARGET := ${BUILDDIR}/host_code +DPU_TARGET := ${BUILDDIR}/dpu_code + +COMMON_INCLUDES := support +HOST_SOURCES := $(wildcard ${HOST_DIR}/*.c) +DPU_SOURCES := $(wildcard ${DPU_DIR}/*.c) + +.PHONY: all clean test + +__dirs := $(shell mkdir -p ${BUILDDIR}) + +COMMON_FLAGS := -Wall -Wextra -g -I${COMMON_INCLUDES} +HOST_FLAGS := ${COMMON_FLAGS} -std=c11 -O3 `dpu-pkg-config --cflags --libs dpu` -DNR_TASKLETS=${NR_TASKLETS} -DNR_DPUS=${NR_DPUS} -DBL=${BL} -D${OP} -D${TYPE} +DPU_FLAGS := ${COMMON_FLAGS} -O2 -flto -DNR_TASKLETS=${NR_TASKLETS} -DBL=${BL} -D${OP} -D${TYPE} + +all: ${HOST_TARGET} ${DPU_TARGET} + +${CONF}: + $(RM) $(call conf_filename,*,*) + touch ${CONF} + +${HOST_TARGET}: ${HOST_SOURCES} ${COMMON_INCLUDES} ${CONF} + $(CC) -o $@ ${HOST_SOURCES} ${HOST_FLAGS} + +${DPU_TARGET}: ${DPU_SOURCES} ${COMMON_INCLUDES} ${CONF} + dpu-upmem-dpurte-clang ${DPU_FLAGS} -o $@ ${DPU_SOURCES} + +clean: + $(RM) -r $(BUILDDIR) + +test: all + ./${HOST_TARGET} diff --git a/Microbenchmarks/Arithmetic-Throughput/dpu/task.c b/Microbenchmarks/Arithmetic-Throughput/dpu/task.c new file mode 100644 index 0000000..2f538ff --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/dpu/task.c @@ -0,0 +1,100 @@ +/* +* Execution of arithmetic operations with multiple tasklets +* +*/ +#include <stdint.h> +#include <stdio.h> +#include <defs.h> +#include <mram.h> +#include <alloc.h> +#include <perfcounter.h> +#include <barrier.h> + +#include "../support/common.h" +#include "../support/cyclecount.h" + +__host dpu_arguments_t DPU_INPUT_ARGUMENTS; +__host dpu_results_t DPU_RESULTS[NR_TASKLETS]; + +// Arithmetic operation +static void update(T *bufferA, T scalar) { + //#pragma unroll + for (unsigned int i = 0; i < BLOCK_SIZE / sizeof(T); i++){ + // WRAM READ + T temp = bufferA[i]; +#ifdef ADD + temp += scalar; // ADD +#elif SUB + temp -= scalar; // SUB +#elif MUL + temp *= scalar; // MUL +#elif DIV + temp /= scalar; // DIV +#endif + // WRAM WRITE + bufferA[i] = temp; + } +} + +// Barrier +BARRIER_INIT(my_barrier, NR_TASKLETS); + +extern int main_kernel1(void); + +int (*kernels[nr_kernels])(void) = {main_kernel1}; + +int main(void) { + // Kernel + return kernels[DPU_INPUT_ARGUMENTS.kernel](); +} + +// main_kernel1 +int main_kernel1() { + unsigned int tasklet_id = me(); +#if PRINT + printf("tasklet_id = %u\n", tasklet_id); +#endif + if (tasklet_id == 0){ // Initialize once the cycle counter + mem_reset(); // Reset the heap + + perfcounter_config(COUNT_CYCLES, true); + } + // Barrier + barrier_wait(&my_barrier); + perfcounter_cycles cycles; + + uint32_t input_size_dpu = DPU_INPUT_ARGUMENTS.size / sizeof(T); + T scalar = (T)input_size_dpu; // Simply use this number as a scalar + + dpu_results_t *result = &DPU_RESULTS[tasklet_id]; + result->cycles = 0; + + // Address of the current processing block in MRAM + uint32_t mram_base_addr_A = (uint32_t)(DPU_MRAM_HEAP_POINTER + (tasklet_id << BLOCK_SIZE_LOG2)); + uint32_t mram_base_addr_B = (uint32_t)(DPU_MRAM_HEAP_POINTER + (tasklet_id << BLOCK_SIZE_LOG2) + input_size_dpu * sizeof(T)); + + // Initialize a local cache to store the MRAM block + T *cache_A = (T *) mem_alloc(BLOCK_SIZE); + + for(unsigned int byte_index = 0; byte_index < input_size_dpu * sizeof(T); byte_index += BLOCK_SIZE * NR_TASKLETS){ + + // Load cache with current MRAM block + mram_read((__mram_ptr void const*)(mram_base_addr_A + byte_index), cache_A, BLOCK_SIZE); + + // Barrier + barrier_wait(&my_barrier); + timer_start(&cycles); // START TIMER + + // Update + update(cache_A, scalar); + + result->cycles += timer_stop(&cycles); // STOP TIMER + // Barrier + barrier_wait(&my_barrier); + + // Write cache to current MRAM block + mram_write(cache_A, (__mram_ptr void*)(mram_base_addr_B + byte_index), BLOCK_SIZE); + } + + return 0; +} diff --git a/Microbenchmarks/Arithmetic-Throughput/host/app.c b/Microbenchmarks/Arithmetic-Throughput/host/app.c new file mode 100644 index 0000000..5933d1d --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/host/app.c @@ -0,0 +1,224 @@ +/** +* app.c +* Arithmetic Throughput Host Application Source File +* +*/ +#include <stdio.h> +#include <stdlib.h> +#include <stdbool.h> +#include <string.h> +#include <dpu.h> +#include <dpu_log.h> +#include <unistd.h> +#include <getopt.h> +#include <assert.h> + +#include "../support/common.h" +#include "../support/timer.h" +#include "../support/params.h" + +// Define the DPU Binary path as DPU_BINARY here +#ifndef DPU_BINARY +#define DPU_BINARY "./bin/dpu_code" +#endif + +// Pointer declaration +static T* A; +static T* B; +static T* C2; + +// Create input arrays +static void read_input(T* A, T* B, unsigned int nr_elements) { + srand(0); + printf("nr_elements\t%u\t", nr_elements); + for (unsigned int i = 0; i < nr_elements; i++) { + A[i] = (T) (rand()); + B[i] = (T) (rand()); + } +} + +// Compute output in the host +static void update_host(T* C, T* A, unsigned int nr_elements) { + for (unsigned int i = 0; i < nr_elements; i++) { +#if ADD + C[i] = A[i] + (nr_elements / NR_DPUS); +#elif SUB + C[i] = A[i] - (nr_elements / NR_DPUS); +#elif MUL + C[i] = A[i] * (nr_elements / NR_DPUS); +#elif DIV + C[i] = A[i] / (nr_elements / NR_DPUS); +#endif + } +} + +// Main of the Host Application +int main(int argc, char **argv) { + + struct Params p = input_params(argc, argv); + + struct dpu_set_t dpu_set, dpu; + uint32_t nr_of_dpus; + + // Allocate DPUs and load binary + DPU_ASSERT(dpu_alloc(NR_DPUS, NULL, &dpu_set)); + DPU_ASSERT(dpu_load(dpu_set, DPU_BINARY, NULL)); + DPU_ASSERT(dpu_get_nr_dpus(dpu_set, &nr_of_dpus)); + printf("Allocated %d DPU(s)\n", nr_of_dpus); + + unsigned int i = 0; + double cc = 0; + double cc_min = 0; + const unsigned int input_size = p.exp == 0 ? p.input_size * nr_of_dpus : p.input_size; + + // Input/output allocation + A = malloc(input_size * sizeof(T)); + B = malloc(input_size * sizeof(T)); + T *bufferA = A; + T *bufferB = B; + C2 = malloc(input_size * sizeof(T)); + + // Create an input file with arbitrary data + read_input(A, B, input_size); + + // Timer declaration + Timer timer; + + printf("NR_TASKLETS\t%d\tBL\t%d\n", NR_TASKLETS, BL); + + // Loop over main kernel + for(int rep = 0; rep < p.n_warmup + p.n_reps; rep++) { + + // Compute output on CPU (performance comparison and verification purposes) + if(rep >= p.n_warmup) + start(&timer, 0, rep - p.n_warmup); + update_host(C2, A, input_size); + if(rep >= p.n_warmup) + stop(&timer, 0); + + printf("Load input data\n"); + if(rep >= p.n_warmup) + start(&timer, 1, rep - p.n_warmup); + // Input arguments + const unsigned int input_size_dpu = input_size / nr_of_dpus; + unsigned int kernel = 0; + dpu_arguments_t input_arguments = {input_size_dpu * sizeof(T), kernel}; + DPU_ASSERT(dpu_copy_to(dpu_set, "DPU_INPUT_ARGUMENTS", 0, (const void *)&input_arguments, sizeof(input_arguments))); + // Copy input arrays + i = 0; + DPU_FOREACH (dpu_set, dpu) { + DPU_ASSERT(dpu_copy_to(dpu, DPU_MRAM_HEAP_POINTER_NAME, 0, bufferA + input_size_dpu * i, input_size_dpu * sizeof(T))); + i++; + } + if(rep >= p.n_warmup) + stop(&timer, 1); + + printf("Run program on DPU(s) \n"); + // Run DPU kernel + if(rep >= p.n_warmup) + start(&timer, 2, rep - p.n_warmup); + DPU_ASSERT(dpu_launch(dpu_set, DPU_SYNCHRONOUS)); + if(rep >= p.n_warmup) + stop(&timer, 2); + +#if PRINT + { + unsigned int each_dpu = 0; + printf("Display DPU Logs\n"); + DPU_FOREACH (dpu_set, dpu) { + printf("DPU#%d:\n", each_dpu); + DPU_ASSERT(dpulog_read_for_dpu(dpu.dpu, stdout)); + each_dpu++; + } + } +#endif + + printf("Retrieve results\n"); + if(rep >= p.n_warmup) + start(&timer, 3, rep - p.n_warmup); + dpu_results_t results[nr_of_dpus]; + i = 0; + DPU_FOREACH (dpu_set, dpu) { + // Copy output array + DPU_ASSERT(dpu_copy_from(dpu, DPU_MRAM_HEAP_POINTER_NAME, input_size_dpu * sizeof(T), bufferB + input_size_dpu * i, input_size_dpu * sizeof(T))); + +#if PERF + results[i].cycles = 0; + // Retrieve tasklet timings + for (unsigned int each_tasklet = 0; each_tasklet < NR_TASKLETS; each_tasklet++) { + dpu_results_t result; + result.cycles = 0; + DPU_ASSERT(dpu_copy_from(dpu, "DPU_RESULTS", each_tasklet * sizeof(dpu_results_t), &result, sizeof(dpu_results_t))); + if (result.cycles > results[i].cycles) + results[i].cycles = result.cycles; + } +#endif + i++; + } + if(rep >= p.n_warmup) + stop(&timer, 3); + +#if PERF + uint64_t max_cycles = 0; + uint64_t min_cycles = 0xFFFFFFFFFFFFFFFF; + // Print performance results + if(rep >= p.n_warmup){ + i = 0; + DPU_FOREACH(dpu_set, dpu) { + if(results[i].cycles > max_cycles) + max_cycles = results[i].cycles; + if(results[i].cycles < min_cycles) + min_cycles = results[i].cycles; + i++; + } + cc += (double)max_cycles; + cc_min += (double)min_cycles; + } +#endif + + } +#ifdef ADD + printf("ADD\n"); +#elif SUB + printf("SUB\n"); +#elif MUL + printf("MUL\n"); +#elif DIV + printf("DIV\n"); +#endif + printf("DPU cycles = %g cc\n", cc / p.n_reps); + + // Print timing results + printf("CPU "); + print(&timer, 0, p.n_reps); + printf("CPU-DPU "); + print(&timer, 1, p.n_reps); + printf("DPU Kernel "); + print(&timer, 2, p.n_reps); + printf("DPU-CPU "); + print(&timer, 3, p.n_reps); + + // Check output + bool status = true; + for (i = 0; i < input_size; i++) { + if(C2[i] != bufferB[i]){ + status = false; +#if PRINT + printf("%d: %u -- %u\n", i, C2[i], bufferB[i]); +#endif + } + } + if (status) { + printf("[" ANSI_COLOR_GREEN "OK" ANSI_COLOR_RESET "] Outputs are equal\n"); + } else { + printf("[" ANSI_COLOR_RED "ERROR" ANSI_COLOR_RESET "] Outputs differ!\n"); + } + + // Deallocation + free(A); + free(B); + free(C2); + DPU_ASSERT(dpu_free(dpu_set)); + + return status ? 0 : -1; +} diff --git a/Microbenchmarks/Arithmetic-Throughput/run.sh b/Microbenchmarks/Arithmetic-Throughput/run.sh new file mode 100755 index 0000000..9200e54 --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/run.sh @@ -0,0 +1,17 @@ +#!/bin/bash + +for i in ADD SUB MUL DIV +do + for j in INT32 FLOAT UINT32 INT64 DOUBLE UINT64 + do + for k in 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 + do + NR_DPUS=1 NR_TASKLETS=$k BL=10 OP=$i TYPE=$j make all + wait + ./bin/host_code -w 0 -e 1 -i 1048576 > profile/${i}_${j}_tl${k}.txt + wait + make clean + wait + done + done +done diff --git a/Microbenchmarks/Arithmetic-Throughput/support/common.h b/Microbenchmarks/Arithmetic-Throughput/support/common.h new file mode 100755 index 0000000..0d1cc8a --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/support/common.h @@ -0,0 +1,48 @@ +#ifndef _COMMON_H_ +#define _COMMON_H_ + +// Structures used by both the host and the dpu to communicate information +typedef struct { + uint32_t size; + enum kernels { + kernel1 = 0, + nr_kernels = 1, + } kernel; +} dpu_arguments_t; + +typedef struct { + uint64_t cycles; +} dpu_results_t; + +// Transfer size between MRAM and WRAM +#ifdef BL +#define BLOCK_SIZE_LOG2 BL +#define BLOCK_SIZE (1 << BLOCK_SIZE_LOG2) +#else +#define BLOCK_SIZE_LOG2 8 +#define BLOCK_SIZE (1 << BLOCK_SIZE_LOG2) +#define BL BLOCK_SIZE_LOG2 +#endif + +// Data type +#ifdef UINT32 +#define T uint32_t +#elif UINT64 +#define T uint64_t +#elif INT32 +#define T int32_t +#elif INT64 +#define T int64_t +#elif FLOAT +#define T float +#elif DOUBLE +#define T double +#endif + +#define PERF 1 // Use perfcounters? +#define PRINT 0 + +#define ANSI_COLOR_RED "\x1b[31m" +#define ANSI_COLOR_GREEN "\x1b[32m" +#define ANSI_COLOR_RESET "\x1b[0m" +#endif diff --git a/Microbenchmarks/Arithmetic-Throughput/support/cyclecount.h b/Microbenchmarks/Arithmetic-Throughput/support/cyclecount.h new file mode 100644 index 0000000..4ef377e --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/support/cyclecount.h @@ -0,0 +1,20 @@ +#include <perfcounter.h> + +// Timer +typedef struct perfcounter_cycles{ + perfcounter_t start; + perfcounter_t end; + perfcounter_t end2; + +}perfcounter_cycles; + +void timer_start(perfcounter_cycles *cycles){ + cycles->start = perfcounter_get(); // START TIMER +} + +uint64_t timer_stop(perfcounter_cycles *cycles){ + cycles->end = perfcounter_get(); // STOP TIMER + cycles->end2 = perfcounter_get(); // STOP TIMER + return(((uint64_t)((uint32_t)(((cycles->end >> 4) - (cycles->start >> 4)) - ((cycles->end2 >> 4) - (cycles->end >> 4))))) << 4); +} + diff --git a/Microbenchmarks/Arithmetic-Throughput/support/params.h b/Microbenchmarks/Arithmetic-Throughput/support/params.h new file mode 100644 index 0000000..4618411 --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/support/params.h @@ -0,0 +1,56 @@ +#ifndef _PARAMS_H_ +#define _PARAMS_H_ + +#include "common.h" + +typedef struct Params { + unsigned int input_size; + int n_warmup; + int n_reps; + int exp; +}Params; + +static void usage() { + fprintf(stderr, + "\nUsage: ./program [options]" + "\n" + "\nGeneral options:" + "\n -h help" + "\n -w <W> # of untimed warmup iterations (default=1)" + "\n -e <E> # of timed repetition iterations (default=3)" + "\n -x <X> Weak (0) or strong (1) scaling (default=0)" + "\n" + "\nBenchmark-specific options:" + "\n -i <I> input size (default=8K elements)" + "\n"); +} + +struct Params input_params(int argc, char **argv) { + struct Params p; + p.input_size = 8 << 10; + p.n_warmup = 1; + p.n_reps = 3; + p.exp = 0; + + int opt; + while((opt = getopt(argc, argv, "hi:w:e:x:")) >= 0) { + switch(opt) { + case 'h': + usage(); + exit(0); + break; + case 'i': p.input_size = atoi(optarg); break; + case 'w': p.n_warmup = atoi(optarg); break; + case 'e': p.n_reps = atoi(optarg); break; + case 'x': p.exp = atoi(optarg); break; + default: + fprintf(stderr, "\nUnrecognized option!\n"); + usage(); + exit(0); + } + } + assert(NR_DPUS > 0 && "Invalid # of dpus!"); + + return p; +} +#endif diff --git a/Microbenchmarks/Arithmetic-Throughput/support/timer.h b/Microbenchmarks/Arithmetic-Throughput/support/timer.h new file mode 100755 index 0000000..eedc385 --- /dev/null +++ b/Microbenchmarks/Arithmetic-Throughput/support/timer.h @@ -0,0 +1,59 @@ +/*
+ * Copyright (c) 2016 University of Cordoba and University of Illinois
+ * All rights reserved.
+ *
+ * Developed by: IMPACT Research Group
+ * University of Cordoba and University of Illinois
+ * http://impact.crhc.illinois.edu/
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a copy
+ * of this software and associated documentation files (the "Software"), to deal
+ * with the Software without restriction, including without limitation the
+ * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
+ * sell copies of the Software, and to permit persons to whom the Software is
+ * furnished to do so, subject to the following conditions:
+ *
+ * > Redistributions of source code must retain the above copyright notice,
+ * this list of conditions and the following disclaimers.
+ * > Redistributions in binary form must reproduce the above copyright
+ * notice, this list of conditions and the following disclaimers in the
+ * documentation and/or other materials provided with the distribution.
+ * > Neither the names of IMPACT Research Group, University of Cordoba,
+ * University of Illinois nor the names of its contributors may be used
+ * to endorse or promote products derived from this Software without
+ * specific prior written permission.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * CONTRIBUTORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS WITH
+ * THE SOFTWARE.
+ *
+ */
+
+#include <sys/time.h>
+
+typedef struct Timer{
+
+ struct timeval startTime[4];
+ struct timeval stopTime[4];
+ double time[4];
+
+}Timer;
+
+void start(Timer *timer, int i, int rep) {
+ if(rep == 0) {
+ timer->time[i] = 0.0;
+ }
+ gettimeofday(&timer->startTime[i], NULL);
+}
+
+void stop(Timer *timer, int i) {
+ gettimeofday(&timer->stopTime[i], NULL);
+ timer->time[i] += (timer->stopTime[i].tv_sec - timer->startTime[i].tv_sec) * 1000000.0 +
+ (timer->stopTime[i].tv_usec - timer->startTime[i].tv_usec);
+}
+
+void print(Timer *timer, int i, int REP) { printf("Time (ms): %f\t", timer->time[i] / (1000 * REP)); }
|