type_crs_add Subroutine

private subroutine type_crs_add(alpha, A, B, C)

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in) :: alpha
type(type_crs), intent(in) :: A
type(type_crs), intent(in) :: B
type(type_crs), intent(inout) :: C

Called by

proc~~type_crs_add~~CalledByGraph proc~type_crs_add type_crs_add interface~add add interface~add->proc~type_crs_add

Source Code

    subroutine type_crs_add(alpha, A, B, C)
        ! C := alpha*A + B
        !
        ! [ATTENTION] Assumes A, B, and C have the exact same sparsity pattern.
        !
        implicit none
        real(real64), intent(in) :: alpha
        type(type_crs), intent(in) :: A
        type(type_crs), intent(in) :: B
        type(type_crs), intent(inout) :: C

        integer(int32) :: i

        !$omp parallel do
        do i = 1, A%nnz
            C%val(i) = alpha * A%val(i) + B%val(i)
        end do
        !$omp end parallel do

    end subroutine type_crs_add