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