Hi, I am developing an asynchronous method for a FEM solver, and need a custom implementation of the VecScatterBegin and VecScatterEnd routines. Since PETSc uses its own limited set of MPI functions, could you tell what would be the best way possible to extend upon it and use for example the non-blocking collectives, igatherv and so on?
I hope the question is specific enough, let me know if anything, I can provide with more information. I would very much appreciate any help, thanks in advance! Best, Tamara
