# Gridap.jl solution at mesh nodes?

**URL:** https://discourse.julialang.org/t/gridap-jl-solution-at-mesh-nodes/60993
**Category:** General Usage
**Tags:** question, gridap
**Created:** [May 12, 2021, 2:41am UTC](https://discourse.julialang.org/t/gridap-jl-solution-at-mesh-nodes/60993 "2021-05-12T02:41:03Z")
**Posts on this page:** 1
**Page:** 1

<div class="post-metadata">

### Author: ![kc104](https://avatars.discourse-cdn.com/v4/letter/k/ed8c4c/32.png) [@kc104](https://discourse.julialang.org/u/kc104)
#### Post date: [May 12, 2021, 2:41am UTC](https://discourse.julialang.org/t/gridap-jl-solution-at-mesh-nodes/60993/1 "2021-05-12T02:41:03Z")

</div>

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:

```julia
> 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!