[论文解读] Documentation for the ratpoints program

Michael Stoll|arXiv (Cornell University)|Mar 21, 2008

Advanced Numerical Analysis Techniques被引用 3

一句话总结

本文提出 ratpoints，一个高度优化的 C 程序，用于计算双曲椭圆曲线上指定高度界限内的所有有理点。它使用二次筛法结合位级并行计算与多精度算术，通过 256 位 AVX 寄存器实现高达 2 倍的性能提升，使其成为枚举曲线上有理点最高效的工具之一。

ABSTRACT

This note explains how to obtain, install, and use the ratpoints program. The program finds rational points up to a specified height on hyperelliptic curves using a highly optimized quadratic sieving algorithm.

研究动机与目标

• 开发一种高效算法，用于在指定高度界限内计算双曲椭圆曲线上所有有理点。
• 通过位操作和移位算术最小化昂贵操作（如整数除法）以优化性能。
• 通过 GMP 库支持高次多项式（最高 100 次）和任意精度系数。
• 利用现代 CPU 指令集（SSE、AVX）进行位级并行操作，加速筛法处理。
• 为数论学家和算术几何研究人员提供一个稳健、可扩展且免费开源的工具。

提出的方法

• 采用二次筛法，通过检查模小素数的局部可解性，高效排除候选分母。
• 使用位数组存储筛法信息，替代较慢的字符数组，实现紧凑且高效的运算。
• 通过 GMP 库实现多精度算术，确保点验证的精确性并处理大系数。
• 基于预期成功率选择素数，以减少不必要的计算。
• 在筛法过程中使用 128 位（SSE）和 256 位（AVX）寄存器并行处理位域，显著提升吞吐量。
• 应用多项式反转等变换，当反转后的曲线能提供更好过滤效果时，提升筛法效率。

实验结果

研究问题

• RQ1 如何使双曲椭圆曲线上有理点的枚举在大高度界限下更加高效？
• RQ2 低级算术与内存访问模式的哪些优化能显著降低有理点搜索的计算成本？
• RQ3 现代 CPU 指令集（SSE、AVX）在数论算法的位级筛法中能实现多大程度的加速？
• RQ4 ratpoints 算法的性能如何随曲线次数和系数大小的增加而扩展？
• RQ5 多项式反转在提升非首一或病态条件曲线筛法有效性方面起到什么作用？

主要发现

• 与之前的 128 位 SSE 版本相比，使用 256 位 AVX 寄存器时，ratpoints 程序的性能最高可提升 2 倍。
• 得益于 AVX2 支持，256 位寄存器的使用带来了显著的性能提升，尤其在具有宽向量单元的现代处理器上表现突出。
• 集成 256 位运算需要在 rp-private.h 中通过宏进行仔细的抽象，以管理不同的字长和指令集。
• 程序现在支持最高 100 次的多项式，次数上限可通过头文件中的 RATPOINTS_MAX_DEGREE 进行配置。
• 新版本（2.2）包含性能改进，如位数组中更快的零值检查，以及在首次筛法阶段对 AVX/SSE 寄存器更高效的使用。
• 通过广泛的测试（包括 testbase 和 testbase2 的回归检查）保持了正确性和鲁棒性，并通过 p-进数和雅克比符号条件处理了非首一首项系数等边界情况。