Commit: 905702c77dea10868c8be2a619325c1f5b9e6695 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revisão: 905702c77dea10868c8be2a619325c1f5b9e6695 (tree)
• Hora: 2008-06-09 19:39:47
• Autor: iselllo
• Commiter: iselllo

Mensagem de Log

I modified the script povray_test.py, which now saves the .pov files corresponding to each snapshot of the system
sequentially. Then I have to post-process them with povray.

Mudança Sumário

• modified: Python-codes/povray_test.py

Diff

diff -r fdb17f24d848 -r 905702c77dea Python-codes/povray_test.py
--- a/Python-codes/povray_test.py	Fri Jun 06 16:07:42 2008 +0000
+++ b/Python-codes/povray_test.py	Mon Jun 09 10:39:47 2008 +0000
@@ -229,19 +229,12 @@
   
 #  pos_list=s.arange(30)*1.
   #gamma = (sqrt(5)-1)/2
-  file = File("my_output.pov")
-  Camera(location=(40,40,100), look_at=(40,40,40)).write(file)
+  file = File()
+  Camera(location=(90,90,90), look_at=(40,40,40)).write(file)
 #   LightSource((40,40,0), color=(0,0,1)).write(file)
 #   LightSource((0,40,40), color=(0,0,1)).write(file)
 #   LightSource((40,0,40), color=(0,0,1)).write(file)
-  LightSource((40,40,80), color=(1,1,1)).write(file)
-  LightSource((0,40,80), color=(1,1,1)).write(file)
-  #LightSource((80,0,80), color=(1,1,1)).write(file)
- # LightSource((0,0,0), color=(1,1,1)).write(file)
-  LightSource((40,0,40), color=(1,1,1)).write(file)
-
-
-
+  LightSource((80,80,80), color=(1,1,1)).write(file)
   #theta = 0.0
   for i in range(len(pos_list)/3):
     #r = i * 0.5  #position vector
@@ -253,19 +246,107 @@
         Finish(
           phong =  0.8, 
           ambient = 0.0,
-          diffuse = 0.8,
-          reflection = 0.25,
+          diffuse = 0.5,
+          reflection = 0.15,
           specular = 0.5
         ),
         Pigment ( color=color ))
         #pigment { color rgb <1,1,0> }
         #Pigment(color=(9,9,9)))
     ).write(file)
+
    
 
 
 
-z=my_spiral()
+def visualize_box(snapshot_label):
+  " Fibonacci spiral "
+
+  box=(5000./0.01)**(1./3.)
+  snapshot="read_pos_%1d"%snapshot_label
+  print "snapshot is, ", snapshot
+  pos_list=p.load(snapshot)
+  #pos_list=s.reshape(pos_list,(1,15000))
+  pos_list=s.remainder(pos_list,box)
+  #pos_list=pos_list[0:30]
+
+  print "pos_list is, ", pos_list
+  print "s.shape(pos_list) is, ", s.shape(pos_list)
+  
+#  pos_list=s.arange(30)*1.
+  #gamma = (sqrt(5)-1)/2
+  #file = File("my_output.pov")
+  snout="snap_out_%05d"%snapshot_label
+  snout = snout + ".pov"
+  print "snout is, ", snout
+  #file = File("snapshot.pov")
+  file = File(snout)
+  #file = File()
+  Camera(location=(90,40,40), look_at=(40,40,40)).write(file)
+#   LightSource((40,40,0), color=(0,0,1)).write(file)
+#   LightSource((0,40,40), color=(0,0,1)).write(file)
+#   LightSource((40,0,40), color=(0,0,1)).write(file)
+  LightSource((85,40,40), color=(2.5,2.5,2.5)).write(file)
+  LightSource((85,45,85), color=(2.5,2.5,2.5)).write(file)
+  LightSource((45,85,85), color=(2.5,2.5,2.5)).write(file)
+  
+  #theta = 0.0
+  for i in range(len(pos_list)/3):
+    #r = i * 0.5  #position vector
+    color = 0,0,1
+    
+    v = [ pos_list[i, 0], pos_list[i, 1], pos_list[i, 2]] #particle actual positions
+    Sphere( v, 0.5,
+      Texture(
+        Finish(
+          phong =  0.8, 
+          ambient = 0.0,
+          diffuse = 0.5,
+          reflection = 0.15,
+          specular = 0.5
+        ),
+        Pigment ( color=color ))
+        #pigment { color rgb <1,1,0> }
+        #Pigment(color=(9,9,9)))
+    ).write(file)
+ 
+
+  
+
+    
+
+#    os.system('enter your bash command here')
+   
 
 
+
+#z=my_spiral()
+
+#snap=20
+
+d=s.arange(0,1500,10)
+
+for i in d :
+  print "i is, ", i
+  #z=visualize_box(snap)
+
+  z=visualize_box(i)
+
+
+
+# label=20
+# snout="snap_out_%05d"%label
+# snout = snout + ".pov"
+# print "snout is, ", snout
+
+# #`basename $i|sed -e s/.gz//`
+
+# os.system('povray snap_out_00020.pov -w1024 -h1024')
+
+# #os.system('povray `$snout`')
+
+
+
+#z=spiral()
+
 print "So far so good"