VerticalEdgeBasis
Evaluate basis functions on left and right edge of biunit Quadrangle.
Interface
- ܀ Interface
- ️܀ See example
- ↢
INTERFACE
MODULE PURE FUNCTION VerticalEdgeBasis_Quadrangle(qe1, qe2, x, y) &
& RESULT(ans)
INTEGER(I4B), INTENT(IN) :: qe1
!! order on left vertical edge (e1), it should be greater than 1
!! It should be greater than 2
INTEGER(I4B), INTENT(IN) :: qe2
!! order on right vertical edge(e2), it should be greater than 1
!! It should be greater than 2
REAL(DFP), INTENT(IN) :: x(:), y(:)
!! point of evaluation
!! these points should be between [-1, 1].
REAL(DFP) :: ans(SIZE(x), qe1 + qe2 - 2)
END FUNCTION VerticalEdgeBasis_Quadrangle
END INTERFACE
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