                  edge = edgeX;
                    elseif dir == 2
                        edge = edgeY;
                    end
            
                    % check if the edge belong to the boundary
                    if norm(edge-ref) < 1.0e-12
                        
                        % Compute the length of the edge
                        dx = edgeX(2) - edgeX(1);
                        dy = edgeY(2) - edgeY(1);
                        l = sqrt(dx*dx + dy*dy);
            
              