The procedure to refine the internal edges is similar to the one for the external edges. Here too routines edgedivide.f and newnodes.f are called. The new nodes are perturbed similarly as for the external edges and reordered in an aleatoric way. Just the insertion of the nodes in cavity.f is more simple:

guido dhondt 2018-12-15