def METHOD_NAME(self):
    base = np.random.random((1, 3))[0]
    # create a random axis. A constant value has been added to avoid [0,0,0]
    axis = np.random.random((1, 3))[0] + 0.01
    plane = BluemiraPlane(base, axis)
    lx = 20
    ly = 10
    bmface = plane.to_face(lx, ly)
    np.testing.assert_almost_equal(bmface.center_of_mass, plane.base)
    np.testing.assert_almost_equal(
        np.array([bmface.length, bmface.area]), np.array([2 * (lx + ly), lx * ly])
    )