N1, N2 = 10, 5 # Number of elements
l1, l2 = 4.0, 2.0 # Mesh size
fs = 20.0 # fontsize
mesh = Mesh()
nodes = mesh.nodes
nodes.add_node(label=1, x=0.0, y=0.0)
nodes.add_node(label=2, x=1.0, y=0.0)
nodes.add_node(label=3, x=0.0, y=1.0)
nodes.add_node(label=4, x=1.5, y=1.0)
nodes.add_node(label=5, x=1.0, y=2.0)
mesh.add_element(label=1, connectivity=[1, 2, 3], space=2)
mesh.add_element(label=2, connectivity=[2, 4, 5, 3], space=2)
centroids = mesh.centroids()
plt.figure(figsize=(8, 3))
plt.gca().set_aspect("equal")
nodes = mesh.nodes
xn, yn, zn = np.array(nodes.x), np.array(nodes.y), np.array(nodes.z) # Nodes coordinates
xe, ye, ze = mesh.get_edges() # Mesh edges
xb, yb, zb = mesh.get_border() # Mesh border
plt.plot(xe, ye, "r-", label="Edges")
plt.plot(xb, yb, "b-", label="Border")
plt.plot(xn, yn, "go", label="Nodes")
plt.xlim([-0.1 * l1, 1.1 * l1])
plt.ylim([-0.1 * l2, 1.1 * l2])
plt.xlabel("$x$", fontsize=fs)
N1, N2 = 10, 5 # Number of elements
l1, l2 = 4., 2. # Mesh size
fs = 20. # fontsize
mesh = Mesh()
nodes = mesh.nodes
nodes.add_node(label=1, x=0., y=0.)
nodes.add_node(label=2, x=1., y=0.)
nodes.add_node(label=3, x=0., y=1.)
nodes.add_node(label=4, x=1.5, y=1.)
nodes.add_node(label=5, x=1., y=2.)
mesh.add_element(label=1, connectivity=[1, 2, 3], space=2)
mesh.add_element(label=2, connectivity=[2, 4, 5, 3], space=2)
centroids = mesh.centroids()
plt.figure(figsize=(8, 3))
plt.gca().set_aspect('equal')
nodes = mesh.nodes
xn, yn, zn = np.array(nodes.x), np.array(
nodes.y), np.array(nodes.z) # Nodes coordinates
xe, ye, ze = mesh.get_edges() # Mesh edges
xb, yb, zb = mesh.get_border() # Mesh border
plt.plot(xe, ye, 'r-', label='Edges')
plt.plot(xb, yb, 'b-', label='Border')
plt.plot(xn, yn, 'go', label='Nodes')
plt.xlim([-.1 * l1, 1.1 * l1])
plt.ylim([-.1 * l2, 1.1 * l2])
