Dear Floriane, On 25/04/16 10:15, Floriane Gidel [RPG] wrote:
Dear all,
Do you have any tips concerning the two questions below ? Thank you and best regards,
I don't have an idea as to why in the CG2 case you don't get identical values everywhere. Although I note that you are solving an equation so you can't expect to get a solution that is closer than the solver tolerance (by default 1e-8) in the general case. As to transferring values from the top surface to the bottom. You can do that like so: f = Function(V) # The nodes topologically at the bottom of the domain: bottom = V.bottom_nodes() # Sim at the top top = V.top_nodes() # Copy from top nodes into bottom nodes f.dat.data_with_halos[bottom] = f.dat.data_with_halos[top] This will only work on extruded meshes (since there there is an easy identification of the "top" and "bottom" of the domain). Cheers, Lawrence