Context Navigation

main.c @ 588

Last change on this file since 588 was 502, checked in by alain, 10 years ago

1) Introduce distributed barriers in the multi-threads applications
(classif) transpose, convol, sort, gameoflife)

2) Introducing support for architectures containing empty clusters
in the mapping of these multi-threaded applications.

3) Removing the "command line arguments" in the sort application
(replaced by the giet_procs_number() system call.

File size: 25.4 KB

Rev	Line
[502]	1	///////////////////////////////////////////////////////////////////////////////////////
[334]	2	// File : main.c (for convol application)
	3	// Date : june 2014
	4	// author : Alain Greiner
[502]	5	///////////////////////////////////////////////////////////////////////////////////////
[444]	6	// This multi-threaded application application implements a 2D convolution product.
	7	// The convolution kernel is [201]*[35] pixels, but it can be factored in two
	8	// independant line and column convolution products.
	9	// It can run on a multi-processors, multi-clusters architecture, with one thread
[502]	10	// per processor.
[444]	11	//
	12	// The (1024 * 1024) pixels image is read from a file (2 bytes per pixel).
[334]	13	//
[502]	14	// - number of clusters containing processors must be power of 2 no larger than 256.
	15	// - number of processors per cluster must be power of 2 no larger than 8.
	16	///////////////////////////////////////////////////////////////////////////////////////
[334]	17
	18	#include "stdio.h"
[353]	19	#include "stdlib.h"
[502]	20	#include "user_barrier.h"
[384]	21	#include "malloc.h"
[334]	22
[502]	23	#define USE_SQT_BARRIER 1
[353]	24	#define VERBOSE 0
	25	#define SUPER_VERBOSE 0
[334]	26
[502]	27	#define X_SIZE_MAX 16
	28	#define Y_SIZE_MAX 16
	29	#define PROCS_MAX 8
	30	#define CLUSTERS_MAX (X_SIZE_MAX * Y_SIZE_MAX)
	31
[353]	32	#define INITIAL_DISPLAY_ENABLE 0
[334]	33	#define FINAL_DISPLAY_ENABLE 1
	34
	35	#define PIXEL_SIZE 2
	36	#define NL 1024
	37	#define NP 1024
	38	#define NB_PIXELS (NP * NL)
	39	#define FRAME_SIZE (NB_PIXELS * PIXEL_SIZE)
	40
	41	#define TA(c,l,p) (A[c][((NP) * (l)) + (p)])
	42	#define TB(c,p,l) (B[c][((NL) * (p)) + (l)])
	43	#define TC(c,l,p) (C[c][((NP) * (l)) + (p)])
	44	#define TD(c,l,p) (D[c][((NP) * (l)) + (p)])
	45	#define TZ(c,l,p) (Z[c][((NP) * (l)) + (p)])
	46
	47	#define max(x,y) ((x) > (y) ? (x) : (y))
	48	#define min(x,y) ((x) < (y) ? (x) : (y))
	49
[372]	50	// global instrumentation counters (cluster_id, lpid]
	51
[502]	52	unsigned int START[CLUSTERS_MAX][PROCS_MAX];
	53	unsigned int H_BEG[CLUSTERS_MAX][PROCS_MAX];
	54	unsigned int H_END[CLUSTERS_MAX][PROCS_MAX];
	55	unsigned int V_BEG[CLUSTERS_MAX][PROCS_MAX];
	56	unsigned int V_END[CLUSTERS_MAX][PROCS_MAX];
	57	unsigned int D_BEG[CLUSTERS_MAX][PROCS_MAX];
	58	unsigned int D_END[CLUSTERS_MAX][PROCS_MAX];
[334]	59
[372]	60	// global synchronization barrier
[334]	61
[502]	62	#if USE_SQT_BARRIER
	63	giet_sqt_barrier_t barrier;
[372]	64	#else
	65	giet_barrier_t barrier;
	66	#endif
[334]	67
[372]	68	volatile unsigned int barrier_init_ok = 0;
	69	volatile unsigned int load_image_ok = 0;
	70	volatile unsigned int instrumentation_ok = 0;
[334]	71
[372]	72	// global pointers on distributed buffers in all clusters
[502]	73	unsigned short * GA[CLUSTERS_MAX];
	74	int * GB[CLUSTERS_MAX];
	75	int * GC[CLUSTERS_MAX];
	76	int * GD[CLUSTERS_MAX];
	77	unsigned char * GZ[CLUSTERS_MAX];
[372]	78
[334]	79	///////////////////////////////////////////
	80	__attribute__ ((constructor)) void main()
	81	///////////////////////////////////////////
	82	{
	83	//////////////////////////////////
	84	// convolution kernel parameters
	85	// The content of this section is
	86	// Philips proprietary information.
	87	///////////////////////////////////
	88
	89	int vnorm = 115;
	90	int vf[35] = { 1, 1, 2, 2, 2,
	91	2, 3, 3, 3, 4,
	92	4, 4, 4, 5, 5,
	93	5, 5, 5, 5, 5,
	94	5, 5, 4, 4, 4,
	95	4, 3, 3, 3, 2,
	96	2, 2, 2, 1, 1 };
	97
	98	int hrange = 100;
	99	int hnorm = 201;
	100
	101	unsigned int date = 0;
	102
	103	int c; // cluster index for loops
	104	int l; // line index for loops
	105	int p; // pixel index for loops
	106	int z; // vertical filter index for loops
	107
[502]	108	// plat-form parameters
	109	unsigned int x_size; // number of clusters in a row
	110	unsigned int y_size; // number of clusters in a column
	111	unsigned int nprocs; // number of processors per cluster
	112
	113	giet_procs_number( &x_size , &y_size , &nprocs );
	114
[432]	115	// processor identifiers
[502]	116	unsigned int x; // x coordinate
	117	unsigned int y; // y coordinate
	118	unsigned int lpid; // local proc/task id
[432]	119	giet_proc_xyp( &x, &y, &lpid );
	120
[502]	121	int file = 0; // file descriptor
	122	unsigned int nclusters = x_size * y_size; // number of clusters
	123	unsigned int cluster_id = (x * y_size) + y; // continuous cluster index
	124	unsigned int task_id = (cluster_id * nprocs) + lpid; // continuous task index
	125	unsigned int ntasks = nclusters * nprocs; // number of tasks
	126	unsigned int frame_size = FRAME_SIZE; // total size (bytes)
	127	unsigned int nblocks = frame_size / 512; // number of blocks/frame
[334]	128
[502]	129	unsigned int lines_per_task = NL / ntasks; // lines per task
	130	unsigned int lines_per_cluster = NL / nclusters; // lines per cluster
	131	unsigned int pixels_per_task = NP / ntasks; // columns per task
	132	unsigned int pixels_per_cluster = NP / nclusters; // columns per cluster
[334]	133
	134	int first, last;
	135
	136	date = giet_proctime();
	137	START[cluster_id][lpid] = date;
	138
[353]	139	#if VERBOSE
	140	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] starts at cycle %d\n", x,y,lpid, date );
	141	#endif
	142
[334]	143	// parameters checking
	144
[502]	145	if ((nprocs != 1) && (nprocs != 2) && (nprocs != 4) && (nprocs != 8))
	146	giet_exit( "[CONVOL ERROR] NB_PROCS_MAX must be 1, 2, 4 or 8\n");
[334]	147
[502]	148	if ((x_size!=1) && (x_size!=2) && (x_size!=4) && (x_size!=8) && (x_size!=16))
	149	giet_exit( "[CONVOL ERROR] x_size must be 1, 2, 4, 8, 16\n");
[334]	150
[502]	151	if ((y_size!=1) && (y_size!=2) && (y_size!=4) && (y_size!=8) && (y_size!=16))
	152	giet_exit( "[CONVOL ERROR] y_size must be 1, 2, 4, 8, 16\n");
[334]	153
	154	if ( NL % nclusters != 0 )
[502]	155	giet_exit( "[CONVOL ERROR] CLUSTERS_MAX must be a divider of NL");
[334]	156
	157	if ( NP % nclusters != 0 )
[502]	158	giet_exit( "[CONVOL ERROR] CLUSTERS_MAX must be a divider of NP");
[334]	159
[372]	160
[334]	161	///////////////////////////////////////////////////////////////////
[372]	162	// task[0][0][0] makes barrier initialisation
[334]	163	///////////////////////////////////////////////////////////////////
[372]	164
[432]	165	if ( (x==0) && (y==0) && (lpid==0) )
[372]	166	{
	167	giet_shr_printf("\n[CONVOL] task[0,0,0] starts barrier init at cycle %d\n"
[502]	168	"- CLUSTERS = %d\n"
	169	"- PROCS = %d\n"
	170	"- TASKS = %d\n"
	171	"- BLOCKS = %d\n",
[444]	172	giet_proctime(), nclusters, nprocs, ntasks, nblocks );
[502]	173	#if USE_SQT_BARRIER
	174	sqt_barrier_init( &barrier, x_size , y_size , nprocs );
[372]	175	#else
	176	barrier_init( &barrier, ntasks );
	177	#endif
[334]	178
[372]	179	giet_shr_printf( "\n[CONVOL] task[0,0,0] completes barrier init at cycle %d\n",
	180	giet_proctime() );
	181
	182	barrier_init_ok = 1;
	183	}
	184	else
	185	{
	186	while ( barrier_init_ok == 0 );
	187	}
	188
	189	///////////////////////////////////////////////////////////////////
	190	// All task[x][y][0] allocate the global buffers in cluster(x,y)
	191	// These buffers mut be sector-aligned.
	192	///////////////////////////////////////////////////////////////////
	193	if ( lpid == 0 )
	194	{
	195
	196	#if VERBOSE
[502]	197	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] enters malloc at cycle %d\n",
	198	x,y,lpid, date );
