Bug fixes and GNumAdj/GAdj/Gmap

I have implemented and tested the creation of GNumAdj/GAdj/Gmap with the agreed-upon contract that:

GNumAdj(n) is the number of neighbors (not including self) for local index n
GAdj(:,n) is an array of global indices of neighbors for local index n
Gmap(n) is the global index of local index n

The above are 1-based.
The above contain no references to inactive elements.