[cig-commits] r14555 - in cs/cigma/trunk/examples/citcomcu: . vtk
luis at geodynamics.org
luis at geodynamics.org
Wed Apr 1 13:24:57 PDT 2009
Author: luis
Date: 2009-04-01 13:24:56 -0700 (Wed, 01 Apr 2009)
New Revision: 14555
Added:
cs/cigma/trunk/examples/citcomcu/steady-state.out
cs/cigma/trunk/examples/citcomcu/steady-state.sh
cs/cigma/trunk/examples/citcomcu/vtk/
cs/cigma/trunk/examples/citcomcu/vtk/.keep
Log:
Driver script used in steady-state example
Added: cs/cigma/trunk/examples/citcomcu/steady-state.out
===================================================================
--- cs/cigma/trunk/examples/citcomcu/steady-state.out (rev 0)
+++ cs/cigma/trunk/examples/citcomcu/steady-state.out 2009-04-01 20:24:56 UTC (rev 14555)
@@ -0,0 +1,145 @@
+luis at stress:~/cigma/examples/citcomcu$ time ./steady-state.sh
++ T64=20000
++ T32=9900
++ T16=4700
++ T8=2100
++ DATA=/home/luis/bm2/cases
++ for v in temperature velocity
++ for n in 32 16 8
++ ta=T64
++ tb=T32
++ cigma compare -v -a /home/luis/bm2/cases/64/case64.20000.pvtr:temperature -b /home/luis/bm2/cases/32/case32.9900.pvtr:temperature -o steady-state.h5:/error_temperature_64_32_20000_9900
+
+Comparing /home/luis/bm2/cases/64/case64.20000.pvtr:temperature and /home/luis/bm2/cases/32/case32.9900.pvtr:temperature
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 3404.5 1.2833 0 1.2833 100%
+
+Summary of comparison:
+ L2 = 0.00494269807472
+ Linf = 0.060186021859
+ volume = 1
+ L2/sqrt(volume) = 0.00494269807472
+ h1 = 0.0270632938683
+ h2 = 0.0541265877365
+
+Creating 'steady-state.h5:/error_temperature_64_32_20000_9900'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values -m /home/luis/bm2/cases/64/case64.0.pvtr -i steady-state.h5:/error_temperature_64_32_20000_9900 -o vtk/steady_state.log_error_temperature_64_32.vtk:log_error_temperature
+Writing residuals to VTK file 'vtk/steady_state.log_error_temperature_64_32.vtk'
++ for n in 32 16 8
++ ta=T64
++ tb=T16
++ cigma compare -v -a /home/luis/bm2/cases/64/case64.20000.pvtr:temperature -b /home/luis/bm2/cases/16/case16.4700.pvtr:temperature -o steady-state.h5:/error_temperature_64_16_20000_4700
+
+Comparing /home/luis/bm2/cases/64/case64.20000.pvtr:temperature and /home/luis/bm2/cases/16/case16.4700.pvtr:temperature
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 3591 1.2167 0 1.2167 100%
+
+Summary of comparison:
+ L2 = 0.0253326544796
+ Linf = 0.294745572254
+ volume = 1
+ L2/sqrt(volume) = 0.0253326544796
+ h1 = 0.0270632938683
+ h2 = 0.108253175473
+
+Updating 'steady-state.h5:/error_temperature_64_16_20000_4700'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values -m /home/luis/bm2/cases/64/case64.0.pvtr -i steady-state.h5:/error_temperature_64_16_20000_4700 -o vtk/steady_state.log_error_temperature_64_16.vtk:log_error_temperature
+Writing residuals to VTK file 'vtk/steady_state.log_error_temperature_64_16.vtk'
++ for n in 32 16 8
++ ta=T64
++ tb=T8
++ cigma compare -v -a /home/luis/bm2/cases/64/case64.20000.pvtr:temperature -b /home/luis/bm2/cases/8/case8.2100.pvtr:temperature -o steady-state.h5:/error_temperature_64_8_20000_2100
+
+Comparing /home/luis/bm2/cases/64/case64.20000.pvtr:temperature and /home/luis/bm2/cases/8/case8.2100.pvtr:temperature
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 3744.9 1.1667 0 1.1667 100%
+
+Summary of comparison:
+ L2 = 0.0852241175851
+ Linf = 0.587950559527
+ volume = 1
+ L2/sqrt(volume) = 0.0852241175851
+ h1 = 0.0270632938683
+ h2 = 0.216506350946
+
+Updating 'steady-state.h5:/error_temperature_64_8_20000_2100'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values -m /home/luis/bm2/cases/64/case64.0.pvtr -i steady-state.h5:/error_temperature_64_8_20000_2100 -o vtk/steady_state.log_error_temperature_64_8.vtk:log_error_temperature
+Writing residuals to VTK file 'vtk/steady_state.log_error_temperature_64_8.vtk'
++ for v in temperature velocity
++ for n in 32 16 8
++ ta=T64
++ tb=T32
++ cigma compare -v -a /home/luis/bm2/cases/64/case64.20000.pvtr:velocity -b /home/luis/bm2/cases/32/case32.9900.pvtr:velocity -o steady-state.h5:/error_velocity_64_32_20000_9900
+
+Comparing /home/luis/bm2/cases/64/case64.20000.pvtr:velocity and /home/luis/bm2/cases/32/case32.9900.pvtr:velocity
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 3196.9 1.3667 0 1.3667 100%
+
+Summary of comparison:
+ L2 = 0.455909906788
+ Linf = 2.23385858757
+ volume = 1
+ L2/sqrt(volume) = 0.455909906788
+ h1 = 0.0270632938683
+ h2 = 0.0541265877365
+
+Updating 'steady-state.h5:/error_velocity_64_32_20000_9900'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values -m /home/luis/bm2/cases/64/case64.0.pvtr -i steady-state.h5:/error_velocity_64_32_20000_9900 -o vtk/steady_state.log_error_velocity_64_32.vtk:log_error_velocity
+Writing residuals to VTK file 'vtk/steady_state.log_error_velocity_64_32.vtk'
++ for n in 32 16 8
++ ta=T64
++ tb=T16
++ cigma compare -v -a /home/luis/bm2/cases/64/case64.20000.pvtr:velocity -b /home/luis/bm2/cases/16/case16.4700.pvtr:velocity -o steady-state.h5:/error_velocity_64_16_20000_4700
+
+Comparing /home/luis/bm2/cases/64/case64.20000.pvtr:velocity and /home/luis/bm2/cases/16/case16.4700.pvtr:velocity
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 3449.3 1.2667 0 1.2667 100%
+
+Summary of comparison:
+ L2 = 2.2028427971
+ Linf = 10.2394865905
+ volume = 1
+ L2/sqrt(volume) = 2.2028427971
+ h1 = 0.0270632938683
+ h2 = 0.108253175473
+
+Updating 'steady-state.h5:/error_velocity_64_16_20000_4700'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values -m /home/luis/bm2/cases/64/case64.0.pvtr -i steady-state.h5:/error_velocity_64_16_20000_4700 -o vtk/steady_state.log_error_velocity_64_16.vtk:log_error_velocity
+Writing residuals to VTK file 'vtk/steady_state.log_error_velocity_64_16.vtk'
++ for n in 32 16 8
++ ta=T64
++ tb=T8
++ cigma compare -v -a /home/luis/bm2/cases/64/case64.20000.pvtr:velocity -b /home/luis/bm2/cases/8/case8.2100.pvtr:velocity -o steady-state.h5:/error_velocity_64_8_20000_2100
+
+Comparing /home/luis/bm2/cases/64/case64.20000.pvtr:velocity and /home/luis/bm2/cases/8/case8.2100.pvtr:velocity
+
+Progress of comparison
+ cells rate mins eta total progress
+ 262144 3640.9 1.2 0 1.2 100%
+
+Summary of comparison:
+ L2 = 11.2325798478
+ Linf = 64.0976236333
+ volume = 1
+ L2/sqrt(volume) = 11.2325798478
+ h1 = 0.0270632938683
+ h2 = 0.216506350946
+
+Updating 'steady-state.h5:/error_velocity_64_8_20000_2100'
++ vtk-residuals --divide-by-sqrt-cell-volumes --output-log-values -m /home/luis/bm2/cases/64/case64.0.pvtr -i steady-state.h5:/error_velocity_64_8_20000_2100 -o vtk/steady_state.log_error_velocity_64_8.vtk:log_error_velocity
+Writing residuals to VTK file 'vtk/steady_state.log_error_velocity_64_8.vtk'
+
+real 8m6.053s
+user 7m53.438s
+sys 0m2.692s
+
Added: cs/cigma/trunk/examples/citcomcu/steady-state.sh
===================================================================
--- cs/cigma/trunk/examples/citcomcu/steady-state.sh (rev 0)
+++ cs/cigma/trunk/examples/citcomcu/steady-state.sh 2009-04-01 20:24:56 UTC (rev 14555)
@@ -0,0 +1,39 @@
+#!/bin/bash
+
+# Enable bash debugging mode
+set -x
+
+#
+# Matching timesteps for all four cases
+#
+
+T64=20000
+T32=9900
+T16=4700
+T8=2100
+
+#
+# Location of benchmark data
+#
+DATA="$HOME/bm2/cases"
+
+#
+# Run the comparisons over all lower resolutions
+#
+for v in temperature velocity; do
+ for n in 32 16 8; do
+ ta="T64"
+ tb="T$n"
+ cigma compare -v \
+ -a "${DATA}/64/case64.${!ta}.pvtr:${v}" \
+ -b "${DATA}/${n}/case${n}.${!tb}.pvtr:${v}" \
+ -o "steady-state.h5:/error_${v}_64_${n}_${!ta}_${!tb}"
+ vtk-residuals \
+ --divide-by-sqrt-cell-volumes \
+ --output-log-values \
+ -m "${DATA}/64/case64.0.pvtr" \
+ -i "steady-state.h5:/error_${v}_64_${n}_${!ta}_${!tb}" \
+ -o "vtk/steady_state.log_error_${v}_64_${n}.vtk:log_error_${v}"
+ done
+done
+
Property changes on: cs/cigma/trunk/examples/citcomcu/steady-state.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/citcomcu/vtk/.keep
===================================================================
More information about the CIG-COMMITS
mailing list