[372]	199	#endif
	200
[384]	201	GA[cluster_id] = remote_malloc( (FRAME_SIZE/nclusters) , x , y );
	202	GB[cluster_id] = remote_malloc( (FRAME_SIZE/nclusters)*2 , x , y );
	203	GC[cluster_id] = remote_malloc( (FRAME_SIZE/nclusters)*2 , x , y );
	204	GD[cluster_id] = remote_malloc( (FRAME_SIZE/nclusters)*2 , x , y );
	205	GZ[cluster_id] = remote_malloc( (FRAME_SIZE/nclusters)/2 , x , y );
[372]	206
	207	#if VERBOSE
	208	giet_shr_printf( "\n[CONVOL] Shared Buffer Virtual Addresses in cluster(%d,%d)\n"
	209	"### GA = %x\n"
	210	"### GB = %x\n"
	211	"### GC = %x\n"
	212	"### GD = %x\n"
	213	"### GZ = %x\n",
	214	x, y,
	215	GA[cluster_id],
	216	GB[cluster_id],
	217	GC[cluster_id],
	218	GD[cluster_id],
	219	GZ[cluster_id] );
	220	#endif
	221	}
	222
[377]	223	///////////////////////////////
[502]	224	#if USE_SQT_BARRIER
	225	sqt_barrier_wait( &barrier );
