Calculate_SOR_preconditioner Function

public function Calculate_SOR_preconditioner(A, omega, alpha) result(L)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in), dimension(:, :) :: A
real(kind=dp), intent(in) :: omega
real(kind=dp), intent(in) :: alpha

Return Value real(kind=dp), dimension(size(A, 1), size(A, 2))


Calls

proc~~calculate_sor_preconditioner~~CallsGraph proc~calculate_sor_preconditioner Calculate_SOR_preconditioner proc~diag Diag proc~calculate_sor_preconditioner->proc~diag

Source Code

    function Calculate_SOR_preconditioner(A, omega, alpha) result(L)
        real(dp), dimension(:, :), intent(in) :: A
        real(dp), intent(in) :: omega, alpha
        real(dp), dimension(size(A, 1), size(A, 2)) :: L
        integer :: N, i

        N = size(A, 1)

        L = 0.d0

        if (any(Diag(A) < TOL_CONVERGENCE_dp)) stop "ERROR :: Zero diagonal in SOR preconditioner"
        do i = 1, size(A, 1)
            L(i, i) = 1.d0 / omega * A(i, i)
            L(i, 1:i - 1) = A(i, 1:i - 1)
        end do

        L = alpha * L

    end function Calculate_SOR_preconditioner