Matlab简单代码mp-quadrature-用于生成通用1D、2D和3D正交规则的多精度算法

1. 引言

在许多数值分析领域中，高阶正交规则（例如Gauss-Legendre，Gauss-Jacobi，Gauss-Lobatto等）的精确计算和列表化至关重要。标准的双精度算术通常仅足以获得14（或更少）个点和权重的精度，因此需要多精度代数库来改善这种情况。尽管用于计算正交规则的标准技术已经有一段时间了，但是某些方法在计算任意精度规则方面比其他方法更好。在这里，我们基于免费提供的GMP，MPFR和GMPFRXX库收集了（希望增长）多种算法，用于生成正交规则。该代码用于将有限元库中的一些一维正交规则制成表格。

2. 安装

要构建库，请键入 ./configure 和 make。您必须同时安装GMP和MPFR库才能构建mp-quadrature库。至少有两个选择：

• 运行包含的 build_gmp_mpfr.sh 脚本。这会将GMP和MPFR从源代码下载、构建并安装到 ./gmp 和 ./mpfr 目录中。然后，配置脚本将自动找到这些脚本。
• 使用以下选项来配置系统的GMP和MPFR安装位置：

--with-gmp-include=/path/to/gmp/include 等参数来指定库的路径。