Fixing a corner case in the volume computation

At the border of the domain an edge has only one neighbour and so the
computation of the normal of this edge was not completed properly.