Hello, I am working on a code in which a DMSTAG object is used to solve a fluid flow problem and I need to gather this flow data on a single process to interact with an existing (serial) library at each timestep of my simulation. After looking around the solution I've tried is:
-use DMStagVecSplitToDMDA to extract vectors of each component of the flow -use DMDACreateNaturalVector and DMDAGlobalToNatural to get the components naturally ordered -use VecScatterCreateToZero to set up and then do the scatter to gather on the single process Unless I'm misunderstanding something this method results in a lot of memory allocation/freeing happening at each step of the evolution and I was wondering if there is a way to directly perform such a scatter from the DMSTAG object without splitting as I'm doing here. Any advice would be much appreciated! Best, Colton Bryant
