Yes, the second argument to

evalexpr,

*xp*, can be an array [n_sdim, n_points] for any points (where points outside the geometry will evaluate to 'nan').

EDIT: Evaluation in mesh nodes is faster, as points in the middle of cells are evaluated with direct finite element interpolation, but "evalexpr" accepts points everywhere.