[377]	226	#else
	227	barrier_wait( &barrier );
	228	#endif
[372]	229
	230	///////////////////////////////////////////////////////////////////
	231	// All tasks initialise in their private stack a copy of the
	232	// arrays of pointers on the shared, distributed buffers.
	233	///////////////////////////////////////////////////////////////////
	234
[502]	235	unsigned short * A[CLUSTERS_MAX];
	236	int * B[CLUSTERS_MAX];
	237	int * C[CLUSTERS_MAX];
	238	int * D[CLUSTERS_MAX];
	239	unsigned char * Z[CLUSTERS_MAX];
[334]	240
	241	for (c = 0; c < nclusters; c++)
	242	{
[372]	243	A[c] = GA[c];
	244	B[c] = GB[c];
	245	C[c] = GC[c];
	246	D[c] = GD[c];
	247	Z[c] = GZ[c];
[334]	248	}
	249
	250	///////////////////////////////////////////////////////////////////////////
[372]	251	// task[0,0,0] open the file containing image, and load it from disk
	252	// to all A[c] buffers (nblocks / nclusters loaded in each cluster).
[334]	253	// Other tasks are waiting on the init_ok condition.
	254	//////////////////////////////////////////////////////////////////////////
[432]	255	if ( (x==0) && (y==0) && (lpid==0) )
[334]	256	{
	257	// open file
	258	file = giet_fat_open("misc/philips_image.raw", 0 );
	259	if ( file < 0 ) giet_exit( "[CONVOL ERROR] task[0,0,0] cannot open"
	260	" file misc/philips_image.raw" );
	261
	262	giet_shr_printf( "\n[CONVOL] task[0,0,0] open file misc/philips_image.raw"
	263	" at cycle %d\n", giet_proctime() );
	264
[502]	265	for ( c = 0 ; c < nclusters ; c++ )
[334]	266	{
[353]	267	giet_shr_printf( "\n[CONVOL] task[0,0,0] starts load "
	268	"for cluster %d at cycle %d\n", c, giet_proctime() );
	269
[334]	270	giet_fat_read( file,
	271	A[c],
	272	nblocks/nclusters,
	273	(nblocks/nclusters)*c );
	274
[353]	275	giet_shr_printf( "\n[CONVOL] task[0,0,0] completes load "
	276	"for cluster %d at cycle %d\n", c, giet_proctime() );
[334]	277	}
[372]	278	load_image_ok = 1;
[334]	279	}
	280	else
	281	{
[372]	282	while ( load_image_ok == 0 );
[334]	283	}
	284
	285	/////////////////////////////////////////////////////////////////////////////
	286	// Optionnal parallel display of the initial image stored in A[c] buffers.
	287	// Eah task displays (NL/ntasks) lines. (one byte per pixel).
	288	/////////////////////////////////////////////////////////////////////////////
	289
	290	if ( INITIAL_DISPLAY_ENABLE )
	291	{
	292
[353]	293	#if VERBOSE
	294	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] starts initial display"
	295	" at cycle %d\n",
	296	x, y, lpid, giet_proctime() );
	297	#endif
	298
