starpu-max/tests/microbenchs/parallel_independent_homogeneous_tasks.c

/* StarPU --- Runtime system for heterogeneous multicore architectures.
 *
 * Copyright (C) 2016  Bérangère Subervie
 *
 * StarPU is free software; you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation; either version 2.1 of the License, or (at
 * your option) any later version.
 *
 * StarPU is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 *
 * See the GNU Lesser General Public License in COPYING.LGPL for more details.
 */

#include <stdbool.h>
#include <starpu.h>
#include "../helper.h"

#define TIME 0.010
#define TASK_COEFFICIENT 100
#define MARGIN 0.05
#define SECONDS_SCALE_COEFFICIENT_TIMING_NOW 1000000

void wait(void *descr[] STARPU_ATTRIBUTE_UNUSED, void *_args){
	starpu_sleep(TIME);
}


double cost_function(struct starpu_task *t, struct starpu_perfmodel_arch *a, unsigned i)
{
	return TIME * 1000000;
}


static struct starpu_perfmodel perf_model =
{
	.type = STARPU_PER_ARCH,
	.arch_cost_function = cost_function,
};

static struct starpu_codelet cl =
{
	.cpu_funcs = { wait },
	.cuda_funcs = { wait }, 
	.opencl_funcs = { wait },
	.cpu_funcs_name = { "wait" },
	.nbuffers = 0,
	.flags = STARPU_CODELET_SIMGRID_EXECUTE,
	.model = &perf_model,
};

int main(int argc, char *argv[]){
	
	int ret;
	
	ret = starpu_initialize(NULL, &argc, &argv); 
	if (ret == -ENODEV) return STARPU_TEST_SKIPPED;
	STARPU_CHECK_RETURN_VALUE(ret, "starpu_init");

	unsigned nb_tasks, nb_workers, i;
	double begin_time, end_time, time_m, time_s, speed_up, expected_speed_up, percentage_expected_speed_up;
	bool check, check_sup;
	
	nb_workers = starpu_worker_get_count_by_type(STARPU_CPU_WORKER) + starpu_worker_get_count_by_type(STARPU_CUDA_WORKER) + starpu_worker_get_count_by_type(STARPU_OPENCL_WORKER);
	nb_tasks = nb_workers*TASK_COEFFICIENT;
	
	begin_time = starpu_timing_now();

	/*execution des tasks*/
	
	for (i=0; i<nb_tasks; i++){
		starpu_task_insert(&cl,0);
	}

	starpu_task_wait_for_all();	
	
	end_time = starpu_timing_now();

	/*on determine si le temps mesure est satisfaisant ou pas*/

	time_m = (end_time - begin_time)/SECONDS_SCALE_COEFFICIENT_TIMING_NOW; //pour ramener en secondes
	time_s = nb_tasks * TIME;
	speed_up = time_s/time_m;
	expected_speed_up =  nb_workers;
	percentage_expected_speed_up = 100 * (speed_up/expected_speed_up);
	check = speed_up >= ((1 - MARGIN) * expected_speed_up);
	check_sup = speed_up <= ((1 + MARGIN) * expected_speed_up);

	printf("measured time = %f seconds\nsequential time = %f seconds\nspeed up = %f\nnumber of workers = %d\nnumber of tasks = %d\nexpected speed up = %f\npercentage of expected speed up = %.2f%%\n", time_m, time_s, speed_up, nb_workers, nb_tasks, expected_speed_up, percentage_expected_speed_up);

	starpu_shutdown();
	
	if (check && check_sup){ //test reussi ou test echoue
		return EXIT_SUCCESS;
	}
	else{
		return EXIT_FAILURE;
	}
}