g.addStructuredSurface([0,1,2,3])
g.addSurface([4,1])
# ---- Create mesh ----------------------------------------------------------
# Element type 3 is quad. (2 is triangle. See user manual for more element types)
elType = 3
# Degrees of freedom per node.
dofsPerNode = 1
coords, edof, dofs, bdofs, elementmarkers = cfm.createGmshMesh(g,
elSizeFactor = 0.05,
elType = elType,
dofsPerNode = dofsPerNode)
# ---- Visualise mesh -------------------------------------------------------
# Draw geometry
cfv.drawGeometry(g)
# Draw mesh
vv.figure()
cfv.drawMesh(coords=coords, edof=edof, dofsPerNode=dofsPerNode, elType=elType, filled=True)
# Enter main loop
g.addStructuredVolume([0,1,2,3,4,5], 0, marker=90)
# ---- Create mesh ----------------------------------------------------------
# Element type 5 is hexahedron. (See user manual for more element types)
elType = 5
# Degrees of freedom per node.
dofsPerNode= 1
# Create mesh
coords, edof, dofs, bdofs, _ = cfm.createGmshMesh(geometry = g, elType = elType,
dofsPerNode = dofsPerNode)
# ---- Visualise mesh -------------------------------------------------------
# Hold Left Mouse button to rotate.
# Hold right mouse button to zoom.
# Hold SHIFT and left mouse button to pan.
# Hold SHIFT and right mouse button to change the field of view.
# Hold Ctrl and left mouse button to roll the camera.
# Draw geometry
cfv.drawGeometry(g, drawPoints=False)
# Draw mesh
