秦飛有些好為人師的搖了搖頭,索性也不看風景了,就這麼饒有興致的盯著李東的螢幕,權當是在旅途中找點樂子了。
此時,李東也沒注意旁邊有個高手在,他只最佳化自己的演算法。
秦飛看著,先是暗自點頭。
「嗯,不錯。」
「知道拋棄原始的尤拉乘積公式,直接上Rienn—Siegelforla。」
「主和項的迴圈寫得很規範,餘項展開的泰勒截斷也算合理。」
「時間複雜度標準的0(t^(1/2)),對於一個普通學生來說,這程式碼底子算得上紮實了。」
秦飛在心裡給出了一個「中規中矩」的評價。
然而,就在秦飛以為李東接下來會順理成章的去寫多執行緒加速程式碼時。
李東卻突然停了下來。
然後就按下了退格鍵,竟然把剛剛寫好的,那段堪稱教科書般標準的黎曼一西格爾主和項迴圈計算程式碼————
刪掉了一大半!
秦飛愣了一下。
「他幹嘛?寫錯撤回了?」
但接下來李東敲出的程式碼,卻讓秦飛徹底看迷糊了。
李東沒有去糾結如何加快單次計算的速度。
他有【程式碼直覺】還有《黎曼絕筆》裡的那種降維數學思維,他現在已經可以初步找到通往現代計算機底層邏輯的完美橋樑了。
他在程式碼頂部,重新引入了一個名為&It;fftw3。h>的標頭檔案。
這是著名的快速傅立葉變換(FFT)開源庫。
然後又在程式碼中定義一個密集的評估網格。
voidltiPointEvaluation Grid(。。。)
他將需要逐個點代入t值去硬算的黎曼Zeta函式,透過構建區域性多項式,強行轉化為了在網格點上的泰勒級數展開。
緊接著,李東呼叫了FFT演算法,讓這些多項式在頻域內進行高速相乘。
秦飛已經把降噪耳取了下來,眼睛看著李東的電腦螢幕,呼吸都變得有些急促了。
「等等————這特麼是在幹什麼?」
秦飛想用自己學過的知識去去解析這套程式碼背後的邏輯。
「算黎曼Zeta函式,你引個快速傅立葉變換(FFT)進來幹嘛?這兩者八竿子打不著啊!」
「還有這個多點求值矩陣?他把泰勒展開放在網格上————」
」?算運式項多的行並個一了包打行強,值求點零有所的間區大個一把,特積卷效高的TFF用利在是他!點單算在是不他————他「麼什了到想像好他,然突
」?銷開算計的怖恐那項和主式公爾格西—曼黎攤分去值求點多用「
。路思的算計點零的現發界學世後被未沒,裡稿手在寫曼黎是就實其這
。度維個一整整了下往又度雜複間時的者後把,上層底學數在至甚,紀世半個一近了早法算演egahnohcS—okzyldO的開公來後比
!啊道知能可不飛秦但
。限下論理的算計點零單是就,度雜複))2/1(^t(0的來帶式公爾格西—曼黎,裡知認的他在
。架框學數個這出跳能可不本,化佳最限極的化程工是只也,法算演S0的快最界學是怕哪
?嗎構架算計個整翻推要想是?麼什幹在底到他,伙傢個這前眼,在現而
」!?能可麼怎這「
。了瘋他覺飛秦
」!的炸級數指呈會對絕差誤斷截的數點浮,搞麼這你?制控麼怎差誤但,度雜複均平低降能然雖上論理,算運陣矩開展項和主把「
。完駁反裡心在飛秦等沒還,而然
!界邊學數的小極個一了在鎖的死死差誤斷截的徑路分積將的完,碼式程波濾值的式函限帶段幾了穿,間格網的值求點多在續繼東李
。雅優很的用都,用複的憶記次一每,走遊的標指行一每
」!現實銷開憶記的低麼這。碼式程的麼這用能可不也,法算演S0是算就————對不「
」!?法算演的來哪是底到這————路思的析分值數代現是不本,輯邏值個這他「
。了呆看底徹飛秦
。了上不跟全完經已他,輯邏法算演的起一在合組但,法語的++C些那識認他然雖
」?吧寫瞎在是會不他「
」?嗎通跑以可的真————寫樣這「
。中疑懷我自的深深了陷飛秦