[334]	299	unsigned int line;
	300	unsigned int offset = lines_per_task * lpid;
	301
	302	for ( l = 0 ; l < lines_per_task ; l++ )
	303	{
	304	line = offset + l;
	305
	306	for ( p = 0 ; p < NP ; p++ )
	307	{
	308	TZ(cluster_id, line, p) = (unsigned char)(TA(cluster_id, line, p) >> 8);
	309	}
	310
[444]	311	giet_fbf_sync_write( NP(l + (task_id lines_per_task) ),
	312	&TZ(cluster_id, line, 0),
	313	NP);
[334]	314	}
	315
[353]	316	#if VERBOSE
	317	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] completes initial display"
	318	" at cycle %d\n",
	319	x, y, lpid, giet_proctime() );
	320	#endif
[334]	321
[377]	322	////////////////////////////
[502]	323	#if USE_SQT_BARRIER
	324	sqt_barrier_wait( &barrier );
[377]	325	#else
[372]	326	barrier_wait( &barrier );
[377]	327	#endif
[334]	328
	329	}
	330
	331	////////////////////////////////////////////////////////
	332	// parallel horizontal filter :
	333	// B <= transpose(FH(A))
	334	// D <= A - FH(A)
	335	// Each task computes (NL/ntasks) lines
	336	// The image must be extended :
	337	// if (z<0) TA(cluster_id,l,z) == TA(cluster_id,l,0)
	338	// if (z>NP-1) TA(cluster_id,l,z) == TA(cluster_id,l,NP-1)
	339	////////////////////////////////////////////////////////
	340
	341	date = giet_proctime();
	342	H_BEG[cluster_id][lpid] = date;
	343
[353]	344	#if VERBOSE
	345	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] starts horizontal filter"
	346	" at cycle %d\n",
	347	x, y, lpid, date );
	348	#else
[432]	349	if ( (x==0) && (y==0) && (lpid==0) )
	350	giet_shr_printf( "\n[CONVOL] task[0,0,0] starts horizontal filter"
	351	" at cycle %d\n", date );
[353]	352	#endif
	353
