“这不可能!”刘凯忍不住出声质疑,“这里的非线性项是强耦合的,怎么可能直接线性化?你这是在乱写!”
徐辰手中的笔没有停,头也不抬地回了一句:
“在切空间上,它是线性的。你没学过微分流形吗?”
“……”
刘凯被这一句话噎得差点背过气去。
微分流形?那是数学系研究生的课,他一个搞计算物理的,哪懂那个?
……
而张乐阳,作为博三的资深学长,虽然也有些跟不上徐辰的跳跃思维,但他能看懂大方向。
越看,他越心惊。
“这……这是蛙跳算法的高阶改进版?!”
“他竟然引入了一个自適应的时间步长,利用局部曲率来动態调整积分精度?!”
“天哪!这个思路……好厉害!既保证了能量守恆,又避开了那些刚性区域的陷阱!”
张乐阳感觉自己的头皮一阵阵发麻。
“他真的是一个大一学生吗?”
……
五分钟后。
徐辰停笔。
草稿纸上,已经写满了一整页的数学公式。
“好了。”
徐辰把纸往旁边一放,双手放在了键盘上。
他转过头,看了一眼刚才还在质疑他的刘凯,嘴角勾起一抹淡淡的笑意,半开玩笑地说道:
“公式我都推导好了,逻辑都在这儿。要不……代码你来写?”
刘凯看著那张写满了希腊字母和微分算子的草稿纸,脸涨得通红。
写?写个屁啊!
他连那个哈密顿方程的变分形式都没看懂,更別说把它翻译成代码了!
他只能尷尬地低下头,装作没听见。
“看来还是得我自己来。”徐辰笑了笑,並没有继续为难他,转过身面对屏幕。
噼里啪啦的键盘敲击声,如同密集的雨点般响起。
徐辰没有看文档,没有查谷歌,甚至连停顿都没有。
那些复杂的张量运算、梯度反向传播、辛积分步骤,就像早已刻在他脑子里一样,行云流水般地流淌到了屏幕上。
defhamiltonian_dynamics(q,p,step_size,num_steps):
。。。
p-=step_size*grad_u(q)2