l}(1:nLinNodes), this.ELEM{el}(1)];

                <span class="comment">% Loop through the edges of the element</span>
                <span class="keyword">for</span> j = 1:nEdges

                    <span class="comment">% coordinates of the edge</span>
                    edgeX = [cX(j) , cX(j+1)];
                    edgeY = [cY(j) , cY(j+1)];

                    <span class="comment">% select the edge</span>
                    <span class="keyword">if</span> dir == 1
                        edge 