Assume that the domain split into two subregions 1 and 2. The exact solution on each region is u1 and u2, respectively. You can draw the exact solution on the domain by using the characteristic functions. For examples,
Code:
# Draw exact solution
char_1 = GridFunction(L2(mesh,order=0))
char_1.Set(CoefficientFunction([1,0]))
char_2 = GridFunction(L2(mesh,order=0))
char_2.Set(CoefficientFunction([0,1]))
u = char_1*u1+char_2*u2
Draw(u, mesh, "uexact")