第329章 需要緩緩
韓柏岩聲音都在發顫：「肖教授，你……你怎麼算出我們沒公開的原始數據的？」
肖宿眨了眨眼。
他不太理解韓柏岩和林瑾為什麼突然露出這種表情。
「推出來的。」他說。
林瑾微微前傾了身體，語氣比剛才又鄭重了幾分：
「你是說，你看到這些模型的輸出結果，就能反推出實驗的原始輸入數據？」
「這不是很明顯嗎？
你們雖然改了一些數值，但是這些數值和整體運算邏輯不符，它們在同一個耦合矩陣里和其他變量的縮放比例根本不一致，只要簡單驗算一下，就能知道原來的數值是什麼了。」
韓柏岩嘴角抽了一下。
簡單驗算？
開兩台伺服器，跑好幾個鐘頭才能還原出原始數據的近似值，這樣的簡單驗算嗎？
那可是幾萬核時超算跑出來的數據啊，涉及幾十個耦合變量，每一個變量之間的縮放關係還都經過了好幾輪疊代。
而肖宿就這麼直接心算出來了？！
韓柏岩覺得自己這把年紀了，什麼大風大浪沒見過，但此刻，還是生出了一股深深的無力感。
緩了好一會兒，他才從這股無力感里掙脫出來。
然後他忽然想到了另一件事。
「肖教授，既然你能根據輸出數據反推出原始數據，那其他實驗室公布的那些實驗數據，你是不是也能看出來他們的實際數值和計算方法？」
肖宿皺了皺眉，不太理解韓柏岩為什麼忽然把話題跳到這裡。
不是在討論飛行器再入段熱流計算模型的修改方案嗎？
怎麼又扯到其他實驗數據了？
但林瑾幾乎是立馬就明白了韓柏岩的用意。
他的眼皮跳了一下，身體微微前傾。
如果肖宿真的有這種能力，能看到別人論文裡被處理過的數據背後藏著什麼，那這件事的價值，可能比他們今天來請教的所有技術問題加起來都要大。
他接過韓柏岩的話頭，把問題又往前推了一步：
「肖教授，你說我們改過的數值跟整體運算邏輯不符，所以一眼就能看出來。
那同樣的道理，其他實驗室在論文裡公布的數據，如果他們被修飾過的數據和計算方法，你是不是也能看出來？」
肖宿點了點頭，簡單地「嗯」了一聲。
他從很早之前就發現這點了。
很多頂刊的論文，順著他們的計算公式往回推，他們報的數值和他們實際實驗得到的數值是根本不一樣的。
有的是把誤差棒畫短了，有的則是把對照組的數據挑著報，還有的是在數據處理過程中悄悄換了一種統計方法，讓P值剛好卡在零點零五以下。
更離譜的是有些論文裡公開的所謂最優計算模型，肖宿在腦子裡跑了一遍，就發現他們連基本的數學工具都用錯了，用錯了還不說，發出來的結果倒是漂漂亮亮的，而大家好像都看不見，那種論文竟然還拿了獎，上了封面，被同行引用了無數次。
從那時候起，肖宿就養成了驗算的習慣。
每看一篇有數據的論文，他都會無意識地順著公式往回推一遍，看看他們報出來的數據和他們的計算框架是不是自洽的。
然後他發現，不自洽的論文比他預想的要多得多。
韓柏岩和林瑾對視了一眼。
那一秒鐘的眼神交換里包含了太多東西，有震驚，有警覺，還有一種獵人忽然發現自己手裡多了一把新槍的興奮。
韓柏岩站起來，聲音壓低了幾分：「肖教授，能不能再請你幫我們看一個數據？」
韓柏岩的表情很認真，好像肖宿今天要是不幫他們，他們是不走了一樣。
「什麼數據？」
「稍等一下，數據我沒帶在身上，我出去打個電話，很快回來。」
他推門出去打電話的時候，書房裡安靜了下來。
肖宿拿起筆，重新翻開了那份被改得密密麻麻的材料。
他把材料翻到飛行器再入段熱流計算那部分，在被圈掉的內容旁邊，重新寫了一個簡潔的計算框架。
比起回答他們的問題，他覺得還是直接改一個模型更簡單一些。
肖宿一旦開始思考問題，渾身的氣場都變得，直接自成一個世界。
林瑾坐在一旁，看著肖宿專心在材料上寫寫畫畫，心情一時有些難言。
書房裡就這麼安靜了幾分鐘，直到韓柏岩推門進來。
「我跟助理說完了，他馬上把材料送過來，大概二十分鐘左右。」
林瑾點點頭。
兩人就這麼安靜坐著，看著肖宿在紙上寫完所有內容。
直到肖宿落筆，林槿才開口問道：
「肖教授，之前拜託你幫忙設計的那個計算模型，現在進展怎麼樣了？」
肖宿站起身，走到旁邊的柜子前，從裡面取出了之前林瑾給他的小行李箱。
「在這裡。」
林瑾眼睛亮了一下，當場打開筆記本，看了看那個模型。
突然，他像是想到了什麼，問道：
「那這個模型，能不能也用在霍爾推進器的磁場優化、超燃衝壓發動機的燃燒室流場分析，還有飛行器的再入段熱流計算上呢？」
肖宿想了想，開口道：
「計算模型從本質上說都是通用的，不管你是做推進器磁場優化還是做發動機流場分析，底層的數學結構都是耦合系統在模式空間裡的譜分解。
不同的應用場景之間，差別的不是數學骨架，而是輸入參數的數量級和耦合關係的複雜度。」
他指了指那台筆記本：「所以從理論上說，這個模型也能算那些問題，只是會比較慢而已。
它當初不是針對那些問題做的，參數接口和輸出格式都不一樣，跑一輪要花不少時間在數據轉換上。」
林瑾追問道：「那有沒有辦法讓它快起來？」
肖宿看了他一眼。
「可以，只要把輸入層的多物理場耦合矩陣做一次預處理，用並行QR分解把耦合矩陣的稀疏結構提取出來，然後在傳熱和流體兩個模塊之間加一層自適應時間步長調整模塊，這樣就可以了。
這兩步做完之後，模型的效率大概能提七八倍，和專門針對那些問題重新寫一個模型差別不大了。」
他說得很輕鬆，但林瑾聽到並行QR分解和自適應時間步長調整模塊這兩個詞的時候，頭都大了。
並行QR分解在計算數學裡不算罕見，但要把一個多物理場強耦合的矩陣，先做並行分解再提取稀疏結構，這個操作的數學複雜度光是想想就讓人頭疼。
更別說自適應時間步長調整了，那可是數值計算領域公認的一個大分支，多少博士生苦熬五年，就為了在這個方向上往前推那麼一點點。
結果到了肖宿嘴裡，聽起來就那麼回事兒。
韓柏岩和林瑾很自然的安靜了下來，沒再繼續追問細節。
他們今天受到的衝擊已經夠多了，需要緩緩。
肖宿眨了眨眼。
他不太理解韓柏岩和林瑾為什麼突然露出這種表情。
「推出來的。」他說。
林瑾微微前傾了身體，語氣比剛才又鄭重了幾分：
「你是說，你看到這些模型的輸出結果，就能反推出實驗的原始輸入數據？」
「這不是很明顯嗎？
你們雖然改了一些數值，但是這些數值和整體運算邏輯不符，它們在同一個耦合矩陣里和其他變量的縮放比例根本不一致，只要簡單驗算一下，就能知道原來的數值是什麼了。」
韓柏岩嘴角抽了一下。
簡單驗算？
開兩台伺服器，跑好幾個鐘頭才能還原出原始數據的近似值，這樣的簡單驗算嗎？
那可是幾萬核時超算跑出來的數據啊，涉及幾十個耦合變量，每一個變量之間的縮放關係還都經過了好幾輪疊代。
而肖宿就這麼直接心算出來了？！
韓柏岩覺得自己這把年紀了，什麼大風大浪沒見過，但此刻，還是生出了一股深深的無力感。
緩了好一會兒，他才從這股無力感里掙脫出來。
然後他忽然想到了另一件事。
「肖教授，既然你能根據輸出數據反推出原始數據，那其他實驗室公布的那些實驗數據，你是不是也能看出來他們的實際數值和計算方法？」
肖宿皺了皺眉，不太理解韓柏岩為什麼忽然把話題跳到這裡。
不是在討論飛行器再入段熱流計算模型的修改方案嗎？
怎麼又扯到其他實驗數據了？
但林瑾幾乎是立馬就明白了韓柏岩的用意。
他的眼皮跳了一下，身體微微前傾。
如果肖宿真的有這種能力，能看到別人論文裡被處理過的數據背後藏著什麼，那這件事的價值，可能比他們今天來請教的所有技術問題加起來都要大。
他接過韓柏岩的話頭，把問題又往前推了一步：
「肖教授，你說我們改過的數值跟整體運算邏輯不符，所以一眼就能看出來。
那同樣的道理，其他實驗室在論文裡公布的數據，如果他們被修飾過的數據和計算方法，你是不是也能看出來？」
肖宿點了點頭，簡單地「嗯」了一聲。
他從很早之前就發現這點了。
很多頂刊的論文，順著他們的計算公式往回推，他們報的數值和他們實際實驗得到的數值是根本不一樣的。
有的是把誤差棒畫短了，有的則是把對照組的數據挑著報，還有的是在數據處理過程中悄悄換了一種統計方法，讓P值剛好卡在零點零五以下。
更離譜的是有些論文裡公開的所謂最優計算模型，肖宿在腦子裡跑了一遍，就發現他們連基本的數學工具都用錯了，用錯了還不說，發出來的結果倒是漂漂亮亮的，而大家好像都看不見，那種論文竟然還拿了獎，上了封面，被同行引用了無數次。
從那時候起，肖宿就養成了驗算的習慣。
每看一篇有數據的論文，他都會無意識地順著公式往回推一遍，看看他們報出來的數據和他們的計算框架是不是自洽的。
然後他發現，不自洽的論文比他預想的要多得多。
韓柏岩和林瑾對視了一眼。
那一秒鐘的眼神交換里包含了太多東西，有震驚，有警覺，還有一種獵人忽然發現自己手裡多了一把新槍的興奮。
韓柏岩站起來，聲音壓低了幾分：「肖教授，能不能再請你幫我們看一個數據？」
韓柏岩的表情很認真，好像肖宿今天要是不幫他們，他們是不走了一樣。
「什麼數據？」
「稍等一下，數據我沒帶在身上，我出去打個電話，很快回來。」
他推門出去打電話的時候，書房裡安靜了下來。
肖宿拿起筆，重新翻開了那份被改得密密麻麻的材料。
他把材料翻到飛行器再入段熱流計算那部分，在被圈掉的內容旁邊，重新寫了一個簡潔的計算框架。
比起回答他們的問題，他覺得還是直接改一個模型更簡單一些。
肖宿一旦開始思考問題，渾身的氣場都變得，直接自成一個世界。
林瑾坐在一旁，看著肖宿專心在材料上寫寫畫畫，心情一時有些難言。
書房裡就這麼安靜了幾分鐘，直到韓柏岩推門進來。
「我跟助理說完了，他馬上把材料送過來，大概二十分鐘左右。」
林瑾點點頭。
兩人就這麼安靜坐著，看著肖宿在紙上寫完所有內容。
直到肖宿落筆，林槿才開口問道：
「肖教授，之前拜託你幫忙設計的那個計算模型，現在進展怎麼樣了？」
肖宿站起身，走到旁邊的柜子前，從裡面取出了之前林瑾給他的小行李箱。
「在這裡。」
林瑾眼睛亮了一下，當場打開筆記本，看了看那個模型。
突然，他像是想到了什麼，問道：
「那這個模型，能不能也用在霍爾推進器的磁場優化、超燃衝壓發動機的燃燒室流場分析，還有飛行器的再入段熱流計算上呢？」
肖宿想了想，開口道：
「計算模型從本質上說都是通用的，不管你是做推進器磁場優化還是做發動機流場分析，底層的數學結構都是耦合系統在模式空間裡的譜分解。
不同的應用場景之間，差別的不是數學骨架，而是輸入參數的數量級和耦合關係的複雜度。」
他指了指那台筆記本：「所以從理論上說，這個模型也能算那些問題，只是會比較慢而已。
它當初不是針對那些問題做的，參數接口和輸出格式都不一樣，跑一輪要花不少時間在數據轉換上。」
林瑾追問道：「那有沒有辦法讓它快起來？」
肖宿看了他一眼。
「可以，只要把輸入層的多物理場耦合矩陣做一次預處理，用並行QR分解把耦合矩陣的稀疏結構提取出來，然後在傳熱和流體兩個模塊之間加一層自適應時間步長調整模塊，這樣就可以了。
這兩步做完之後，模型的效率大概能提七八倍，和專門針對那些問題重新寫一個模型差別不大了。」
他說得很輕鬆，但林瑾聽到並行QR分解和自適應時間步長調整模塊這兩個詞的時候，頭都大了。
並行QR分解在計算數學裡不算罕見，但要把一個多物理場強耦合的矩陣，先做並行分解再提取稀疏結構，這個操作的數學複雜度光是想想就讓人頭疼。
更別說自適應時間步長調整了，那可是數值計算領域公認的一個大分支，多少博士生苦熬五年，就為了在這個方向上往前推那麼一點點。
結果到了肖宿嘴裡，聽起來就那麼回事兒。
韓柏岩和林瑾很自然的安靜了下來，沒再繼續追問細節。
他們今天受到的衝擊已經夠多了，需要緩緩。