Dennis Gläser authored Feb 09, 2019 and Timo Koch committed Feb 28, 2019

This way we can formulate the flux dependent on both the matrix as
well as the facet permeability. This significantly reduces the condition
number of the resulting matrix and thus improves newton convergence.