How to get components of vector solution in Gridap

Specific Domains Numerics
,
1

Dear Julia users,

I’m using Gridap, and I’m learning from the tutorials. I have some questions:

a) How can we recover x- or y-components from a numerical solution uh, as the solution to linear elasticity (like in the tutorials). Note that

typeof(uh)
Gridap.FESpaces.SingleFieldFEFunction{true}

there are get_free_values(uh) and get_dirichlet_values(uh) throwing Array{Float64,1} values, but I’m puzzled about getting something like ux=uh[1], uy=uh[2] (sorry for the incorrect syntax in the latter, it is just a way to express it).

b) Does anyone tried to solve anisotropic elasticity using Gridap? Some advice?

Thanks.

1 Like