块预处理orthomin方法的程序实现方法

？在现代油藏数值模拟中,由于全隐式方法、井的隐式处理方法、混合网格嵌套、局部网格加密、热采问题等造成传导率差别很大,从而使方程组ax=b中的矩阵a严重病态和严重的非对称.国外一些文献已证明许多传统的迭代法,像sor、sip和adip等依赖于矩阵性质和迭代参数的求解方法,对此类问题已经无效,而预处理orthomin方法是目前解决不对称问题有效的迭代方法之一.