My program boundary condition is also u = 0 on the boundary; The default setting of tgv is tgv = 1e30, but in fact, according to the obtained finite element solution, the value on the boundary is not equal to 0.
