AN EFFICIENT NUMERICAL IMPLEMENTATION FOR LOCALLY EMBEDDED STRONG DISCONTINUITIES