Changeset 8331 for trunk/anuga_work/development/mem_time_tests/parameters/timelen-over-timestep/main.py

Timestamp:

Feb 10, 2012, 10:34:08 AM (12 years ago)

Author:

pittj

Message:

update for the new experiments

File:

• trunk/anuga_work/development/mem_time_tests/parameters/timelen-over-timestep/main.py (modified) (3 diffs)

trunk/anuga_work/development/mem_time_tests/parameters/timelen-over-timestep/main.py

@@ -5 +5 @@
 import anuga
 import time
-import random
 import subprocess
 import csv
@@ -19 +18 @@
 scenariodir = add_directories(home, ["data","mem_time_test", "parameters",
                                      "timelen-over-timestep"])
-file = 'ex1.csv'
-file_path = os.path.join(scenariodir, file)
 meta = 'metalog.csv'
 meta_path = os.path.join(scenariodir, meta)
@@ -26 +23 @@
 final_path = os.path.join(scenariodir, final)
 
-#create each file
-firstex1 = open(file_path, 'wb')
-spamWriter = csv.writer(firstex1)
-spamWriter.writerow(['Time Length(s)','Time Step (s)' , 'Time Taken','Space Used'])
+# main loops that give the time length (m) and the time step (n)
+for m in range(10,1000,200):
 
-# main loops that give the time length (m) and the time step (n)
-for m in range(1,1000,50):
-
-    for n in range(1,1000,5):
+    for n in range(10,1000,100):
         n = n/100.0 #adjust it so its small enough to see the nature of the relationship
-
-        x = time.time()#time it
         subprocess.call(['python2.5', 'ex1.py',str(m),str(n)])#run it
-        y = time.time()#time it
-        spamWriter.writerow([m,n ,(y-x),'b'])#write the results
-
 print 'Done'
 
-analyse_log(scenariodir, os.path.join(scenariodir,'metalog.csv')) #get the memory usage from the log files
+#get the important data for the experiments from the anuga experiments
+analyse_log(scenariodir, os.path.join(scenariodir,'metalog.csv')) 
 
-#close ex1.csv so we can read from it for a different csv reader object
-firstex1.close()
-
-#merge the metalog useful memory info and the ex1 recorded info into one csv file named final
-ex1 = csv.reader(open(file_path,'rb'))
+#open files to read from and write to
 metalog = csv.reader(open(meta_path,'rb'))
 final = csv.writer(open(final_path,'wb'))
 
+#list stores the index of the values requird
+indexlist = []
 
-for row in ex1:
-    d = metalog.next()
-    final.writerow([row[0],row[1],d[0],d[3],d[6],d[9],d[12]])#,d[15],d[18]])
+#read in the first row
+firstrow = metalog.next()
+
+#get the indices of the values we want, so that the data can be condensed
+indexlist.append(firstrow.index("beforetime"))
+indexlist.append(firstrow.index("aftertime"))
+indexlist.append(firstrow.index("beforesimulationmemory"))
+indexlist.append(firstrow.index("aftermeshmemory"))
+indexlist.append(firstrow.index("afterinitialconditionsmemory"))
+indexlist.append(firstrow.index("afterboundarymemory"))
+indexlist.append(firstrow.index("aftersimulationmemory"))
+indexlist.append(firstrow.index("variable1"))
+
+#not all experiments have 3 variables, but some do
+try:
+    indexlist.append(firstrow.index("variable2"))
+except ValueError:
+    indexlist.append(firstrow.index("beforetime"))
+try:
+    indexlist.append(firstrow.index("variable3"))
+except ValueError:
+    indexlist.append(firstrow.index("aftertime"))
+
+#write the header for the final csv
+final.writerow(["TimeTaken(s)",firstrow[(indexlist[2])],
+                firstrow[(indexlist[3])],firstrow[(indexlist[4])],
+                firstrow[(indexlist[5])],firstrow[(indexlist[6])],
+                firstrow[(indexlist[7])],firstrow[(indexlist[8])],
+                firstrow[(indexlist[9])]])
+
+#write the data for each column in the final csv
+for row in metalog:
+
+    #manipulate the beginning and end time to get the time taken
+    begin = time.strptime(row[(indexlist[0])],'%Y%m%d_%H%M%S')
+    end = time.strptime(row[(indexlist[1])],'%Y%m%d_%H%M%S')
+    taken = time.mktime(end) - time.mktime(begin)
+
+    #write to file
+    final.writerow([str(taken),row[(indexlist[2])],
+                    row[(indexlist[3])],row[(indexlist[4])],
+                    row[(indexlist[5])],row[(indexlist[6])],
+                    row[(indexlist[7])],row[(indexlist[8])],
+                    row[(indexlist[9])]])
+
+