Refactoring is possible if new Matrix class functions applied
For few matrix elements, NumPy and OpenFOAM yielded few-decimal point differences
//Info<< "Mx - b:" << nl << (Mx - source) << endl; affects the subsequent program states despite encapsulation. Memory leak might be the reason, will be checked