lpolyn.f

来自「网络带宽测试工具」· F 代码 · 共 34 行

      Subroutine lpolyn( n, nel, m, indx, rowp, matvals, u, v, gamma )
! ----------------------------------------------------------------------
! --- Left polynomial preconditioner of order 'm'.
!     
! ----------------------------------------------------------------------
      Use         numerics
      Use         floptime
      Implicit    None

      Integer  :: n, nel, m
      Integer  :: indx(nel), rowp(n+1)
      Real(l_) :: matvals(nel)
      Real(l_) :: u(n), v(n)
      Real(l_) :: gamma(m+1)

      Real(l_) :: w(n)
      Integer  :: i, il
! ----------------------------------------------------------------------
!omp parallel do
      Do il = 1, n
         v(il) = gamma(m+1)*u(il)
      End Do
      Do i = 1, m
         Call spmxv( n, nel, indx, rowp, matvals, v, w )
!omp parallel do
         Do il = 1, n
            v(il) = v(il) - w(il)
            v(il) = v(il) + gamma(m-i+1)*u(il)
         End Do
      End Do
      flops = flops + n*( 1 + 3*m )
! ----------------------------------------------------------------------
      End Subroutine lpolyn