def test_quadrilateral():
quadpy.quadrilateral.integrate(
lambda x: numpy.exp(x[0]),
quadpy.quadrilateral.rectangle_points([0.0, 1.0], [0.0, 1.0]),
quadpy.quadrilateral.Stroud('C2 5-4')
)
val = quadpy.quadrilateral.integrate(
lambda x: [numpy.exp(x[0]), numpy.exp(x[1])],
numpy.stack([
quadpy.quadrilateral.rectangle_points([0.0, 1.0], [0.0, 1.0]),
quadpy.quadrilateral.rectangle_points([0.0, 1.0], [0.0, 1.0]),
quadpy.quadrilateral.rectangle_points([0.0, 1.0], [0.0, 1.0]),
quadpy.quadrilateral.rectangle_points([0.0, 1.0], [0.0, 1.0]),
quadpy.quadrilateral.rectangle_points([0.0, 1.0], [0.0, 1.0]),
], axis=-2),
quadpy.quadrilateral.Stroud('C2 3-1')
)
assert val.shape == (2, 5)
return
评论列表
文章目录