Function test_approximate2D

tests/test_nurbs.py:347–389

(lam, penalty)

@mark.parametrize(("penalty", "lam"), [(3, 1e-6), (4, 1e-6), (5, 1e-6), (2, 0)])
def test_approximate2D(lam, penalty):

    t = torus(5, 1).face()

    # double periodic surface
    us = np.linspace(0, 1, endpoint=False)
    vs = np.linspace(0, 1, endpoint=False)

    pts = np.array(
        [t.positionAt(u * 2 * np.pi, v * 2 * np.pi).toTuple() for v in vs for u in us]
    )

    surf = approximate2D(
        pts,
        us[None, :].repeat(len(vs), 0).ravel(),
        vs[:, None].repeat(len(us), 1).ravel(),
        3,
        3,
        15,
        15,
        uperiodic=True,
        vperiodic=True,
        penalty=penalty,
        lam=lam,
    )

    f = surf.face()

    # general sanity checks
    assert f.isValid()
    assert f.Area() == approx(t.Area(), rel=1e-3)

    # check the stability of the parameters
    us0 = us[None, :].repeat(len(vs), 0).ravel()
    vs0 = vs[:, None].repeat(len(us), 1).ravel()

    us2, vs2 = f.params(array2vec(pts))

    delta_u = np.array(us2 - us0)
    delta_v = np.array(vs2 - vs0)

    assert np.allclose(np.where(delta_u >= 1, delta_u, 0) % 1, 0)
    assert np.allclose(np.where(delta_v >= 1, delta_v, 0) % 1, 0)


EDGES = (