Hi Randall, Thanks for providing a pointer to the DMDAGetRay function!
After looking at its implementation, I came up with a solution that creates a natural ordered slice vector on the same subset of processors as the DMDA ordered slice vector (by scattering from the DMDA order slice to a natural ordered slice by using an AO associated with a temporary 2D DMDA object that lives only on the subset of ranks where the slice vector lives). I've attached the code for the same should it be of interest to anyone who reads this. -- Sajid Ali (he/him) | PhD Candidate Applied Physics Northwestern University s-sajid-ali.github.io
ex_slice_nat.c
Description: Binary data
