Gridap.jl solution at mesh nodes?

I’m using Gridap.jl for some PDE work and would like to access the solution at the mesh nodes + the mesh coordinates in the form of a standard julia array (the end goal is to interpolate the solution at a new set of coordinates which don’t align with the mesh).

To clarify I’m hoping for something like the following:

> uh = solve(..)
> u_array = get_array(uh)
> x = get_coordinates_array(trian)
> size(u_array)
> size(x)

Any help would be appreciated. (As a temporary solution I’ve been writing to a VTK file and then reading the output to an array which isn’t ideal).

PS: Really enjoying working with a PDE solver built from the ground up in Julia, big thanks to the creators!