ConnorMallon

Hey @amartinhuertas @fverdugo, seeing as multifield autodiff is fixed it would be nice to also have this functionality. The first problem I am having is in the `change_argument` function: ```julia...

Hey @oriolcg, that problem can be solved by adding the triangulation that you are looping over at that point as an argument to the `get_cell_dof_values` function (see [LOC of fork](https://github.com/gridap/Gridap.jl/commit/a5fd7153e92c421451162908e4ce1e9c2b325fe7#diff-89bb8c41412f11b17822062405e79a1b992e7a1cc02c083b95b15de6c964429bR39)...

hey @amartinhuertas, have you had a chance to take a look at this?

that does not solve the issue for me. The output with LD_DEBUG=libs is : ``` process. 11218: find library=imobiledevice.so [0]; searching 11218: search cache=/etc/ld.so.cache 11218: search path=/lib/x86_64-linux-gnu:/usr/lib/x86_64-linux-gnu:/lib:/usr/lib (system search path)...

Hey @fverdugo! The problem still remains after increasing `degree`.

Hey @oriolcg, I ended up not needing this functionality so I don't have a work-around. I am still happy however to help this get implemented.

Hey @oriolcg , have you ever had a problem similar to this when solving for flow problems in a channel ? If I solve the Stokes problem in a channel...

problem is related to the linear solver. change the ls or add terms in the bottom right block to avoid the issue (e.g. `h^2 (∇p ⋅ ∇q)` )

Is it possible to tag these isolated volumes? I mean to somehow know which sections of the domain are part of a set ISOLATED (red) or the set NOT ISOLATED...