1. 引言

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

2. 安装

要构建库,请键入 ./configuremake。您必须同时安装GMPMPFR库才能构建mp-quadrature库。至少有两个选择:

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

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