[cig-commits] r14559 - cs/cigma/trunk/examples/gale
luis at geodynamics.org
luis at geodynamics.org
Wed Apr 1 13:25:26 PDT 2009
Author: luis
Date: 2009-04-01 13:25:23 -0700 (Wed, 01 Apr 2009)
New Revision: 14559
Added:
cs/cigma/trunk/examples/gale/02-compare-fields.out
cs/cigma/trunk/examples/gale/02-compare-fields.sh
cs/cigma/trunk/examples/gale/03-visualize-fields.out
cs/cigma/trunk/examples/gale/03-visualize-fields.sh
Modified:
cs/cigma/trunk/examples/gale/03-visualize-analytic.sh
Log:
Comparisons between the Gale results, for inclusion benchmark
Added: cs/cigma/trunk/examples/gale/02-compare-fields.out
===================================================================
--- cs/cigma/trunk/examples/gale/02-compare-fields.out (rev 0)
+++ cs/cigma/trunk/examples/gale/02-compare-fields.out 2009-04-01 20:25:23 UTC (rev 14559)
@@ -0,0 +1,95 @@
+luis at stress:~/cigma/examples/gale$ time ./02-compare-fields.sh
++ DATA=/home/luis/bm3
++ pressure=PressureField
++ velocity=VelocityField
++ for v in pressure velocity
++ for n in 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:PressureField
++ b=/home/luis/bm3/256_8/fields.00000.pvts:PressureField
++ e=inclusion.h5:error_pressure_512_256
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:PressureField /home/luis/bm3/256_8/fields.00000.pvts:PressureField -o inclusion.h5:error_pressure_512_256 -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:PressureField and /home/luis/bm3/256_8/fields.00000.pvts:PressureField
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 5461.3 0.8 0 0.8 100%
+
+Summary of comparison:
+ L2 = 0.0126291594164
+ Linf = 0.257430935978
+ volume = 64
+ L2/sqrt(volume) = 0.00157864492705
+ h1 = 0.0220970869121
+ h2 = 0.0441941738242
+
+Updating 'inclusion.h5:error_pressure_512_256'
++ for n in 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:PressureField
++ b=/home/luis/bm3/128_8/fields.00000.pvts:PressureField
++ e=inclusion.h5:error_pressure_512_128
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:PressureField /home/luis/bm3/128_8/fields.00000.pvts:PressureField -o inclusion.h5:error_pressure_512_128 -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:PressureField and /home/luis/bm3/128_8/fields.00000.pvts:PressureField
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 5698.8 0.76667 0 0.76667 100%
+
+Summary of comparison:
+ L2 = 0.0325962809388
+ Linf = 0.534620275095
+ volume = 64
+ L2/sqrt(volume) = 0.00407453511735
+ h1 = 0.0220970869121
+ h2 = 0.0883883476483
+
+Updating 'inclusion.h5:error_pressure_512_128'
++ for v in pressure velocity
++ for n in 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:VelocityField
++ b=/home/luis/bm3/256_8/fields.00000.pvts:VelocityField
++ e=inclusion.h5:error_velocity_512_256
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:VelocityField /home/luis/bm3/256_8/fields.00000.pvts:VelocityField -o inclusion.h5:error_velocity_512_256 -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:VelocityField and /home/luis/bm3/256_8/fields.00000.pvts:VelocityField
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 5461.3 0.8 0 0.8 100%
+
+Summary of comparison:
+ L2 = 0.000180938107628
+ Linf = 0.00280286609087
+ volume = 64
+ L2/sqrt(volume) = 2.26172634535e-05
+ h1 = 0.0220970869121
+ h2 = 0.0441941738242
+
+Updating 'inclusion.h5:error_velocity_512_256'
++ for n in 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:VelocityField
++ b=/home/luis/bm3/128_8/fields.00000.pvts:VelocityField
++ e=inclusion.h5:error_velocity_512_128
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:VelocityField /home/luis/bm3/128_8/fields.00000.pvts:VelocityField -o inclusion.h5:error_velocity_512_128 -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:VelocityField and /home/luis/bm3/128_8/fields.00000.pvts:VelocityField
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 5461.3 0.8 0 0.8 100%
+
+Summary of comparison:
+ L2 = 0.000833054185043
+ Linf = 0.00580947167519
+ volume = 64
+ L2/sqrt(volume) = 0.00010413177313
+ h1 = 0.0220970869121
+ h2 = 0.0883883476483
+
+Updating 'inclusion.h5:error_velocity_512_128'
+
+real 4m5.794s
+user 4m1.367s
+sys 0m3.368s
+
Added: cs/cigma/trunk/examples/gale/02-compare-fields.sh
===================================================================
--- cs/cigma/trunk/examples/gale/02-compare-fields.sh (rev 0)
+++ cs/cigma/trunk/examples/gale/02-compare-fields.sh 2009-04-01 20:25:23 UTC (rev 14559)
@@ -0,0 +1,36 @@
+#!/bin/bash
+
+# Enable bash debugging
+set -x
+
+#
+# Location of data files
+#
+DATA="$HOME/bm3"
+
+
+#
+# Names of the fields in the VTK .pvts files
+#
+pressure=PressureField
+velocity=VelocityField
+
+#
+# Compare solutions against highest resolution field
+#
+for v in pressure velocity; do
+ for n in 256 128; do
+
+ # The first field (whose discretization is also used for the integration mesh)
+ a="${DATA}/512_8/fields.00000.pvts:${!v}"
+
+ # The second field
+ b="${DATA}/${n}_8/fields.00000.pvts:${!v}"
+
+ # Output HDF5 array for storing the comparison
+ e="inclusion.h5:error_${v}_512_${n}"
+
+ cigma compare ${a} ${b} -o ${e} -v
+
+ done
+done
Property changes on: cs/cigma/trunk/examples/gale/02-compare-fields.sh
___________________________________________________________________
Name: svn:executable
+ *
Modified: cs/cigma/trunk/examples/gale/03-visualize-analytic.sh
===================================================================
--- cs/cigma/trunk/examples/gale/03-visualize-analytic.sh 2009-04-01 20:25:16 UTC (rev 14558)
+++ cs/cigma/trunk/examples/gale/03-visualize-analytic.sh 2009-04-01 20:25:23 UTC (rev 14559)
@@ -14,17 +14,17 @@
m="${DATA}/512_8/fields.00000.pvts"
#m="${n}_8/fields.00000.pvts"
- # Location of comparison array we want
+ # Location of the comparison data
e="inclusion.h5:error_${v}_${n}"
- # Write out the normalized errors
+ # Write out the errors, normalized by cell size
vtk-residuals \
--divide-by-sqrt-cell-volumes \
--mesh="${m}" \
--input="${e}" \
--output="vtk/error_${v}_${n}.vtk:error_${v}"
- # Write out normalized errors using log-scale
+ # Same as above, but using a log-scale
vtk-residuals \
--divide-by-sqrt-cell-volumes \
--output-log-values \
Added: cs/cigma/trunk/examples/gale/03-visualize-fields.out
===================================================================
--- cs/cigma/trunk/examples/gale/03-visualize-fields.out (rev 0)
+++ cs/cigma/trunk/examples/gale/03-visualize-fields.out 2009-04-01 20:25:23 UTC (rev 14559)
@@ -0,0 +1,37 @@
+luis at stress:~/cigma/examples/gale$ time ./03-visualize-fields.sh
++ DATA=/home/luis/bm3
++ for v in pressure velocity
++ for n in 256 128
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ e=inclusion.h5:error_pressure_512_256
++ vtk-residuals --divide-by-sqrt-cell-volumes --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_pressure_512_256 --output=vtk/error_pressure_512_256.vtk:error_pressure
+Writing residuals to VTK file 'vtk/error_pressure_512_256.vtk'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_pressure_512_256 --output=vtk/log_error_pressure_512_256.vtk:log_error_pressure
+Writing residuals to VTK file 'vtk/log_error_pressure_512_256.vtk'
++ for n in 256 128
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ e=inclusion.h5:error_pressure_512_128
++ vtk-residuals --divide-by-sqrt-cell-volumes --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_pressure_512_128 --output=vtk/error_pressure_512_128.vtk:error_pressure
+Writing residuals to VTK file 'vtk/error_pressure_512_128.vtk'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_pressure_512_128 --output=vtk/log_error_pressure_512_128.vtk:log_error_pressure
+Writing residuals to VTK file 'vtk/log_error_pressure_512_128.vtk'
++ for v in pressure velocity
++ for n in 256 128
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ e=inclusion.h5:error_velocity_512_256
++ vtk-residuals --divide-by-sqrt-cell-volumes --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_velocity_512_256 --output=vtk/error_velocity_512_256.vtk:error_velocity
+Writing residuals to VTK file 'vtk/error_velocity_512_256.vtk'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_velocity_512_256 --output=vtk/log_error_velocity_512_256.vtk:log_error_velocity
+Writing residuals to VTK file 'vtk/log_error_velocity_512_256.vtk'
++ for n in 256 128
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ e=inclusion.h5:error_velocity_512_128
++ vtk-residuals --divide-by-sqrt-cell-volumes --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_velocity_512_128 --output=vtk/error_velocity_512_128.vtk:error_velocity
+Writing residuals to VTK file 'vtk/error_velocity_512_128.vtk'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values --mesh=/home/luis/bm3/512_8/fields.00000.pvts --input=inclusion.h5:error_velocity_512_128 --output=vtk/log_error_velocity_512_128.vtk:log_error_velocity
+Writing residuals to VTK file 'vtk/log_error_velocity_512_128.vtk'
+
+real 1m40.068s
+user 1m33.974s
+sys 0m5.592s
+
Added: cs/cigma/trunk/examples/gale/03-visualize-fields.sh
===================================================================
--- cs/cigma/trunk/examples/gale/03-visualize-fields.sh (rev 0)
+++ cs/cigma/trunk/examples/gale/03-visualize-fields.sh 2009-04-01 20:25:23 UTC (rev 14559)
@@ -0,0 +1,35 @@
+#!/bin/bash
+
+# Enable bash debugging
+set -x
+
+# Location of data
+DATA="$HOME/bm3"
+
+for v in pressure velocity; do
+ for n in 256 128; do
+
+ # Since the first field was always the same, we also used
+ # the same mesh for all the comparisons
+ m="${DATA}/512_8/fields.00000.pvts"
+
+ # Location of the comparison data
+ e="inclusion.h5:error_${v}_512_${n}"
+
+ # Write out the errors, normalized by the cell size
+ vtk-residuals \
+ --divide-by-sqrt-cell-volumes \
+ --mesh="${m}" \
+ --input="${e}" \
+ --output="vtk/error_${v}_512_${n}.vtk:error_${v}"
+
+ # Same as above, but using a log-scale
+ vtk-residuals \
+ --divide-by-sqrt-cell-volumes \
+ --output-log-values \
+ --mesh="${m}" \
+ --input="${e}" \
+ --output="vtk/log_error_${v}_512_${n}.vtk:log_error_${v}"
+
+ done
+done
Property changes on: cs/cigma/trunk/examples/gale/03-visualize-fields.sh
___________________________________________________________________
Name: svn:executable
+ *
More information about the CIG-COMMITS
mailing list