[cig-commits] r8634 -
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/figs
brad at geodynamics.org
brad at geodynamics.org
Mon Dec 10 12:59:16 PST 2007
Author: brad
Date: 2007-12-10 12:59:16 -0800 (Mon, 10 Dec 2007)
New Revision: 8634
Modified:
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/figs/plot_error.py
Log:
Improved viz.
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/figs/plot_error.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/figs/plot_error.py 2007-12-10 20:58:55 UTC (rev 8633)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/figs/plot_error.py 2007-12-10 20:59:16 UTC (rev 8634)
@@ -13,7 +13,7 @@
from plot_geometry import PlotScene
shape = "tet4"
-res = 250
+res = 1000
showSlice = True
class PlotError(PlotScene):
@@ -33,24 +33,30 @@
script.add_source(VTKDataSource(data=data))
script.engine.current_object.name = "Error"
+ if showSlice:
+ slice = ScalarCutPlane()
+ script.add_module(slice)
+ slice.actor.property.opacity = 0.5
+ slice.implicit_plane.origin = (12.0, 12.0, -12.0)
+ slice.implicit_plane.normal = (0, -1.0, 0.0)
+
threshold = Threshold()
script.add_filter(threshold)
- threshold.lower_threshold = -3.5
+ threshold.lower_threshold = -3.0
surf = Surface()
script.add_filter(surf)
if showSlice:
- surf.actor.property.opacity = 0.2
- slice = ScalarCutPlane()
- script.add_module(slice)
- slice.implicit_plane.origin = (12.0, 12.0, -12.0)
- slice.implicit_plane.normal = (0, -1.0, 0.0)
+ surf.actor.property.opacity = 0.3
- colorbar = script.engine.current_object.module_manager.scalar_lut_manager
+ for obj in [slice, surf]:
+ colorbar = obj.module_manager.scalar_lut_manager
+ colorbar.data_range = (threshold.lower_threshold, -2.0)
+ colorbar.lut_mode = "hot"
+ colorbar.reverse_lut = True
colorbar.show_scalar_bar = True
- colorbar.data_range = (-3.5, -2.0)
- colorbar.number_of_labels = 7
- colorbar.scalar_bar.label_format = "%-4.2f"
+ colorbar.number_of_labels = 5
+ colorbar.scalar_bar.label_format = "%-3.1f"
w,h = colorbar.scalar_bar.position2
colorbar.scalar_bar.position2 = (w, 0.1)
@@ -82,11 +88,14 @@
assert(spaceDim == 3)
assert(ncorners == 4)
cellType = tvtk.Tetra().cell_type
+ elif shape == "hex8":
+ assert(spaceDim == 3)
+ assert(ncorners == 8)
+ cellType = tvtk.Hexahedron().cell_type
else:
raise ValueError("Unknown shape '%s'." % shape)
errorLog10 = numpy.log10(error)
- disp = numpy.reshape(disp, (ncells, spaceDim))
data = tvtk.UnstructuredGrid()
data.points = vertices
More information about the cig-commits
mailing list