们正在做大量电镜照片的二维傅立叶运算处理。
祁旻找了个没人的位置坐下，看到旁边一个戴着森海塞尔耳机的小伙子正在调试代码。
她不禁好奇地瞥了一眼，发现他正在写的是一个RNA结构预测工具。其实RNA结构预测已经是研究比较充分的领域了，现有的算法各有千秋，然而对于某类特定的RNA，或许需要根据其已知的信息而设计合适的特殊算法。
而这段代码看起来就像是这个小伙子自己写的，因为祁旻已经看他修修改改快十分钟了，输出的预测数组仍然是一堆负数。
碰巧祁旻原先写过RNA结构预测的代码，此时顺着他修改看了一遍，很快就找到了问题。
“嗯……这个遍历方式好像不太对吧。”祁旻不好意思直接对他说，只是坐在旁边嘀咕道。
一开始这小伙子戴着耳机装没听见，压根儿没搭理她。
祁旻看他这么调真是费劲，于是忍不住说道：“对于最近邻居算法(i，n-i)应该从i大于等于2开始，i是0和1的情况要用氢键最多原理……”
“是‘能量最低’原理。”调代码的小伙子摘下耳机，转身对她说道，“但是初测值的实验方式不同，这两种算法压根儿不能通用。”
“谁说的？第一对氢键赋值根本不影响。你先试试吧。”祁旻对他说道。
小伙子将信将疑地加上了一段儿IF，从另外一个单纯使用能量最低原理的算法里复制出了相关代码。
他再运行之后，输出的矩阵里果然不全是负数了。
但很明显结果也不对。
负数还是太多了，导致最大值才只有不到4，而很明显用于测试的那段30个碱基的线性RNA不可能只形成一个碱基配对。
“这又是怎么回事儿？”小伙子这回直接对她问道。
祁旻凑到屏幕前仔细看了看，而后一边敲着键盘一边说道：“你这个对比最大值的嵌套FOR语句里没有加FLAG，只要加上就好了……”
代码再运行了一遍，这回输出的矩阵终于看起来正常了。
小伙子松了口气儿，刚要从椅子上站起来，祁旻却说道：“哎，这还是不对。你看这段RNA的已知结构，这里面有三个LOOP，但你的算法算出来有四个……哦，我知道了，还得加一个防止无效配对的IF，连续两个碱基以下的配对不应该存在。”
她把这个IF语句加进去之后，最终画出的RNA二级结构终于
本章未完，请点击下一页继续阅读！