But why Kokkos Models?

kokkos.cpp

#pragma omp parallel for
for(int i=0;i<lots;++i){
  // insert psychics here
}

__global__ void cuda(){
  int thread = threadIdx.x +
    blockdim.x * blockIdx.x;
  if(thread < lots){
    // insert psychics here
  }
}

Kokkos::parallel_for("do_psychics",lots, 
KOKKOS_LAMBDA(int i){
  // insert psychics here
});