addind gray matter mask to repo
bug fix: changed np.min to min
BUG FIX: unscaled dual residual norm was returned, did not impact on convergence but on adaptive rho