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)