[334]	354	// l = absolute line index / p = absolute pixel index
	355	// first & last define which lines are handled by a given task
	356
	357	first = task_id * lines_per_task;
	358	last = first + lines_per_task;
	359
	360	for (l = first; l < last; l++)
	361	{
	362	// src_c and src_l are the cluster index and the line index for A & D
	363	int src_c = l / lines_per_cluster;
	364	int src_l = l % lines_per_cluster;
	365
	366	// We use the specific values of the horizontal ep-filter for optimisation:
	367	// sum(p) = sum(p-1) + TA[p+hrange] - TA[p-hrange-1]
	368	// To minimize the number of tests, the loop on pixels is split in three domains
	369
	370	int sum_p = (hrange + 2) * TA(src_c, src_l, 0);
	371	for (z = 1; z < hrange; z++)
	372	{
	373	sum_p = sum_p + TA(src_c, src_l, z);
	374	}
	375
	376	// first domain : from 0 to hrange
	377	for (p = 0; p < hrange + 1; p++)
	378	{
	379	// dst_c and dst_p are the cluster index and the pixel index for B
	380	int dst_c = p / pixels_per_cluster;
	381	int dst_p = p % pixels_per_cluster;
	382	sum_p = sum_p + (int) TA(src_c, src_l, p + hrange) - (int) TA(src_c, src_l, 0);
	383	TB(dst_c, dst_p, l) = sum_p / hnorm;
	384	TD(src_c, src_l, p) = (int) TA(src_c, src_l, p) - sum_p / hnorm;
	385	}
	386	// second domain : from (hrange+1) to (NP-hrange-1)
	387	for (p = hrange + 1; p < NP - hrange; p++)
	388	{
	389	// dst_c and dst_p are the cluster index and the pixel index for B
	390	int dst_c = p / pixels_per_cluster;
	391	int dst_p = p % pixels_per_cluster;
	392	sum_p = sum_p + (int) TA(src_c, src_l, p + hrange)
	393	- (int) TA(src_c, src_l, p - hrange - 1);
	394	TB(dst_c, dst_p, l) = sum_p / hnorm;
	395	TD(src_c, src_l, p) = (int) TA(src_c, src_l, p) - sum_p / hnorm;
	396	}
	397	// third domain : from (NP-hrange) to (NP-1)
	398	for (p = NP - hrange; p < NP; p++)
	399	{
	400	// dst_c and dst_p are the cluster index and the pixel index for B
	401	int dst_c = p / pixels_per_cluster;
	402	int dst_p = p % pixels_per_cluster;
	403	sum_p = sum_p + (int) TA(src_c, src_l, NP - 1)
	404	- (int) TA(src_c, src_l, p - hrange - 1);
	405	TB(dst_c, dst_p, l) = sum_p / hnorm;
	406	TD(src_c, src_l, p) = (int) TA(src_c, src_l, p) - sum_p / hnorm;
	407	}
	408
[353]	409	#if SUPER_VERBOSE
	410	giet_shr_printf(" - line %d computed at cycle %d\n", l, giet_proctime() );
	411	#endif
	412
[334]	413	}
	414
	415	date = giet_proctime();
	416	H_END[cluster_id][lpid] = date;
	417
[353]	418	#if VERBOSE
	419	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] completes horizontal filter"
	420	" at cycle %d\n",
	421	x, y, lpid, date );
	422	#else
[432]	423	if ( (x==0) && (y==0) && (lpid==0) )
	424	giet_shr_printf( "\n[CONVOL] task[0,0,0] completes horizontal filter"
	425	" at cycle %d\n", date );
[353]	426	#endif
	427
[377]	428	/////////////////////////////
[502]	429	#if USE_SQT_BARRIER
	430	sqt_barrier_wait( &barrier );
[377]	431	#else
	432	barrier_wait( &barrier );
	433	#endif
[334]	434
[372]	435
[353]	436	///////////////////////////////////////////////////////////////
[334]	437	// parallel vertical filter :
	438	// C <= transpose(FV(B))
	439	// Each task computes (NP/ntasks) columns
	440	// The image must be extended :
	441	// if (l<0) TB(cluster_id,p,l) == TB(cluster_id,p,0)
	442	// if (l>NL-1) TB(cluster_id,p,l) == TB(cluster_id,p,NL-1)
[353]	443	///////////////////////////////////////////////////////////////
[334]	444
	445	date = giet_proctime();
	446	V_BEG[cluster_id][lpid] = date;
	447
[353]	448	#if VERBOSE
	449	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] starts vertical filter"
	450	" at cycle %d\n",
	451	x, y, lpid, date );
	452	#else
[432]	453	if ( (x==0) && (y==0) && (lpid==0) )
	454	giet_shr_printf( "\n[CONVOL] task[0,0,0] starts vertical filter"
	455	" at cycle %d\n", date );
