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)