made many last changes, including hidden reference arguments in Moperator, etc. code compiles in both C and C++, but does not converge in C