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.