def METHOD_NAME():
    """Build and test lagrangian ds"""
    q0[...] = [1, 2, 3]
    v0[...] = [4, 5, 6]
    mass = np.asarray(np.diag([1, 2, 3]), dtype=np.float64)
    ds = sk.LagrangianDS(q0, v0, mass)
    ec = ds.computeKineticEnergy()
    assert ec == 87.0
    assert ds.dimension() == ndof
    assert np.allclose(ds.mass(), mass)