Initial version

Committed on the Free edition of March Hare Software CVSNT Server. Upgrade to CVS Suite for more features and support: http://march-hare.com/cvsnt/ git-svn-id: http://moon:8086/svn/vhdl/trunk@790 cc03376c-175c-47c8-b038-4cd826a8556b
2010-03-13 14:13:53 +00:00
parent ab41a4f80f
commit 191af92d47
63 changed files with 17582 additions and 0 deletions
@@ -0,0 +1,54 @@
+// by Erik Wrenholt
+#include <stdio.h>
+#include <sys/time.h>
+
+#define BAILOUT 16
+#define MAX_ITERATIONS 1000
+
+int mandelbrot(double x, double y)
+{
+	double cr = y - 0.5;
+	double ci = x;
+	double zi = 0.0;
+	double zr = 0.0;
+	int i = 0;
+
+	while(1) {
+		i ++;
+		double temp = zr * zi;
+		double zr2 = zr * zr;
+		double zi2 = zi * zi;
+		zr = zr2 - zi2 + cr;
+		zi = temp + temp + ci;
+		if (zi2 + zr2 > BAILOUT)
+			return i;
+		if (i > MAX_ITERATIONS)
+			return 0;
+	}
+	
+}
+
+int main (int argc, const char * argv[]) {
+	struct timeval aTv;
+	gettimeofday(&aTv, NULL);
+	long init_time = aTv.tv_sec;
+	long init_usec = aTv.tv_usec;
+
+	int x,y;
+	for (y = -39; y < 39; y++) {
+		printf("\n");
+		for (x = -39; x < 39; x++) {
+			int i = mandelbrot(x/40.0, y/40.0);
+			if (i==0)
+				printf("*");
+			else
+				printf(" ");
+		}	
+	}
+	printf ("\n");
+	
+	gettimeofday(&aTv,NULL);
+	double query_time = (aTv.tv_sec - init_time) + (double)(aTv.tv_usec - init_usec)/1000000.0;	
+	printf ("C Elapsed %0.2f\n", query_time);
+    return 0;
+}