12 LOGICAL,
OPTIONAL,
INTENT(IN) :: edge_stab
14 IF (present(edge_stab))
THEN
15 mesh%edge_stab=edge_stab
17 mesh%edge_stab=.false.
20 nws =
SIZE(mesh%jjs,1)
21 IF (nw==3 .AND. nws==2)
THEN
23 ELSE IF (nw==6 .AND. nws==3)
THEN
27 ELSE IF (nw==0 .AND. nws==0)
THEN
33 ELSE IF (nw==4 .AND. nws==3)
THEN
34 WRITE(*,*)
' Finite element not yet programmed '
37 ELSE IF (nw==10 .AND. nws==6)
THEN
38 WRITE(*,*)
' Finite element not yet programmed '
42 WRITE(*,*)
' Finite element not yet programmed '
subroutine, public gauss_points_2d_p2(mesh)
subroutine, public gauss_points_2d_p1(mesh)
Prepares gauss points for P1 approximation in 2D.
subroutine, public gen_gauss(mesh, edge_stab)