Added MC73 ordering routine
- Needs the F90 library and modules that I could not convert
Updated to new KSPSetOperators() interface