Hi,
Following my question and Dokken’s answer in another topic, it looks like you are missing a pair of parenthesis in
uD.interpolate(fem.Expression(ufl.real(x0), V.element.interpolation_points, dtype=np.complex128))
Just write V.element.interpolation_points()
and it should work.