Fixing du_dot_du for non-nodal FEs (closes #3784)

The computation of du_dot_du assumed nodal FE and was putting the
computed value into each quad-point. This did not work for non-nodal FEs
like monomials. Those require the "function" to be projected into the
space.  Fortunately, all of our time scheme have du_dot_du as a scalar
value, thus, we do not need the solution vector and we can store it as a
single value.

Also, adding a test case that will check the Jacobian by FDing it (on
monomial FEs).
Add checking to getParam<>()
refs #3934
Fixup objects that were retrieving before checking
closes #3934
Removing unused object from Modules
Merge pull request #3941 from permcody/check_param_3934

Check param 3934
Merge pull request #3948 from andrsd/fix-du-dot-du

Fixing du_dot_du for non-nodal FEs (closes #3784)