Attachment 29206 Details for Bug 23700 – example cuda source file

example cuda source file

kernel.cu (text/plain), 2.49 KB, created by Huseyin Temucin on 2009-09-04 03:01:11 UTC

(hide)

Description:

Filename:

MIME Type:

Creator: Huseyin Temucin

Created: 2009-09-04 03:01:11 UTC

Size: 2.49 KB

patch

obsolete

>// includes, system
>#include <stdlib.h>
>#include <stdio.h>
>#include <string.h>
>#include <math.h>
>
>// includes, project
>#include <cutil.h>
>
>#define lastError() printf("%s\n", cudaGetErrorString(cudaGetLastError()))
>#define NUM_THREADS 32
>
>//GPU KERNEL
>__global__ void process_array(float *array, int size){
>	int ind = blockIdx.x * blockDim.x + threadIdx.x;
>
>	if(ind < size)
>		array[ind] = pow(array[ind], (float)2) * 123.4 /44.3;
>}
>
>void write_file(int id, float *array, int size, char type){
>	FILE *fp;
>	char file_name[20];
>
>	if(type == 0)
>		strcpy(file_name, "device.txt");
>	else
>		strcpy(file_name, "host.txt");
>
>	fp = fopen(file_name, "a+");
>	fprintf(fp, "TEST : %d:\n", id);	
>
>	for(unsigned int i=0; i<size; i++)
>		fprintf(fp, "%d: %.2f\n",i, array[i]);
>	fclose(fp);
>}
>
>void init_array(float *array, int size){
>	for( unsigned int i = 0; i < size; ++i) {
>		 array[i] = (float)i;
>	}
>}
>
>void makeTest(int id, int ARRAY_SIZE){
>	float *h_array;
>	float *d_array;
>	size_t mem_size = ARRAY_SIZE * sizeof(float);
>	unsigned int timer;
>
>	printf("TEST SIZE : %d\n-------------------\n",ARRAY_SIZE);
>
>	h_array = (float *)malloc(mem_size);
>
>	init_array(h_array, ARRAY_SIZE);
>
>//	DEVICE PROCESS!!
>	CUDA_SAFE_CALL( cudaMalloc( (void**) &d_array, mem_size));
>    	// copy host memory to device
>	CUDA_SAFE_CALL( cudaMemcpy( d_array, h_array, mem_size,cudaMemcpyHostToDevice) );
>
>	int block_size = (ARRAY_SIZE/NUM_THREADS) + (ARRAY_SIZE % NUM_THREADS != 0);
>	dim3 grid(block_size, 1, 1);
>	printf("Block Size : %d\n", block_size);
>	cutCreateTimer(&timer);
>	cutStartTimer(timer);
>	process_array<<<block_size, NUM_THREADS>>>(d_array, 32);
>	cudaThreadSynchronize();
>
>	cutStopTimer(timer);
>	printf("Processing Time of Device : %.3f sec.\n", cutGetTimerValue(timer));
>	cutDeleteTimer(timer);
>	
>	CUDA_SAFE_CALL(cudaMemcpy( h_array, d_array, mem_size,cudaMemcpyDeviceToHost) );
>	write_file(id, h_array, ARRAY_SIZE, 0);
>
>//	HOST PROCESS!!
>	init_array(h_array, ARRAY_SIZE);
>
>	cutCreateTimer(&timer);
>	cutStartTimer(timer);
>	for(unsigned int i=0; i< ARRAY_SIZE; i++)
>		h_array[i] = pow(h_array[i], 2) * 123.4 /44.3;
>
>	cutStopTimer(timer);
>	printf("\nProcessing Time of Host : %.3f sec.\n", cutGetTimerValue(timer));
>	cutDeleteTimer(timer);
>	write_file(id, h_array, ARRAY_SIZE, 1);
>	
>	CUDA_SAFE_CALL( cudaFree(d_array) );
>//	free(&h_array);
>	printf("------------------------------------\n");
>}
>
>
>int main(){
>	int size = 1000;
>
>	makeTest(1, size);
>	makeTest(2, size*2);
>	makeTest(3, size*3);
>	makeTest(4,  size*4);
>	makeTest(5, size*5);
>	makeTest(6, size*10);
>
>	return 0;
>}

Actions: View

Attachments on bug 23700: 29206 | 32903

Use of freedesktop.org services, including Bugzilla, is subject to our Code of Conduct. How we collect and use information is described in our Privacy Policy.