Hello,

I have a LinearImplicitSystem with a mesh imported from an Exodus file.
When assembling the system matrix and rhs, I proceed according to your
examples iterating over the elements with the active_local_elements
iterator.

I want to do some calculations in between time steps using the solution
vector and iterating over the nodes. When running in parallel, depending on
the iterator, I obtain inconsistent results. Specifically, the nodes that
are shared between processors don't have numerical value assigned to them
which results in a strange region that propagates error (pic attached when
running with 2 proc).

I have tried using different iterators: local_nodes, active_local_nodes,
pid_nodes, and iterating over the list obtained
from dof_map.get_send_list(). The results are inconsistent and never match
the conditions when running in 1 processor. Is there any way I can access
the list of nodes that comprise that boundary? Should I be using some
function from the Parallel communicator?

Thanks!
Best,
Viviana.

 --
*Viviana Palacio Betancur*
PhD Candidate || Graduate Research Assistant || de Pablo group
<https://pme.uchicago.edu/de_pablo_lab/>
Office: ERC 351 || Email: vpalac...@uchicago.edu || PME @ UChicago
Phone: +1 312 889 2434

_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to