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) (100,) > size(x) (100,2)
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!