[cig-commits] r9141 - cs/benchmark/cigma/trunk/src/tests

luis at geodynamics.org luis at geodynamics.org
Fri Jan 25 09:09:15 PST 2008


Author: luis
Date: 2008-01-25 09:09:15 -0800 (Fri, 25 Jan 2008)
New Revision: 9141

Added:
   cs/benchmark/cigma/trunk/src/tests/strikeslip-compare.py
Log:
Python script for running all cases of the pylith strikeslip benchmark
(16 comparisons involving two resolutions, two element types)


Added: cs/benchmark/cigma/trunk/src/tests/strikeslip-compare.py
===================================================================
--- cs/benchmark/cigma/trunk/src/tests/strikeslip-compare.py	                        (rev 0)
+++ cs/benchmark/cigma/trunk/src/tests/strikeslip-compare.py	2008-01-25 17:09:15 UTC (rev 9141)
@@ -0,0 +1,53 @@
+#!/usr/bin/env python
+
+
+import os
+from time import time
+
+
+elements = ['tet4', 'hex8']
+resolutions = ['1000m', '0500m']
+
+cases = []
+for e in elements:
+    for r in resolutions:
+        cases.append('%s_%s' % (e,r))
+
+
+pathstr = './strikeslip_%s_t0.vtk:displacements_t0'
+outputstr = './residuals%d%d-%s-%s.vtk'
+
+comparisons = []
+for i,a in enumerate(cases):
+    for j,b in enumerate(cases):
+        input1 = pathstr % a
+        input2 = pathstr % b
+        output = outputstr % (i,j,a,b)
+        comparison = dict(first=input1, second=input2, output=output, case=(i,j))
+        comparisons.append(comparison)
+
+
+commandstr = ' '.join([
+    '../cigma', 'compare',
+    '--output=%(output)s',
+    '--first=%(first)s',
+    '--second=%(second)s',
+    '--verbose',
+    '--output-frequency=1000',
+])
+
+t0 = time()
+for comparison in comparisons:
+    i,j = comparison['case']
+    print "Case %d%d" % (i,j)
+    print "Comparing %s with %s" % (cases[i], cases[j])
+    command = commandstr % comparison
+    print command
+    os.system(command)
+    print "-" * 78
+t1 = time()
+
+dt = (t1 - t0)/60.0
+
+print "Total: %lf mins" % dt
+


Property changes on: cs/benchmark/cigma/trunk/src/tests/strikeslip-compare.py
___________________________________________________________________
Name: svn:executable
   + *



More information about the cig-commits mailing list