令为三角形,令上的平滑函数。TTffTT
我们可以使用中点求积,其中的中点。∫fdx≈|T|⋅f(xM)∫fdx≈|T|⋅f(xM)xMxMTT
你能给我提供(参考)单纯形的高阶公式吗?
拥有最有效的规则是否重要?您是否要求所有权重都是正数和/或所有点都在该区域内?两个通用资源:
如果您正在寻找一些实现(在 Fortran 中),那么我可以推荐这个页面。在那里,您可以找到针对不同几何形状(如三角形、单纯形等)的多种积分规则。当然,您还可以在那里找到对原始论文的参考。
经典的三角形论文是 Dunavant,Symmetrical Gaussian Quadrature Rules,1985 的论文。您的规则最高为 20 阶。可以在此处找到实现。