Substitute own GaussAlgorithm implementation by Eigen version

Enable vectorial output for secondary variables.

this is already done at start of iteration

The least squares extrapolation proceeds component-wise and produces
vector-valued fields.

Now time, current solution and the process DOF table are passed to the
extrapolation routine, enabling more versatile nodal output
possiblities.

[Mat/Fracture] support 3d in MC model