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()