Replace numpy by par loop in construction of map from facets to 4 dofs on facet

Final step in replacing the numpy code by par loops in the construction of redirection maps in the BDFM1 lumped mass matrix.