MSP问题解法的并行化研究.docVIP

MSP问题解法的并行化研究 　　摘 要：为提高MSP问题的多项式时间算法ZH算法的计算速度，使其能够进行更大规模多级图的测试，本文对ZH算法进行了性能分析与并行的可行性评估，针对ZH算法中循环体较多的特点，分别在巨型机和普通PC机上进行MSP问题求解算法的首次并行化实践，随之对并行化算法提出优化方法，在实验中取得了较高的加速比。 　　关键词：MSP；并行化；巨型机 　　中图分类号：TP301.6 文献标识码：A 　　1 引 言 　　NP完全问题（NP complete）是指对任意的NP问题，都能在多项式时间内转化成为的问题，NPC问题一直以来都是计算科学领域的经典问题，具有很高的理论价值和应用价值，也吸引了很多学者的关注。文献[1]中提出了一个新的NP问题――多级图的简单路径问题（MSP问题），在文献中作者提出了MSP问题的一个多项式时间算法ZH算法，并且将哈密顿图问题归结到MSP问题上，从而证明了MSP问题是一个NP完全问题。大量的证明工作和测试显示，ZH算法是一个多项式时间的算法。但是由于该算法的复杂性的多项式的阶较高，导致运算需要耗费大量的时间和空间，普通的PC机难以支持较大规模多级图的测试。 　　有鉴于此，本文考虑利用超级计算机的强大计算能力对ZH算法进行并行化加速，并且在PC机上进行ZH算法的并行化研究。本文余下的本分按照如下的结构安排：第二节介绍MSP问题以及并行化研究的相关工作；第三节中本文提出了一种一般化的多项式时间算法的并行化思路，并对MSP问题的解决算法进行了并行化，之后提出了ZH算法并行化方案的优化方法；第四节分两个部分展示了ZH算法在超级计算机和普通PC机上的并行化实验结果，并对实验结果进行了分析；最后一节对本文工作进行了总结，并针对可能改进的方面提出了新的展望。 　　2 相关工作 　　2.1 MSP问题及相关工作 　　在文献[1]中，作者提出了一个新的NP完全问题，即MSP问题。该问题可简述如下： 　　对任意给定的一个单源单汇的多级图G，其中除源点外的其它所有顶点都附加有一个由图的边的构成的子集，问G中是否存在自源点到汇点的一条路径，使得对于路径上除源点外的所有顶点v，v附加的边集都包含该路径自源点至v的部分。 　　图1 一个多级图的例子 　　图1给出了一个多级图的例子。文献[1]找出了求解MSP问题的多项式算法ZH算法并给予了证明。MSP问题具有强大的表达能力，作者提出了一种解决MSP问题的算法，之后将图的撕裂或压缩的等价变形与归纳法进行了结合，证明了ZH算法是一个多项式时间算法。ZH算法主要是通过循环的方式对多级图中每一级上每个节点的边集和每条边的可达路径进行修改和约束，这个过程中多级图的边集信息将会进行收敛，直至整个多级图的信息不再发生变化。这些针对边集和路径所做的约束保证了多级图中的简单路径存在与否。ZH算法中的约束是一种复杂的耦合关系，某一个顶点的边集或者某一条边的可达路径都是依赖于多级图中其他节点的信息的，而这种约束反过来也会影响到多级图中其他节点和边信息的约束。 　　近几年，樊硕等人[4]将更多的NPC问题归结到了MSP问题上，并给予了证明，比如SAT问题、子集和问题等，扩展了MSP问题算法的应用空间。我们将由归结得到的SAT问题的多项式时间解法MSP-SAT进行了一定规模的测试，取得了较好的效果，但是受限于PC机的运算速度，没有进行更大规模的测试。 　　2.2 并行化方法介绍与相关工作 　　最近十年，中国的超级计算机取得了长足的发展与进步，超级计算机有着广泛的应用价值，已经与国计民生密不可分，其应用领域包括气候预测、生命科学、核物理模拟等等。文献[6]对高性能计算机的并行计算进行了理论方面的研究。[5]则将超算应用到了雷达信号的处理上，[7]利用超级计算机解决了卫星的重力测量问题，由此可见，超级计算机正在被应用到越来越重要的领域。在普通计算机或者是超级计算机的一个节点内，openMP库是进行算法并行化必须要考虑的工具，openMP是对并行算法的抽象，为程序员提供了简便地方法对程序进行并行化改造，openMP库提供的抽象使得程序员能够很大程度上忽略底层的硬件信息，从而投入更多地精力到并行算法本身，而不必过多关心具体实现。同时，openMP库能够部分地解决并行算法中的负载平衡，线程冲突等问题，为算法的并行化带来了很多便利。例如文献[12]中，作者利用openMP平台对语音识别算法进行了并行化，取得了可观的效率改进；白丽泽[10]等人则基于openMP平台改进了分子动力学算法。在文献[13]中，作者利用openMP对SAT问题进行了并行化改造。 　　3 MSP问题的求解算法的并行化研究 　　3.1 多项式时间算法的并行化思路 　　任务的并行是物理世