Fix J Integral adjacent to free surfaces ref #3898

If a crack intersects a free surface, we need to zero out the q function
on that surface for all crack front nodes except the one on that surface.
Add an option to list a set of free surfaces intersected by the crack
to do this.
Merge pull request #3899 from bwspenc/jint_fix

Fix J Integral adjacent to free surfaces
Merge commit '8e1b518bb080b2a9e80e0b170286f89900b42df2'
Fix the comparison in PetscJacobianTester (closes #3912).

It was missing the type cast, so the result was wrong even if the
numbers did not satisfy the relation.