Coefficients in contact 3D suddenly infinity

[jorgensd]

CI started failing May 17th:

  coverage run --rcfile=.coveragerc python/demos/demo_contact_3D.py --gmsh --theta 0 --timing

which has coefficients equal to inf.

Setting theta to be small (0.000001) resolves it.