def test_lobatto():
f = lambda x: 12 * x ** 11
F = lambda x: x ** 12
exact = F(1) - F(0)
x, w = lobatto_quad(7)
x, w = map_pts_wts(x, w, 0, 1)
assert(len(x) == 7)
est = np.sum(f(x) * w)
np.testing.assert_almost_equal(exact, est)
def test_lobatto():
f = lambda x: 12 * x**11
F = lambda x: x**12
exact = F(1) - F(0)
x, w = lobatto_quad(7)
x, w = map_pts_wts(x, w, 0, 1)
assert (len(x) == 7)
est = np.sum(f(x) * w)
np.testing.assert_almost_equal(exact, est)
def gll_basis(degree):
""" A basis from the Gauss-Lobatto-Lagrange nodes """
x, w = lobatto_quad(degree + 1)
nodes, w = map_pts_wts(x, w, 0, 1)
return basis_from_nodes(nodes)
def lobatto(N):
x, w = qc.lobatto_quad(N)
x, w = qc.map_pts_wts(x, w, 0.0, 1.0)
return QuadratureInfo(0.0, x, w)
