Groundwater Flow Process

This is based on #521 and must not be merged before.

The commits (starting with [[PL] Add GlobalSetup as template to GW process.](https://github.com/endJunction/ogs/commit/3bb51a7976b678798dcad7d18fff75d95a1a2e51)) fills in the initialize() and the solve() functions in the GroundwaterFlowProcess. The creation of the process is documented in the commit messages.

One think needs special attention because of incompatibility between the MathLib::DenseMatrix and Eigen matrices used: see the [[ML] Remove size check, because eigen has no getNRows()](https://github.com/endJunction/ogs/commit/aa7c84541900c375c165da2cf2a9df5fa8b8d72b) commit for details.

For those who wonder where the boundary conditions are applied, please don't: these are going into another PR.
As a preliminary implementation,  ``_hydraulic_conductivity`` could be an identity matrix multiplied with a constant number of conductivity for isotropic material. For assembly of the Laplacian, we can consider an option of
 - Isotropic conductivity: take scalar number like what presented here.
- Anisotropic conductivity: take tensor/matrix