To compile and run CUDA code, you need a CUDA-enabled GPU, appropriate drivers, the CUDA toolkit, and the nvcc compiler. Unlike other libraries, you need special hardware and compiler to enable CUDA. Because the hardware is always evolving, older GPUs do not support all the capabilities of newer ones. In order to use CUDA with deal.II, you will need your GPU to have compute capability 3.5 or higher. Independently from the GPU itself, you also need a version of CUDA recent enough. deal.II supports CUDA 8.0 and higher. Finally to be able to configure deal.II, you will need CMake 3.9 or higher.
To configure deal.II with CUDA use the following option:
-DDEAL_II_WITH_CUDA=ONDepending on you system, this may be enough to get CUDA to work. If you are using CUDA 8 with gcc 5.4, you will need to turn off support for C++14:
-DDEAL_II_WITH_CXX14=OFFIf you are using CUDA 9 or CUDA 10, you will need to turn off support for C++17 similarly. By default, we try to detect the compute capability of your device but you can easily set your own CUDA flags:
-DDEAL_II_CUDA_FLAGS="-arch=sm_60"
-DDEAL_II_CUDA_FLAGS_DEBUG
and
-DDEAL_II_CUDA_FLAGS_RELEASE
are also available if you want
a finer control on the CUDA flags. The CUDA compiler and the
CUDA toolkit root directory can be set using
-DDEAL_II_CUDA_COMPILER
and
-DDEAL_II_CUDA_TOOLKIT_ROOT_DIR
.
Finally, the CUDA host compiler is the same as the C++ compiler
by default, but can be changed using CUDA flags as well.
Several MPI implementations are able to perform MPI operations with data located in device memory directly without the need to copy to CPU memory explicitly first. This feature is commonly known as "CUDA-aware MPI". In case deal.II is compiled with MPI support and the MPI implementation supports this feature, you can tell deal.II to use it by configuring with
-DDEAL_II_WITH_CUDA=ON -DDEAL_II_WITH_MPI=ON -DDEAL_II_MPI_WITH_CUDA_SUPPORT=ONNote, that there is no check that detects if the MPI implementation really is CUDA-aware. Activating this flag for incompatible MPI libraries will lead to segmentation faults in MPI calls.
Using CUDA in combination with architecture-specific C++ compiler flags
like -march=native
is known to be fragile and there might be
compatibility issues with other libraries, e.g. using CUDA 10.1 with
-DDEAL_II_WITH_THREADS=ON
and
-DDEAL_II_CXX_FLAGS=-march=native
results in compile time
errors like:
/usr/lib/gcc/x86_64-linux-gnu/7/include/avx512fintrin.h(11265): error: identifier "__builtin_ia32_scalefsd_round" is undefined /usr/lib/gcc/x86_64-linux-gnu/7/include/avx512fintrin.h(11274): error: identifier "__builtin_ia32_scalefss_round" is undefinedSince vectorization in VectorizedArray is disabled when compiling with CUDA support anyway, it is recommended to drop the compile flag in that case.