A C++ implementation of the "double obstacle phase field algorithm" for forced curvature flow