第四章：超算中心#

第(2/3)页

体系上进行操作的“湿实验”很少，没有什么有害试剂和微生物的污染，因此用不着把实验台跟办公桌分离开。而同样与生物信息学特征相对应的是，实验室里每个人的办公桌上都立着至少两台屏幕和大得异常的主机，房间里空调开得很大，祁旻从其中一台没有人但仍然常亮的屏幕上看到他们正在做大量电镜照片的二维傅立叶运算处理。

　　祁旻找了个没人的位置坐下，看到旁边一个戴着森海塞尔耳机的小伙子正在调试代码。

　　她不禁好奇地瞥了一眼，发现他正在写的是一个RNA结构预测工具。其实RNA结构预测已经是研究比较充分的领域了，现有的算法各有千秋，然而对于某类特定的RNA，或许需要根据其已知的信息而设计合适的特殊算法。

　　而这段代码看起来就像是这个小伙子自己写的，因为祁旻已经看他修修改改快十分钟了，输出的预测数组仍然是一堆负数。

　　碰巧祁旻原先写过RNA结构预测的代码，此时顺着他修改看了一遍，很快就找到了问题。

　　“嗯……这个遍历方式好像不太对吧。”祁旻不好意思直接对他说，只是坐在旁边嘀咕道。

　　一开始这小伙子戴着耳机装没听见，压根儿没搭理她。

　　祁旻看他这么调真是费劲，于是忍不住说道：“对于最近邻居算法(i，n-i)应该从i大于等于2开始，i是0和1的情况要用氢键最多原理……”

　　“是‘能量最低’原理。”调代码的小伙子摘下耳机，转身对她说道，“但是初测值的实验方式不同，这两种算法压根儿不能通用。”

　　“谁说的？第一对氢键赋值根本不影响。你先试试吧。”祁旻对他说道。

　　小伙子将信将疑地加上了一段儿IF，从另外一个单纯使用能量最低原理的算法里复制出了相关代码。

　　他再运行之后，输出的矩阵里果然不全是负数了。

　　但很明显结果也不对。

　　负数还是太多了，导致最大值才只有不到4，而很明显用于测试的那段30个碱基的线性RNA不可能只形成一个碱基配对。

　　“这又是怎么回事儿？”小伙子这回直接对她问道。

　　

(本章未完,请翻页)
记住手机版网址：m.lw00.net