Handle unused boundary nodes correctly
For "unused" boundary nodes the special value -1 is used. Mapping
this through isBoundaryNode[-1] would result in an out-of-bounds
memory access and invalid values in the vertices_c_style array.
Closes: #54 (closed), #55 (closed)