Bernd Flemisch authored 4 years ago and Timo Koch committed 4 years ago

So far, the `LinearPDESolver` (re)assembles the matrix in its
`solve` function unconditionally. This might be unnecessary if the
matrix did not change compared to the previous call to `solve`.
Make a first attempt to avoid the unnecessary calculation by
adding a corresponding control flag.

Feature/box scvf num outside scvs

See merge request !2469

Fix/pnm gridgeo

See merge request !2482

[md] Add multidomain convergence writer

See merge request !2475

Feature/pore network model

See merge request !2208

Timo Koch authored 4 years ago and Kilian Weishaupt committed 4 years ago

The standard library distributions use implementation defined algorithms
which make them produce different random numbers on different library version
(different compiler versions or systems).
We use here now our own implementation which should have a high enough
quality for this purpose and it makes the tests reproducible.

[box] In elementsolution.hh use Dune::ReservedVector instead of...

Closes #989

See merge request !2479

[box] In elementsolution.hh use Dune::ReservedVector instead of Dune::Blockvector to have a fixed maximum size at compile time.

Feature/upwind interface

See merge request !2414

Fix/intersecting entities

See merge request !2274