本文整理汇总了Python中hermes2d.Mesh.plot方法的典型用法代码示例。如果您正苦于以下问题：Python Mesh.plot方法的具体用法？Python Mesh.plot怎么用？Python Mesh.plot使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类hermes2d.Mesh的用法示例。
在下文中一共展示了Mesh.plot方法的5个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: RefSystem
# 需要导入模块: from hermes2d import Mesh [as 别名]
# 或者: from hermes2d.Mesh import plot [as 别名]
# Assemble and solve the fine mesh problem
rs = RefSystem(ls)
rs.assemble()
rs.solve_system(sln_fine)
# Either solve on coarse mesh or project the fine mesh solution
# on the coarse mesh.
if SOLVE_ON_COARSE_MESH:
ls.assemble()
ls.solve_system(sln_coarse)
else:
ls.project_global(sln_fine, sln_coarse)
# View the solution and mesh
sview.show(sln_coarse);
mesh.plot(space=space)
# Calculate error estimate wrt. fine mesh solution
hp = H1Adapt(ls)
hp.set_solutions([sln_coarse], [sln_fine])
err_est = hp.calc_error() * 100
print("Error estimate: %d" % err_est)
# If err_est too large, adapt the mesh
if (err_est < ERR_STOP):
done = True
else:
done = hp.adapt(selector, THRESHOLD, STRATEGY, MESH_REGULARITY)
if (ls.get_num_dofs() >= NDOF_STOP):
done = True
示例2: get_example_mesh
# 需要导入模块: from hermes2d import Mesh [as 别名]
# 或者: from hermes2d.Mesh import plot [as 别名]
#! /usr/bin/env python
# This example shows how to load a mesh, perform various types
# of "manual" element refinements.
# Import modules
from hermes2d import Mesh, MeshView
from hermes2d.examples import get_example_mesh
# Load the mesh file
domain_mesh = get_example_mesh()
mesh = Mesh()
mesh.load(domain_mesh)
# Perform some sample initial refinements
mesh.refine_all_elements(); # Refines all elements.
#mesh.refine_towards_vertex(3, 4); # Refines mesh towards vertex #3 (4x).
#mesh.refine_towards_boundary(2, 4); # Refines all elements along boundary 2 (4x).
# Display the mesh
mesh.plot()
示例3: H1Space
# 需要导入模块: from hermes2d import Mesh [as 别名]
# 或者: from hermes2d.Mesh import plot [as 别名]
# Create the x- and y- displacement space using the default H1 shapeset
xdisp = H1Space(mesh, P_INIT)
ydisp = H1Space(mesh, P_INIT)
set_bc(xdisp, ydisp)
# Initialize the weak formulation
wf = WeakForm(2)
set_forms(wf)
# Initialize the linear system.
ls = LinSystem(wf)
ls.set_spaces(xdisp, ydisp)
# Assemble and solve the matrix problem
xsln = Solution()
ysln = Solution()
ls.assemble()
ls.solve_system(xsln, ysln, lib="scipy")
# Visualize the solution
view = ScalarView("Von Mises stress [Pa]", 50, 50, 1200, 600)
E = float(200e9)
nu = 0.3
l = (E * nu) / ((1 + nu) * (1 - 2*nu))
mu = E / (2*(1 + nu))
stress = VonMisesFilter(xsln, ysln, mu, l)
view.show(stress)
# Visualize the mesh
mesh.plot(space=xdisp)
示例4: Mesh
# 需要导入模块: from hermes2d import Mesh [as 别名]
# 或者: from hermes2d.Mesh import plot [as 别名]
#! /usr/bin/env python
import sys
from hermes2d import Mesh
m = Mesh()
m.load(sys.argv[1])
m.plot(lib="mpl", method="simple")
示例5: RefSystem
# 需要导入模块: from hermes2d import Mesh [as 别名]
# 或者: from hermes2d.Mesh import plot [as 别名]
rs = RefSystem(ls)
rs.assemble()
rs.solve_system(u_sln_fine, v_sln_fine, lib="scipy")
# Either solve on coarse mesh or project the fine mesh solution
# on the coarse mesh.
if SOLVE_ON_COARSE_MESH:
ls.assemble()
ls.solve_system(u_sln_coarse, v_sln_coarse, lib="scipy")
else:
ls.project_global()
# View the solution and meshes
uview.show(u_sln_coarse)
vview.show(v_sln_coarse)
umesh.plot(space=uspace)
vmesh.plot(space=vspace)
# Calculate element errors and total error estimate
hp = H1Adapt(ls)
hp.set_solutions([u_sln_coarse, v_sln_coarse], [u_sln_fine, v_sln_fine]);
set_hp_forms(hp)
err_est = hp.calc_error() * 100
print("Error estimate: %s" % err_est)
# If err_est too large, adapt the mesh
if err_est < ERR_STOP:
done = True
else:
MULTI = False if MULTI == True else True