A High Order Method for Determining the Edges in the Gradient of a Function. Commun. Comput. Phys. 2018, 5 (2-4), 694-711.