I’m dealing with the 3d Elastodynamics problem in the examples, where we have a 3D displacement vector field. Obvioulsly, if we have `N`

points in the mesh, then the solution vector `u`

will be of length `3N`

.

My question is about the ordering in the following solution vector accessed via:

```
u.vector().get_local(vertex_to_dof_map(self.V))
```

my guess is that the vector is ordered as `uh = [uh_x, uh_y, uh_z]`

, where each `uh_x`

is of length `N`

. So that if we access the `i-th`

item in the coordinates array, then the associated solution vector shoud be `[uh_x[i], uh_y[i], uh_z[i]]`

, is that reasonable?