FEniCS:自定义正交规则

计算科学 正交 芬尼克斯
2021-12-06 03:29:35

对于我的 PDE 在二维三角网格上的反应项的数值积分,我想使用方案 其中是三角形节点。

TΦ(x)|T|3(Φ(x0)+Φ(x1)+Φ(x2))
xi

assemble在ing 系统时如何指定这个?

2个回答

实际上,您没有为汇编程序指定正交规则,而是为形成编译器指定正交规则。您可以ffc/quadrature_schemes.py根据需要进行调整,也可以调整生成的代码。在这两种情况下,您只需提供正交点和权重。

您可以通过设置调整正交度

parameters["form_compiler_parameters"]["quadrature_degree"] = q

这将改变正交次数,使得 q 次数的多项式被精确积分。

FEniCS 将尝试自动检测必要的正交度,因此通常不需要。但是,在某些情况下,您可能希望通过设置较低的 q 值来进行欠积分。