How to get nemrical solution value for a given point?

General
2

If you are using dolfin version 2018.1.0, then

u = Function(V)
print(u.vector().get_local())

is supported.
If your point does not match a mesh vertex, you need to call the function eval function, i.e.
print(u(Point(0.25))) to get the solution at that point.

2 Likes