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