[353]	456	#endif
	457
[334]	458	// l = absolute line index / p = absolute pixel index
	459	// first & last define which pixels are handled by a given task
	460
	461	first = task_id * pixels_per_task;
	462	last = first + pixels_per_task;
	463
	464	for (p = first; p < last; p++)
	465	{
	466	// src_c and src_p are the cluster index and the pixel index for B
	467	int src_c = p / pixels_per_cluster;
	468	int src_p = p % pixels_per_cluster;
	469
	470	int sum_l;
	471
	472	// We use the specific values of the vertical ep-filter
	473	// To minimize the number of tests, the NL lines are split in three domains
	474
	475	// first domain : explicit computation for the first 18 values
	476	for (l = 0; l < 18; l++)
	477	{
	478	// dst_c and dst_l are the cluster index and the line index for C
	479	int dst_c = l / lines_per_cluster;
	480	int dst_l = l % lines_per_cluster;
	481
	482	for (z = 0, sum_l = 0; z < 35; z++)
	483	{
	484	sum_l = sum_l + vf[z] * TB(src_c, src_p, max(l - 17 + z,0) );
	485	}
	486	TC(dst_c, dst_l, p) = sum_l / vnorm;
	487	}
	488	// second domain
	489	for (l = 18; l < NL - 17; l++)
	490	{
	491	// dst_c and dst_l are the cluster index and the line index for C
	492	int dst_c = l / lines_per_cluster;
	493	int dst_l = l % lines_per_cluster;
	494
	495	sum_l = sum_l + TB(src_c, src_p, l + 4)
	496	+ TB(src_c, src_p, l + 8)
	497	+ TB(src_c, src_p, l + 11)
	498	+ TB(src_c, src_p, l + 15)
	499	+ TB(src_c, src_p, l + 17)
	500	- TB(src_c, src_p, l - 5)
	501	- TB(src_c, src_p, l - 9)
	502	- TB(src_c, src_p, l - 12)
	503	- TB(src_c, src_p, l - 16)
	504	- TB(src_c, src_p, l - 18);
	505
	506	TC(dst_c, dst_l, p) = sum_l / vnorm;
	507	}
	508	// third domain
	509	for (l = NL - 17; l < NL; l++)
	510	{
	511	// dst_c and dst_l are the cluster index and the line index for C
	512	int dst_c = l / lines_per_cluster;
	513	int dst_l = l % lines_per_cluster;
	514
	515	sum_l = sum_l + TB(src_c, src_p, min(l + 4, NL - 1))
	516	+ TB(src_c, src_p, min(l + 8, NL - 1))
	517	+ TB(src_c, src_p, min(l + 11, NL - 1))
	518	+ TB(src_c, src_p, min(l + 15, NL - 1))
	519	+ TB(src_c, src_p, min(l + 17, NL - 1))
	520	- TB(src_c, src_p, l - 5)
	521	- TB(src_c, src_p, l - 9)
	522	- TB(src_c, src_p, l - 12)
	523	- TB(src_c, src_p, l - 16)
	524	- TB(src_c, src_p, l - 18);
	525
	526	TC(dst_c, dst_l, p) = sum_l / vnorm;
	527	}
	528
[353]	529	#if SUPER_VERBOSE
	530	giet_shr_printf(" - column %d computed at cycle %d\n", p, giet_proctime());
	531	#endif
	532
[334]	533	}
	534
	535	date = giet_proctime();
	536	V_END[cluster_id][lpid] = date;
	537
[353]	538	#if VERBOSE
	539	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] completes vertical filter"
	540	" at cycle %d\n",
	541	x, y, lpid, date );
	542	#else
[432]	543	if ( (x==0) && (y==0) && (lpid==0) )
	544	giet_shr_printf( "\n[CONVOL] task[0,0,0] completes vertical filter"
	545	" at cycle %d\n", date );
[353]	546	#endif
	547
[377]	548	////////////////////////////
[502]	549	#if USE_SQT_BARRIER
	550	sqt_barrier_wait( &barrier );
[377]	551	#else
	552	barrier_wait( &barrier );
	553	#endif
