B循环

采用易于并行化的循环归约算法,对块三对角矩阵进行最小填充因子分解。分解块的存储允许将逆运算应用到分解时可能未知的多个右侧。块行数的可伸缩性是通过循环缩减实现的,而块大小的可伸缩性是使用用于块矩阵操作的多线程例程(OpenMP、GotoBLAS)实现的。这种双重可伸缩性是这个新解算器的一个值得注意的特性,以及它能够有效地处理任意(非2的幂次)块行和处理器数量。与目前最先进的并行稀疏解算器进行了比较。预计这种新的求解器将允许许多物理应用程序优化使用当前超级计算机上的并行资源。举例说明了求解器在磁流体力学(MHD)和高温聚变等离子体三维平衡求解器中的应用。