getAlpha1D Module Function

  1. smod_legendre_1D.f90
  2. smod_legendre_1D
  3. getAlpha1D

function getAlpha1D(N) result(alpha)

Arguments

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

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


Calls

proc~~getalpha1d~~CallsGraph proc~getalpha1d getAlpha1D proc~linsolve_quick linsolve_quick proc~getalpha1d->proc~linsolve_quick proc~eye eye proc~getalpha1d->proc~eye proc~linsolve linsolve proc~linsolve_quick->proc~linsolve dgetrf dgetrf proc~linsolve->dgetrf dgetrs dgetrs proc~linsolve->dgetrs

Contents

Source Code

getAlpha1D

Source Code

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

    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)    :: x

      x = getx(N)

      do ii = 1,N
        A(ii,:) = getArow1D(N, x(ii))
      enddo

      return
    end function getA
  end function getAlpha1D