对于我的 PDE 在二维三角网格上的反应项的数值积分,我想使用方案 其中是三角形节点。
assemble
在ing 系统时如何指定这个?
对于我的 PDE 在二维三角网格上的反应项的数值积分,我想使用方案 其中是三角形节点。
assemble
在ing 系统时如何指定这个?
实际上,您没有为汇编程序指定正交规则,而是为形成编译器指定正交规则。您可以ffc/quadrature_schemes.py
根据需要进行调整,也可以调整生成的代码。在这两种情况下,您只需提供正交点和权重。
您可以通过设置调整正交度
parameters["form_compiler_parameters"]["quadrature_degree"] = q
这将改变正交次数,使得 q 次数的多项式被精确积分。
FEniCS 将尝试自动检测必要的正交度,因此通常不需要。但是,在某些情况下,您可能希望通过设置较低的 q 值来进行欠积分。