第四章:超算中心#-《类脑体》
第(2/3)页
她不禁好奇地瞥了一眼,发现他正在写的是一个RNA结构预测工具。其实RNA结构预测已经是研究比较充分的领域了,现有的算法各有千秋,然而对于某类特定的RNA,或许需要根据其已知的信息而设计合适的特殊算法。
而这段代码看起来就像是这个小伙子自己写的,因为祁旻已经看他修修改改快十分钟了,输出的预测数组仍然是一堆负数。
碰巧祁旻原先写过RNA结构预测的代码,此时顺着他修改看了一遍,很快就找到了问题。
“嗯……这个遍历方式好像不太对吧。”祁旻不好意思直接对他说,只是坐在旁边嘀咕道。
一开始这小伙子戴着耳机装没听见,压根儿没搭理她。
祁旻看他这么调真是费劲,于是忍不住说道:“对于最近邻居算法(i,n-i)应该从i大于等于2开始,i是0和1的情况要用氢键最多原理……”
“是‘能量最低’原理。”调代码的小伙子摘下耳机,转身对她说道,“但是初测值的实验方式不同,这两种算法压根儿不能通用。”
“谁说的?第一对氢键赋值根本不影响。你先试试吧。”祁旻对他说道。
小伙子将信将疑地加上了一段儿IF,从另外一个单纯使用能量最低原理的算法里复制出了相关代码。
他再运行之后,输出的矩阵里果然不全是负数了。
但很明显结果也不对。
负数还是太多了,导致最大值才只有不到4,而很明显用于测试的那段30个碱基的线性RNA不可能只形成一个碱基配对。
“这又是怎么回事儿?”小伙子这回直接对她问道。
祁旻凑到屏幕前仔细看了看,而后一边敲着键盘一边说道:“你这个对比最大值的嵌套FOR语句里没有加FLAG,只要加上就好了……”
代码再运行了一遍,这回输出的矩阵终于看起来正常了。
第(2/3)页