47 implicit none
48
49
50
51 integer, intent(in) :: nelem
52 integer, intent(in) :: connectivity(max_nb_nodes_per_elt,nelem)
53 integer, intent(in) :: nb_nodes_mini(nelem)
54 integer, intent(in) :: e1
55 integer, intent(in) :: e2
56
57
58
59 logical :: bool
60 integer :: count
61 integer :: i,j
62
63
64
65 bool = .false.
66 count = 0
67
68 i = 1
69 j = 1
70 do while (i <= max_nb_nodes_per_elt .and. j <= max_nb_nodes_per_elt)
71 if(connectivity(i,e1) == 0 .or. connectivity(j,e2) == 0) then
72 exit
73 elseif (connectivity(i,e1) == connectivity(j,e2)) then
74 count = count + 1
75 i = i + 1
76 j = j + 1
77 else if (connectivity(i,e1) > connectivity(j,e2)) then
78 i = i + 1
79 else
80 j = j + 1
81 end if
82 end do
83
84 if (count >= nb_nodes_mini(e1) .or. count >= nb_nodes_mini(e2)) then
85 bool = .true.
86 end if