第88章 Zeta函数非平凡零点的计算残卷(2 / 3)
“必须要用斯特林公式进行渐近展开,然后分离实部和虚部,接著还要计算超越数π
和自然对数i的高精度小数值,最后还要做三角函数的泰勒级数展开。”
“而且这还只是算一个点!”
“为了捕捉到符號改变的瞬间,还必须在t=14到t=15之间密集取点。”
“每一次取点,都要重复上面那一长串极其噁心的纯四则运算——”
“哪怕中间有一个小数点进位算错,前面的力气全部白费!”
於是不服气的李东服气了。
不过他毕竟是群主,还是挺不要脸的——
“不对啊,我跟他较什么劲?”
“我现在是21世纪啊!我有计算机啊!让人脑去干计算器的活,这不是纯纯的脑瘫吗?
“,他立刻打开那台联想笔记本电脑,调出python的pycharm界面。
他凭著脑子里的基础编程知识,直接把最基础的欧拉—麦克劳林求和公式翻译成了代码。
为了保证精度,他还调用了python的高精度十进位库,强制保留了25位有效数字,通过对t密集取点的暴力扫描方式寻找零点。
“走起!”
屏幕上开始滚动出零点的数值,100个——————500个————1000个————
但很快,李东就发现了问题。
笔记本的散热风扇开始发疯狂转动,键盘表面也开始烫手。
屏幕上代码的运行速度渐渐的慢了下来。
他写的代码完全没有做內存优化,每一次计算的所有中间变量、歷史取点数据都被一股脑的塞进了列表里,没有任何释放机制。
高精度运算带来的海量临时对象疯狂堆积,windows任务管理器里,python进程的內存占用率从4g一下子就来到了14g!
当零点计算推进到4120个时,画面彻底卡死。
滑鼠变成了圈圈,怎么点都没反应。
紧接著屏幕变成了蓝色。
【终止代码:memory—management】
李东傻傻的看著屏幕。
內心直呼“傻唄联想。”
完全忘记了自己的代码也跟粪一样。
不过这也让李东更加的震撼。
“就算我的电脑不行,代码也有待优化,可是这毕竟是距离黎曼150多年的硅基產物呀。”
“黎曼是怎么硬生生手算出前1104个零点的?” ↑返回顶部↑
和自然对数i的高精度小数值,最后还要做三角函数的泰勒级数展开。”
“而且这还只是算一个点!”
“为了捕捉到符號改变的瞬间,还必须在t=14到t=15之间密集取点。”
“每一次取点,都要重复上面那一长串极其噁心的纯四则运算——”
“哪怕中间有一个小数点进位算错,前面的力气全部白费!”
於是不服气的李东服气了。
不过他毕竟是群主,还是挺不要脸的——
“不对啊,我跟他较什么劲?”
“我现在是21世纪啊!我有计算机啊!让人脑去干计算器的活,这不是纯纯的脑瘫吗?
“,他立刻打开那台联想笔记本电脑,调出python的pycharm界面。
他凭著脑子里的基础编程知识,直接把最基础的欧拉—麦克劳林求和公式翻译成了代码。
为了保证精度,他还调用了python的高精度十进位库,强制保留了25位有效数字,通过对t密集取点的暴力扫描方式寻找零点。
“走起!”
屏幕上开始滚动出零点的数值,100个——————500个————1000个————
但很快,李东就发现了问题。
笔记本的散热风扇开始发疯狂转动,键盘表面也开始烫手。
屏幕上代码的运行速度渐渐的慢了下来。
他写的代码完全没有做內存优化,每一次计算的所有中间变量、歷史取点数据都被一股脑的塞进了列表里,没有任何释放机制。
高精度运算带来的海量临时对象疯狂堆积,windows任务管理器里,python进程的內存占用率从4g一下子就来到了14g!
当零点计算推进到4120个时,画面彻底卡死。
滑鼠变成了圈圈,怎么点都没反应。
紧接著屏幕变成了蓝色。
【终止代码:memory—management】
李东傻傻的看著屏幕。
內心直呼“傻唄联想。”
完全忘记了自己的代码也跟粪一样。
不过这也让李东更加的震撼。
“就算我的电脑不行,代码也有待优化,可是这毕竟是距离黎曼150多年的硅基產物呀。”
“黎曼是怎么硬生生手算出前1104个零点的?” ↑返回顶部↑