#include "../../interface/sycldevice.hpp" #include // for MPI sycl device awareness #include // SIGSEGV #include #include #if PetscDefined(USE_LOG) PETSC_INTERN PetscErrorCode PetscLogInitialize(void); #else #define PetscLogInitialize() 0 #endif namespace Petsc { // definition for static std::array SyclDevice::devices_array_ = {}; SyclDevice::SyclDeviceInternal** SyclDevice::devices_ = &SyclDevice::devices_array_[1]; int SyclDevice::defaultDevice_ = PETSC_SYCL_DEVICE_NONE; bool SyclDevice::initialized_ = false; static std::jmp_buf MPISyclAwareJumpBuffer; static bool MPISyclAwareJumpBufferSet; // internal "impls" class for SyclDevice. Each instance represents a single sycl device class PETSC_NODISCARD SyclDevice::SyclDeviceInternal { const int id_; // -1 for the host device; 0 and up for gpu devices bool devInitialized_; const sycl::device syclDevice_; public: // default constructor SyclDeviceInternal(int id) noexcept : id_(id),devInitialized_(false),syclDevice_(chooseSYCLDevice_(id)){} int id() const {return id_;} bool initialized() const {return devInitialized_;} PETSC_NODISCARD PetscErrorCode initialize() noexcept { PetscFunctionBegin; if (devInitialized_) PetscFunctionReturn(0); if (syclDevice_.is_gpu() && use_gpu_aware_mpi) { if (!isMPISyclAware_()) { (*PetscErrorPrintf)("PETSc is configured with sycl support, but your MPI is not aware of sycl GPU devices. For better performance, please use a sycl GPU-aware MPI.\n"); (*PetscErrorPrintf)("If you do not care, add option -use_gpu_aware_mpi 0. To not see the message again, add the option to your .petscrc, OR add it to the env var PETSC_OPTIONS.\n"); PETSCABORT(PETSC_COMM_SELF,PETSC_ERR_LIB); } } devInitialized_ = true; PetscFunctionReturn(0); } PETSC_NODISCARD PetscErrorCode view(PetscViewer viewer) const noexcept { PetscErrorCode ierr; MPI_Comm comm; PetscMPIInt rank; PetscBool iascii; PetscFunctionBegin; if (PetscUnlikely(!devInitialized_)) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_COR,"Device %d being viewed before it was initialized or configured",id_); ierr = PetscObjectTypeCompare(reinterpret_cast(viewer),PETSCVIEWERASCII,&iascii);CHKERRQ(ierr); ierr = PetscObjectGetComm(reinterpret_cast(viewer),&comm);CHKERRQ(ierr); if (iascii) { ierr = MPI_Comm_rank(comm,&rank);CHKERRMPI(ierr); ierr = PetscViewerASCIIPushSynchronized(viewer);CHKERRQ(ierr); ierr = PetscViewerASCIISynchronizedPrintf(viewer,"[%d] device: %s\n",rank,syclDevice_.get_info().c_str());CHKERRQ(ierr); // flush the assignment information ierr = PetscViewerFlush(viewer);CHKERRQ(ierr); ierr = PetscViewerASCIIPushTab(viewer);CHKERRQ(ierr); ierr = PetscViewerASCIISynchronizedPrintf(viewer,"-> Device vendor: %s\n",syclDevice_.get_info().c_str());CHKERRQ(ierr); ierr = PetscViewerASCIIPopTab(viewer);CHKERRQ(ierr); ierr = PetscViewerFlush(viewer);CHKERRQ(ierr); ierr = PetscViewerASCIIPopSynchronized(viewer);CHKERRQ(ierr); } PetscFunctionReturn(0); } private: static sycl::device chooseSYCLDevice_(int id) { if (id == PETSC_SYCL_DEVICE_HOST) { return sycl::device(sycl::host_selector()); } else { return sycl::device::get_devices(sycl::info::device_type::gpu)[id]; } } // Is the underlying MPI aware of sycl (GPU) devices? bool isMPISyclAware_() noexcept { PetscErrorCode ierr; const int bufSize = 2; const int hbuf[bufSize] = {1,0}; int *dbuf = nullptr; bool awareness = false; const auto SyclSignalHandler = [](int signal, void *ptr) -> PetscErrorCode { if ((signal == SIGSEGV) && MPISyclAwareJumpBufferSet) std::longjmp(MPISyclAwareJumpBuffer,1); return PetscSignalHandlerDefault(signal,ptr); }; PetscFunctionBegin; auto Q = sycl::queue(syclDevice_); dbuf = sycl::malloc_device(bufSize,Q); Q.memcpy(dbuf,hbuf,sizeof(int)*bufSize).wait(); ierr = PetscPushSignalHandler(SyclSignalHandler,nullptr);CHKERRABORT(PETSC_COMM_SELF,ierr); MPISyclAwareJumpBufferSet = true; if (setjmp(MPISyclAwareJumpBuffer)) { // if a segv was triggered in the MPI_Allreduce below, it is very likely due to MPI not being GPU-aware awareness = false; PetscStackPop; } else if (!MPI_Allreduce(dbuf,dbuf+1,1,MPI_INT,MPI_SUM,PETSC_COMM_SELF)) awareness = true; MPISyclAwareJumpBufferSet = false; ierr = PetscPopSignalHandler();CHKERRABORT(PETSC_COMM_SELF,ierr); sycl::free(dbuf,Q); PetscFunctionReturn(awareness); } }; PetscErrorCode SyclDevice::initialize(MPI_Comm comm, PetscInt *defaultDeviceId, PetscDeviceInitType *defaultInitType) noexcept { PetscInt initType = *defaultInitType,id = *defaultDeviceId; PetscBool view = PETSC_FALSE,flg; PetscInt ngpus; PetscErrorCode ierr; PetscFunctionBegin; if (initialized_) PetscFunctionReturn(0); initialized_ = true; ierr = PetscRegisterFinalize(finalize_);CHKERRQ(ierr); ierr = PetscOptionsBegin(comm,nullptr,"PetscDevice SYCL Options","Sys");CHKERRQ(ierr); ierr = PetscOptionsEList("-device_enable_sycl","How (or whether) to initialize a device","SyclDevice::initialize()",PetscDeviceInitTypes,3,PetscDeviceInitTypes[initType],&initType,nullptr);CHKERRQ(ierr); ierr = PetscOptionsRangeInt("-device_select_sycl","Which sycl device to use? Pass -2 for host, PETSC_DECIDE (-1) to let PETSc decide, 0 and up for GPUs","PetscDeviceCreate",id,&id,nullptr,-2,std::numeric_limits::max());CHKERRQ(ierr); ierr = PetscOptionsBool("-device_view_sycl","Display device information and assignments (forces eager initialization)",nullptr,view,&view,&flg);CHKERRQ(ierr); ierr = PetscOptionsEnd();CHKERRQ(ierr); // post-process the options and lay the groundwork for initialization if needs be std::vector gpu_devices = sycl::device::get_devices(sycl::info::device_type::gpu); ngpus = static_cast(gpu_devices.size()); if (PetscUnlikely(ngpus == 0 && id >= 0)) SETERRQ1(comm,PETSC_ERR_USER_INPUT,"You specified a sycl gpu device with -device_select_sycl %d but there is no GPU", (int)id); if (PetscUnlikely(ngpus > 0 && id >= ngpus)) SETERRQ2(comm,PETSC_ERR_USER_INPUT,"You specified a sycl gpu device with -device_select_sycl %d but there are only %d GPU", (int)id, (int)ngpus); if (initType == PETSC_DEVICE_INIT_NONE) id = PETSC_SYCL_DEVICE_NONE; /* user wants to disable all sycl devices */ else { ierr = PetscDeviceCheckDeviceCount_Internal(ngpus);CHKERRQ(ierr); if (id == PETSC_DECIDE) { /* petsc will choose a GPU device if any, otherwise a CPU device */ if (ngpus) { PetscMPIInt rank; ierr = MPI_Comm_rank(comm,&rank);CHKERRMPI(ierr); id = rank % ngpus; } else id = PETSC_SYCL_DEVICE_HOST; } view = static_cast(view && flg); if (view) initType = PETSC_DEVICE_INIT_EAGER; } if (id == -2) id = PETSC_SYCL_DEVICE_HOST; // user passed in '-device_select_sycl -2'. We transform it into canonical form defaultDevice_ = static_cast(id); if (PetscUnlikely(initType == PETSC_DEVICE_INIT_EAGER && id == PETSC_SYCL_DEVICE_NONE)) SETERRQ(comm,PETSC_ERR_USER_INPUT,"Cannot eagerly initialize sycl devices as you disabled them by -device_enable_sycl none"); if (initType == PETSC_DEVICE_INIT_EAGER) { devices_[defaultDevice_] = new SyclDeviceInternal(defaultDevice_); ierr = devices_[defaultDevice_]->initialize();CHKERRQ(ierr); if (view) { PetscViewer viewer; ierr = PetscLogInitialize();CHKERRQ(ierr); ierr = PetscViewerASCIIGetStdout(comm,&viewer);CHKERRQ(ierr); ierr = devices_[defaultDevice_]->view(viewer);CHKERRQ(ierr); } } // record the results of the initialization *defaultInitType = static_cast(initType); *defaultDeviceId = id; PetscFunctionReturn(0); } PetscErrorCode SyclDevice::finalize_() noexcept { PetscFunctionBegin; if (!initialized_) PetscFunctionReturn(0); for (auto&& devPtr : devices_array_) delete devPtr; defaultDevice_ = PETSC_SYCL_DEVICE_NONE; // disabled by default initialized_ = false; PetscFunctionReturn(0); } PetscErrorCode SyclDevice::getDevice(PetscDevice device, PetscInt id) const noexcept { PetscErrorCode ierr; PetscFunctionBegin; if (PetscUnlikely(defaultDevice_ == PETSC_SYCL_DEVICE_NONE)) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_WRONGSTATE,"Trying to retrieve a SYCL PetscDevice when it has been disabled"); if (id == PETSC_DECIDE) id = defaultDevice_; if ((id < PETSC_SYCL_DEVICE_HOST) || (id-PETSC_SYCL_DEVICE_HOST >= PETSC_DEVICE_MAX_DEVICES)) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Only supports %zu number of devices but trying to get device with id %" PetscInt_FMT,devices_array_.size(),id); if (devices_[id]) { if (id != devices_[id]->id()) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_PLIB,"Entry %" PetscInt_FMT " contains device with mismatching id %" PetscInt_FMT,id,devices_[id]->id()); } else devices_[id] = new SyclDeviceInternal(id); ierr = devices_[id]->initialize();CHKERRQ(ierr); device->deviceId = devices_[id]->id(); // technically id = devices_[id]->id_ here device->ops->createcontext = create_; device->ops->configure = this->configureDevice; device->ops->view = this->viewDevice; PetscFunctionReturn(0); } PetscErrorCode SyclDevice::configureDevice(PetscDevice device) noexcept { PetscFunctionBegin; // Nothing for now PetscFunctionReturn(0); } PetscErrorCode SyclDevice::viewDevice(PetscDevice device, PetscViewer viewer) noexcept { PetscErrorCode ierr; PetscFunctionBegin; ierr = devices_[device->deviceId]->view(viewer);CHKERRQ(ierr); PetscFunctionReturn(0); } } // namespace Petsc