Added lots of debug information and fixed errors in matrix forming code
(missing absolute values around locel and lnods mostly.
Fixed incorrect nxtk update for interior midpoint nodes in GETCOLIDX.
Added additional solution analysis to help with debugging.
Added even more debugging output!