Only compute derivatives if computing Jacobians
Merge pull request #104 from pbauman/boundary_condition_fixes

Only compute derivatives if computing Jacobians