Commit: a9047d613054dcc50dc8f5fe228dfcefd4922d77 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revisão: a9047d613054dcc50dc8f5fe228dfcefd4922d77 (tree)
• Hora: 2008-11-12 23:11:03
• Autor: iselllo
• Commiter: iselllo

Mensagem de Log

I generalized the code in order to be able to take copies of 1D arrays
as well.

Mudança Sumário

• modified: Python-codes/replica.py

Diff

diff -r 87b3de796b4e -r a9047d613054 Python-codes/replica.py
--- a/Python-codes/replica.py	Wed Nov 12 13:54:51 2008 +0000
+++ b/Python-codes/replica.py	Wed Nov 12 14:11:03 2008 +0000
@@ -8,20 +8,35 @@
     #temp=pos_arr
     dim=s.shape(pos_arr)
     n_row=dim[0]
-    n_col=dim[1]
-    print "n_row and n_col are, ", n_row, n_col
-    
-    copy=s.zeros(((n_row*N_rep),n_col))
+    if (len(dim)>1):
+        n_col=dim[1]
+        print "n_row and n_col are, ", n_row, n_col
+    else:
+        print "n_row is, ", n_row
+
+
+    if (len(dim)>1):
+        copy=s.zeros(((n_row*N_rep),n_col))
+    else:
+        copy=s.zeros((n_row*N_rep))
+        
+
+
+        
     for i in xrange(N_rep):
-        copy[(i*n_row):((i+1)*n_row),:]=pos_arr
+        if (len(dim)>1):
+            copy[(i*n_row):((i+1)*n_row),:]=pos_arr
+        else:
+            copy[(i*n_row):((i+1)*n_row)]=pos_arr
+            
 
     return copy
 
-cluster=p.load("initial_folded_larger_box")
+cluster=p.load("eta_0")
 
-my_copy=replica(cluster,1000)
+my_copy=replica(cluster,200)
 
-p.save("copy_of_cluster50_1000", my_copy)
+p.save("eta_0_200", my_copy)
     
 print "So far so good"