overhauled MirrorOpt to load first and second thresholds (will need to generate them), fixed bug that prevented MoperatorJacobian preallocation from working correctly for MPI single process (np =1)