[334]	554
	555	////////////////////////////////////////////////////////////////
	556	// Optional parallel display of the final image Z <= D + C
	557	// Eah task displays (NL/ntasks) lines. (one byte per pixel).
	558	////////////////////////////////////////////////////////////////
	559
	560	if ( FINAL_DISPLAY_ENABLE )
	561	{
	562	date = giet_proctime();
	563	D_BEG[cluster_id][lpid] = date;
	564
[353]	565	#if VERBOSE
	566	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] starts final display"
	567	" at cycle %d\n",
	568	x, y, lpid, date);
	569	#else
[432]	570	if ( (x==0) && (y==0) && (lpid==0) )
	571	giet_shr_printf( "\n[CONVOL] task[0,0,0] starts final display"
	572	" at cycle %d\n", date );
[353]	573	#endif
	574
[334]	575	unsigned int line;
	576	unsigned int offset = lines_per_task * lpid;
	577
	578	for ( l = 0 ; l < lines_per_task ; l++ )
	579	{
	580	line = offset + l;
	581
	582	for ( p = 0 ; p < NP ; p++ )
	583	{
	584	TZ(cluster_id, line, p) =
	585	(unsigned char)( (TD(cluster_id, line, p) +
	586	TC(cluster_id, line, p) ) >> 8 );
	587	}
	588
[444]	589	giet_fbf_sync_write( NP(l + (task_id lines_per_task) ),
	590	&TZ(cluster_id, line, 0),
	591	NP);
[334]	592	}
	593
	594	date = giet_proctime();
	595	D_END[cluster_id][lpid] = date;
[353]	596
	597	#if VERBOSE
	598	giet_shr_printf( "\n[CONVOL] task[%d,%d,%d] completes final display"
	599	" at cycle %d\n",
	600	x, y, lpid, date);
	601	#else
[432]	602	if ( (x==0) && (y==0) && (lpid==0) )
	603	giet_shr_printf( "\n[CONVOL] task[0,0,0] completes final display"
	604	" at cycle %d\n", date );
[353]	605	#endif
[334]	606
[377]	607	//////////////////////////////
[502]	608	#if USE_SQT_BARRIER
	609	sqt_barrier_wait( &barrier );
[377]	610	#else
	611	barrier_wait( &barrier );
	612	#endif
[372]	613
[334]	614	}
	615
	616	/////////////////////////////////////////////////////////
	617	// Task[0,0,0] makes the instrumentation
	618	/////////////////////////////////////////////////////////
	619
