Is there a way to retrieve global array when running in parallel
4 weeks ago by
Community: FEniCS Project
25 days ago by
comm = mesh.mpi_comm() arr_local = ... # some array that is different on each rank arrs = comm.gather(arr_local) if comm.rank == 0: # You have a list of arrays from each rank # arrs[rank] should be equal to arr_local on # the given rank else: assert arrs is None
You will probably want to send both vertex coordinates and vertex values. You can even send a tuple of
(coords, values)if you want (i think ...)
Please login to add an answer/comment or follow this question.