Shifted van-Genuchten curve added.  Fixes #3350

This is for 2-phase flow.  It may make the Jacobian nonsingular in
2phase models when one phase disappears and residual saturation is
zero.  Time will tell whether this is useful.

Included are doco and tests
Shifted van-Genuchten curve added.  Fixes #3350

This is for 2-phase flow.  It makes the Jacobian nonsingular
in 2phase models when one phase disappears and the residual
saturation is zero.  In limited testing this appears to give
good convergence characteristics.

Included are doco and tests.

Also, I have changed RichardsSeff::dseff and d2seff so they
accept references to the answers, rather than providing the
answers as return values.  According to jwpeterson this might
speed up computations by not having to do extra copying.
Since I have so many RichardsSeff type of classes, this required
a lot of minor changes
How 'bout this fix, eh?
closes #3388
Merge pull request #3390 from permcody/max_parallel

How 'bout this fix, eh?
Merge commit '7bdc29057de3ceff0df67da8235e3df95eb0f0a1'
Merge pull request #3353 from WilkAndy/shiftedVG_3350

Shifted van-Genuchten curve added.  Fixes #3350
Merge commit 'ba7c8e5b8f40d4211f28c02c561f5a4cde3576d9'
Doco for VGshifted added into user.tex.

Fixes #3393

This document had fallen a few weeks behind the code
so was substantially revised to with my current thoughts.
Merge pull request #3395 from WilkAndy/vgshifted_doco_3393

Doco for VGshifted added into user.tex.
Merge commit '8bf28ba7f5ca53ae82a36c35754bc1dcb0f0198c'