Solution adaptivity using a triangular mesh