Added a subroutine to modify the cmc matrix using this: (D^-0.5 CMC_petsc D^-0.5) D^0.5 X = D^-0.5 rhs_p. The obtained matrix has a lower condition number but in the cases I have run so far the results are worse. However, it might be useful in the future...
Merge branch 'fluidity-core/MultiFluids_Dev' of github.com:jrper/MultiFluidFluidity into fluidity-core/MultiFluids_Dev