def test_scheme(scheme, tol):
degree = check_degree(
lambda poly: quadpy.e1r.integrate(poly, scheme),
lambda k: math.factorial(k[0]),
1,
scheme.degree + 1,
tol=tol
)
assert degree == scheme.degree, \
'Observed: {} expected: {}'.format(degree, scheme.degree)
return
评论列表
文章目录