getAlpha2D Module Function

function getAlpha2D(N) result(alpha)

Arguments

Type IntentOptional AttributesName
integer, intent(in) :: N

Return Value real(kind=wp), dimension(N,N)


Calls

proc~~getalpha2d~~CallsGraph proc~getalpha2d getAlpha2D proc~linsolve_quick linsolve_quick proc~getalpha2d->proc~linsolve_quick proc~eye eye proc~getalpha2d->proc~eye proc~linsolve linsolve proc~linsolve_quick->proc~linsolve dgetrf dgetrf proc~linsolve->dgetrf dgetrs dgetrs proc~linsolve->dgetrs

Contents

Source Code


Source Code

module function getAlpha2D(N) result(alpha)
    integer, intent(in)       :: N
    real(wp), dimension(N,N)  :: alpha

    integer                   :: ii
    real(wp), dimension(N,N)  :: A, B

    A = getA(N)
    B = eye(N)

    call linsolve_quick(N, A, N, B, alpha)

    return
contains
    function getA(N) result(A)
        integer, intent(in)       :: N
        real(wp), dimension(N,N)  :: A

        integer                   :: ii
        real(wp), dimension(N,2)  :: xy

        xy = getXY_2D(N)

        do ii = 1,N
            A(ii,:) = getArow(N, xy(ii,1), xy(ii,2))
        enddo

        return
    end function getA
end function getAlpha2D