Triangle

Quadrature points on unit triangle

PROGRAM main
USE easifemBase
IMPLICIT NONE
TYPE(QuadraturePoint_) :: obj
TYPE(ReferenceTriangle_) :: refelem
INTEGER(I4B) :: order
refelem = ReferenceTriangle(nsd=2_I4B)
order = 4_I4B
CALL initiate(obj=obj, &
  & refelem=refelem, &
  & order=order,  &
  & quadratureType=GaussLegendre)
CALL display(obj, "ans: ")
END PROGRAM main

note

In the case of triangle quadratureType is ineffective, that is, we do not use quadratureType.

See results

x1	9.15762E-02	0.81685	9.15762E-02	0.44595	0.44595	0.1081
x2	9.15762E-02	9.15762E-02	0.81685	0.1081	0.44595	0.44595
w	5.49759E-02	5.49759E-02	5.49759E-02	0.11169	0.11169	0.11169

Quadrature points on biunit triangle

PROGRAM main
USE easifemBase
IMPLICIT NONE
TYPE(QuadraturePoint_) :: obj
TYPE(ReferenceTriangle_) :: refelem
INTEGER(I4B) :: order
refelem = ReferenceTriangle(nsd=2_I4B, xij=RefCoord_Triangle("BIUNIT"))
order = 4_I4B
CALL Initiate(obj=obj, &
  & refelem=refelem, &
  & order=order,  &
  & quadratureType=GaussLegendre)
CALL Display(mdencode(obj) , "")
END PROGRAM main

See results

x1	-0.81685	0.6337	-0.81685	-0.1081	-0.1081	-0.78379
x2	-0.81685	-0.81685	0.6337	-0.78379	-0.1081	-0.1081
w	0.2199	0.2199	0.2199	0.44676	0.44676	0.44676