Feature/multidomain with soldep params

Numerical differentation together with solution-dependent spatial parameters didn't work properly, because the flux variables cache was not reset adequately after the deflection of primary variables. This introduces the necessary resets of the flux variables caches such that at the beginning of each derivative computation it is in its undeflected state.