Interpolation nodes in 2D

import matplotlib.pyplot as pt
import numpy as np
import modepy as mp

nodes = mp.warp_and_blend_nodes(2, 10)
pt.plot(nodes[0], nodes[1], "x")

tri = np.array([(-1, -1), (1, -1), (-1, 1), (-1, -1)]).T
pt.plot(nodes[0], nodes[1], "x")
pt.plot(tri[0], tri[1], "-b")
pt.gca().set_aspect("equal")