[432]	620	if ( (x==0) && (y==0) && (lpid==0) )
[334]	621	{
	622	date = giet_proctime();
[353]	623	giet_shr_printf("\n[CONVOL] task[0,0,0] starts instrumentation"
	624	" at cycle %d\n\n", date );
[334]	625
	626	int cc, pp;
	627
	628	unsigned int min_start = 0xFFFFFFFF;
	629	unsigned int max_start = 0;
	630
	631	unsigned int min_h_beg = 0xFFFFFFFF;
	632	unsigned int max_h_beg = 0;
	633
	634	unsigned int min_h_end = 0xFFFFFFFF;
	635	unsigned int max_h_end = 0;
	636
	637	unsigned int min_v_beg = 0xFFFFFFFF;
	638	unsigned int max_v_beg = 0;
	639
	640	unsigned int min_v_end = 0xFFFFFFFF;
	641	unsigned int max_v_end = 0;
	642
	643	unsigned int min_d_beg = 0xFFFFFFFF;
	644	unsigned int max_d_beg = 0;
	645
	646	unsigned int min_d_end = 0xFFFFFFFF;
	647	unsigned int max_d_end = 0;
	648
	649	for (cc = 0; cc < nclusters; cc++)
	650	{
	651	for (pp = 0; pp < nprocs; pp++ )
	652	{
	653	if (START[cc][pp] < min_start) min_start = START[cc][pp];
	654	if (START[cc][pp] > max_start) max_start = START[cc][pp];
	655
	656	if (H_BEG[cc][pp] < min_h_beg) min_h_beg = H_BEG[cc][pp];
	657	if (H_BEG[cc][pp] > max_h_beg) max_h_beg = H_BEG[cc][pp];
	658
	659	if (H_END[cc][pp] < min_h_end) min_h_end = H_END[cc][pp];
	660	if (H_END[cc][pp] > max_h_end) max_h_end = H_END[cc][pp];
	661
	662	if (V_BEG[cc][pp] < min_v_beg) min_v_beg = V_BEG[cc][pp];
	663	if (V_BEG[cc][pp] > max_v_beg) max_v_beg = V_BEG[cc][pp];
	664
	665	if (V_END[cc][pp] < min_v_end) min_v_end = V_END[cc][pp];
	666	if (V_END[cc][pp] > max_v_end) max_v_end = V_END[cc][pp];
	667
	668	if (D_BEG[cc][pp] < min_d_beg) min_d_beg = D_BEG[cc][pp];
	669	if (D_BEG[cc][pp] > max_d_beg) max_d_beg = D_BEG[cc][pp];
	670
	671	if (D_END[cc][pp] < min_d_end) min_d_end = D_END[cc][pp];
	672	if (D_END[cc][pp] > max_d_end) max_d_end = D_END[cc][pp];
	673	}
	674	}
	675
	676	giet_shr_printf(" - START : min = %d / max = %d / med = %d / delta = %d\n",
	677	min_start, max_start, (min_start+max_start)/2, max_start-min_start);
	678
	679	giet_shr_printf(" - H_BEG : min = %d / max = %d / med = %d / delta = %d\n",
	680	min_h_beg, max_h_beg, (min_h_beg+max_h_beg)/2, max_h_beg-min_h_beg);
	681
[353]	682	giet_shr_printf(" - H_END : min = %d / max = %d / med = %d / delta = %d\n",
[334]	683	min_h_end, max_h_end, (min_h_end+max_h_end)/2, max_h_end-min_h_end);
	684
	685	giet_shr_printf(" - V_BEG : min = %d / max = %d / med = %d / delta = %d\n",
	686	min_v_beg, max_v_beg, (min_v_beg+max_v_beg)/2, max_v_beg-min_v_beg);
	687
[353]	688	giet_shr_printf(" - V_END : min = %d / max = %d / med = %d / delta = %d\n",
[334]	689	min_v_end, max_v_end, (min_v_end+max_v_end)/2, max_v_end-min_v_end);
	690
	691	giet_shr_printf(" - D_BEG : min = %d / max = %d / med = %d / delta = %d\n",
	692	min_d_beg, max_d_beg, (min_d_beg+max_d_beg)/2, max_d_beg-min_d_beg);
	693
[353]	694	giet_shr_printf(" - D_END : min = %d / max = %d / med = %d / delta = %d\n",
[334]	695	min_d_end, max_d_end, (min_d_end+max_d_end)/2, max_d_end-min_d_end);
	696
	697	giet_shr_printf( "\n General Scenario (Kcycles for each step)\n" );
[353]	698	giet_shr_printf( " - BOOT OS = %d\n", (min_start )/1000 );
	699	giet_shr_printf( " - LOAD IMAGE = %d\n", (min_h_beg - min_start)/1000 );
	700	giet_shr_printf( " - H_FILTER = %d\n", (max_h_end - min_h_beg)/1000 );
	701	giet_shr_printf( " - BARRIER HORI/VERT = %d\n", (min_v_beg - max_h_end)/1000 );
	702	giet_shr_printf( " - V_FILTER = %d\n", (max_v_end - min_v_beg)/1000 );
	703	giet_shr_printf( " - BARRIER VERT/DISP = %d\n", (min_d_beg - max_v_end)/1000 );
	704	giet_shr_printf( " - DISPLAY = %d\n", (max_d_end - min_d_beg)/1000 );
[334]	705
[372]	706	instrumentation_ok = 1;
[334]	707	}
	708	else
	709	{
[372]	710	while ( instrumentation_ok == 0 );
[334]	711	}
	712
	713	giet_exit( "completed");
	714
	715	} // end main()
	716
	717	// Local Variables:
	718	// tab-width: 3
	719	// c-basic-offset: 3
	720	// c-file-offsets:((innamespace . 0)(inline-open . 0))
	721	// indent-tabs-mode: nil
	722	// End:
	723
	724	// vim: filetype=cpp:expandtab:shiftwidth=3:tabstop=3:softtabstop=3
	725
	726

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: soft/giet_vm/applications/convol/main.c @ 588

Download in other formats: