HorizontalEdgeBasis

Evaluate basis functions on left and right edge of biunit Quadrangle.

Interface

܀ Interface

INTERFACE
  MODULE PURE FUNCTION HorizontalEdgeBasis_Quadrangle(pe3, pe4, x, y) &
    & RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: pe3
    !! order on bottom vertical edge (e3), it should be greater than 1
    INTEGER(I4B), INTENT(IN) :: pe4
    !! order on top vertical edge(e4), it should be greater than 1
    REAL(DFP), INTENT(IN) :: x(:), y(:)
    !! point of evaluation
    REAL(DFP) :: ans(SIZE(x), pe3 + pe4 - 2)
  END FUNCTION HorizontalEdgeBasis_Quadrangle
END INTERFACE

️܀ See example

program main
  use easifembase
  use easifemClasses
  implicit none
  real(dfp), allocatable :: xij(:,:), avec(:)
  integer(i4b) :: ii, jj, cnt, n
  real(dfp), allocatable :: ans(:,:)
  integer(i4b) :: order, pe1, pe2, pe3
  type( VTKPlot_ ) :: aplot
  character(len=*), parameter :: fname="./results/"

  n = 51
  call reallocate(avec, n)
  call reallocate(xij, 2, int((n+1)*(n+2)/2))
  avec= linspace(0.0_DFP, 1.0_DFP, n)
  cnt=0
  do ii = 1, n
    do jj = 1, n-ii+1
      cnt=cnt+1
      xij(1,cnt) = avec(ii)
      xij(2,cnt) = avec(jj)
    end do
  end do

  order=3; pe1=3; pe2=3; pe3=3

  ans = EdgeBasis_Triangle(&
    & pe1=pe1, pe2=pe2, pe3=pe3, &
    & xij=xij, refTriangle="UNIT")

  do ii  = 1, size(ans,2)
    call aplot%scatter3D(x=xij(1,:), y=xij(2, :), z=ans(:,ii), &
      & filename=fname//"EdgeBasis( " // &
        & ", pe1="//tostring(pe1) // &
        & ", pe2="//tostring(pe2) // &
        & ", pe3="//tostring(pe3) // &
        & " )" // tostring(ii) // &
        & ".vtp", &
      & label="P")
  end do

end program main

↢