thermalOutput = output.outputTypes['scalar field'].clone()
thermalOutput.plugParameter(temperature, "field")
rscl = outputClones.RescaleOutput.clone()
rscl.plugInput(thermalOutput)
rscl.plugParameter(0.0, "minimum")
rscl.plugParameter(0.95, "maximum")
import drawmesh
import psoutput
import color
psdriver = psoutput.PSoutput(psfilename)
psdriver.set_colormap(color.SpectralMap())
drawmesh.draw_mesh(mesh, psdriver,
[positionSum, outputClones.posOutput],
valueOutput=rscl,
drawedges=0)
## psgfx.psoutput(mesh, psfilename, psgfx.PSgeometry(),
## [positionSum,outputClones.posOutput],
## valueOutput=colormap,
## drawedges=1)
else:
for node in mesh.funcnode_iterator():
pos = node.position()
print displacement.value(mesh, node, 0), \
displacement.value(mesh, node, 1), \
temperature.value(mesh, node, 0)
except ErrSetupError, s:
print "Caught ErrSetupError."
thermalOutput = output.outputTypes['scalar field'].clone()
thermalOutput.plugParameter(temperature, "field")
rscl = outputClones.RescaleOutput.clone()
rscl.plugInput(thermalOutput)
rscl.plugParameter(0.0, "minimum")
rscl.plugParameter(0.95, "maximum")
import drawmesh
import psoutput
import color
psdriver = psoutput.PSoutput(psfilename)
psdriver.set_colormap(color.SpectralMap())
drawmesh.draw_mesh(mesh, psdriver,
[positionSum, outputClones.posOutput],
valueOutput=rscl,
drawedges=0)
## psgfx.psoutput(mesh, psfilename, psgfx.PSgeometry(),
## [positionSum,outputClones.posOutput],
## valueOutput=colormap,
## drawedges=1)
else:
for node in mesh.funcnode_iterator():
pos = node.position()
print displacement.value(mesh, node, 0), \
displacement.value(mesh, node, 1), \
temperature.value(mesh, node, 0)
except ErrSetupError, s:
print "Caught ErrSetupError."
print "BC's defined."
trace_enable()
mesh.equilibrate(maxiters, tolerance, krylov_dim)
# trace_disable()
if psfilename:
deformOutput = outputClones.VectorFieldOutput.clone()
deformOutput.plugParameter(displacement, "field")
positionSum = outputClones.VFieldSumOutput.clone()
positionSum.plugInput(deformOutput)
positionSum.plugInput(outputClones.posOutput)
import drawmesh
import psoutput
psdriver = psoutput.PSoutput(psfilename)
drawmesh.draw_mesh(mesh, psdriver,
[positionSum, outputClones.posOutput])
else:
# trace_disable()
for node in mesh.funcnode_iterator():
pos = node.position()
print displacement.value(mesh, node, 0),\
displacement.value(mesh, node, 1)
except ErrSetupError, s:
print "Caught ErrSetupError."
print "Message string: " + s.message()
print "BC's defined."
trace_enable()
mesh.equilibrate(maxiters, tolerance, krylov_dim)
# trace_disable()
if psfilename:
deformOutput = outputClones.VectorFieldOutput.clone()
deformOutput.plugParameter(displacement,"field")
positionSum = outputClones.VFieldSumOutput.clone()
positionSum.plugInput(deformOutput)
positionSum.plugInput(outputClones.posOutput)
import drawmesh
import psoutput
psdriver = psoutput.PSoutput(psfilename)
drawmesh.draw_mesh(mesh, psdriver,
[positionSum, outputClones.posOutput])
else:
# trace_disable()
for node in mesh.funcnode_iterator():
pos = node.position()
print displacement.value(mesh, node, 0),\
displacement.value(mesh, node, 1)
except ErrSetupError, s:
print "Caught ErrSetupError."
print "Message string: " + s.message()
