[cig-commits] r14557 - cs/cigma/trunk/examples/gale
luis at geodynamics.org
luis at geodynamics.org
Wed Apr 1 13:25:10 PDT 2009
Author: luis
Date: 2009-04-01 13:25:09 -0700 (Wed, 01 Apr 2009)
New Revision: 14557
Added:
cs/cigma/trunk/examples/gale/02-compare-analytic.out
cs/cigma/trunk/examples/gale/02-compare-analytic.sh
cs/cigma/trunk/examples/gale/03-visualize-analytic.sh
Log:
Driver script for inclusion analytic comparison
Added: cs/cigma/trunk/examples/gale/02-compare-analytic.out
===================================================================
--- cs/cigma/trunk/examples/gale/02-compare-analytic.out (rev 0)
+++ cs/cigma/trunk/examples/gale/02-compare-analytic.out 2009-04-01 20:25:09 UTC (rev 14557)
@@ -0,0 +1,137 @@
+luis at stress:~/cigma/examples/gale$ time ./02-compare-analytic.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=bm.circular_inclusion.pressure
++ e=inclusion.h5:/error_pressure_512
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:PressureField bm.circular_inclusion.pressure -m /home/luis/bm3/512_8/fields.00000.pvts -o inclusion.h5:/error_pressure_512 -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:PressureField and bm.circular_inclusion.pressure
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 10923 0.4 0 0.4 100%
+
+Summary of comparison:
+ L2 = 0.0196281510417
+ Linf = 0.55285126605
+ volume = 64
+ L2/sqrt(volume) = 0.00245351888022
+ h1 = 0.0220970869121
+
+Creating 'inclusion.h5:/error_pressure_512'
++ for n in 512 256 128
++ a=/home/luis/bm3/256_8/fields.00000.pvts:PressureField
++ b=bm.circular_inclusion.pressure
++ e=inclusion.h5:/error_pressure_256
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ cigma compare /home/luis/bm3/256_8/fields.00000.pvts:PressureField bm.circular_inclusion.pressure -m /home/luis/bm3/512_8/fields.00000.pvts -o inclusion.h5:/error_pressure_256 -v
+
+Comparing /home/luis/bm3/256_8/fields.00000.pvts:PressureField and bm.circular_inclusion.pressure
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 10923 0.4 0 0.4 100%
+
+Summary of comparison:
+ L2 = 0.0272562597333
+ Linf = 0.58406362347
+ volume = 64
+ L2/sqrt(volume) = 0.00340703246666
+ h1 = 0.0441941738242
+
+Updating 'inclusion.h5:/error_pressure_256'
++ for n in 512 256 128
++ a=/home/luis/bm3/128_8/fields.00000.pvts:PressureField
++ b=bm.circular_inclusion.pressure
++ e=inclusion.h5:/error_pressure_128
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ cigma compare /home/luis/bm3/128_8/fields.00000.pvts:PressureField bm.circular_inclusion.pressure -m /home/luis/bm3/512_8/fields.00000.pvts -o inclusion.h5:/error_pressure_128 -v
+
+Comparing /home/luis/bm3/128_8/fields.00000.pvts:PressureField and bm.circular_inclusion.pressure
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 11398 0.38333 0 0.38333 100%
+
+Summary of comparison:
+ L2 = 0.0410194426894
+ Linf = 0.87391931751
+ volume = 64
+ L2/sqrt(volume) = 0.00512743033617
+ h1 = 0.0883883476483
+
+Updating 'inclusion.h5:/error_pressure_128'
++ for v in pressure velocity
++ for n in 512 256 128
++ a=/home/luis/bm3/512_8/fields.00000.pvts:VelocityField
++ b=bm.circular_inclusion.velocity
++ e=inclusion.h5:/error_velocity_512
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ cigma compare /home/luis/bm3/512_8/fields.00000.pvts:VelocityField bm.circular_inclusion.velocity -m /home/luis/bm3/512_8/fields.00000.pvts -o inclusion.h5:/error_velocity_512 -v
+
+Comparing /home/luis/bm3/512_8/fields.00000.pvts:VelocityField and bm.circular_inclusion.velocity
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 10923 0.4 0 0.4 100%
+
+Summary of comparison:
+ L2 = 0.0110382139752
+ Linf = 0.0363051712905
+ volume = 64
+ L2/sqrt(volume) = 0.0013797767469
+ h1 = 0.0220970869121
+
+Updating 'inclusion.h5:/error_velocity_512'
++ for n in 512 256 128
++ a=/home/luis/bm3/256_8/fields.00000.pvts:VelocityField
++ b=bm.circular_inclusion.velocity
++ e=inclusion.h5:/error_velocity_256
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ cigma compare /home/luis/bm3/256_8/fields.00000.pvts:VelocityField bm.circular_inclusion.velocity -m /home/luis/bm3/512_8/fields.00000.pvts -o inclusion.h5:/error_velocity_256 -v
+
+Comparing /home/luis/bm3/256_8/fields.00000.pvts:VelocityField and bm.circular_inclusion.velocity
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 11398 0.38333 0 0.38333 100%
+
+Summary of comparison:
+ L2 = 0.011072234192
+ Linf = 0.0361963869941
+ volume = 64
+ L2/sqrt(volume) = 0.001384029274
+ h1 = 0.0441941738242
+
+Updating 'inclusion.h5:/error_velocity_256'
++ for n in 512 256 128
++ a=/home/luis/bm3/128_8/fields.00000.pvts:VelocityField
++ b=bm.circular_inclusion.velocity
++ e=inclusion.h5:/error_velocity_128
++ m=/home/luis/bm3/512_8/fields.00000.pvts
++ cigma compare /home/luis/bm3/128_8/fields.00000.pvts:VelocityField bm.circular_inclusion.velocity -m /home/luis/bm3/512_8/fields.00000.pvts -o inclusion.h5:/error_velocity_128 -v
+
+Comparing /home/luis/bm3/128_8/fields.00000.pvts:VelocityField and bm.circular_inclusion.velocity
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 10923 0.4 0 0.4 100%
+
+Summary of comparison:
+ L2 = 0.011538895333
+ Linf = 0.0383012086686
+ volume = 64
+ L2/sqrt(volume) = 0.00144236191663
+ h1 = 0.0883883476483
+
+Updating 'inclusion.h5:/error_velocity_128'
+
+real 4m3.575s
+user 3m53.499s
+sys 0m5.744s
+
Added: cs/cigma/trunk/examples/gale/02-compare-analytic.sh
===================================================================
--- cs/cigma/trunk/examples/gale/02-compare-analytic.sh (rev 0)
+++ cs/cigma/trunk/examples/gale/02-compare-analytic.sh 2009-04-01 20:25:09 UTC (rev 14557)
@@ -0,0 +1,40 @@
+#!/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 fields against corresponding analytical solution
+#
+for v in pressure velocity; do
+ for n in 512 256 128; do
+
+ # Path to the field we want
+ a="${DATA}/${n}_8/fields.00000.pvts:${!v}"
+
+ # Name of analytic function
+ b="bm.circular_inclusion.${v}"
+
+ # Name of output array for storing errors for the $v field
+ e="inclusion.h5:/error_${v}_${n}"
+
+ # Use the same mesh for all comparisons (not necessary,
+ # but useful when extracting selections in paraview).
+ m="${DATA}/512_8/fields.00000.pvts"
+
+ cigma compare ${a} ${b} -m ${m} -o ${e} -v
+
+ done
+done
+
Property changes on: cs/cigma/trunk/examples/gale/02-compare-analytic.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/gale/03-visualize-analytic.sh
===================================================================
--- cs/cigma/trunk/examples/gale/03-visualize-analytic.sh (rev 0)
+++ cs/cigma/trunk/examples/gale/03-visualize-analytic.sh 2009-04-01 20:25:09 UTC (rev 14557)
@@ -0,0 +1,36 @@
+#!/bin/bash
+
+# Enable bash debugging
+set -x
+
+# Location of data
+DATA="$HOME/bm3"
+
+# Create visualizations of the field comparisons
+for v in pressure velocity; do
+ for n in 512 256 128; do
+
+ # We used the same mesh for all comparisons
+ m="512_8/fields.00000.pvts"
+ #m="${n}_8/fields.00000.pvts"
+
+ # Location of comparison array we want
+ e="inclusion.h5:error_${v}_${n}"
+
+ # Write out the normalized errors
+ 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
+ vtk-residuals \
+ --divide-by-sqrt-cell-volumes \
+ --output-log-values \
+ --mesh="${m}" \
+ --input="${e}" \
+ --output="vtk/log_error_${v}_${n}.vtk:log_error_${v}"
+
+ done
+done
Property changes on: cs/cigma/trunk/examples/gale/03-visualize-analytic.sh
___________________________________________________________________
Name: svn:executable
+ *
More information about the CIG-COMMITS
mailing list