我想知道哪种方法是数值积分勒让德多项式的最佳方法。我正在构建一个不连续的 Galerkin CFD 代码,其中勒让德多项式用作光谱投影中的基函数。
在以前的代码中,我习惯于使用 Gauss-Legendre 数值求积,其中只有允许获得度多项式的精确积分的点. 现在,由于这样的交点是勒让德多项式的根,我需要使用积分点-th Legendre 多项式(否则多项式在其根处求值,所有高阶积分为零)。还有其他有效的方法吗?
非常感谢您的帮助。
我想知道哪种方法是数值积分勒让德多项式的最佳方法。我正在构建一个不连续的 Galerkin CFD 代码,其中勒让德多项式用作光谱投影中的基函数。
在以前的代码中,我习惯于使用 Gauss-Legendre 数值求积,其中只有允许获得度多项式的精确积分的点. 现在,由于这样的交点是勒让德多项式的根,我需要使用积分点-th Legendre 多项式(否则多项式在其根处求值,所有高阶积分为零)。还有其他有效的方法吗?
非常感谢您的帮助。
一般而言,您为有限元空间使用哪个基对于求积无关紧要。例如,如果您使用次数多项式(无论是拉格朗日基、勒让德基、单项式还是其他任何东西)那么质量矩阵的被积函数具有多项式次数,你需要一个高斯正交点正确整合。这是因为这样的高斯求积可以积分任何次数的多项式确切地说,无论您实际上如何根据基数表示该多项式。