The background is that in the case of global caching, we need to pass
the grid variables as non-const in order to be able to deflect them
during numeric differentiation. It would be nice if we could find a way
around it.

[gridmanager] Add interface to query for existence of grid data

See merge request !2534

Feature/timestep methods

See merge request !2296

Based on initial draft by Timo Koch.

[md][newton] Fixup conflict in solutionChanged between !2291 and !2527

See merge request !2533

[propertysystem] Some cleanup

See merge request !2512

[newton] Fixup solution reset in line search update after !2291

See merge request !2532

Feature/newton linesearch

See merge request !2507

Bump package versions.

See merge request !2518

Fix/multidomain newton residual

See merge request !2527