I'm trying to understand how VecScatters work with GPU-native Kokkos Vecs. 

Specifically, I'm interested in what will happen in code like in 
src/vec/vec/tests/ex22.c, 

        ierr = VecScatterCreate(x,is1,y,is2,&ctx);CHKERRQ(ierr);

(from https://gitlab.com/petsc/petsc/-/blob/master/src/vec/vec/tests/ex22.c#L44 
<https://gitlab.com/petsc/petsc/-/blob/master/src/vec/vec/tests/ex22.c#L44>)

Here, x and y can be set to type KOKKOS using -vec_type kokkos at the command 
line. But is1 and is2 are (I think), always
CPU/host data. Assuming that the scatter itself can happen on the GPU, the 
indices must make it to the device somehow - are they copied there when the 
scatter is created? Is there a way to create the scatter using indices already 
on the GPU (Maybe using SF more directly)?

Reply via email to