xref: /petsc/src/sys/objects/kokkos/kinit.kokkos.cxx (revision b9d03b0c270a89a1d7975e414e2df3c3436115f7)
1*b9d03b0cSStefano Zampini #include <petscdevice.h>
2c2b86a48SJunchao Zhang #include <petsc/private/petscimpl.h>
3c2b86a48SJunchao Zhang #include <Kokkos_Core.hpp>
4c2b86a48SJunchao Zhang 
545639126SStefano Zampini PetscBool PetscKokkosInitialized = PETSC_FALSE;
645639126SStefano Zampini 
7c2b86a48SJunchao Zhang PetscErrorCode PetscKokkosFinalize_Private(void)
8c2b86a48SJunchao Zhang {
9c2b86a48SJunchao Zhang   PetscFunctionBegin;
10c2b86a48SJunchao Zhang   Kokkos::finalize();
11c2b86a48SJunchao Zhang   PetscFunctionReturn(0);
12c2b86a48SJunchao Zhang }
13c2b86a48SJunchao Zhang 
14c2b86a48SJunchao Zhang PetscErrorCode PetscKokkosIsInitialized_Private(PetscBool *isInitialized)
15c2b86a48SJunchao Zhang {
16c2b86a48SJunchao Zhang   PetscFunctionBegin;
17c2b86a48SJunchao Zhang   *isInitialized = Kokkos::is_initialized() ? PETSC_TRUE : PETSC_FALSE;
18c2b86a48SJunchao Zhang   PetscFunctionReturn(0);
19c2b86a48SJunchao Zhang }
20375e5adfSJunchao Zhang 
2134766dafSJunchao Zhang /* Initialize Kokkos if not yet */
22375e5adfSJunchao Zhang PetscErrorCode PetscKokkosInitializeCheck(void)
23375e5adfSJunchao Zhang {
2455c8bb78SStefano Zampini #if defined(KOKKOS_ENABLE_CUDA) || defined(KOKKOS_ENABLE_HIP)
25375e5adfSJunchao Zhang   PetscErrorCode        ierr;
2655c8bb78SStefano Zampini #endif
2734766dafSJunchao Zhang   Kokkos::InitArguments args;
2834766dafSJunchao Zhang   int                   devId = -1;
29375e5adfSJunchao Zhang 
30375e5adfSJunchao Zhang   PetscFunctionBegin;
3134766dafSJunchao Zhang   if (!Kokkos::is_initialized()) {
32375e5adfSJunchao Zhang    #if defined(KOKKOS_ENABLE_CUDA)
33*b9d03b0cSStefano Zampini     cudaError_t cerr;
34*b9d03b0cSStefano Zampini 
35375e5adfSJunchao Zhang     ierr = PetscCUDAInitializeCheck();CHKERRQ(ierr);
36*b9d03b0cSStefano Zampini     cerr = cudaGetDevice(&devId);CHKERRCUDA(cerr);
3734766dafSJunchao Zhang    #elif defined(KOKKOS_ENABLE_HIP) /* Kokkos does not support CUDA and HIP at the same time */
38*b9d03b0cSStefano Zampini     hipError_t herr;
39*b9d03b0cSStefano Zampini 
40375e5adfSJunchao Zhang     ierr = PetscHIPInitializeCheck();CHKERRQ(ierr);
41*b9d03b0cSStefano Zampini     herr = hipGetDevice(&devId);CHKERRHIP(herr);
42375e5adfSJunchao Zhang    #endif
4334766dafSJunchao Zhang     args.device_id   = devId;
4434766dafSJunchao Zhang     Kokkos::initialize(args);
4559e55d94SJunchao Zhang     PetscBeganKokkos = PETSC_TRUE;
4659e55d94SJunchao Zhang   }
4745639126SStefano Zampini   PetscKokkosInitialized = PETSC_TRUE;
48375e5adfSJunchao Zhang   PetscFunctionReturn(0);
49375e5adfSJunchao Zhang }
50