🔨 Added scripts for automatic running with custom sizes

Makefile

@@ -29,3 +29,4 @@ clean:
 	@rm -f distributed-fractal
 	@rm -f fractal-images/distrib-fractal.jpg
 	@rm -f gpu-fractal
+	@rm -f fractal-images/cuda-fractal.jpg

code/cuda-fractal.cu

@@ -98,6 +98,12 @@ static void run(double xmin, double xmax, double ymin, double ymax, int width, i
 }
 
 int main(int argc, char** argv){
+
+	FILE *out_file = fopen("code/data/cuda.txt", "a");
+    if (out_file == NULL) {   
+		printf("Error! Could not open file\n"); 
+        exit(-1);
+    } 
     cudaError_t err = cudaSuccess;
 
     int width = 300;
@@ -109,11 +115,28 @@ int main(int argc, char** argv){
 	double ymin=-1.0;
 	double ymax= 1.0;
 
+
+	if (argc > 1)
+    	height = atoi(argv[1]);
+	if(argc>2){
+		width = atoi(argv[2]);
+	}
+	if(argc>3){
+		max_iter = atoi(argv[3]);
+	}
+	clock_t begin = clock();
+
     char *result = NULL;
     err = cudaMalloc(&result, width*height*sizeof(char));
     checkErr(err, "Failed to allocate result memory on gpu");
 
     run(xmin, xmax, ymin, ymax, width, height, max_iter, result);
+	clock_t end = clock();
+double time_spent = (double)(end - begin) / CLOCKS_PER_SEC;
+
+	fprintf(out_file, "%f & ", time_spent); // write to file 
+
+	printf("time took for execution of sequential algorithm: %f\n", time_spent);
 
     cudaFree(result);
 	return 0;

code/data/cuda.txt

@@ -0,0 +1 @@
+0.097496 & 0.130999 & 0.249084 & 0.608844 & 

code/data/seq.txt

@@ -0,0 +1 @@
+0.072927 & 0.170139 & 0.918021 & 11.324259 & 2.699279 & 10.882355 & 42.348594 & 168.265071 & 

code/distributed-fractal.c

@@ -35,8 +35,8 @@ void compute_image( double xmin, double xmax, double ymin, double ymax, int maxi
 	double xstep = (xmax-xmin) / (width-1);
 	double ystep = (ymax-ymin) / (height-1);
 
-    #pragma omp parallel shared(result, maxiter, start, end) private(i,iter)
-    #pragma omp for schedule(runtime)
+//    #pragma omp parallel shared(result, maxiter, start, end) private(i,iter)
+    #pragma omp parallel for shared(result, maxiter, start, end) private(i,iter)
     for (i = start; i < end; i++) {
 
 		double x = xmin + (i%width)*xstep;
@@ -94,6 +94,13 @@ void run(int maxiter, int current_processor, int processors_amount, int width, i
 
 int  main( int argc, char **argv ){
 
+
+	FILE *out_file = fopen("code/data/distributed.txt", "a");
+    if (out_file == NULL) {   
+		printf("Error! Could not open file\n"); 
+        exit(-1);
+    } 
+
 	double xmin=-1.5;
 	double xmax= 0.5;
 	double ymin=-1.0;
@@ -103,6 +110,14 @@ int  main( int argc, char **argv ){
     int height = 1200;
 	int maxiter=5000;
 
+	if (argc > 1)
+    	height = atoi(argv[1]);
+	if(argc>2){
+		width = atoi(argv[2]);
+	}
+	if(argc>3){
+		maxiter = atoi(argv[3]);
+	}
 
     int size, rank;
 

code/seq-fractal.c

@@ -44,22 +44,37 @@ void compute_image( double xmin, double xmax, double ymin, double ymax, int maxi
 		}
 
 	}
-     stbi_write_jpg("fractal-images/fractal.jpg", width, height, 1, buffer, 200);
 
+    stbi_write_jpg("fractal-images/fractal.jpg", width, height, 1, buffer, 200);
+
 }
 
 int main( int argc, char *argv[] )
 {
-
+	FILE *out_file = fopen("code/data/seq.txt", "a");
+    if (out_file == NULL) {   
+		printf("Error! Could not open file\n"); 
+        exit(-1);
+    } 
+
 	double xmin=-1.5;
 	double xmax= 0.5;
 	double ymin=-1.0;
 	double ymax= 1.0;
 
-    int width = 1000;
-    int height = 1000;
+    int width = 128;
+    int height = 128;
 
-	int maxiter=1000;
+	int maxiter=100;
+
+	if (argc > 1)
+    	height = atoi(argv[1]);
+	if(argc>2){
+		width = atoi(argv[2]);
+	}
+	if(argc>3){
+		maxiter = atoi(argv[3]);
+	}
 
 	printf("Timer started\n");
 
@@ -68,6 +83,8 @@ int main( int argc, char *argv[] )
 	clock_t end = clock();
 
 	double time_spent = (double)(end - begin) / CLOCKS_PER_SEC;
+
+	fprintf(out_file, "%f & ", time_spent); // write to file 
 
 	printf("time took for execution of sequential algorithm: %f\n", time_spent);
 

code/shared-fractal.c

@@ -86,6 +86,11 @@ void compute_image_opt( double xmin, double xmax, double ymin, double ymax, int
 
 int main( int argc, char *argv[] )
 {
+	FILE *out_file = fopen("code/data/shared.txt", "a");
+    if (out_file == NULL) {   
+		printf("Error! Could not open file\n"); 
+        exit(-1);
+    } 
 	double xmin=-1.5;
 	double xmax= 0.5;
 	double ymin=-1.0;
@@ -95,13 +100,18 @@ int main( int argc, char *argv[] )
     int height = 1200;
 	int threadct = omp_get_max_threads();
 	int maxiter=5000;
-
+	
 	if (argc > 1)
-    	threadct = atoi(argv[1]);
+    	height = atoi(argv[1]);
 	if(argc>2){
-		maxiter = atoi(argv[2]);
+		width = atoi(argv[2]);
+	}
+	if(argc>3){
+		maxiter = atoi(argv[3]);
 	}
 
+
+
 	printf("Coordinates: %lf %lf %lf %lf\n",xmin,xmax,ymin,ymax);
 	printf("Timer started\n");
 

run_cuda_fractals.sh

@@ -0,0 +1,13 @@
+#!/bin/bash
+make cuda
+
+./gpu-fractal 128 128 5000
+
+./gpu-fractal 256 256 5000
+
+./gpu-fractal 512 512 5000
+
+./gpu-fractal 1024 1024 5000
+
+
+

run_distributed_fractals.sh

@@ -0,0 +1,13 @@
+#!/bin/bash
+make distributed-fractal
+
+mpiexec ./distributed-fractal 128 128 5000
+
+mpiexec ./distributed-fractal 256 256 5000
+
+mpiexec ./distributed-fractal 512 512 5000
+
+mpiexec ./distributed-fractal 1024 1024 5000
+
+
+

run_seq_fractals.sh

@@ -0,0 +1,11 @@
+#!/bin/bash
+make seq-fractal
+
+./seq-fractal 128 128 5000
+
+./seq-fractal 256 256 5000
+
+./seq-fractal 512 512 5000
+
+./seq-fractal 1024 1024 5000
+

run_shared_fractals.sh

@@ -0,0 +1,13 @@
+#!/bin/bash
+make shared-fractal
+
+./shared-fractal 128 128 5000
+
+./shared-fractal 256 256 5000
+
+./shared-fractal 512 512 5000
+
+./shared-fractal 1024 1024 5000
+
+
+