From: hackbard <hackbard@sage.physik.uni-augsburg.de>
Date: Mon, 13 Aug 2007 16:32:45 +0000 (+0200)
Subject: visualization stuff
X-Git-Url: https://hackdaworld.org/gitweb/?a=commitdiff_plain;h=e67bb08cb64bf21a6d4e7e8275007bdc5e769154;p=physik%2Fposic.git

visualization stuff
---

diff --git a/ppm2avi b/ppm2avi
index 6987317..3680b80 100755
--- a/ppm2avi
+++ b/ppm2avi
@@ -2,4 +2,7 @@
 #	convert $i $1/$(basename $i ppm)png
 #done
 
-mencoder mf://$1/*.png -ovc copy -o $1/md.avi >/dev/null 2>&1
+OUTFILE=md.avi
+[ -n "$2" ] && OUTFILE=$2
+
+mencoder mf://$1/*.png -ovc copy -o $1/$OUTFILE >/dev/null 2>&1
diff --git a/sic.c b/sic.c
index bf33e16..be75aad 100644
--- a/sic.c
+++ b/sic.c
@@ -26,8 +26,8 @@
 #define R_C		1.0
 #define T_C		10.0
 #define LCNT		5
-#define PRERUN		200
-#define POSTRUN		500
+#define PRERUN		500
+#define POSTRUN		1000
 
 typedef struct s_hp {
 	int a_count;	/* atom count */
@@ -361,8 +361,8 @@ int main(int argc,char **argv) {
 	moldyn_set_log(&md,LOG_TOTAL_ENERGY,1);
 	moldyn_set_log(&md,LOG_TEMPERATURE,1);
 	moldyn_set_log(&md,LOG_PRESSURE,1);
-	moldyn_set_log(&md,VISUAL_STEP,100);
-	moldyn_set_log(&md,SAVE_STEP,100);
+	moldyn_set_log(&md,VISUAL_STEP,10);
+	moldyn_set_log(&md,SAVE_STEP,10);
 	moldyn_set_log(&md,CREATE_REPORT,0);
 
 	/*
diff --git a/visualize b/visualize
index f24e7b7..6255e89 100755
--- a/visualize
+++ b/visualize
@@ -5,6 +5,18 @@
 # author: frank.zirkelbach@physik.uni-augsburg.de
 #
 
+# help function
+draw_cyl() {
+	cat >> temp.pov <<-EOF
+cylinder {
+<$1, $3, $2>, <$4, $6, $5>, 0.05
+pigment { color White }
+}
+EOF
+}
+
+	
+
 directory="doesnt_exist____for_sure"
 width="1024"
 height="768"
@@ -45,7 +57,7 @@ while [ "$1" ]; do
 				exit 1;;
 	esac
 done
-		
+
 if [ ! -d $directory ] ; then
 	echo "no valid directory"
 	exit 1
@@ -76,8 +88,8 @@ EOF
 		export x0 y0 z0 x1 y1 z1 radius
 		cat $file | grep -v '#' | awk '\
 		BEGIN {
-			x0=ENVIRON["x0"]; y0==ENVIRON["y0"]; z0==ENVIRON["z0"];
-			x1=ENVIRON["x1"]; y0==ENVIRON["y1"]; z0==ENVIRON["z1"];
+			x0=ENVIRON["x0"]; y0=ENVIRON["y0"]; z0=ENVIRON["z0"];
+			x1=ENVIRON["x1"]; y1=ENVIRON["y1"]; z1=ENVIRON["z1"];
 			radius=ENVIRON["radius"];
 		}
 		{
@@ -85,7 +97,7 @@ EOF
 			   ($2<=x1)&&($3<=y1)&&($4<=z1)) {
 				print "sphere { <"$2","$4","$3">, "radius" ";
 				print "texture { pigment { color "$5" } ";
-				print "finish { phong 1, metallic } } }";
+				print "finish { phong 1 metallic } } }";
 			}
 		}' >> temp.pov
 	else
@@ -115,6 +127,22 @@ pigment { color White }
 }
 EOF
 	done
+	else
+		# manually drawing the 3x4 boundaries ...
+		draw_cyl $x0 $y0 $z0 $x1 $y0 $z0
+		draw_cyl $x0 $y0 $z0 $x0 $y1 $z0
+		draw_cyl $x1 $y1 $z0 $x1 $y0 $z0
+		draw_cyl $x0 $y1 $z0 $x1 $y1 $z0
+
+		draw_cyl $x0 $y0 $z1 $x1 $y0 $z1
+		draw_cyl $x0 $y0 $z1 $x0 $y1 $z1
+		draw_cyl $x1 $y1 $z1 $x1 $y0 $z1
+		draw_cyl $x0 $y1 $z1 $x1 $y1 $z1
+
+		draw_cyl $x0 $y0 $z1 $x0 $y0 $z0
+		draw_cyl $x0 $y1 $z1 $x0 $y1 $z0
+		draw_cyl $x1 $y0 $z1 $x1 $y0 $z0
+		draw_cyl $x1 $y1 $z1 $x1 $y1 $z0
 	fi
 
 	# add camera and light source