[cig-commits] r14560 - cs/cigma/trunk/examples/gale

[luis at geodynamics.org]

Author: luis
Date: 2009-04-01 13:32:06 -0700 (Wed, 01 Apr 2009)
New Revision: 14560

Added:
   cs/cigma/trunk/examples/gale/calculate-norms.out
   cs/cigma/trunk/examples/gale/calculate-norms.sh
Log:
Driver script for calculating the norms of the
gale results for the inclusion problem

Added: cs/cigma/trunk/examples/gale/calculate-norms.out
===================================================================
--- cs/cigma/trunk/examples/gale/calculate-norms.out	                        (rev 0)
+++ cs/cigma/trunk/examples/gale/calculate-norms.out	2009-04-01 20:32:06 UTC (rev 14560)
@@ -0,0 +1,119 @@
+luis at stress:~/cigma/examples/gale$ time ./calculate-norms.sh 
++ DATA=/home/luis/bm3
++ pressure=PressureField
++ velocity=VelocityField
++ for v in pressure velocity
++ for n in 512 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:PressureField
++ b=zero
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:PressureField zero -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:PressureField and zero
+
+Progress of comparison
+    cells      rate      mins       eta     total   progress
+   262144     10923       0.4         0       0.4       100%                     
+
+Summary of comparison: 
+  L2 = 0.0797157166649
+  Linf = 0.877584640981
+  volume = 64
+  L2/sqrt(volume) = 0.00996446458312
+  h1 = 0.0220970869121
+
++ for n in 512 256 128
++ a=/home/luis/bm3/256_8/fields.00000.pvts:PressureField
++ b=zero
++ cigma compare /home/luis/bm3/256_8/fields.00000.pvts:PressureField zero -v
+
+Comparing /home/luis/bm3/256_8/fields.00000.pvts:PressureField and zero
+
+Progress of comparison
+    cells      rate      mins       eta     total   progress
+    65536     10923       0.1         0       0.1       100%                     
+
+Summary of comparison: 
+  L2 = 0.075451883375
+  Linf = 0.669648916386
+  volume = 64
+  L2/sqrt(volume) = 0.00943148542188
+  h1 = 0.0441941738242
+
++ for n in 512 256 128
++ a=/home/luis/bm3/128_8/fields.00000.pvts:PressureField
++ b=zero
++ cigma compare /home/luis/bm3/128_8/fields.00000.pvts:PressureField zero -v
+
+Comparing /home/luis/bm3/128_8/fields.00000.pvts:PressureField and zero
+
+Progress of comparison
+    cells      rate      mins       eta     total   progress
+    16384      8192  0.033333         0  0.033333       100%                    
+
+Summary of comparison: 
+  L2 = 0.0622282671898
+  Linf = 0.393081675605
+  volume = 64
+  L2/sqrt(volume) = 0.00777853339873
+  h1 = 0.0883883476483
+
++ for v in pressure velocity
++ for n in 512 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:VelocityField
++ b=zero
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:VelocityField zero -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:VelocityField and zero
+
+Progress of comparison
+    cells      rate      mins       eta     total   progress
+   262144     10486   0.41667         0   0.41667       100%                     
+
+Summary of comparison: 
+  L2 = 52.2557823347
+  Linf = 11.3090388351
+  volume = 64
+  L2/sqrt(volume) = 6.53197279183
+  h1 = 0.0220970869121
+
++ for n in 512 256 128
++ a=/home/luis/bm3/256_8/fields.00000.pvts:VelocityField
++ b=zero
++ cigma compare /home/luis/bm3/256_8/fields.00000.pvts:VelocityField zero -v
+
+Comparing /home/luis/bm3/256_8/fields.00000.pvts:VelocityField and zero
+
+Progress of comparison
+    cells      rate      mins       eta     total   progress
+    65536     10923       0.1         0       0.1       100%                     
+
+Summary of comparison: 
+  L2 = 52.255782307
+  Linf = 11.3043691711
+  volume = 64
+  L2/sqrt(volume) = 6.53197278838
+  h1 = 0.0441941738242
+
++ for n in 512 256 128
++ a=/home/luis/bm3/128_8/fields.00000.pvts:VelocityField
++ b=zero
++ cigma compare /home/luis/bm3/128_8/fields.00000.pvts:VelocityField zero -v
+
+Comparing /home/luis/bm3/128_8/fields.00000.pvts:VelocityField and zero
+
+Progress of comparison
+    cells      rate      mins       eta     total   progress
+    16384     16384  0.016667         0  0.016667       100%                    
+
+Summary of comparison: 
+  L2 = 52.2557822561
+  Linf = 11.2950298425
+  volume = 64
+  L2/sqrt(volume) = 6.53197278201
+  h1 = 0.0883883476483
+
+
+real	1m36.221s
+user	1m34.134s
+sys	0m1.680s
+

Added: cs/cigma/trunk/examples/gale/calculate-norms.sh
===================================================================
--- cs/cigma/trunk/examples/gale/calculate-norms.sh	                        (rev 0)
+++ cs/cigma/trunk/examples/gale/calculate-norms.sh	2009-04-01 20:32:06 UTC (rev 14560)
@@ -0,0 +1,33 @@
+#!/bin/bash
+
+# Enable bash debugging
+set -x
+
+#
+# Location of data files
+#
+DATA="$HOME/bm3"
+
+#
+# Names of the fields
+#
+pressure=PressureField
+velocity=VelocityField
+
+#
+# Calculate norm of each field.
+#
+for v in pressure velocity; do
+    for n in 512 256 128; do
+
+        # Location of the field we want
+        a="${DATA}/${n}_8/fields.00000.pvts:${!v}"
+
+        # Function which always evaluates to 0
+        b="zero"
+
+        cigma compare ${a} ${b} -v
+
+    done
+done
+


Property changes on: cs/cigma/trunk/examples/gale/calculate-norms.sh
___________________________________________________________________
Name: svn:executable
   + *