Fix a bug when splitting parcels

Since H is defined on local stereographic coordinate system, different
parcels have different coordinate system. After the parcel is split, H
of new parcels should not constructed directly from old U, V of old H,
but from old S (shrinked) and a major axis vertex.
Found a bug related to shape function and remapping weights

Split the parcel into two parcels, but the original parcel is reserved.
Comment out the parcel merging code for later testing.
Comment out the skeleton reset code for later testing.
The initial parcel shape size affects the solid rotation test results.
Finish the preliminary parcel splitting and merging algorithm

After a long time struggling, the scheme finally got on the right track!
Use ConfigManager in geomtk to set parameters