开放在网格上……”

突然，他好像想到了什么

“他……他不是在算单点！他是在利用fft的高效卷积特性，把一个大区间内的所有零点求值，强行打包成了一个并行的多项式运算？”

“用多点求值去分摊黎曼-西格尔公式主和项那恐怖的计算开销？”

这其实就是黎曼写在手稿里，没未被后世学界发现的零点计算的思路。

比后来公开的odlyzko–schonhage算法早了近一个半世纪，甚至在数学底层上，把后者的时间复杂度又往下压了整整一个维度。

但秦飞不可能知道啊！

在他的认知里，黎曼-西格尔公式带来的o(t(1/2))复杂度，就是单零点计算的理论下限。

哪怕是学界最快的os算法，也只是工程化的极限优化，根本不可能跳出这个数学框架。

而现在，眼前这个家伙，他到底在干什么？是想要推翻整个计算架构吗？

“这怎么可能？！”

秦飞感觉他疯了。

“把主和项展开成矩阵运算，理论上虽然能降低平均复杂度，但误差怎么控制？你这么搞，浮点数的截断误差绝对会呈指数级爆炸的！”

然而，还没等秦飞在心里反驳完。

李东继续在多点求值的网格间，穿插了几段带限函数的插值滤波代码，完美的将积分路径的截断误差死死的锁在了一个极小的数学边界内！

每一行指针的游走，每一次内存的复用，都用的很优雅。

“不对……就算是os算法，也不可能用这么少的代码、这么低的内存开销实现！”

“他这个插值逻辑，根本不是现代数值分析的思路……这到底是哪来的算法？！”

秦飞彻底看呆了。

虽然他认识那些c++的语法，但组合在一起的算法逻辑，他已经完全跟不上了。

“他不会是在瞎写吧？”

“这样写……真的可以跑通吗？”

秦飞陷入了深深的自我怀疑中。