[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