Writing a constraint for a jump

Dear colleagues;
I have trouble creating a jump over an edge for my field, \psi, a P2 element.
I am solving, where q and \mu are known, and r is a multiplier:
Screenshot from 2022-04-23 21-45-11
Over a square, with an extracted circle from the center as follow:
Screenshot from 2022-04-23 22-06-26

I want a difference of \psi over the blue dashed line: i.e., \psi|_{below}=\psi|_{above}+2 \pi for instance, while \nabla \psi remains continuos on both sides (without cutting or any alterations to the domain).

What terms some I add to the

solve Mahdi( \psi , v ) = int2d(Th)( r*(dx(\psi)*dx(v) + dy(\psi)dy(v)) ) + int2d(Th)( ( \mu x-rqx )dx(v) + ( \mu y-rqy )*dy(v) )

Is there such a thing available in FreeFem for a P2 element?

Please accept my apologies for my weak knowledge.
Thanks for your resourceful thoughts in advance.

Writing as:
solve Mahdi(\psi , v )= int2d(Th)( r*(dx(\psi)*dx(v) + dy(\psi)dy(v)) ) + int2d(Th)( (\mu x-rqx )dx(v) + ( \mu y-rqy )dy(v) )
+
intalledges(Th)((label==1)
v * (2 *pi)*jump(\psi));

Where the label for the blue dashed line is 1, it didn’t solve the problem.