Reformulate penalty factor
Description
The penalty factor \beta / h
currently depends on the length of the adjacent fiber segment parts at the interface. We discussed that it should only depend on the local cell geometry
Proposal
Reformulate the penalty factor to only depend on the local geometry:
h = \min\left(|T_1|, |T_2|\right)/ |F|