CUDA C Programming Guide
CUDA C Programming Guide
Design Guide
CHANGES FROM VERSION 7.5
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 2
TABLE OF CONTENTS
Chapter 1. Introduction 1
1. From Graphics Processing to General Purpose Parallel Computing 1
2. CUDA®: A General-Purpose Parallel Computing Platform and Programming Model 3
3. A Scalable Programming Model 4
4. Document Structure 6
Chapter 2. Programming Model 8
1. Kernels 8
2. Thread Hierarchy 9
3. Memory Hierarchy 11
4. Heterogeneous Programming 13
5. Compute Capability 15
Chapter 3. Programming Interface 16
1. Compilation with NVCC 16
1. Compilation Workflow 17
1. Offline Compilation 17
2. Just-in-Time Compilation 17
2. Binary Compatibility 17
3. PTX Compatibility 18
4. Application Compatibility 18
5. C/C++ Compatibility 19
6. 64-Bit Compatibility 19
2. CUDA C Runtime 19
1. Initialization 20
2. Device Memory 20
3. Shared Memory 23
4. Page-Locked Host Memory 28
1. Portable Memory 29
2. Write-Combining Memory 29
3. Mapped Memory 29
5. Asynchronous Concurrent Execution 30
1. Concurrent Execution between Host and Device 31
2. Concurrent Kernel Execution 31
3. Overlap of Data Transfer and Kernel Execution 31
4. Concurrent Data Transfers 32
5. Streams 32
6. Events 36
7. Synchronous Calls 36
6. Multi-Device System 37
1. Device Enumeration 37
2. Device Selection 37
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 3
3. Stream and Event Behavior 37
4. Peer-to-Peer Memory Access 38
5. Peer-to-Peer Memory Copy 38
7. Unified Virtual Address Space 39
8. Interprocess Communication 40
9. Error Checking 40
10. Call Stack 41
11. Texture and Surface Memory 41
1. Texture Memory 41
2. Surface Memory 51
3. CUDA Arrays 55
4. Read/Write Coherency 55
12. Graphics Interoperability 55
1. OpenGL Interoperability 56
2. Direct3D Interoperability 58
3. SLI Interoperability 64
3. Versioning and Compatibility 65
4. Compute Modes 66
5. Mode Switches 67
6. Tesla Compute Cluster Mode for Windows 67
Chapter 4. Hardware Implementation 68
1. SIMT Architecture 68
2. Hardware Multithreading 70
Chapter 5. Performance Guidelines 71
1. Overall Performance Optimization Strategies 71
2. Maximize Utilization 71
1. Application Level 71
2. Device Level 72
3. Multiprocessor Level 72
1. Occupancy Calculator 74
3. Maximize Memory Throughput 76
1. Data Transfer between Host and Device 77
2. Device Memory Accesses 78
4. Maximize Instruction Throughput 82
1. Arithmetic Instructions 82
2. Control Flow Instructions 86
3. Synchronization Instruction 87
Appendix A. CUDA-Enabled GPUs 88
Appendix B. C Language Extensions 89
5. Function Type Qualifiers 89
1. device 89
2. global 89
3. host 89
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 4
B.1.4. noinline and forceinline 90
2. Variable Type Qualifiers 90
1. device 90
2. constant 91
3. shared 91
4. managed 92
5. restrict 92
3. Built-in Vector Types 93
1. char, short, int, long, longlong, float, double 93
2. dim3 94
4. Built-in Variables 95
1. gridDim 95
2. blockIdx 95
3. blockDim 95
4. threadIdx 95
5. warpSize 95
5. Memory Fence Functions 95
6. Synchronization Functions 98
7. Mathematical Functions 99
8. Texture Functions 99
1. Texture Object API 100
1. tex1Dfetch() 100
B.8.1.2. tex1D() 100
3. tex1DLod() 100
4. tex1DGrad() 100
B.8.1.5. tex2D() 100
6. tex2DLod() 100
7. tex2DGrad() 101
B.8.1.8. tex3D() 101
9. tex3DLod() 101
10. tex3DGrad() 101
11. tex1DLayered() 101
12. tex1DLayeredLod() 101
13. tex1DLayeredGrad() 102
14. tex2DLayered() 102
15. tex2DLayeredLod() 102
16. tex2DLayeredGrad() 102
17. texCubemap() 102
18. texCubemapLod() 102
19. texCubemapLayered() 103
20. texCubemapLayeredLod() 103
21. tex2Dgather() 103
B.8.2. Texture Reference API 104
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 5
B.8.2.1. tex1Dfetch() 104
B.8.2.2. tex1D() 104
3. tex1DLod() 105
4. tex1DGrad() 105
B.8.2.5. tex2D() 105
6. tex2DLod() 105
7. tex2DGrad() 105
B.8.2.8. tex3D() 106
9. tex3DLod() 106
10. tex3DGrad() 106
11. tex1DLayered() 106
12. tex1DLayeredLod() 107
13. tex1DLayeredGrad() 107
14. tex2DLayered() 107
15. tex2DLayeredLod() 107
16. tex2DLayeredGrad() 108
17. texCubemap() 108
18. texCubemapLod() 108
19. texCubemapLayered() 108
20. texCubemapLayeredLod() 108
21. tex2Dgather() 109
9. Surface Functions 109
1. Surface Object API 109
1. surf1Dread() 109
2. surf1Dwrite 109
3. surf2Dread() 110
4. surf2Dwrite() 110
5. surf3Dread() 110
6. surf3Dwrite() 110
7. surf1DLayeredread() 111
8. surf1DLayeredwrite() 111
9. surf2DLayeredread() 111
10. surf2DLayeredwrite() 111
11. surfCubemapread() 112
12. surfCubemapwrite() 112
13. surfCubemapLayeredread() 112
14. surfCubemapLayeredwrite() 112
2. Surface Reference API 113
1. surf1Dread() 113
2. surf1Dwrite 113
3. surf2Dread() 113
4. surf2Dwrite() 113
5. surf3Dread() 114
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 6
6. surf3Dwrite() 114
7. surf1DLayeredread() 114
8. surf1DLayeredwrite() 114
9. surf2DLayeredread() 115
10. surf2DLayeredwrite() 115
11. surfCubemapread() 115
12. surfCubemapwrite() 115
13. surfCubemapLayeredread() 116
14. surfCubemapLayeredwrite() 116
10. Read-Only Data Cache Load Function 116
11. Time Function 116
12. Atomic Functions 117
1. Arithmetic Functions 118
1. atomicAdd() 118
2. atomicSub() 118
3. atomicExch() 119
4. atomicMin() 119
5. atomicMax() 119
6. atomicInc() 119
7. atomicDec() 120
8. atomicCAS() 120
2. Bitwise Functions 120
1. atomicAnd() 120
2. atomicOr() 120
3. atomicXor() 121
13. Warp Vote Functions 121
14. Warp Shuffle Functions 122
1. Synopsis 122
2. Description 122
3. Return Value 123
4. Notes 123
5. Examples 124
1. Broadcast of a single value across a warp 124
2. Inclusive plus-scan across sub-partitions of 8 threads 124
3. Reduction across a warp 125
15. Profiler Counter Function 125
16. Assertion 125
17. Formatted Output 126
1. Format Specifiers 127
2. Limitations 127
3. Associated Host-Side API 128
4. Examples 129
18. Dynamic Global Memory Allocation and Operations 130
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 7
1. Heap Memory Allocation 130
2. Interoperability with Host Memory API 131
3. Examples 131
1. Per Thread Allocation 131
2. Per Thread Block Allocation 132
3. Allocation Persisting Between Kernel Launches 133
19. Execution Configuration 134
20. Launch Bounds 134
21. #pragma unroll 137
22. SIMD Video Instructions 137
Appendix C. CUDA Dynamic Parallelism 139
1. Introduction 139
1. Overview 139
2. Glossary 139
2. Execution Environment and Memory Model 140
1. Execution Environment 140
1. Parent and Child Grids 140
2. Scope of CUDA Primitives 141
3. Synchronization 141
4. Streams and Events 141
5. Ordering and Concurrency 142
6. Device Management 142
2. Memory Model 142
1. Coherence and Consistency 143
3. Programming Interface 145
1. CUDA C/C++ Reference 145
1. Device-Side Kernel Launch 145
2. Streams 146
3. Events 147
4. Synchronization 147
5. Device Management 147
6. Memory Declarations 148
7. API Errors and Launch Failures 149
8. API Reference 150
2. Device-side Launch from PTX 151
1. Kernel Launch APIs 151
2. Parameter Buffer Layout 153
3. Toolkit Support for Dynamic Parallelism 153
1. Including Device Runtime API in CUDA Code 153
2. Compiling and Linking 154
4. Programming Guidelines 154
1. Basics 154
2. Performance 155
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 8
1. Synchronization 155
2. Dynamic-parallelism-enabled Kernel Overhead 155
3. Implementation Restrictions and Limitations 156
1. Runtime 156
Appendix D. Mathematical Functions 159
1. Standard Functions 159
2. Intrinsic Functions 167
Appendix E. C/C++ Language Support 170
1. C++11 Language Features 170
2. Restrictions 173
1. Host Compiler Extensions 173
2. Preprocessor Symbols 173
1. CUDA_ARCH 173
3. Qualifiers 175
1. Device Memory Qualifiers 175
2. managed Qualifier 176
3. Volatile Qualifier 177
4. Pointers 178
5. Operators 178
1. Assignment Operator 178
2. Address Operator 178
6. Run Time Type Information (RTTI) 178
7. Exception Handling 178
8. Standard Library 178
9. Functions 179
1. External Linkage 179
2. Compiler generated functions 179
3. Function Parameters 179
4. Static Variables within Function 180
5. Function Pointers 180
6. Function Recursion 181
10. Classes 181
1. Data Members 181
2. Function Members 181
3. Virtual Functions 181
4. Virtual Base Classes 181
5. Anonymous Unions 181
6. Windows-Specific 181
11. Templates 182
12. Trigraphs and Digraphs 182
13. Const-qualified variables 182
14. C++11 Features 183
1. Lambda Expressions 183
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 9
2. std::initializer_list 184
3. Rvalue references 185
4. Constexpr functions and function templates 185
5. Constexpr variables 185
6. Inline namespaces 186
7. thread_local 187
8. global functions and function templates 187
9. device / constant / shared variables 189
3. Polymorphic Function Wrappers 189
4. Experimental Feature: Extended Lambdas 191
1. Extended Lambda Type Traits 193
2. Extended Lambda Restrictions 194
3. Notes on host device lambdas 200
4. *this Capture By Value 200
5. Additional Notes 203
5. Code Samples 205
1. Data Aggregation Class 205
2. Derived Class 205
3. Class Template 206
4. Function Template 206
5. Functor Class 207
Appendix F. Texture Fetching 208
1. Nearest-Point Sampling 208
2. Linear Filtering 209
3. Table Lookup 210
Appendix G. Compute Capabilities 212
1. Features and Technical Specifications 212
2. Floating-Point Standard 216
3. Compute Capability 2.x 217
1. Architecture 217
2. Global Memory 218
3. Shared Memory 219
4. Constant Memory 220
4. Compute Capability 3.x 220
1. Architecture 220
2. Global Memory 222
3. Shared Memory 223
5. Compute Capability 5.x 224
1. Architecture 224
2. Global Memory 225
3. Shared Memory 225
6. Compute Capability 6.x 229
1. Architecture 229
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 10
2. Global Memory 229
3. Shared Memory 229
Appendix H. Driver API 230
1. Context 233
2. Module 234
3. Kernel Execution 235
4. Interoperability between Runtime and Driver APIs 237
Appendix I. CUDA Environment Variables 238
Appendix J. Unified Memory Programming 241
5. Unified Memory Introduction 241
1. Simplifying GPU Programming 242
2. Data Migration and Coherency 243
3. GPU Memory Oversubscription 244
4. Multi-GPU Support 244
5. System Requirements 245
6. Programming Model 245
1. Managed Memory Opt In 245
1. Explicit Allocation Using cudaMallocManaged() 245
2. Global-Scope Managed Variables Using managed 246
2. Coherency and Concurrency 247
1. GPU Exclusive Access To Managed Memory 247
2. Explicit Synchronization and Logical GPU Activity 248
3. Managing Data Visibility and Concurrent CPU + GPU Access with Streams 249
4. Stream Association Examples 250
5. Stream Attach With Multithreaded Host Programs 251
6. Advanced Topic: Modular Programs and Data Access Constraints 252
7. Memcpy()/Memset() Behavior With Managed Memory 253
3. Language Integration 253
1. Host Program Errors with managed Variables 254
4. Querying Unified Memory Support 255
1. Device Properties 255
2. Pointer Attributes 255
5. Advanced Topics 255
1. Managed Memory with Multi-GPU Programs on pre-6.x Architectures 255
2. Using fork() with Managed Memory 256
7. Performance Tuning 256
1. Data Prefetching 257
2. Data Usage Hints 258
3. Querying Usage Attributes 259
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 11
LIST OF FIGURES
Figure 1 Floating-Point Operations per Second for the CPU and GPU 1
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 12
LIST OF TABLES
Table 6 Single-Precision Mathematical Standard Library Functions with Maximum ULP Error 159
Table 7 Double-Precision Mathematical Standard Library Functions with Maximum ULP Error... 163
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 13
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 14
Chapter 1.
INTRODUCTION
Figure 1 Floating-Point Operations per Second for the CPU and GPU
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 1
Introduction
ALU ALU
Cache
DRAM DRAM
CPU GPU
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 2
Introduction
there is a lower requirement for sophisticated flow control, and because it is executed on
many data elements and has high arithmetic intensity, the memory access latency can be
hidden with calculations instead of big data caches.
Data-parallel processing maps data elements to parallel processing threads. Many
applications that process large data sets can use a data-parallel programming model
to speed up the computations. In 3D rendering, large sets of pixels and vertices are
mapped to parallel threads. Similarly, image and media processing applications such as
post-processing of rendered images, video encoding and decoding, image scaling, stereo
vision, and pattern recognition can map image blocks and pixels to parallel processing
threads. In fact, many algorithms outside the field of image rendering and processing
are accelerated by data-parallel processing, from general signal processing or physics
simulation to computational finance or computational biology.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 3
Introduction
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 4
Introduction
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 5
Introduction
SM 0 SM 1 SM 0 SM 1 SM 2 SM 3
Block 4 Block 5
Block 6 Block 7
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 6
Introduction
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 7
Chapter 2.
PROGRAMMING MODEL
This chapter introduces the main concepts behind the CUDA programming model by
outlining how they are exposed in C. An extensive description of CUDA C is given in
Programming Interface.
Full code for the vector addition example used in this chapter and the next can be found
in the vectorAdd CUDA sample.
2.1. Kernels
CUDA C extends C by allowing the programmer to define C functions, called kernels,
that, when called, are executed N times in parallel by N different CUDA threads, as
opposed to only once like regular C functions.
A kernel is defined using the global declaration specifier and the number of
CUDA threads that execute that kernel for a given kernel call is specified using a new
<<<...>>> execution configuration syntax (see C Language Extensions). Each thread
that executes the kernel is given a unique thread ID that is accessible within the kernel
through the built-in threadIdx variable.
As an illustration, the following sample code adds two vectors A and B of size N and
stores the result into vector C:
// Kernel definition
global void VecAdd(float* A, float* B, float* C)
{
int i = threadIdx.x;
C[i] = A[i] + B[i];
}
int main()
{
...
// Kernel invocation with N threads
VecAdd<<<1, N>>>(A, B, C);
...
}
Here, each of the N threads that execute VecAdd() performs one pair-wiseaddition.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 8
Programming Model
int main()
{
...
// Kernel invocation with one block of N * N * 1 threads
int numBlocks = 1;
dim3 threadsPerBlock(N, N);
MatAdd<<<numBlocks, threadsPerBlock>>>(A, B, C);
...
}
There is a limit to the number of threads per block, since all threads of a block are
expected to reside on the same processor core and must share the limited memory
resources of that core. On current GPUs, a thread block may contain up to 1024 threads.
However, a kernel can be executed by multiple equally-shaped thread blocks, so that the
total number of threads is equal to the number of threads per block times the number of
blocks.
Blocks are organized into a one-dimensional, two-dimensional, or three-dimensional
grid of thread blocks as illustrated by Figure 6. The number of thread blocks in a grid is
usually dictated by the size of the data being processed or the number of processors in
the system, which it can greatly exceed.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 9
Programming Model
Grid
Block ( 1, 1)
Thr ead (0, 0) Thr ead (1, 0) Thr ead (2, 0) Thr ead (3, 0)
Thr ead (0, 1) Thr ead (1, 1) Thr ead (2, 1) Thr ead (3, 1)
Thr ead (0, 2) Thr ead (1, 2) Thr ead (2, 2) Thr ead (3, 2)
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 10
Programming Model
Extending the previous MatAdd() example to handle multiple blocks, the codebecomes
as follows.
// Kernel definition
global void MatAdd(float A[N][N], float B[N][N],
float C[N][N])
{
int i = blockIdx.x * blockDim.x + threadIdx.x;
int j = blockIdx.y * blockDim.y + threadIdx.y;
if (i < N && j < N)
C[i][j] = A[i][j] + B[i][j];
}
int main()
{
...
// Kernel invocation
dim3 threadsPerBlock(16, 16);
dim3 numBlocks(N / threadsPerBlock.x, N / threadsPerBlock.y);
MatAdd<<<numBlocks, threadsPerBlock>>>(A, B, C);
...
}
A thread block size of 16x16 (256 threads), although arbitrary in this case, is a common
choice. The grid is created with enough blocks to have one thread per matrix element
as before. For simplicity, this example assumes that the number of threads per grid in
each dimension is evenly divisible by the number of threads per block in that dimension,
although that need not be the case.
Thread blocks are required to execute independently: It must be possible to execute
them in any order, in parallel or in series. This independence requirement allows thread
blocks to be scheduled in any order across any number of cores as illustrated by Figure
5, enabling programmers to write code that scales with the number of cores.
Threads within a block can cooperate by sharing data through some shared memory and
by synchronizing their execution to coordinate memory accesses. More precisely, one
can specify synchronization points in the kernel by calling the syncthreads()
intrinsic function; syncthreads() acts as a barrier at which all threads in the block
must wait before any is allowed to proceed. Shared Memory gives an example of using
shared memory.
For efficient cooperation, the shared memory is expected to be a low-latency memory
near each processor core (much like an L1 cache) and syncthreads() is expected to
be lightweight.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 11
Programming Model
memory also offers different addressing modes, as well as data filtering, for some
specific data formats (see Texture and Surface Memory).
The global, constant, and texture memory spaces are persistent across kernel launches
by the same application.
Thr ead
Per- t hread local
mem ory
Thread Block
Per- block shared
mem ory
Grid 0
Grid 1
Global m em ory
Block ( 0, 0) Block ( 1, 0)
Block ( 0, 1) Block ( 1, 1)
Block ( 0, 2) Block ( 1, 2)
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 12
Programming Model
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 13
Programming Model
C Program
Sequential
Execution
Host
Serial code
Host
Serial code
Device
Parallel kernel Grid 1
Kernel1 < < < > > > ()
Block ( 0, 0) Block ( 1, 0)
Block ( 0, 1) Block ( 1, 1)
Block ( 0, 2) Block ( 1, 2)
Serial code executes on the host while parallel code executes on the device.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 14
Programming Model
The compute capability version of a particular GPU should not be confused with the
CUDA version (e.g., CUDA 5.5, CUDA 6, CUDA 6.5), which is the version of the CUDA
software platform. The CUDA platform is used by application developers to create
applications that run on many generations of GPU architectures, including future GPU
architectures yet to be invented. While new versions of the CUDA platform often add
native support for a new GPU architecture by supporting the compute capability
version of that architecture, new versions of the CUDA platform typically also include
software features that are independent of hardware generation.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 15
Chapter 3.
PROGRAMMING INTERFACE
CUDA C provides a simple path for users familiar with the C programming language to
easily write programs for execution by the device.
It consists of a minimal set of extensions to the C language and a runtime library.
The core language extensions have been introduced in Programming Model. They allow
programmers to define a kernel as a C function and use some new syntax to specify the
grid and block dimension each time the function is called. A complete description of all
extensions can be found in C Language Extensions. Any source file that contains some of
these extensions must be compiled with nvcc as outlined in Compilation with NVCC.
The runtime is introduced in Compilation Workflow. It provides C functions that
execute on the host to allocate and deallocate device memory, transfer data between host
memory and device memory, manage systems with multiple devices, etc. A complete
description of the runtime can be found in the CUDA reference manual.
The runtime is built on top of a lower-level C API, the CUDA driver API, which is
also accessible by the application. The driver API provides an additional level of
control by exposing lower-level concepts such as CUDA contexts - the analogue of host
processes for the device - and CUDA modules - the analogue of dynamically loaded
libraries for the device. Most applications do not use the driver API as they do not
need this additional level of control and when using the runtime, context and module
management are implicit, resulting in more concise code. The driver API is introduced
in Driver API and fully described in the reference manual.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 16
Programming Interface
1. Compilation Workflow
1. Offline Compilation
Source files compiled with nvcc can include a mix of host code (i.e., code that executes
on the host) and device code (i.e., code that executes on the device). nvcc's basic
workflow consists in separating device code from host code and then:
‣ compiling the device code into an assembly form (PTX code) and/or binary form
(cubin object),
‣ and modifying the host code by replacing the <<<...>>> syntax introduced in
Kernels (and described in more details in Execution Configuration) by the necessary
CUDA C runtime function calls to load and launch each compiled kernel from the
PTX code and/or cubin object.
The modified host code is output either as C code that is left to be compiled using
another tool or as object code directly by letting nvcc invoke the host compiler during
the last compilation stage.
Applications can then:
‣ Either link to the compiled host code (this is the most common case),
‣ Or ignore the modified host code (if any) and use the CUDA driver API (see Driver
API) to load and execute the PTX code or cubin object.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 17
Programming Interface
compatibility is guaranteed from one minor revision to the next one, but not from one
minor revision to the previous one or across major revisions. In other words, a cubin
object generated for compute capability X.y will only execute on devices of compute
capability X.z where z≥y.
3. PTX Compatibility
Some PTX instructions are only supported on devices of higher compute capabilities.
For example, warp shuffle instructions are only supported on devices of compute
capability 3.0 and above. The -arch compiler option specifies the compute capability
that is assumed when compiling C to PTX code. So, code that contains warp shuffle, for
example, must be compiled with -arch=sm_30 (or higher).
PTX code produced for some specific compute capability can always be compiled to
binary code of greater or equal compute capability.
4. Application Compatibility
To execute code on devices of specific compute capability, an application must load
binary or PTX code that is compatible with this compute capability as described in
Binary Compatibility and PTX Compatibility. In particular, to be able to execute code on
future architectures with higher compute capability (for which no binary code can be
generated yet), an application must load PTX code that will be just-in-time compiled for
these devices (see Just-in-Time Compilation).
Which PTX and binary code gets embedded in a CUDA C application is controlled by
the -arch and -code compiler options or the -gencode compiler option as detailed in
the nvcc user manual. For example,
nvcc x.cu
-gencode arch=compute_20,code=sm_20
-gencode arch=compute_30,code=sm_30
-gencode arch=compute_35,code=\'compute_35,sm_35\'
embeds binary code compatible with compute capability 2.0 and 3.0 (first and second
-gencode options) and PTX and binary code compatible with compute capability3.5
(third -gencode option).
Host code is generated to automatically select at runtime the most appropriate code to
load and execute, which, in the above example, will be:
‣ 2.0 binary code for devices with compute capability 2.0 and 2.1,
‣ 3.0 binary code for devices with compute capability 3.0,
‣ 3.5 binary code for devices with compute capability 3.5 and 3.7,
‣ PTX code which is compiled to binary code at runtime for devices with compute
capability 5.0 and higher.
x.cu can have an optimized code path that uses warp shuffle operations, for example,
which are only supported in devices of compute capability 3.0 and higher. The
CUDA_ARCH macro can be used to differentiate various code paths based on
compute capability. It is only defined for device code. When compiling with -
arch=compute_35 for example, CUDA_ARCH is equal to 350.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 18
Programming Interface
Applications using the driver API must compile code to separate files and explicitly load
and execute the most appropriate file at runtime.
The nvcc user manual lists various shorthand for the -arch, -code, and -gencode
compiler options. For example, -arch=sm_35 is a shorthand for -arch=compute_35 -
code=compute_35,sm_35 (which is the same as -gencode arch=compute_35,code=
\'compute_35,sm_35\').
5. C/C++ Compatibility
The front end of the compiler processes CUDA source files according to C++ syntax
rules. Full C++ is supported for the host code. However, only a subset of C++ is fully
supported for the device code as described in C/C++ Language Support.
6. 64-Bit Compatibility
The 64-bit version of nvcc compiles device code in 64-bit mode (i.e., pointers are 64-bit).
Device code compiled in 64-bit mode is only supported with host code compiled in 64-
bit mode.
Similarly, the 32-bit version of nvcc compiles device code in 32-bit mode and device
code compiled in 32-bit mode is only supported with host code compiled in 32-bit mode.
The 32-bit version of nvcc can compile device code in 64-bit mode also using the -m64
compiler option.
The 64-bit version of nvcc can compile device code in 32-bit mode also using the -m32
compiler option.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 19
Programming Interface
Multi-Device System shows how the programming model extends to a system with
multiple devices attached to the same host.
Error Checking describes how to properly check the errors generated by the runtime.
Call Stack mentions the runtime functions used to manage the CUDA C call stack.
Texture and Surface Memory presents the texture and surface memory spaces that
provide another way to access device memory; they also expose a subset of the GPU
texturing hardware.
Graphics Interoperability introduces the various functions the runtime provides to
interoperate with the two main graphics APIs, OpenGL and Direct3D.
1. Initialization
There is no explicit initialization function for the runtime; it initializes the first time a
runtime function is called (more specifically any function other than functions from the
device and version management sections of the reference manual). One needs to keep
this in mind when timing runtime function calls and when interpreting the error code
from the first call into the runtime.
During initialization, the runtime creates a CUDA context for each device in the system
(see Context for more details on CUDA contexts). This context is the primary context for
this device and it is shared among all the host threads of the application. As part of this
context creation, the device code is just-in-time compiled if necessary (see Just-in-Time
Compilation) and loaded into device memory. This all happens under the hood and the
runtime does not expose the primary context to the application.
When a host thread calls cudaDeviceReset(), this destroys the primary context of the
device the host thread currently operates on (i.e., the current device as defined in Device
Selection). The next runtime function call made by any host thread that has this device
as current will create a new primary context for this device.
2. Device Memory
As mentioned in Heterogeneous Programming, the CUDA programming model
assumes a system composed of a host and a device, each with their own separate
memory. Kernels operate out of device memory, so the runtime provides functions to
allocate, deallocate, and copy device memory, as well as transfer data between host
memory and device memory.
Device memory can be allocated either as linear memory or as CUDA arrays.
CUDA arrays are opaque memory layouts optimized for texture fetching. They are
described in Texture and Surface Memory.
Linear memory exists on the device in a 40-bit address space, so separately allocated
entities can reference one another via pointers, for example, in a binary tree.
Linear memory is typically allocated using cudaMalloc() and freed using cudaFree()
and data transfer between host memory and device memory are typically done using
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 20
Programming Interface
cudaMemcpy(). In the vector addition code sample of Kernels, the vectors need to be
copied from host memory to device memory:
// Device code
global void VecAdd(float* A, float* B, float* C, int N)
{
int i = blockDim.x * blockIdx.x + threadIdx.x;
if (i < N)
C[i] = A[i] + B[i];
}
// Host code
int main()
{
int N = ...;
size_t size = N * sizeof(float);
// Invoke kernel
int threadsPerBlock = 256;
int blocksPerGrid =
(N + threadsPerBlock - 1) / threadsPerBlock;
VecAdd<<<blocksPerGrid, threadsPerBlock>>>(d_A, d_B, d_C, N);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 21
Programming Interface
following code sample allocates a width x height 2D array of floating-point values and
shows how to loop over the array elements in device code:
// Host code
int width = 64, height = 64;
float* devPtr;
size_t pitch;
cudaMallocPitch(&devPtr, &pitch,
width * sizeof(float), height);
MyKernel<<<100, 512>>>(devPtr, pitch, width, height);
// Device code
global void MyKernel(float* devPtr,
{ size_t pitch, int width, int height)
The following code sample allocates a width x height x depth 3D array of floating-
point values and shows how to loop over the array elements in device code:
// Host code
int width = 64, height = 64, depth = 64;
cudaExtent extent = make_cudaExtent(width * sizeof(float),
height, depth);
cudaPitchedPtr devPitchedPtr;
cudaMalloc3D(&devPitchedPtr, extent);
MyKernel<<<100, 512>>>(devPitchedPtr, width, height, depth);
// Device code
global void MyKernel(cudaPitchedPtr devPitchedPtr,
{ int width, int height, int depth)
The reference manual lists all the various functions used to copy memory between
linear memory allocated with cudaMalloc(), linear memory allocated with
cudaMallocPitch() or cudaMalloc3D(), CUDA arrays, and memory allocated for
variables declared in global or constant memory space.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 22
Programming Interface
The following code sample illustrates various ways of accessing global variables via the
runtime API:
constant float constData[256];
float data[256];
cudaMemcpyToSymbol(constData, data, sizeof(data));
cudaMemcpyFromSymbol(data, constData, sizeof(data));
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 23
Programming Interface
// Invoke kernel
dim3 dimBlock(BLOCK_SIZE, BLOCK_SIZE);
dim3 dimGrid(B.width / dimBlock.x, A.height / dimBlock.y);
MatMulKernel<<<dimGrid, dimBlock>>>(d_A, d_B, d_C);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 24
Programming Interface
B. w idt h-1
0 col
B. height
0
A C
A.height
row
A.height -1
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 25
Programming Interface
By blocking the computation this way, we take advantage of fast shared memory and
save a lot of global memory bandwidth since A is only read (B.width / block_size) times
from global memory and B is read (A.height / block_size) times.
The Matrix type from the previous code sample is augmented with a stride field, so that
sub-matrices can be efficiently represented with the same type. device functions are
used to get and set elements and build any sub-matrix from a matrix.
// Matrices are stored in row-major order:
// M(row, col) = *(M.elements + row * M.stride + col)
typedef struct {
int width;
int height;
int stride;
float* elements;
} Matrix;
float value)
{
A.elements[row * A.stride + col] = value;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 26
Programming Interface
cudaMalloc(&d_B.elements, size);
cudaMemcpy(d_B.elements, B.elements, size,
cudaMemcpyHostToDevice);
// Invoke kernel
dim3 dimBlock(BLOCK_SIZE, BLOCK_SIZE);
dim3 dimGrid(B.width / dimBlock.x, A.height / dimBlock.y);
MatMulKernel<<<dimGrid, dimBlock>>>(d_A, d_B, d_C);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 27
Programming Interface
BLOCK _SI ZE
B
B. height
BLOCK _SI ZE
BLOCK_SI ZE- 1
A C
0 col
0
BLOCK _SI ZE
Csub
blockRow
A.height
row
BLOCK_SI ZE-1
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 28
Programming Interface
1. Portable Memory
A block of page-locked memory can be used in conjunction with any device in the
system (see Multi-Device System for more details on multi-device systems), but by
default, the benefits of using page-locked memory described above are only available in
conjunction with the device that was current when the block was allocated (and with all
devices sharing the same unified address space, if any, as described in Unified Virtual
Address Space). To make these advantages available to all devices, the block needs to be
allocated by passing the flag cudaHostAllocPortable to cudaHostAlloc() or page-
locked by passing the flag cudaHostRegisterPortable to cudaHostRegister().
2. Write-Combining Memory
By default page-locked host memory is allocated as cacheable. It can optionally be
allocated as write-combining instead by passing flag cudaHostAllocWriteCombined
to cudaHostAlloc(). Write-combining memory frees up the host's L1 and L2 cache
resources, making more cache available to the rest of the application. In addition, write-
combining memory is not snooped during transfers across the PCI Express bus, which
can improve transfer performance by up to 40%.
Reading from write-combining memory from the host is prohibitively slow, so write-
combining memory should in general be used for memory that the host only writes to.
3. Mapped Memory
A block of page-locked host memory can also be mapped into the address space
of the device by passing flag cudaHostAllocMapped to cudaHostAlloc() or by
passing flag cudaHostRegisterMapped to cudaHostRegister(). Such a block
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 29
Programming Interface
has therefore in general two addresses: one in host memory that is returned by
cudaHostAlloc() or malloc(), and one in device memory that can be retrieved using
cudaHostGetDevicePointer() and then used to access the block from within a
kernel. The only exception is for pointers allocated with cudaHostAlloc() and whena
unified address space is used for the host and the device as mentioned in Unified Virtual
Address Space.
Accessing host memory directly from within a kernel has several advantages:
‣ There is no need to allocate a block in device memory and copy data between this
block and the block in host memory; data transfers are implicitly performed as
needed by the kernel;
‣ There is no need to use streams (see Concurrent Data Transfers) to overlap data
transfers with kernel execution; the kernel-originated data transfers automatically
overlap with kernel execution.
Since mapped page-locked memory is shared between host and device however, the
application must synchronize memory accesses using streams or events (see
Asynchronous Concurrent Execution) to avoid any potential read-after-write, write-
after-read, or write-after-write hazards.
To be able to retrieve the device pointer to any mapped page-locked memory, page-
locked memory mapping must be enabled by calling cudaSetDeviceFlags() with
the cudaDeviceMapHost flag before any other CUDA call is performed. Otherwise,
cudaHostGetDevicePointer() will return an error.
cudaHostGetDevicePointer() also returns an error if the device does not support
mapped page-locked host memory. Applications may query this capability by checking
the canMapHostMemory device property (see Device Enumeration), which is equal to 1
for devices that support mapped page-locked host memory.
Note that atomic functions (see Atomic Functions) operating on mapped page-locked
memory are not atomic from the point of view of the host or other devices.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 30
Programming Interface
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 31
Programming Interface
device copies are initiated using the standard memory copy functions with destination
and source addresses residing on the same device.
5. Streams
Applications manage the concurrent operations described above through streams. A
stream is a sequence of commands (possibly issued by different host threads) that
execute in order. Different streams, on the other hand, may execute their commands out
of order with respect to one another or concurrently; this behavior is not guaranteed and
should therefore not be relied upon for correctness (e.g., inter-kernel communication is
undefined).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 32
Programming Interface
2. Default Stream
Kernel launches and host <-> device memory copies that do not specify any stream
parameter, or equivalently that set the stream parameter to zero, are issued to the default
stream. They are therefore executed in order.
For code that is compiled using the --default-stream per-thread compilation flag
(or that defines the CUDA_API_PER_THREAD_DEFAULT_STREAM macro before including
CUDA headers (cuda.h and cuda_runtime.h)), the default stream is a regular stream
and each host thread has its own default stream.
For code that is compiled using the --default-stream legacy compilation flag, the
default stream is a special stream called the NULL stream and each device has a single
NULL stream used for all host threads. The NULL stream is special as it causes implicit
synchronization as described in Implicit Synchronization.
For code that is compiled without specifying a --default-stream compilation flag, --
default-stream legacy is assumed as the default.
3. Explicit Synchronization
There are various ways to explicitly synchronize streams with each other.
cudaDeviceSynchronize() waits until all preceding commands in all streams of all
host threads have completed.
cudaStreamSynchronize()takes a stream as a parameter and waits until all preceding
commands in the given stream have completed. It can be used to synchronize the host
with a specific stream, allowing other streams to continue executing on the device.
cudaStreamWaitEvent()takes a stream and an event as parameters (see Events for
a description of events)and makes all the commands added to the given stream after
the call to cudaStreamWaitEvent()delay their execution until the given event has
completed. The stream can be 0, in which case all the commands added to any stream
after the call to cudaStreamWaitEvent()wait on the event.
cudaStreamQuery()provides applications with a way to know if all preceding
commands in a stream have completed.
To avoid unnecessary slowdowns, all these synchronization functions are usually best
used for timing purposes or to isolate a launch or memory copy that is failing.
4. Implicit Synchronization
Two commands from different streams cannot run concurrently if any one of the
following operations is issued in-between them by the host thread:
‣ a page-locked host memory allocation,
‣ a device memory allocation,
‣ a device memory set,
‣ a memory copy between two addresses to the same device memory,
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 33
Programming Interface
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 34
Programming Interface
the second kernel launch is issued to stream[1] after the memory copy from device
to host is issued to stream[0], so it is blocked until the first kernel launch issued to
stream[0] is complete as per Implicit Synchronization. If the code is rewritten as
above, the kernel executions overlap (assuming the device supports concurrent kernel
execution) since the second kernel launch is issued to stream[1] before the memory copy
from device to host is issued to stream[0]. In that case however, the memory copy from
device to host issued to stream[0] only overlaps with the last thread blocks of the kernel
launch issued to stream[1] as per Implicit Synchronization, which can represent only a
small portion of the total execution time of the kernel.
3.2.5.5.6. Callbacks
The runtime provides a way to insert a callback at any point into a stream via
cudaStreamAddCallback(). A callback is a function that is executed on the host once
all commands issued to the stream before the callback have completed. Callbacks in
stream 0 are executed once all preceding tasks and commands issued in all streams
before the callback have completed.
The following code sample adds the callback function MyCallback to each of two
streams after issuing a host-to-device memory copy, a kernel launch and a device-to-host
memory copy into each stream. The callback will begin execution on the host after each
of the device-to-host memory copies completes.
The commands that are issued in a stream (or all commands issued to any stream if the
callback is issued to stream 0) after a callback do not start executing before the callback
has completed. The last parameter of cudaStreamAddCallback() is reserved for future
use.
A callback must not make CUDA API calls (directly or indirectly), as it might end up
waiting on itself if it makes such a call leading to a deadlock.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 35
Programming Interface
The following code sample obtains the allowable range of priorities for the current
device, and creates streams with the highest and lowest available priorities
// get the range of stream priorities for this device
int priority_high, priority_low;
cudaDeviceGetStreamPriorityRange(&priority_low, &priority_high);
// create streams with highest and lowest available priorities
cudaStream_t st_high, st_low;
cudaStreamCreateWithPriority(&st_high, cudaStreamNonBlocking, priority_high);
cudaStreamCreateWithPriority(&st_low, cudaStreamNonBlocking, priority_low);
6. Events
The runtime also provides a way to closely monitor the device's progress, as well as
perform accurate timing, by letting the application asynchronously record events at
any point in the program and query when these events are completed. An event has
completed when all tasks - or optionally, all commands in a given stream - preceding the
event have completed. Events in stream zero are completed after all preceding tasks and
commands in all streams are completed.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 36
Programming Interface
6. Multi-Device System
1. Device Enumeration
A host system can have multiple devices. The following code sample shows how to
enumerate these devices, query their properties, and determine the number of CUDA-
enabled devices.
int deviceCount;
cudaGetDeviceCount(&deviceCount);
int device;
for (device = 0; device < deviceCount; ++device) {
cudaDeviceProp deviceProp;
cudaGetDeviceProperties(&deviceProp, device);
printf("Device %d has compute capability %d.%d.\n",
device, deviceProp.major, deviceProp.minor);
}
A memory copy will succeed even if it is issued to a stream that is not associated to the
current device.
cudaEventRecord() will fail if the input event and input stream are associated to
different devices.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 37
Programming Interface
cudaEventElapsedTime() will fail if the two input events are associated to different
devices.
cudaEventSynchronize() and cudaEventQuery() will succeed even if the input
event is associated to a device that is different from the current device.
cudaStreamWaitEvent() will succeed even if the input stream and input event are
associated to different devices. cudaStreamWaitEvent() can therefore be used to
synchronize multiple devices with each other.
Each device has its own default stream (see Default Stream), so commands issued to
the default stream of a device may execute out of order or concurrently with respect to
commands issued to the default stream of any other device.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 38
Programming Interface
A copy (in the implicit NULL stream) between the memories of two different devices:
‣ does not start until all commands previously issued to either device have completed
and
‣ runs to completion before any commands (see Asynchronous Concurrent Execution)
issued after the copy to either device can start.
Consistent with the normal behavior of streams, an asynchronous copy between the
memories of two devices may overlap with copies or kernels in another stream.
Note that if peer-to-peer access is enabled between two devices via
cudaDeviceEnablePeerAccess() as described in Peer-to-Peer Memory Access, peer-
to-peer memory copy between these two devices no longer needs to be staged through
the host and is therefore faster.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 39
Programming Interface
Applications may query if the unified address space is used for a particular device by
checking that the unifiedAddressing device property (see Device Enumeration) is
equal to 1.
8. Interprocess Communication
Any device memory pointer or event handle created by a host thread can be directly
referenced by any other thread within the same process. It is not valid outside this
process however, and therefore cannot be directly referenced by threads belonging to a
different process.
To share device memory pointers and events across processes, an application must use
the Inter Process Communication API, which is described in detail in the reference
manual. The IPC API is only supported for 64-bit processes on Linux and for devices of
compute capability 2.0 and higher.
Using this API, an application can get the IPC handle for a given device memory
pointer using cudaIpcGetMemHandle(), pass it to another process using
standard IPC mechanisms (e.g., interprocess shared memory or files), and use
cudaIpcOpenMemHandle() to retrieve a device pointer from the IPC handle that is a
valid pointer within this other process. Event handles can be shared using similar entry
points.
An example of using the IPC API is where a single master process generates a batch
of input data, making the data available to multiple slave processes without requiring
regeneration or copying.
9. Error Checking
All runtime functions return an error code, but for an asynchronous function (see
Asynchronous Concurrent Execution), this error code cannot possibly report any of the
asynchronous errors that could occur on the device since the function returns before the
device has completed the task; the error code only reports errors that occur on the host
prior to executing the task, typically related to parameter validation; if an asynchronous
error occurs, it will be reported by some subsequent unrelated runtime function call.
The only way to check for asynchronous errors just after some asynchronous
function call is therefore to synchronize just after the call by calling
cudaDeviceSynchronize() (or by using any other synchronization mechanisms
described in Asynchronous Concurrent Execution) and checking the error code returned
by cudaDeviceSynchronize().
The runtime maintains an error variable for each host thread that is initialized to
cudaSuccess and is overwritten by the error code every time an error occurs (be it
a parameter validation error or an asynchronous error). cudaPeekAtLastError()
returns this variable. cudaGetLastError() returns this variable and resets it to
cudaSuccess.
Kernel launches do not return any error code, so cudaPeekAtLastError() or
cudaGetLastError() must be called just after the kernel launch to retrieve any
pre-launch errors. To ensure that any error returned by cudaPeekAtLastError()
or cudaGetLastError() does not originate from calls prior to the kernel launch,
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 40
Programming Interface
one has to make sure that the runtime error variable is set to cudaSuccess just before
the kernel launch, for example, by calling cudaGetLastError() just before the
kernel launch. Kernel launches are asynchronous, so to check for asynchronous
errors, the application must synchronize in-between the kernel launch and the call to
cudaPeekAtLastError() or cudaGetLastError().
Note that cudaErrorNotReady that may be returned by cudaStreamQuery() and
cudaEventQuery() is not considered an error and is therefore not reported by
cudaPeekAtLastError() or cudaGetLastError().
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 41
Programming Interface
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 42
Programming Interface
Cubemap Textures and Cubemap Layered Textures describe a special type of texture,
the cubemap texture.
Texture Gather describes a special texture fetch, texture gather.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 43
Programming Interface
The following code sample applies some simple transformation kernel to a texture.
// Simple transformation kernel
global void transformKernel(float* output,
cudaTextureObject_t texObj,
int width, int height,
float theta)
{
// Calculate normalized texture coordinates
unsigned int x = blockIdx.x * blockDim.x + threadIdx.x;
unsigned int y = blockIdx.y * blockDim.y + threadIdx.y;
float u = x / (float)width;
float v = y / (float)height;
// Transform coordinates
u -= 0.5f;
v -= 0.5f;
float tu = u * cosf(theta) - v * sinf(theta) + 0.5f;
float tv = v * cosf(theta) + u * sinf(theta) + 0.5f;
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 44
Programming Interface
// Host code
int main()
{
// Allocate CUDA array in device memory
cudaChannelFormatDesc channelDesc =
cudaCreateChannelDesc(32, 0, 0, 0,
cudaChannelFormatKindFloat);
cudaArray* cuArray;
cudaMallocArray(&cuArray, &channelDesc, width, height);
// Specify texture
struct cudaResourceDesc resDesc;
memset(&resDesc, 0, sizeof(resDesc));
resDesc.resType = cudaResourceTypeArray;
resDesc.res.array.array = cuArray;
// Invoke kernel
dim3 dimBlock(16, 16);
dim3 dimGrid((width + dimBlock.x - 1) / dimBlock.x,
(height + dimBlock.y - 1) / dimBlock.y);
transformKernel<<<dimGrid, dimBlock>>>(output,
texObj, width, height,
angle);
return 0;
}
3.2.11.1.2. Texture Reference API
Some of the attributes of a texture reference are immutable and must be known at
compile time; they are specified when declaring the texture reference. A texture
reference is declared at file scope as a variable of type texture:
texture<DataType, Type, ReadMode> texRef;
where:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 45
Programming Interface
where x, y, z, and w are equal to the number of bits of each component of the
returned value and f is:
‣ cudaChannelFormatKindSigned if these components are of signed integer
type,
‣ cudaChannelFormatKindUnsigned if they are of unsigned integer type,
‣ cudaChannelFormatKindFloat if they are of floating point type.
‣ See reference manual for sRGB, maxAnisotropy, mipmapFilterMode,
mipmapLevelBias, minMipmapLevelClamp, and maxMipmapLevelClamp.
normalized, addressMode, and filterMode may be directly modified in host code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 46
Programming Interface
Before a kernel can use a texture reference to read from texture memory, the
texture reference must be bound to a texture using cudaBindTexture() or
cudaBindTexture2D() for linear memory, or cudaBindTextureToArray() for CUDA
arrays. cudaUnbindTexture() is used to unbind a texture reference. Once a texture
reference has been unbound, it can be safely rebound to another array, even if kernels
that use the previously bound texture have not completed. It is recommended to allocate
two-dimensional textures in linear memory using cudaMallocPitch() and use the
pitch returned by cudaMallocPitch() as input parameter to cudaBindTexture2D().
The following code samples bind a 2D texture reference to linear memory pointed to by
devPtr:
‣ Using the low-level API:
texture<float, cudaTextureType2D,
cudaReadModeElementType> texRef;
textureReference* texRefPtr;
cudaGetTextureReference(&texRefPtr, &texRef);
cudaChannelFormatDesc channelDesc =
cudaCreateChannelDesc<float>();
size_t offset;
cudaBindTexture2D(&offset, texRefPtr, devPtr, &channelDesc,
width, height, pitch);
‣ Using the high-level API:
texture<float, cudaTextureType2D,
cudaReadModeElementType> texRef;
cudaChannelFormatDesc channelDesc =
cudaCreateChannelDesc<float>();
size_t offset;
cudaBindTexture2D(&offset, texRef, devPtr, channelDesc,
width, height, pitch);
The following code samples bind a 2D texture reference to a CUDA array cuArray:
‣ Using the low-level API:
texture<float, cudaTextureType2D,
cudaReadModeElementType> texRef;
textureReference* texRefPtr;
cudaGetTextureReference(&texRefPtr, &texRef);
cudaChannelFormatDesc channelDesc;
cudaGetChannelDesc(&channelDesc, cuArray);
cudaBindTextureToArray(texRef, cuArray, &channelDesc);
‣ Using the high-level API:
texture<float, cudaTextureType2D,
cudaReadModeElementType> texRef;
cudaBindTextureToArray(texRef, cuArray);
The format specified when binding a texture to a texture reference must match the
parameters specified when declaring the texture reference; otherwise, the results of
texture fetches are undefined.
There is a limit to the number of textures that can be bound to a kernel as specified in
Table 13.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 47
Programming Interface
The following code sample applies some simple transformation kernel to a texture.
// 2D float texture
texture<float, cudaTextureType2D, cudaReadModeElementType> texRef;
float u = x / (float)width;
float v = y / (float)height;
// Transform coordinates
u -= 0.5f;
v -= 0.5f;
float tu = u * cosf(theta) - v * sinf(theta) + 0.5f;
float tv = v * cosf(theta) + u * sinf(theta) + 0.5f;
// Host code
int main()
{
// Allocate CUDA array in device memory
cudaChannelFormatDesc channelDesc =
cudaCreateChannelDesc(32, 0, 0, 0,
cudaChannelFormatKindFloat);
cudaArray* cuArray;
cudaMallocArray(&cuArray, &channelDesc, width, height);
// Invoke kernel
dim3 dimBlock(16, 16);
dim3 dimGrid((width + dimBlock.x - 1) / dimBlock.x,
(height + dimBlock.y - 1) / dimBlock.y);
transformKernel<<<dimGrid, dimBlock>>>(output, width, height,
angle);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 48
Programming Interface
4. Layered Textures
A one-dimensional or two-dimensional layered texture (also known as texture array in
Direct3D and array texture in OpenGL) is a texture made up of a sequence of layers, all of
which are regular textures of same dimensionality, size, and data type.
A one-dimensional layered texture is addressed using an integer index and a floating-
point texture coordinate; the index denotes a layer within the sequence and the
coordinate addresses a texel within that layer. A two-dimensional layered texture is
addressed using an integer index and two floating-point texture coordinates; the index
denotes a layer within the sequence and the coordinates address a texel within that layer.
A layered texture can only be a CUDA array by calling cudaMalloc3DArray() with the
cudaArrayLayered flag (and a height of zero for one-dimensional layered texture).
Layered textures are fetched using the device functions described in tex1DLayered(),
tex1DLayered(), tex2DLayered(), and tex2DLayered(). Texture filtering (see Texture
Fetching) is done only within a layer, not across layers.
Layered textures are only supported on devices of compute capability 2.0 and higher.
5. Cubemap Textures
A cubemap texture is a special type of two-dimensional layered texture that has six layers
representing the faces of a cube:
‣ The width of a layer is equal to its height.
‣ The cubemap is addressed using three texture coordinates x, y, and z that are
interpreted as a direction vector emanating from the center of the cube and pointing
to one face of the cube and a texel within the layer corresponding to that face. More
specifically, the face is selected by the coordinate with largest magnitude m and the
corresponding layer is addressed using coordinates (s/m+1)/2 and (t/m+1)/2 where s
and t are defined in Table 1.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 49
Programming Interface
A layered texture can only be a CUDA array by calling cudaMalloc3DArray() with the
cudaArrayCubemap flag.
Cubemap textures are fetched using the device function described in texCubemap() and
texCubemap().
Cubemap textures are only supported on devices of compute capability 2.0 and higher.
7. Texture Gather
Texture gather is a special texture fetch that is available for two-dimensional textures
only. It is performed by the tex2Dgather() function, which has the same parameters
as tex2D(), plus an additional comp parameter equal to 0, 1, 2, or 3 (see tex2Dgather()
and tex2Dgather()). It returns four 32-bit numbers that correspond to the value of the
component comp of each of the four texels that would have been used for bilinear
filtering during a regular texture fetch. For example, if these texels are of values
(253, 20, 31, 255), (250, 25, 29, 254), (249, 16, 37, 253), (251, 22, 30, 250), and comp is 2,
tex2Dgather() returns (31, 29, 37, 30).
Texture gather is only supported for CUDA arrays created with the
cudaArrayTextureGather flag and of width and height less than the maximum
specified in Table 13 for texture gather, which is smaller than for regular texture fetch.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 50
Programming Interface
Texture gather is only supported on devices of compute capability 2.0 and higher.
2. Surface Memory
For devices of compute capability 2.0 and higher, a CUDA array (described in Cubemap
Surfaces), created with the cudaArraySurfaceLoadStore flag, can be read and written
via a surface object or surface reference using the functions described in Surface Functions.
Table 13 lists the maximum surface width, height, and depth depending on the compute
capability of the device.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 51
Programming Interface
The following code sample applies some simple transformation kernel to a texture.
// Simple copy kernel
global void copyKernel(cudaSurfaceObject_t inputSurfObj,
cudaSurfaceObject_t outputSurfObj,
int width, int height)
{
// Calculate surface coordinates
unsigned int x = blockIdx.x * blockDim.x + threadIdx.x;
unsigned int y = blockIdx.y * blockDim.y + threadIdx.y;
if (x < width && y < height) {
uchar4 data;
// Read from input surface
surf2Dread(&data, inputSurfObj, x * 4, y);
// Write to output surface
surf2Dwrite(data, outputSurfObj, x * 4, y);
}
}
// Host code
int main()
{
// Allocate CUDA arrays in device memory
cudaChannelFormatDesc channelDesc =
cudaCreateChannelDesc(8, 8, 8, 8,
cudaChannelFormatKindUnsigned);
cudaArray* cuInputArray;
cudaMallocArray(&cuInputArray, &channelDesc, width, height,
cudaArraySurfaceLoadStore);
cudaArray* cuOutputArray;
cudaMallocArray(&cuOutputArray, &channelDesc, width, height,
cudaArraySurfaceLoadStore);
// Specify surface
struct cudaResourceDesc resDesc;
memset(&resDesc, 0, sizeof(resDesc));
resDesc.resType = cudaResourceTypeArray;
// Invoke kernel
dim3 dimBlock(16, 16);
dim3 dimGrid((width + dimBlock.x - 1) / dimBlock.x,
(height + dimBlock.y - 1) / dimBlock.y);
copyKernel<<<dimGrid, dimBlock>>>(inputSurfObj,
outputSurfObj,
width, height);
return 0;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 52
Programming Interface
where Type specifies the type of the surface reference and is equal to
cudaSurfaceType1D, cudaSurfaceType2D, cudaSurfaceType3D,
cudaSurfaceTypeCubemap, cudaSurfaceType1DLayered,
cudaSurfaceType2DLayered, or cudaSurfaceTypeCubemapLayered; Type is an
optional argument which defaults to cudaSurfaceType1D. A surface reference can only
be declared as a static global variable and cannot be passed as an argument to a function.
Before a kernel can use a surface reference to access a CUDA array, the surface reference
must be bound to the CUDA array using cudaBindSurfaceToArray().
The following code samples bind a surface reference to a CUDA array cuArray:
‣ Using the low-level API:
surface<void, cudaSurfaceType2D> surfRef;
surfaceReference* surfRefPtr;
cudaGetSurfaceReference(&surfRefPtr, "surfRef");
cudaChannelFormatDesc channelDesc;
cudaGetChannelDesc(&channelDesc, cuArray);
cudaBindSurfaceToArray(surfRef, cuArray, &channelDesc);
‣ Using the high-level API:
surface<void, cudaSurfaceType2D> surfRef;
cudaBindSurfaceToArray(surfRef, cuArray);
A CUDA array must be read and written using surface functions of matching
dimensionality and type and via a surface reference of matching dimensionality;
otherwise, the results of reading and writing the CUDA array are undefined.
Unlike texture memory, surface memory uses byte addressing. This means that
the x-coordinate used to access a texture element via texture functions needs to be
multiplied by the byte size of the element to access the same element via a surface
function. For example, the element at texture coordinate x of a one-dimensional
floating-point CUDA array bound to a texture reference texRef and a surfacereference
surfRef is read using tex1d(texRef, x) via texRef, but surf1Dread(surfRef,
4*x) via surfRef. Similarly, the element at texture coordinate x and y of a two-
dimensional floating-point CUDA array bound to a texture reference texRef and a
surface reference surfRef is accessed using tex2d(texRef, x, y) via texRef, but
surf2Dread(surfRef, 4*x, y) via surfRef (the byte offset of the y-coordinate is
internally calculated from the underlying line pitch of the CUDA array).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 53
Programming Interface
The following code sample applies some simple transformation kernel to a texture.
// 2D surfaces
surface<void, 2> inputSurfRef;
surface<void, 2> outputSurfRef;
// Host code
int main()
{
// Allocate CUDA arrays in device memory
cudaChannelFormatDesc channelDesc =
cudaCreateChannelDesc(8, 8, 8, 8,
cudaChannelFormatKindUnsigned);
cudaArray* cuInputArray;
cudaMallocArray(&cuInputArray, &channelDesc, width, height,
cudaArraySurfaceLoadStore);
cudaArray* cuOutputArray;
cudaMallocArray(&cuOutputArray, &channelDesc, width, height,
cudaArraySurfaceLoadStore);
// Invoke kernel
dim3 dimBlock(16, 16);
dim3 dimGrid((width + dimBlock.x - 1) / dimBlock.x,
(height + dimBlock.y - 1) / dimBlock.y);
copyKernel<<<dimGrid, dimBlock>>>(width, height);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 54
Programming Interface
3. CUDA Arrays
CUDA arrays are opaque memory layouts optimized for texture fetching. They are one
dimensional, two dimensional, or three-dimensional and composed of elements, each of
which has 1, 2 or 4 components that may be signed or unsigned 8-, 16-, or 32-bit integers,
16-bit floats, or 32-bit floats. CUDA arrays are only accessible by kernels through texture
fetching as described in Texture Memory or surface reading and writing as described in
Surface Memory.
4. Read/Write Coherency
The texture and surface memory is cached (see Device Memory Accesses) and within
the same kernel call, the cache is not kept coherent with respect to global memory
writes and surface memory writes, so any texture fetch or surface read to an address
that has been written to via a global write or a surface write in the same kernel call
returns undefined data. In other words, a thread can safely read some texture or surface
memory location only if this memory location has been updated by a previous kernel
call or memory copy, but not if it has been previously updated by the same thread or
another thread from the same kernel call.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 55
Programming Interface
A mapped resource can be read from or written to by kernels using the device memory
address returned by cudaGraphicsResourceGetMappedPointer() for buffers and
cudaGraphicsSubResourceGetMappedArray() for CUDA arrays.
Accessing a resource through OpenGL, Direct3D, or another CUDA context while
it is mapped produces undefined results. OpenGL Interoperability and Direct3D
Interoperability give specifics for each graphics API and some code samples. SLI
Interoperability gives specifics for when the system is in SLI mode.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 56
Programming Interface
The following code sample uses a kernel to dynamically modify a 2D width x height
grid of vertices stored in a vertex buffer object:
GLuint positionsVBO;
struct cudaGraphicsResource* positionsVBO_CUDA;
int main()
{
// Initialize OpenGL and GLUT for device 0
// and make the OpenGL context current
...
glutDisplayFunc(display);
...
}
void display()
{
// Map buffer object for writing from CUDA
float4* positions;
cudaGraphicsMapResources(1, &positionsVBO_CUDA, 0);
size_t num_bytes;
cudaGraphicsResourceGetMappedPointer((void**)&positions,
&num_bytes,
positionsVBO_CUDA));
// Execute kernel
dim3 dimBlock(16, 16, 1);
dim3 dimGrid(width / dimBlock.x, height / dimBlock.y, 1);
createVertices<<<dimGrid, dimBlock>>>(positions, time,
width, height);
// Swap buffers
glutSwapBuffers();
glutPostRedisplay();
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 57
Programming Interface
void deleteVBO()
{
cudaGraphicsUnregisterResource(positionsVBO_CUDA);
glDeleteBuffers(1, &positionsVBO);
}
// Calculate uv coordinates
float u = x / (float)width;
float v = y / (float)height;
u = u * 2.0f - 1.0f;
v = v * 2.0f - 1.0f;
// Write positions
positions[y * width + x] = make_float4(u, w, v, 1.0f);
}
On Windows and for Quadro GPUs, cudaWGLGetDevice() can be used to retrieve the
CUDA device associated to the handle returned by wglEnumGpusNV(). Quadro GPUs
offer higher performance OpenGL interoperability than GeForce and Tesla GPUs in a
multi-GPU configuration where OpenGL rendering is performed on the Quadro GPU
and CUDA computations are performed on other GPUs in the system.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 58
Programming Interface
int main()
{
int dev;
// Initialize Direct3D
D3D = Direct3DCreate9Ex(D3D_SDK_VERSION);
// Create device
...
D3D->CreateDeviceEx(adapter, D3DDEVTYPE_HAL, hWnd,
D3DCREATE_HARDWARE_VERTEXPROCESSING,
¶ms, NULL, &device);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 59
Programming Interface
void Render()
{
// Map vertex buffer for writing from CUDA
float4* positions;
cudaGraphicsMapResources(1, &positionsVB_CUDA, 0);
size_t num_bytes;
cudaGraphicsResourceGetMappedPointer((void**)&positions,
&num_bytes,
positionsVB_CUDA));
// Execute kernel
dim3 dimBlock(16, 16, 1);
dim3 dimGrid(width / dimBlock.x, height / dimBlock.y, 1);
createVertices<<<dimGrid, dimBlock>>>(positions, time,
width, height);
void releaseVB()
{
cudaGraphicsUnregisterResource(positionsVB_CUDA);
positionsVB->Release();
}
// Calculate uv coordinates
float u = x / (float)width;
float v = y / (float)height;
u = u * 2.0f - 1.0f;
v = v * 2.0f - 1.0f;
// Write positions
positions[y * width + x] =
make_float4(u, w, v, int_as_float(0xff00ff00));
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 60
Programming Interface
int main()
{
int dev;
// Get a CUDA-enabled adapter
IDXGIFactory* factory;
CreateDXGIFactory( uuidof(IDXGIFactory), (void**)&factory);
IDXGIAdapter* adapter = 0;
for (unsigned int i = 0; !adapter; ++i) {
if (FAILED(factory->EnumAdapters(i, &adapter))
break;
if (cudaD3D10GetDevice(&dev, adapter) == cudaSuccess)
break;
adapter->Release();
}
factory->Release();
cudaGraphicsResourceSetMapFlags(positionsVB_CUDA,
cudaGraphicsMapFlagsWriteDiscard);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 61
Programming Interface
void Render()
{
// Map vertex buffer for writing from CUDA
float4* positions;
cudaGraphicsMapResources(1, &positionsVB_CUDA, 0);
size_t num_bytes;
cudaGraphicsResourceGetMappedPointer((void**)&positions,
&num_bytes,
positionsVB_CUDA));
// Execute kernel
dim3 dimBlock(16, 16, 1);
dim3 dimGrid(width / dimBlock.x, height / dimBlock.y, 1);
createVertices<<<dimGrid, dimBlock>>>(positions, time,
width, height);
void releaseVB()
{
cudaGraphicsUnregisterResource(positionsVB_CUDA);
positionsVB->Release();
}
// Calculate uv coordinates
float u = x / (float)width;
float v = y / (float)height;
u = u * 2.0f - 1.0f;
v = v * 2.0f - 1.0f;
// Write positions
positions[y * width + x] =
make_float4(u, w, v, int_as_float(0xff00ff00));
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 62
Programming Interface
int main()
{
int dev;
// Get a CUDA-enabled adapter
IDXGIFactory* factory;
CreateDXGIFactory( uuidof(IDXGIFactory), (void**)&factory);
IDXGIAdapter* adapter = 0;
for (unsigned int i = 0; !adapter; ++i) {
if (FAILED(factory->EnumAdapters(i, &adapter))
break;
if (cudaD3D11GetDevice(&dev, adapter) == cudaSuccess)
break;
adapter->Release();
}
factory->Release();
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 63
Programming Interface
void Render()
{
// Map vertex buffer for writing from CUDA
float4* positions;
cudaGraphicsMapResources(1, &positionsVB_CUDA, 0);
size_t num_bytes;
cudaGraphicsResourceGetMappedPointer((void**)&positions,
&num_bytes,
positionsVB_CUDA));
// Execute kernel
dim3 dimBlock(16, 16, 1);
dim3 dimGrid(width / dimBlock.x, height / dimBlock.y, 1);
createVertices<<<dimGrid, dimBlock>>>(positions, time,
width, height);
void releaseVB()
{
cudaGraphicsUnregisterResource(positionsVB_CUDA);
positionsVB->Release();
}
// Calculate uv coordinates
float u = x / (float)width;
float v = y / (float)height;
u = u * 2.0f - 1.0f;
v = v * 2.0f - 1.0f;
// Write positions
positions[y * width + x] =
make_float4(u, w, v, int_as_float(0xff00ff00));
}
3.2.12.3. SLI Interoperability
In a system with multiple GPUs, all CUDA-enabled GPUs are accessible via the CUDA
driver and runtime as separate devices. There are however special considerations as
described below when the system is in SLI mode.
First, an allocation in one CUDA device on one GPU will consume memory on other
GPUs that are part of the SLI configuration of the Direct3D or OpenGL device. Because
of this, allocations may fail earlier than otherwise expected.
Second, applications should create multiple CUDA contexts, one for each GPU in the SLI
configuration. While this is not a strict requirement, it avoids unnecessary data transfers
between devices. The application can use the cudaD3D[9|10|11]GetDevices() for
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 64
Programming Interface
Direct3D and cudaGLGetDevices() for OpenGL set of calls to identify the CUDA
device handle(s) for the device(s) that are performing the rendering in the current
and next frame. Given this information the application will typically choose the
appropriate device and map Direct3D or OpenGL resources to the CUDA device
returned by cudaD3D[9|10|11]GetDevices() or cudaGLGetDevices() when the
deviceList parameter is set to cudaD3D[9|10|11]DeviceListCurrentFrame or
cudaGLDeviceListCurrentFrame.
Please note that resource returned from cudaGraphicsD9D[9|10|
11]RegisterResource and cudaGraphicsGLRegister[Buffer|Image] must be
only used on device the registration happened. Therefore on SLI configurations when
data for different frames is computed on different CUDA devices it is necessary to
register the resources for each separatly.
See Direct3D Interoperability and OpenGL Interoperability for details on how the
CUDA runtime interoperate with Direct3D and OpenGL, respectively.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 65
Programming Interface
Compatible Incompatible
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 66
Programming Interface
5. Mode Switches
GPUs that have a display output dedicate some DRAM memory to the so-called primary
surface, which is used to refresh the display device whose output is viewed by the user.
When users initiate a mode switch of the display by changing the resolution or bit depth
of the display (using NVIDIA control panel or the Display control panel on Windows),
the amount of memory needed for the primary surface changes. For example, if the user
changes the display resolution from 1280x1024x32-bit to 1600x1200x32-bit, the system
must dedicate 7.68 MB to the primary surface rather than 5.24 MB. (Full-screen graphics
applications running with anti-aliasing enabled may require much more display
memory for the primary surface.) On Windows, other events that may initiate display
mode switches include launching a full-screen DirectX application, hitting Alt
+Tab to task switch away from a full-screen DirectX application, or hitting Ctrl+Alt+Del
to lock the computer.
If a mode switch increases the amount of memory needed for the primary surface, the
system may have to cannibalize memory allocations dedicated to CUDA applications.
Therefore, a mode switch results in any call to the CUDA runtime to fail and return an
invalid context error.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 67
Chapter 4.
HARDWARE IMPLEMENTATION
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 68
Hardware Implementation
way a block is partitioned into warps is always the same; each warp contains threads
of consecutive, increasing thread IDs with the first warp containing thread 0. Thread
Hierarchy describes how thread IDs relate to thread indices in the block.
A warp executes one common instruction at a time, so full efficiency is realized when
all 32 threads of a warp agree on their execution path. If threads of a warp diverge via a
data-dependent conditional branch, the warp serially executes each branch path taken,
disabling threads that are not on that path, and when all paths complete, the threads
converge back to the same execution path. Branch divergence occurs only within a
warp; different warps execute independently regardless of whether they are executing
common or disjoint code paths.
The SIMT architecture is akin to SIMD (Single Instruction, Multiple Data) vector
organizations in that a single instruction controls multiple processing elements. A key
difference is that SIMD vector organizations expose the SIMD width to the software,
whereas SIMT instructions specify the execution and branching behavior of a single
thread. In contrast with SIMD vector machines, SIMT enables programmers to write
thread-level parallel code for independent, scalar threads, as well as data-parallel code
for coordinated threads. For the purposes of correctness, the programmer can essentially
ignore the SIMT behavior; however, substantial performance improvements can be
realized by taking care that the code seldom requires threads in a warp to diverge. In
practice, this is analogous to the role of cache lines in traditional code: Cache line size
can be safely ignored when designing for correctness but must be considered in the code
structure when designing for peak performance. Vector architectures, on the other hand,
require the software to coalesce loads into vectors and manage divergence manually.
Notes
The threads of a warp that are on that warp's current execution path are called the active
threads, whereas threads not on the current path are inactive (disabled). Threads can be
inactive because they have exited earlier than other threads of their warp, or because
they are on a different branch path than the branch path currently executed by the warp,
or because they are the last threads of a block whose number of threads is not a multiple
of the warp size.
If a non-atomic instruction executed by a warp writes to the same location in global or
shared memory for more than one of the threads of the warp, the number of serialized
writes that occur to that location varies depending on the compute capability of the
device (see Compute Capability 2.x, Compute Capability 3.x, Compute Capability 5.x,
and Compute Capability 6.x), and which thread performs the final write is undefined.
If an atomic instruction executed by a warp reads, modifies, and writes to the same
location in global memory for more than one of the threads of the warp, each read/
modify/write to that location occurs and they are all serialized, but the order in which
they occur is undefined.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 69
Hardware Implementation
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 70
Chapter 5.
PERFORMANCE GUIDELINES
2. Maximize Utilization
To maximize utilization the application should be structured in a way that it exposes
as much parallelism as possible and efficiently maps this parallelism to the various
components of the system to keep them busy most of the time.
1. Application Level
At a high level, the application should maximize parallel execution between the host, the
devices, and the bus connecting the host to the devices, by using asynchronous functions
calls and streams as described in Asynchronous Concurrent Execution. It should assign
to each processor the type of work it does best: serial workloads to the host; parallel
workloads to the devices.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 71
Performance Guidelines
For the parallel workloads, at points in the algorithm where parallelism is broken
because some threads need to synchronize in order to share data with each other,
there are two cases: Either these threads belong to the same block, in which case they
should use syncthreads() and share data through shared memory within the same
kernel invocation, or they belong to different blocks, in which case they must share data
through global memory using two separate kernel invocations, one for writing to and
one for reading from global memory. The second case is much less optimal since it adds
the overhead of extra kernel invocations and global memory traffic. Its occurrence
should therefore be minimized by mapping the algorithm to the CUDA programming
model in such a way that the computations that require inter-thread communication are
performed within a single thread block as much as possible.
2. Device Level
At a lower level, the application should maximize parallel execution between the
multiprocessors of a device.
Multiple kernels can execute concurrently on a device, so maximum utilization can
also be achieved by using streams to enable enough kernels to execute concurrently as
described in Asynchronous Concurrent Execution.
3. Multiprocessor Level
At an even lower level, the application should maximize parallel execution between the
various functional units within a multiprocessor.
As described in Hardware Multithreading, a GPU multiprocessor relies on thread-
level parallelism to maximize utilization of its functional units. Utilization is therefore
directly linked to the number of resident warps. At every instruction issue time, a warp
scheduler selects a warp that is ready to execute its next instruction, if any, and issues
the instruction to the active threads of the warp. The number of clock cycles it takes for
a warp to be ready to execute its next instruction is called the latency, and full utilization
is achieved when all warp schedulers always have some instruction to issue for some
warp at every clock cycle during that latency period, or in other words, when latency is
completely "hidden". The number of instructions required to hide a latency of L clock
cycles depends on the respective throughputs of these instructions (see Arithmetic
Instructions for the throughputs of various arithmetic instructions); assuming maximum
throughput for all instructions, it is:
‣ L for devices of compute capability 2.0 since a multiprocessor issues one instruction
per warp over two clock cycles for two warps at a time, as mentioned in Compute
Capability 2.x,
‣ 2L for devices of compute capability 2.1 since a multiprocessor issues a pair of
instructions per warp over two clock cycles for two warps at a time, as mentioned in
Compute Capability 2.x,
‣ 8L for devices of compute capability 3.x since a multiprocessor issues a pair of
instructions per warp over one clock cycle for four warps at a time, as mentioned in
Compute Capability 3.x.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 72
Performance Guidelines
For devices of compute capability 2.0, the two instructions issued every other cycle are
for two different warps. For devices of compute capability 2.1, the four instructions
issued every other cycle are two pairs for two different warps, each pair being for the
same warp.
For devices of compute capability 3.x, the eight instructions issued every cycle are four
pairs for four different warps, each pair being for the same warp.
The most common reason a warp is not ready to execute its next instruction is that the
instruction's input operands are not available yet.
If all input operands are registers, latency is caused by register dependencies, i.e., some
of the input operands are written by some previous instruction(s) whose execution has
not completed yet. In the case of a back-to-back register dependency (i.e., some input
operand is written by the previous instruction), the latency is equal to the execution
time of the previous instruction and the warp schedulers must schedule instructions for
different warps during that time. Execution time varies depending on the instruction,
but it is typically about 22 clock cycles for devices of compute capability 2.x and about
11 clock cycles for devices of compute capability 3.x, which translates to 22 warps for
devices of compute capability 2.x and 44 warps for devices of compute capability 3.x
and higher (still assuming that warps execute instructions with maximum throughput,
otherwise fewer warps are needed). For devices of compute capability 2.1 and higher,
this is also assuming enough instruction-level parallelism so that schedulers are always
able to issue pairs of instructions for each warp.
If some input operand resides in off-chip memory, the latency is much higher: 400
to 800 clock cycles for devices of compute capability 2.x and about 200 to 400 clock
cycles for devices of compute capability 3.x. The number of warps required to keep the
warp schedulers busy during such high latency periods depends on the kernel code
and its degree of instruction-level parallelism. In general, more warps are required
if the ratio of the number of instructions with no off-chip memory operands (i.e.,
arithmetic instructions most of the time) to the number of instructions with off-chip
memory operands is low (this ratio is commonly called the arithmetic intensity of the
program). For example, assume this ratio is 30, also assume the latencies are 600 cycles
on devices of compute capability 2.x and 300 cycles on devices of compute capability
3.x. Then about 20 warps are required for devices of compute capability 2.x and about
40 for devices of compute capability 3.x (with the same assumptions as in the previous
paragraph).
Another reason a warp is not ready to execute its next instruction is that it is waiting
at some memory fence (Memory Fence Functions) or synchronization point (Memory
Fence Functions). A synchronization point can force the multiprocessor to idle as
more and more warps wait for other warps in the same block to complete execution of
instructions prior to the synchronization point. Having multiple resident blocks per
multiprocessor can help reduce idling in this case, as warps from different blocks do not
need to wait for each other at synchronization points.
The number of blocks and warps residing on each multiprocessor for a given kernel
call depends on the execution configuration of the call (Execution Configuration), the
memory resources of the multiprocessor, and the resource requirements of the kernel as
described in Hardware Multithreading. Register and shared memory usage are reported
by the compiler when compiling with the -ptxas-options=-v option.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 73
Performance Guidelines
The total amount of shared memory required for a block is equal to the sum of the
amount of statically allocated shared memory and the amount of dynamically allocated
shared memory.
The number of registers used by a kernel can have a significant impact on the number
of resident warps. For example, for devices of compute capability 2.x, if a kernel uses
32 registers and each block has 512 threads and requires very little shared memory,
then two blocks (i.e., 32 warps) can reside on the multiprocessor since they require
2x512x32 registers, which exactly matches the number of registers available on the
multiprocessor. But as soon as the kernel uses one more register, only one block (i.e.,
16 warps) can be resident since two blocks would require 2x512x17 registers, which are
more registers than are available on the multiprocessor. Therefore, the compiler attempts
to minimize register usage while keeping register spilling (see Device Memory Accesses)
and the number of instructions to a minimum. Register usage can be controlled using the
maxrregcount compiler option or launch bounds as described in LaunchBounds.
Each double variable and each long long variable uses two registers.
The effect of execution configuration on performance for a given kernel call generally
depends on the kernel code. Experimentation is therefore recommended. Applications
can also parameterize execution configurations based on register file size and shared
memory size, which depends on the compute capability of the device, as well as on the
number of multiprocessors and memory bandwidth of the device, all of which can be
queried using the runtime (see reference manual).
The number of threads per block should be chosen as a multiple of the warp size to
avoid wasting computing resources with under-populated warps as much as possible.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 74
Performance Guidelines
The following code sample calculates the occupancy of MyKernel. It then reports the
occupancy level with the ratio between concurrent warps versus maximum warps per
multiprocessor.
// Device code
global void MyKernel(int *d, int *a, int *b)
{
int idx = threadIdx.x + blockIdx.x * blockDim.x;
d[idx] = a[idx] * b[idx];
}
// Host code
int main()
{
int numBlocks; // Occupancy in terms of active blocks
int blockSize = 32;
cudaGetDevice(&device);
cudaGetDeviceProperties(&prop, device);
cudaOccupancyMaxActiveBlocksPerMultiprocessor(
&numBlocks,
MyKernel,
blockSize,
0);
std::cout << "Occupancy: " << (double)activeWarps / maxWarps * 100 << "%" <<
std::endl;
return 0;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 75
Performance Guidelines
// Host code
int launchMyKernel(int *array, int arrayCount)
{
int blockSize; // The launch configurator returned block size
int minGridSize; // The minimum grid size needed to achieve the
// maximum occupancy for a full device
// launch
int gridSize; // The actual grid size needed, based on input
// size
cudaOccupancyMaxPotentialBlockSize(
&minGridSize,
&blockSize,
(void*)MyKernel,
0,
arrayCount);
return 0;
}
The CUDA Toolkit also provides a self-documenting, standalone occupancy calculator
and launch configurator implementation in <CUDA_Toolkit_Path>/include/
cuda_occupancy.h for any use cases that cannot depend on the CUDA software stack.
A spreadsheet version of the occupancy calculator is also provided. The spreadsheet
version is particularly useful as a learning tool that visualizes the impact of changes
to the parameters that affect occupancy (block size, registers per thread, and shared
memory per thread).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 76
Performance Guidelines
on devices of compute capability 2.x and 3.x, L2 cache available on devices of compute
capability 2.x and higher, texture cache and constant cache available on all devices).
Shared memory is equivalent to a user-managed cache: The application explicitly
allocates and accesses it. As illustrated in CUDA C Runtime, a typical programming
pattern is to stage data coming from device memory into shared memory; in other
words, to have each thread of a block:
‣ Load data from device memory to shared memory,
‣ Synchronize with all the other threads of the block so that each thread can safely
read shared memory locations that were populated by different threads,
‣ Process the data in shared memory,
‣ Synchronize again if necessary to make sure that shared memory has been updated
with the results,
‣ Write the results back to device memory.
For some applications (e.g., for which global memory access patterns are data-
dependent), a traditional hardware-managed cache is more appropriate to exploit data
locality. As mentioned in Compute Capability 2.x and Compute Capability 3.x, for
devices of compute capability 2.x and 3.x, the same on-chip memory is used for both L1
and shared memory, and how much of it is dedicated to L1 versus shared memory is
configurable for each kernel call.
The throughput of memory accesses by a kernel can vary by an order of magnitude
depending on access pattern for each type of memory. The next step in maximizing
memory throughput is therefore to organize memory accesses as optimally as possible
based on the optimal memory access patterns described in Device Memory Accesses.
This optimization is especially important for global memory accesses as global memory
bandwidth is low, so non-optimal global memory accesses have a higher impact on
performance.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 77
Performance Guidelines
are and that the mapped memory is read or written only once, using mapped page-
locked memory instead of explicit copies between device and host memory can be a win
for performance.
On integrated systems where device memory and host memory are physically the same,
any copy between host and device memory is superfluous and mapped page-locked
memory should be used instead. Applications may query a device is integrated by
checking that the integrated device property (see Device Enumeration) is equal to 1.
Global Memory
Global memory resides in device memory and device memory is accessed via 32-, 64-, or
128-byte memory transactions. These memory transactions must be naturally aligned:
Only the 32-, 64-, or 128-byte segments of device memory that are aligned to their size
(i.e., whose first address is a multiple of their size) can be read or written by memory
transactions.
When a warp executes an instruction that accesses global memory, it coalesces the
memory accesses of the threads within the warp into one or more of these memory
transactions depending on the size of the word accessed by each thread and the
distribution of the memory addresses across the threads. In general, the more
transactions are necessary, the more unused words are transferred in addition to the
words accessed by the threads, reducing the instruction throughput accordingly. For
example, if a 32-byte memory transaction is generated for each thread's 4-byte access,
throughput is divided by 8.
How many transactions are necessary and how much throughput is ultimately affected
varies with the compute capability of the device. Compute Capability 2.x, Compute
Capability 3.x, Compute Capability 5.x and Compute Capability 6.x give more details on
how global memory accesses are handled for various compute capabilities.
To maximize global memory throughput, it is therefore important to maximize
coalescing by:
‣ Following the most optimal access patterns based on Compute Capability 2.x,
Compute Capability 3.x , Compute Capability 5.x and Compute Capability 6.x,
‣ Using data types that meet the size and alignment requirement detailed in Device
Memory Accesses,
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 78
Performance Guidelines
‣ Padding data in some cases, for example, when accessing a two-dimensional array
as described in Device Memory Accesses.
or
struct align (16) {
float x;
float y;
float z;
};
Any address of a variable residing in global memory or returned by one of the memory
allocation routines from the driver or runtime API is always aligned to at least 256 bytes.
Reading non-naturally aligned 8-byte or 16-byte words produces incorrect results (off by
a few words), so special care must be taken to maintain alignment of the starting address
of any value or array of values of these types. A typical case where this might be easily
overlooked is when using some custom global memory allocation scheme, whereby the
allocations of multiple arrays (with multiple calls to cudaMalloc() or cuMemAlloc())
is replaced by the allocation of a single large block of memory partitioned into multiple
arrays, in which case the starting address of each array is offset from the block's starting
address.
Two-Dimensional Arrays
A common global memory access pattern is when each thread of index (tx,ty) uses the
following address to access one element of a 2D array of width width, located at address
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 79
Performance Guidelines
BaseAddress of type type* (where type meets the requirement described in Maximize
Utilization):
BaseAddress + width * ty + tx
For these accesses to be fully coalesced, both the width of the thread block and the width
of the array must be a multiple of the warp size.
In particular, this means that an array whose width is not a multiple of this size will be
accessed much more efficiently if it is actually allocated with a width rounded up to the
closest multiple of this size and its rows padded accordingly. The cudaMallocPitch()
and cuMemAllocPitch() functions and associated memory copy functions described in
the reference manual enable programmers to write non-hardware-dependent code to
allocate arrays that conform to these constraints.
Local Memory
Local memory accesses only occur for some automatic variables as mentioned in
Variable Type Qualifiers. Automatic variables that the compiler is likely to place in local
memory are:
‣ Arrays for which it cannot determine that they are indexed with constant quantities,
‣ Large structures or arrays that would consume too much register space,
‣ Any variable if the kernel uses more registers than available (this is also known as
register spilling).
Inspection of the PTX assembly code (obtained by compiling with the -ptx or-
keep option) will tell if a variable has been placed in local memory during the first
compilation phases as it will be declared using the .local mnemonic and accessed
using the ld.local and st.local mnemonics. Even if it has not, subsequent
compilation phases might still decide otherwise though if they find it consumes too
much register space for the targeted architecture: Inspection of the cubin object using
cuobjdump will tell if this is the case. Also, the compiler reports total localmemory
usage per kernel (lmem) when compiling with the --ptxas-options=-v option. Note
that some mathematical functions have implementation paths that might access local
memory.
The local memory space resides in device memory, so local memory accesses have same
high latency and low bandwidth as global memory accesses and are subject to the same
requirements for memory coalescing as described in Device Memory Accesses.
Local memory is however organized such that consecutive 32-bit words are accessed
by consecutive thread IDs. Accesses are therefore fully coalesced as long as all threads
in a warp access the same relative address (e.g., same index in an array variable, same
member in a structure variable).
On devices of compute capability 2.x and 3.x, local memory accesses are always cached
in L1 and L2 in the same way as global memory accesses (see Compute Capability 2.x
and Compute Capability 3.x).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 80
Performance Guidelines
On devices of compute capability 5.x and 6.x, local memory accesses are always cached
in L2 in the same way as global memory accesses (see Compute Capability 5.x and
Compute Capability 6.x).
Shared Memory
Because it is on-chip, shared memory has much higher bandwidth and much lower
latency than local or global memory.
To achieve high bandwidth, shared memory is divided into equally-sized memory
modules, called banks, which can be accessed simultaneously. Any memory read or
write request made of n addresses that fall in n distinct memory banks can therefore be
serviced simultaneously, yielding an overall bandwidth that is n times as high as the
bandwidth of a single module.
However, if two addresses of a memory request fall in the same memory bank, there is a
bank conflict and the access has to be serialized. The hardware splits a memory request
with bank conflicts into as many separate conflict-free requests as necessary, decreasing
throughput by a factor equal to the number of separate memory requests. If the number
of separate memory requests is n, the initial memory request is said to cause n-way bank
conflicts.
To get maximum performance, it is therefore important to understand how memory
addresses map to memory banks in order to schedule the memory requests so as
to minimize bank conflicts. This is described in Compute Capability 2.x, Compute
Capability 3.x, Compute Capability 5.x, and Compute Capability 6.x for devices of
compute capability 2.x, 3.x, 5.x and 6.x, respectively.
Constant Memory
The constant memory space resides in device memory and is cached in the constant
cache mentioned in Compute Capability 2.x.
A request is then split into as many separate requests as there are different memory
addresses in the initial request, decreasing throughput by a factor equal to the number
of separate requests.
The resulting requests are then serviced at the throughput of the constant cache in case
of a cache hit, or at the throughput of device memory otherwise.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 81
Performance Guidelines
Also, it is designed for streaming fetches with a constant latency; a cache hit reduces
DRAM bandwidth demand but not fetch latency.
Reading device memory through texture or surface fetching present some benefits that
can make it an advantageous alternative to reading device memory from global or
constant memory:
‣ If the memory reads do not follow the access patterns that global or constant
memory reads must follow to get good performance, higher bandwidth can be
achieved providing that there is locality in the texture fetches or surface reads;
‣ Addressing calculations are performed outside the kernel by dedicated units;
‣ Packed data may be broadcast to separate variables in a single operation;
‣ 8-bit and 16-bit integer input data may be optionally converted to 32 bit floating-
point values in the range [0.0, 1.0] or [-1.0, 1.0] (see Texture Memory).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 82
Performance Guidelines
Compute Capability
2.0 2.1 3.0, 3.5, 5.0, 5.3 6.0 6.1 6.2
3.2 3.7 5.2
16-bit floating-
point add,
multiply, N/A N/A N/A N/A N/A 256 128 2 256
multiply-add
32-bit floating-
point add,
multiply, 32 48 192 192 128 128 64 128 128
multiply-add
64-bit floating-
point add,
multiply, 161 4 8 642 4 4 32 4 4
multiply-add
32-bit floating-
point reciprocal,
reciprocal
square root,
base-2 logarithm
( log2f), base 2 4 8 32 32 32 32 16 32 32
exponential
(exp2f), sine
( sinf), cosine (
cosf)
32-bit integer
add, extended-
precision add,
subtract, 32 48 160 160 128 128 64 128 128
extended-
precision
subtract
32-bit integer
multiply,
multiply-add, Multiple Multiple Multiple Multiple Multiple
extended- 16 16 32 32
instruct. instruct. instruct. instruct. instruct.
precision
multiply-add
24-bit integer Multiple Multiple Multiple Multiple Multiple Multiple Multiple Multiple Multiple
multiply instruct. instruct. instruct. instruct. instruct. instruct. instruct. instruct. instruct.
( [u]mul24)
1
4 for GeForce GPUs
2
8 for GeForce GPUs
3
32 for GeForce GPUs
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 83
Performance Guidelines
Compute Capability
2.0 2.1 3.0, 3.5, 5.0, 5.3 6.0 6.1 6.2
3.2 3.7 5.2
32-bit integer bit
reverse, bit field 16 16 32 32 64 64 32 64 64
extract/insert
32-bit bitwise 32 48 160 160 128 128 64 128 128
AND, OR, XOR
count of leading
zeros, most Multiple Multiple Multiple Multiple Multiple
significant non- 16 16 32 32
instruct. instruct. instruct. instruct. instruct.
sign bit
population count 16 16 32 32 32 32 16 32 32
warp shuffle N/A N/A 32 32 32 32 32 32 32
sum of absolute 16 16 32 32 64 64 32 64 64
difference
SIMD video Multiple Multiple Multiple Multiple Multiple
instructions N/A N/A 160 160 instruct. instruct. instruct. instruct. instruct.
vabsdiff2
Other instructions and functions are implemented on top of the native instructions.
The implementation may be different for devices of different compute capabilities, and
the number of native instructions after compilation may fluctuate with every compiler
version. For complicated functions, there can be multiple code paths depending on
input. cuobjdump can be used to inspect a particular implementation in a cubin object.
The implementation of some functions are readily available on the CUDA header files
(math_functions.h, device_functions.h, ...).
In general, code compiled with -ftz=true (denormalized numbers are flushed tozero)
tends to have higher performance than code compiled with -ftz=false. Similarly,
4
4 for GeForce GPUs
5
8 for GeForce GPUs
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 84
Performance Guidelines
code compiled with -prec div=false (less precise division) tends to have higher
performance code than code compiled with -prec div=true, and code compiled
with -prec-sqrt=false (less precise square root) tends to have higherperformance
than code compiled with -prec-sqrt=true. The nvcc user manual describes these
compilation flags in more details.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 85
Performance Guidelines
Integer Arithmetic
Integer division and modulo operation are costly as they compiler to up to 20
instructions. They can be replaced with bitwise operations in some cases: If n is a power
of 2, (i/n) is equivalent to (i>>log2(n)) and (i%n) is equivalent to (i&(n-1)); the
compiler will perform these conversions if n is literal.
brev and popc map to a single instruction and brevll and popcll to a few
instructions.
[u]mul24 are legacy intrinsic functions that have no longer any reason to be used.
Type Conversion
Sometimes, the compiler must insert conversion instructions, introducing additional
execution cycles. This is the case for:
‣ Functions operating on variables of type char or short whose operands generally
need to be converted to int,
‣ Double-precision floating-point constants (i.e., those constants defined without
any type suffix) used as input to single-precision floating-point computations (as
mandated by C/C++ standards).
This last case can be avoided by using single-precision floating-point constants, defined
with an f suffix such as 3.141592653589793f, 1.0f, 0.5f.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 86
Performance Guidelines
Sometimes, the compiler may unroll loops or it may optimize out if or switch
statements by using branch predication instead, as detailed below. In these cases, no
warp can ever diverge. The programmer can also control loop unrolling using the
#pragma unroll directive (see #pragma unroll).
When using branch predication none of the instructions whose execution depends on
the controlling condition gets skipped. Instead, each of them is associated with a per-
thread condition code or predicate that is set to true or false based on the controlling
condition and although each of these instructions gets scheduled for execution, only
the instructions with a true predicate are actually executed. Instructions with a false
predicate do not write results, and also do not evaluate addresses or read operands.
The compiler replaces a branch instruction with predicated instructions only if the
number of instructions controlled by the branch condition is less or equal to a certain
threshold: If the compiler determines that the condition is likely to produce many
divergent warps, this threshold is 7, otherwise it is 4.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 87
Appendix A.
CUDA-ENABLED GPUS
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 88
Appendix B.
C LANGUAGE EXTENSIONS
1. device
The device qualifier declares a function that is:
B.1.2. global
The global qualifier declares a function as being a kernel. Such a function is:
B.1.3. host
The host qualifier declares a function that is:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 89
C Language Extensions
1. device
The device qualifier declares a variable that resides on the device.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 90
C Language Extensions
At most one of the other type qualifiers defined in the next two sections may be used
together with device to further specify which memory space the variable belongs
to. If none of them is present, the variable:
B.2.2. constant
The constant qualifier, optionally used together with device , declares a
variable that:
B.2.3. shared
The shared qualifier, optionally used together with device , declares a
variable that:
the size of the array is determined at launch time (see Execution Configuration). All
variables declared in this fashion, start at the same address in memory, so that the layout
of the variables in the array must be explicitly managed through offsets. For example, if
one wants the equivalent of
short array0[128];
float array1[64];
int array2[256];
in dynamically allocated shared memory, one could declare and initialize the arrays the
following way:
extern shared float array[];
device void func() // device or global function
{ =
short* array0 = (short*)array;
float* array1 = (float*)&array0[128];
int* array2 (int*)&array1[64];
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 91
C Language Extensions
Note that pointers need to be aligned to the type they point to, so the following code, for
example, does not work since array1 is not aligned to 4 bytes.
extern shared float array[];
device void func() // device or global function
{ =
short* array0 = (short*)array;
float* array1 (float*)&array0[127];
}
Alignment requirements for the built-in vector types are listed in Table 3.
B.2.4. managed
The managed qualifier, optionally used together with device , declares a
variable that:
‣ Can be referenced from both device and host code, e.g., its address can be taken or it
can be read or written directly from a device or host function.
‣ Has the lifetime of an application.
See managed Qualifier for more details.
B.2.5. restrict
nvcc supports restricted pointers via the restrict keyword.
Restricted pointers were introduced in C99 to alleviate the aliasing problem that exists in
C-type languages, and which inhibits all kind of optimization from code re-ordering to
common sub-expression elimination.
Here is an example subject to the aliasing issue, where use of restricted pointer can help
the compiler to reduce the number of instructions:
void foo(const float* a,
const float* b,
float* c)
{
c[0] = a[0] * b[0];
c[1] = a[0] * b[0];
c[2] = a[0] * b[0] * a[1];
c[3] = a[0] * a[1];
c[4] = a[0] * b[0];
c[5] = b[0];
...
}
In C-type languages, the pointers a, b, and c may be aliased, so any write through c
could modify elements of a or b. This means that to guarantee functional correctness, the
compiler cannot load a[0] and b[0] into registers, multiply them, and store the result
to both c[0] and c[1], because the results would differ from the abstract execution
model if, say, a[0] is really the same location as c[0]. So the compiler cannot take
advantage of the common sub-expression. Likewise, the compiler cannot just reorder the
computation of c[4] into the proximity of the computation of c[0] and c[1] because
the preceding write to c[3] could change the inputs to the computation of c[4].
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 92
C Language Extensions
By making a, b, and c restricted pointers, the programmer asserts to the compiler that
the pointers are in fact not aliased, which in this case means writes through c would
never overwrite elements of a or b. This changes the function prototype as follows:
void foo(const float* restrict a,
const float* restrict b,
float* restrict c);
Note that all pointer arguments need to be made restricted for the compiler optimizer
to derive any benefit. With the restrict keywords added, the compiler can now
reorder and do common sub-expression elimination at will, while retaining functionality
identical with the abstract execution model:
void foo(const float* restrict a,
const float* restrict b,
float* restrict c)
{
float t0 = a[0];
float t1 = b[0];
float t2 = t0 * t2;
float t3 = a[1];
c[0] = t2;
c[1] = t2;
c[4] = t2;
c[2] = t2 * t3;
c[3] = t0 * t3;
c[5] = t1;
...
}
The effects here are a reduced number of memory accesses and reduced number of
computations. This is balanced by an increase in register pressure due to "cached" loads
and common sub-expressions.
Since register pressure is a critical issue in many CUDA codes, use of restricted pointers
can have negative performance impact on CUDA code, due to reduced occupancy.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 93
C Language Extensions
char1, uchar1 1
char2, uchar2 2
char3, uchar3 1
char4, uchar4 4
short1, ushort1 2
short2, ushort2 4
short3, ushort3 2
short4, ushort4 8
int1, uint1 4
int2, uint2 8
int3, uint3 4
int4, uint4 16
long4, ulong4 16
longlong1, ulonglong1 8
longlong2, ulonglong2 16
float1 4
float2 8
float3 4
float4 16
double1 8
double2 16
B.3.2. dim3
This type is an integer vector type based on uint3 that is used to specify dimensions.
When defining a variable of type dim3, any component left unspecified is initialized to 1.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 94
C Language Extensions
4. Built-in Variables
Built-in variables specify the grid and block dimensions and the block and thread
indices. They are only valid within functions that are executed on the device.
1. gridDim
This variable is of type dim3 (see dim3) and contains the dimensions of thegrid.
2. blockIdx
This variable is of type uint3 (see char, short, int, long, longlong, float, double) and
contains the block index within the grid.
3. blockDim
This variable is of type dim3 (see dim3) and contains the dimensions of theblock.
4. threadIdx
This variable is of type uint3 (see char, short, int, long, longlong, float, double ) and
contains the thread index within the block.
5. warpSize
This variable is of type int and contains the warp size in threads (see SIMT Architecture
for the definition of a warp).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 95
C Language Extensions
‣ A equal to 1 and B equal to 2 (thread 1's writes to X and Y happen after thread 2's
read of X and Y),
‣ A equal to 10 and B equal to 2 (thread 1's write to X happens before thread 2's read of
X and thread 1's write to Y happens after thread 2's read of Y),
‣ A equal to 10 and B equal to 20 (thread 1's writes to X and Y happen before thread 2's
read of X and Y),
Memory fence functions can be used to enforce some ordering.
void threadfence_block();
ensures that:
‣ All writes to shared and global memory made by the calling thread before the call
to threadfence_block() are observed by all threads in the block of the calling
thread as occurring before all writes to shared memory and global memory made by
the calling thread after the call to threadfence_block();
‣ All reads from shared memory and global memory made by the calling thread
before the call to threadfence_block() are performed before all reads from
shared memory and global memory made by the calling thread after the call to
threadfence_block().
void threadfence();
acts as threadfence_block() for all threads in the block of the calling thread and
also ensures that no writes to global memory made by the calling thread after the call
to threadfence() are observed by any thread in the device as occurring before any
write to global memory made by the calling thread before the call to threadfence().
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 96
C Language Extensions
Note that for this ordering guarantee to be true, the observing threads must truly
observe global memory and not cached versions of it; this is ensured by using the
volatile keyword as detailed in Volatile Qualifier.
void threadfence_system();
acts as threadfence_block() for all threads in the block of the calling thread and
also ensures that:
‣ All writes to global memory, page-locked host memory, and the memory of a peer
device made by the calling thread before the call to threadfence_system() are
observed by all threads in the device, host threads, and all threads in peer devices
as occurring before all writes to global memory, page-locked host memory, and the
memory of a peer device made by the calling thread after the call to
threadfence_system().
‣ All reads from shared memory, global memory, page-locked host memory,
and the memory of a peer device made by the calling thread before the call to
threadfence_system() are performed before all reads from shared memory,
global memory, page-locked host memory, and the memory of a peer device made
by the calling thread after the call to threadfence_system().
threadfence_system() is only supported by devices of compute capability 2.xand
higher.
In the previous code sample, inserting a fence function call between X = 10; and Y
= 20; and between int A = X; and int B = Y; would ensure that for thread 1, A
will always be equal to 10 if B is equal to 20. If thread 1 and 2 belong to the sameblock,
it is enough to use threadfence_block(). If thread 1 and 2 do not belong to the
same block, threadfence() must be used if they are CUDA threads from the same
device and threadfence_system() must be used if they are CUDA threads from
two different devices.
A common use case is when threads consume some data produced by other threads as
illustrated by the following code sample of a kernel that computes the sum of an array
of N numbers in one call. Each block first sums a subset of the array and stores the result
in global memory. When all blocks are done, the last block done reads each of these
partial sums from global memory and sums them to obtain the final result. In order to
determine which block is finished last, each block atomically increments a counter to
signal that it is done with computing and storing its partial sum (see Atomic Functions
about atomic functions). The last block is the one that receives the counter value equal
to gridDim.x-1. If no fence is placed between storing the partial sum and incrementing
the counter, the counter might increment before the partial sum is stored and therefore,
might reach gridDim.x-1 and let the last block start reading partial sums before they
have been actually updated in memory.
Memory fence functions only affect the ordering of memory operations by a thread;
they do not ensure that these memory operations are visible to other threads (like
syncthreads() does for threads within a block (see Synchronization Functions)).In
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 97
C Language Extensions
the code sample below, the visibility of memory operations on the result variableis
ensured by declaring it as volatile (see Volatile Qualifier).
if (threadIdx.x == 0) {
if (isLastBlockDone) {
if (threadIdx.x == 0) {
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 98
C Language Extensions
waits until all threads in the thread block have reached this point and all global and
shared memory accesses made by these threads prior to syncthreads() are visible
to all threads in the block.
syncthreads() is used to coordinate communication between the threads of the
same block. When some threads within a block access the same addresses in shared
or global memory, there are potential read-after-write, write-after-read, or write-after-
write hazards for some of these memory accesses. These data hazards can be avoided by
synchronizing threads in-between these accesses.
syncthreads() is allowed in conditional code but only if the conditional evaluates
identically across the entire thread block, otherwise the code execution is likely to hang
or produce unintended side effects.
Devices of compute capability 2.x and higher support three variations of
syncthreads() described below.
int syncthreads_count(int predicate);
is identical to syncthreads() with the additional feature that it evaluates predicate for
all threads of the block and returns non-zero if and only if predicate evaluates to non-
zero for all of them.
int syncthreads_or(int predicate);
is identical to syncthreads() with the additional feature that it evaluates predicate for
all threads of the block and returns non-zero if and only if predicate evaluates to non-
zero for any of them.
7. Mathematical Functions
The reference manual lists all C/C++ standard library mathematical functions that are
supported in device code and all intrinsic functions that are only supported in device
code.
Mathematical Functions provides accuracy information for some of these functions
when relevant.
8. Texture Functions
Texture objects are described in Texture Object API
Texture references are described in Texture Reference API
Texture fetching is described in Texture Fetching.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 99
C Language Extensions
fetches from the region of linear memory specified by the one-dimensional texture
object texObj using integer texture coordinate x. tex1Dfetch() only works with non-
normalized coordinates, so only the border and clamp addressing modes are supported.
It does not perform any texture filtering. For integer types, it may optionally promote
the integer to single-precision floating point.
B.8.1.2. tex1D()
template<class T>
T tex1D(cudaTextureObject_t texObj, float x);
fetches from the CUDA array specified by the one-dimensional texture object texObj
using texture coordinate x.
B.8.1.3. tex1DLod()
template<class T>
T tex1DLod(cudaTextureObject_t texObj, float x, float level);
fetches from the CUDA array specified by the one-dimensional texture object texObj
using texture coordinate x at the level-of-detail level.
B.8.1.4. tex1DGrad()
template<class T>
T tex1DGrad(cudaTextureObject_t texObj, float x, float dx, float dy);
fetches from the CUDA array specified by the one-dimensional texture object texObj
using texture coordinate x. The level-of-detail is derived from the X-gradient dx andY-
gradient dy.
B.8.1.5. tex2D()
template<class T>
T tex2D(cudaTextureObject_t texObj, float x, float y);
fetches from the CUDA array or the region of linear memory specified by the two-
dimensional texture object texObj using texture coordinate (x,y).
B.8.1.6. tex2DLod()
template<class T>
tex2DLod(cudaTextureObject_t texObj, float x, float y, float level);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 100
C Language Extensions
fetches from the CUDA array or the region of linear memory specified by the two-
dimensional texture object texObj using texture coordinate (x,y) at level-of-detail
level.
B.8.1.7. tex2DGrad()
template<class T>
T tex2DGrad(cudaTextureObject_t texObj, float x, float y,
float2 dx, float2 dy);
fetches from the CUDA array specified by the two-dimensional texture object texObj
using texture coordinate (x,y). The level-of-detail is derived from the dx and dy
gradients.
B.8.1.8. tex3D()
template<class T>
T tex3D(cudaTextureObject_t texObj, float x, float y, float z);
fetches from the CUDA array specified by the three-dimensional texture object texObj
using texture coordinate (x,y,z).
B.8.1.9. tex3DLod()
template<class T>
T tex3DLod(cudaTextureObject_t texObj, float x, float y, float z, float level);
fetches from the CUDA array or the region of linear memory specified by the three-
dimensional texture object texObj using texture coordinate (x,y,z) at level-of-detail
level.
B.8.1.10. tex3DGrad()
template<class T>
T tex3DGrad(cudaTextureObject_t texObj, float x, float y, float z,
float4 dx, float4 dy);
fetches from the CUDA array specified by the three-dimensional texture object texObj
using texture coordinate (x,y,z) at a level-of-detail derived from the X and Y gradients
dx and dy.
B.8.1.11. tex1DLayered()
template<class T>
T tex1DLayered(cudaTextureObject_t texObj, float x, int layer);
fetches from the CUDA array specified by the one-dimensional texture object texObj
using texture coordinate x and index layer, as described in Layered Textures
B.8.1.12. tex1DLayeredLod()
template<class T>
T tex1DLayeredLod(cudaTextureObject_t texObj, float x, int layer, float level);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 101
C Language Extensions
fetches from the CUDA array specified by the one-dimensional layered texture at layer
layer using texture coordinate x and level-of-detail level.
B.8.1.13. tex1DLayeredGrad()
template<class T>
T tex1DLayeredGrad(cudaTextureObject_t texObj, float x, int layer,
float dx, float dy);
fetches from the CUDA array specified by the one-dimensional layered texture at layer
layer using texture coordinate x and a level-of-detail derived from the dx and dy
gradients.
B.8.1.14. tex2DLayered()
template<class T>
T tex2DLayered(cudaTextureObject_t texObj,
float x, float y, int layer);
fetches from the CUDA array specified by the two-dimensional texture object texObj
using texture coordinate (x,y) and index layer, as described in Layered Textures.
B.8.1.15. tex2DLayeredLod()
template<class T>
T tex2DLayeredLod(cudaTextureObject_t texObj, float x, float y, int layer,
float level);
fetches from the CUDA array specified by the two-dimensional layered texture at layer
layer using texture coordinate (x,y).
B.8.1.16. tex2DLayeredGrad()
template<class T>
T tex2DLayeredGrad(cudaTextureObject_t texObj, float x, float y, int layer,
float2 dx, float2 dy);
fetches from the CUDA array specified by the two-dimensional layered texture at layer
layer using texture coordinate (x,y) and a level-of-detail derived from the dx and dy
X and Y gradients.
B.8.1.17. texCubemap()
template<class T>
T texCubemap(cudaTextureObject_t texObj, float x, float y, float z);
fetches the CUDA array specified by the three-dimensional texture object texObj using
texture coordinate (x,y,z), as described in Cubemap Textures.
B.8.1.18. texCubemapLod()
template<class T>
T texCubemapLod(cudaTextureObject_t texObj, float x, float, y, float z,
float level);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 102
C Language Extensions
fetches from the CUDA array specified by the three-dimensional texture object texObj
using texture coordinate (x,y,z) as described inCubemap Textures. The level-of-detail
used is given by level.
B.8.1.19. texCubemapLayered()
template<class T>
T texCubemapLayered(cudaTextureObject_t texObj,
float x, float y, float z, int layer);
fetches from the CUDA array specified by the cubemap layered texture object texObj
using texture coordinates (x,y,z), and index layer, as described in Cubemap Layered
Textures.
B.8.1.20. texCubemapLayeredLod()
template<class T>
T texCubemapLayeredLod(cudaTextureObject_t texObj, float x, float y, float z,
int layer, float level);
fetches from the CUDA array specified by the cubemap layered texture object texObj
using texture coordinate (x,y,z) and index layer, as described in Cubemap Layered
Textures, at level-of-detail level level.
B.8.1.21. tex2Dgather()
template<class T>
T tex2Dgather(cudaTextureObject_t texObj,
float x, float y, int comp = 0);
fetches from the CUDA array specified by the 2D texture object texObj using texture
coordinates x and y and the comp parameter as described in Texture Gather.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 103
C Language Extensions
float tex1Dfetch(
texture<unsigned char, cudaTextureType1D,
cudaReadModeNormalizedFloat> texRef,
int x);
float tex1Dfetch(
texture<signed char, cudaTextureType1D,
cudaReadModeNormalizedFloat> texRef,
int x);
float tex1Dfetch(
texture<unsigned short, cudaTextureType1D,
cudaReadModeNormalizedFloat> texRef,
int x);
float tex1Dfetch(
texture<signed short, cudaTextureType1D,
cudaReadModeNormalizedFloat> texRef,
int x);
fetches from the region of linear memory bound to the one-dimensional texture
reference texRef using integer texture coordinate x. tex1Dfetch() only works with
non-normalized coordinates, so only the border and clamp addressing modes are
supported. It does not perform any texture filtering. For integer types, it may optionally
promote the integer to single-precision floating point.
Besides the functions shown above, 2-, and 4-tuples are supported; for example:
float4 tex1Dfetch(
texture<uchar4, cudaTextureType1D,
cudaReadModeNormalizedFloat> texRef,
int x);
fetches from the region of linear memory bound to texture reference texRef using
texture coordinate x.
B.8.2.2. tex1D()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex1D(texture<DataType, cudaTextureType1D, readMode> texRef,
float x);
fetches from the CUDA array bound to the one-dimensional texture reference texRef
using texture coordinate x. Type is equal to DataType except when readMode is equal
to cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
equal to the matching floating-point type.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 104
C Language Extensions
B.8.2.3. tex1DLod()
template<class DataType, enum
cudaTextureReadMode readMode>
Type tex1DLod(texture<DataType, cudaTextureType1D, readMode> texRef, float x,
float level);
fetches from the CUDA array bound to the one-dimensional texture reference texRef
using texture coordinate x. The level-of-detail is given by level. Type is the same as
DataType except when readMode is cudaReadModeNormalizedFloat (see Texture
Reference API), in which case Type is the corresponding floating-pointtype.
B.8.2.4. tex1DGrad()
template<class DataType, enum
cudaTextureReadMode readMode>
Type tex1DGrad(texture<DataType, cudaTextureType1D, readMode> texRef, float x,
float dx, float dy);
fetches from the CUDA array bound to the one-dimensional texture reference
texRef using texture coordinate x. The level-of-detail is derived from the dx and
dy X- and Y-gradients. Type is the same as DataType except when readMode is
cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
the corresponding floating-point type.
B.8.2.5. tex2D()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex2D(texture<DataType, cudaTextureType2D, readMode> texRef,
float x, float y);
fetches from the CUDA array or the region of linear memory bound to the two-
dimensional texture reference texRef using texture coordinates x and y. Type is equal
to DataType except when readMode is equal to cudaReadModeNormalizedFloat (see
Texture Reference API), in which case Type is equal to the matching floating-pointtype.
B.8.2.6. tex2DLod()
template<class DataType, enum
cudaTextureReadMode readMode>
Type tex2DLod(texture<DataType, cudaTextureType2D, readMode> texRef,
float x, float y, float level);
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y). The level-of-detail is given by level. Type is the same
as DataType except when readMode is cudaReadModeNormalizedFloat (see Texture
Reference API), in which case Type is the corresponding floating-pointtype.
B.8.2.7. tex2DGrad()
template<class DataType, enum
cudaTextureReadMode readMode>
Type tex2DGrad(texture<DataType, cudaTextureType2D, readMode> texRef,
float x, float y, float2 dx, float2 dy);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 105
C Language Extensions
fetches from the CUDA array bound to the two-dimensional texture reference
texRef using texture coordinate (x,y). The level-of-detail is derived from the dx
and dy X- and Y-gradients. Type is the same as DataType except when readMode is
cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
the corresponding floating-point type.
B.8.2.8. tex3D()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex3D(texture<DataType, cudaTextureType3D, readMode> texRef,
float x, float y, float z);
fetches from the CUDA array bound to the three-dimensional texture reference texRef
using texture coordinates x, y, and z. Type is equal to DataType except when readMode
is equal to cudaReadModeNormalizedFloat (see Texture Reference API), in which case
Type is equal to the matching floating-point type.
B.8.2.9. tex3DLod()
template<class DataType, enum
cudaTextureReadMode readMode>
Type tex3DLod(texture<DataType, cudaTextureType3D, readMode> texRef,
float x, float y, float z, float level);
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y,z). The level-of-detail is given by level. Type is the
same as DataType except when readMode is cudaReadModeNormalizedFloat (see
Texture Reference API), in which case Type is the corresponding floating-pointtype.
B.8.2.10. tex3DGrad()
template<class DataType, enum
cudaTextureReadMode readMode>
Type tex3DGrad(texture<DataType, cudaTextureType3D, readMode> texRef,
float x, float y, float z, float4 dx, float4 dy);
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y,z). The level-of-detail is derived from the dx and
dy X- and Y-gradients. Type is the same as DataType except when readMode is
cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
the corresponding floating-point type.
B.8.2.11. tex1DLayered()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex1DLayered(
texture<DataType, cudaTextureType1DLayered, readMode> texRef,
float x, int layer);
fetches from the CUDA array bound to the one-dimensional layered texture
reference texRef using texture coordinate x and index layer, as described in
Layered Textures. Type is equal to DataType except when readMode is equal to
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 106
C Language Extensions
B.8.2.12. tex1DLayeredLod()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex1DLayeredLod(texture<DataType, cudaTextureType1D, readMode> texRef,
float x, int layer, float level);
fetches from the CUDA array bound to the one-dimensional texture reference texRef
using texture coordinate x and index layer as described in Layered Textures. The level-
of-detail is given by level. Type is the same as DataType except when readMode is
cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
the corresponding floating-point type.
B.8.2.13. tex1DLayeredGrad()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex1DLayeredGrad(texture<DataType, cudaTextureType1D, readMode> texRef,
float x, int layer, float dx, float dy);
fetches from the CUDA array bound to the one-dimensional texture reference texRef
using texture coordinate x and index layer as described in Layered Textures. The
level-of-detail is derived from the dx and dy X- and Y-gradients. Type is the same as
DataType except when readMode is cudaReadModeNormalizedFloat (see Texture
Reference API), in which case Type is the corresponding floating-pointtype.
B.8.2.14. tex2DLayered()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex2DLayered(
texture<DataType, cudaTextureType2DLayered, readMode> texRef,
float x, float y, int layer);
fetches from the CUDA array bound to the two-dimensional layered texture
reference texRef using texture coordinates x and y, and index layer, as described
in Texture Memory. Type is equal to DataType except when readMode is equal to
cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
equal to the matching floating-point type.
B.8.2.15. tex2DLayeredLod()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex2DLayeredLod(texture<DataType, cudaTextureType2D, readMode> texRef,
float x, float y, int layer, float level);
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y) and index layer as described in Layered Textures. The
level-of-detail is given by level. Type is the same as DataType except when readMode
is cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
the corresponding floating-point type.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 107
C Language Extensions
B.8.2.16. tex2DLayeredGrad()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex2DLayeredGrad(texture<DataType, cudaTextureType2D, readMode> texRef,
float x, float y, int layer, float2 dx, float2 dy);
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y) and index layer as described in Layered Textures. The
level-of-detail is derived from the dx and dy X- and Y-gradients. Type is the same as
DataType except when readMode is cudaReadModeNormalizedFloat (see Texture
Reference API), in which case Type is the corresponding floating-pointtype.
B.8.2.17. texCubemap()
template<class DataType, enum cudaTextureReadMode readMode>
Type texCubemap(
texture<DataType, cudaTextureTypeCubemap, readMode> texRef,
float x, float y, float z);
fetches from the CUDA array bound to the cubemap texture reference texRef using
texture coordinates x, y, and z, as described in Cubemap Textures. Type is equal to
DataType except when readMode is equal to cudaReadModeNormalizedFloat (see
Texture Reference API), in which case Type is equal to the matching floating-pointtype.
B.8.2.18. texCubemapLod()
template<class DataType, enum cudaTextureReadMode readMode>
Type texCubemapLod(texture<DataType, cudaTextureType3D, readMode> texRef,
float x, float y, float z, float level);
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y,z). The level-of-detail is given by level. Type is the
same as DataType except when readMode is cudaReadModeNormalizedFloat (see
Texture Reference API), in which case Type is the corresponding floating-pointtype.
B.8.2.19. texCubemapLayered()
template<class DataType, enum cudaTextureReadMode readMode>
Type texCubemapLayered(
texture<DataType, cudaTextureTypeCubemapLayered, readMode> texRef,
float x, float y, float z, int layer);
fetches from the CUDA array bound to the cubemap layered texture reference texRef
using texture coordinates x, y, and z, and index layer, as described in Cubemap
Layered Textures. Type is equal to DataType except when readMode is equal to
cudaReadModeNormalizedFloat (see Texture Reference API), in which case Type is
equal to the matching floating-point type.
B.8.2.20. texCubemapLayeredLod()
template<class DataType, enum cudaTextureReadMode readMode>
Type texCubemapLayeredLod(texture<DataType, cudaTextureType3D, readMode> texRef,
float x, float y, float z, int layer, float level);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 108
C Language Extensions
fetches from the CUDA array bound to the two-dimensional texture reference texRef
using texture coordinate (x,y,z) and index layer as described in Layered Textures.
The level-of-detail is given by level. Type is the same as DataType except when
readMode is cudaReadModeNormalizedFloat (see Texture Reference API), in which
case Type is the corresponding floating-point type.
B.8.2.21. tex2Dgather()
template<class DataType, enum cudaTextureReadMode readMode>
Type tex2Dgather(
texture<DataType, cudaTextureType2D, readMode> texRef,
float x, float y, int comp = 0);
fetches from the CUDA array bound to the 2D texture reference texRef using texture
coordinates x and y and the comp parameter as described in Texture Gather. Type is a 4-
component vector type. It is based on the base type of DataType except when readMode
is equal to cudaReadModeNormalizedFloat (see Texture Reference API), in which case
it is always float4.
9. Surface Functions
Surface functions are only supported by devices of compute capability 2.0 and higher.
Surface objects are described in described in Surface Object API
reads the CUDA array specified by the one-dimensional surface object surfObj using
coordinate x.
B.9.1.2. surf1Dwrite
template<class T>
void surf1Dwrite(T data,
cudaSurfaceObject_t surfObj,
int x,
boundaryMode = cudaBoundaryModeTrap);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 109
C Language Extensions
writes value data to the CUDA array specified by the one-dimensional surface object
surfObj at coordinate x.
B.9.1.3. surf2Dread()
template<class T>
T surf2Dread(cudaSurfaceObject_t surfObj,
int x, int y,
boundaryMode = cudaBoundaryModeTrap);
template<class T>
void surf2Dread(T* data,
cudaSurfaceObject_t surfObj,
int x, int y,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array specified by the two-dimensional surface object surfObj using
coordinates x and y.
B.9.1.4. surf2Dwrite()
template<class T>
void surf2Dwrite(T data,
cudaSurfaceObject_t surfObj,
int x, int y,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array specified by the two-dimensional surface object
surfObj at coordinate x and y.
B.9.1.5. surf3Dread()
template<class T>
T surf3Dread(cudaSurfaceObject_t surfObj,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
template<class T>
void surf3Dread(T* data,
cudaSurfaceObject_t surfObj,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array specified by the three-dimensional surface object surfObj using
coordinates x, y, and z.
B.9.1.6. surf3Dwrite()
template<class T>
void surf3Dwrite(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array specified by the three-dimensional object surfObj
at coordinate x, y, and z.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 110
C Language Extensions
B.9.1.7. surf1DLayeredread()
template<class T>
T surf1DLayeredread(
cudaSurfaceObject_t surfObj,
int x, int layer,
boundaryMode = cudaBoundaryModeTrap);
template<class T>
void surf1DLayeredread(T data,
cudaSurfaceObject_t surfObj,
int x, int layer,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array specified by the one-dimensional layered surface object surfObj
using coordinate x and index layer.
B.9.1.8. surf1DLayeredwrite()
template<class Type>
void surf1DLayeredwrite(T data,
cudaSurfaceObject_t surfObj,
int x, int layer,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array specified by the two-dimensional layered surface
object surfObj at coordinate x and index layer.
B.9.1.9. surf2DLayeredread()
template<class T>
T surf2DLayeredread(
cudaSurfaceObject_t surfObj,
int x, int y, int layer,
boundaryMode = cudaBoundaryModeTrap);
template<class T>
void surf2DLayeredread(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int layer,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array specified by the two-dimensional layered surface object surfObj
using coordinate x and y, and index layer.
B.9.1.10. surf2DLayeredwrite()
template<class T>
void surf2DLayeredwrite(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int layer,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array specified by the one-dimensional layered surface
object surfObj at coordinate x and y, and index layer.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 111
C Language Extensions
B.9.1.11. surfCubemapread()
template<class T>
T surfCubemapread(
cudaSurfaceObject_t surfObj,
int x, int y, int face,
boundaryMode = cudaBoundaryModeTrap);
template<class T>
void surfCubemapread(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int face,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array specified by the cubemap surface object surfObj using
coordinate x and y, and face index face.
B.9.1.12. surfCubemapwrite()
template<class T>
void surfCubemapwrite(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int face,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array specified by the cubemap object surfObj at
coordinate x and y, and face index face.
B.9.1.13. surfCubemapLayeredread()
template<class T>
T surfCubemapLayeredread(
cudaSurfaceObject_t surfObj,
int x, int y, int layerFace,
boundaryMode = cudaBoundaryModeTrap);
template<class T>
void surfCubemapLayeredread(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int layerFace,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array specified by the cubemap layered surface object surfObj using
coordinate x and y, and index layerFace.
B.9.1.14. surfCubemapLayeredwrite()
template<class T>
void surfCubemapLayeredwrite(T data,
cudaSurfaceObject_t surfObj,
int x, int y, int layerFace,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array specified by the cubemap layered object surfObj at
coordinate x and y, and index layerFace.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 112
C Language Extensions
B.9.2.2. surf1Dwrite
template<class Type>
void surf1Dwrite(Type data,
surface<void, cudaSurfaceType1D> surfRef,
int x,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the one-dimensional surfacereference
surfRef at coordinate x.
B.9.2.3. surf2Dread()
template<class Type>
Type surf2Dread(surface<void, cudaSurfaceType2D> surfRef,
int x, int y,
boundaryMode = cudaBoundaryModeTrap);
template<class Type>
void surf2Dread(Type* data,
surface<void, cudaSurfaceType2D> surfRef,
int x, int y,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array bound to the two-dimensional surface reference surfRef using
coordinates x and y.
B.9.2.4. surf2Dwrite()
template<class Type>
void surf3Dwrite(Type data,
surface<void, cudaSurfaceType3D> surfRef,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the two-dimensional surfacereference
surfRef at coordinate x and y.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 113
C Language Extensions
B.9.2.5. surf3Dread()
template<class Type>
Type surf3Dread(surface<void, cudaSurfaceType3D> surfRef,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
template<class Type>
void surf3Dread(Type* data,
surface<void, cudaSurfaceType3D> surfRef,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array bound to the three-dimensional surface reference surfRef using
coordinates x, y, and z.
B.9.2.6. surf3Dwrite()
template<class Type>
void surf3Dwrite(Type data,
surface<void, cudaSurfaceType3D> surfRef,
int x, int y, int z,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the three-dimensional surface reference
surfRef at coordinate x, y, and z.
B.9.2.7. surf1DLayeredread()
template<class Type>
Type surf1DLayeredread(
surface<void, cudaSurfaceType1DLayered> surfRef,
int x, int layer,
boundaryMode = cudaBoundaryModeTrap);
template<class Type>
void surf1DLayeredread(Type data,
surface<void, cudaSurfaceType1DLayered> surfRef,
int x, int layer,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array bound to the one-dimensional layered surface reference surfRef
using coordinate x and index layer.
B.9.2.8. surf1DLayeredwrite()
template<class Type>
void surf1DLayeredwrite(Type data,
surface<void, cudaSurfaceType1DLayered> surfRef,
int x, int layer,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the two-dimensional layered surface
reference surfRef at coordinate x and index layer.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 114
C Language Extensions
B.9.2.9. surf2DLayeredread()
template<class Type>
Type surf2DLayeredread(
surface<void, cudaSurfaceType2DLayered> surfRef,
int x, int y, int layer,
boundaryMode = cudaBoundaryModeTrap);
template<class Type>
void surf2DLayeredread(Type data,
surface<void, cudaSurfaceType2DLayered> surfRef,
int x, int y, int layer,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array bound to the two-dimensional layered surface reference surfRef
using coordinate x and y, and index layer.
B.9.2.10. surf2DLayeredwrite()
template<class Type>
void surf2DLayeredwrite(Type data,
surface<void, cudaSurfaceType2DLayered> surfRef,
int x, int y, int layer,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the one-dimensional layered surface
reference surfRef at coordinate x and y, and index layer.
B.9.2.11. surfCubemapread()
template<class Type>
Type surfCubemapread(
surface<void, cudaSurfaceTypeCubemap> surfRef,
int x, int y, int face,
boundaryMode = cudaBoundaryModeTrap);
template<class Type>
void surfCubemapread(Type data,
surface<void, cudaSurfaceTypeCubemap> surfRef,
int x, int y, int face,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array bound to the cubemap surface reference surfRef using
coordinate x and y, and face index face.
B.9.2.12. surfCubemapwrite()
template<class Type>
void surfCubemapwrite(Type data,
surface<void, cudaSurfaceTypeCubemap> surfRef,
int x, int y, int face,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the cubemap reference surfRef at
coordinate x and y, and face index face.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 115
C Language Extensions
B.9.2.13. surfCubemapLayeredread()
template<class Type>
Type surfCubemapLayeredread(
surface<void, cudaSurfaceTypeCubemapLayered> surfRef,
int x, int y, int layerFace,
boundaryMode = cudaBoundaryModeTrap);
template<class Type>
void surfCubemapLayeredread(Type data,
surface<void, cudaSurfaceTypeCubemapLayered> surfRef,
int x, int y, int layerFace,
boundaryMode = cudaBoundaryModeTrap);
reads the CUDA array bound to the cubemap layered surface reference surfRef using
coordinate x and y, and index layerFace.
B.9.2.14. surfCubemapLayeredwrite()
template<class Type>
void surfCubemapLayeredwrite(Type data,
surface<void, cudaSurfaceTypeCubemapLayered> surfRef,
int x, int y, int layerFace,
boundaryMode = cudaBoundaryModeTrap);
writes value data to the CUDA array bound to the cubemap layered reference surfRef
at coordinate x and y, and index layerFace.
returns the data of type T located at address address, where T is char, short, int,
long long unsigned char, unsigned short, unsigned int, unsigned long
long, int2, int4, uint2, uint4, float, float2, float4, double, or double2. The
operation is cached in the read-only data cache (see Global Memory).
when executed in device code, returns the value of a per-multiprocessor counter that is
incremented every clock cycle. Sampling this counter at the beginning and at the end of
a kernel, taking the difference of the two samples, and recording the result per thread
provides a measure for each thread of the number of clock cycles taken by the device to
completely execute the thread, but not of the number of clock cycles the device actually
spent executing thread instructions. The former number is greater than the latter since
threads are time sliced.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 116
C Language Extensions
void foo() {
int *addr;
cudaMallocManaged(&addr, 4);
*addr = 0;
mykernel<<<...>>>(addr);
sync_fetch_and_add(addr, 10); // CPU atomic operation
}
The new scoped versions of atomics are available for all atomics listed below only for
compute capabilities 6.x.
Note that any atomic operation can be implemented based on atomicCAS() (Compare
And Swap). For example, atomicAdd() for double-precision floating-point numbers
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 117
C Language Extensions
is not available on devices with compute capability lower than 6.0 but it can be
implemented as follows:
#if CUDA_ARCH < 600
device double atomicAdd(double* address, double val)
{
unsigned long long int* address_as_ull =
(unsigned long long int*)address;
unsigned long long int old = *address_as_ull, assumed;
do {
assumed = old;
old = atomicCAS(address_as_ull, assumed,
double_as_longlong(val +
longlong_as_double(assumed)));
// Note: uses integer comparison to avoid hang in case of NaN (since NaN !=
NaN)
} while (assumed != old);
return longlong_as_double(old);
}
#endif
1. Arithmetic Functions
1. atomicAdd()
int atomicAdd(int* address, int val);
unsigned int atomicAdd(unsigned int* address,
unsigned int val);
unsigned long long int atomicAdd(unsigned long long int* address,
unsigned long long int val);
float atomicAdd(float* address, float val);
double atomicAdd(double* address, double val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes (old + val), and stores the result back to memory at the same
address. These three operations are performed in one atomic transaction. The function
returns old.
The 32-bit floating-point version of atomicAdd() is only supported by devicesof
compute capability 2.x and higher.
The 64-bit floating-point version of atomicAdd() is only supported by devicesof
compute capability 6.x and higher.
B.12.1.2. atomicSub()
int atomicSub(int* address, int val);
unsigned int atomicSub(unsigned int* address,
unsigned int val);
reads the 32-bit word old located at the address address in global or shared memory,
computes (old - val), and stores the result back to memory at the same address.
These three operations are performed in one atomic transaction. The function returns
old.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 118
C Language Extensions
B.12.1.3. atomicExch()
int atomicExch(int* address, int val);
unsigned int atomicExch(unsigned int* address,
unsigned int val);
unsigned long long int atomicExch(unsigned long long int* address,
unsigned long long int val);
float atomicExch(float* address, float val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory and stores val back to memory at the same address. These two operationsare
performed in one atomic transaction. The function returns old.
B.12.1.4. atomicMin()
int atomicMin(int* address, int val);
unsigned int atomicMin(unsigned int* address,
unsigned int val);
unsigned long long int atomicMin(unsigned long long int* address,
unsigned long long int val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes the minimum of old and val, and stores the result back to memory
at the same address. These three operations are performed in one atomic transaction.
The function returns old.
The 64-bit version of atomicMin() is only supported by devices of computecapability
3.5 and higher.
B.12.1.5. atomicMax()
int atomicMax(int* address, int val);
unsigned int atomicMax(unsigned int* address,
unsigned int val);
unsigned long long int atomicMax(unsigned long long int* address,
unsigned long long int val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes the maximum of old and val, and stores the result back to memory
at the same address. These three operations are performed in one atomic transaction.
The function returns old.
The 64-bit version of atomicMax() is only supported by devices of computecapability
3.5 and higher.
B.12.1.6. atomicInc()
unsigned int atomicInc(unsigned int* address,
unsigned int val);
reads the 32-bit word old located at the address address in global or shared memory,
computes ((old >= val) ? 0 : (old+1)), and stores the result back to memory at
the same address. These three operations are performed in one atomic transaction. The
function returns old.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 119
C Language Extensions
B.12.1.7. atomicDec()
unsigned int atomicDec(unsigned int* address,
unsigned int val);
reads the 32-bit word old located at the address address in global or shared memory,
computes (((old == 0) | (old > val)) ? val : (old-1) ), and stores the
result back to memory at the same address. These three operations are performed in one
atomic transaction. The function returns old.
B.12.1.8. atomicCAS()
int atomicCAS(int* address, int compare, int val);
unsigned int atomicCAS(unsigned int* address,
unsigned int compare,
unsigned int val);
unsigned long long int atomicCAS(unsigned long long int* address,
unsigned long long int compare,
unsigned long long int val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes (old == compare ? val : old) , and stores the result back to
memory at the same address. These three operations are performed in one atomic
transaction. The function returns old (Compare And Swap).
2. Bitwise Functions
1. atomicAnd()
int atomicAnd(int* address, int val);
unsigned int atomicAnd(unsigned int* address,
unsigned int val);
unsigned long long int atomicAnd(unsigned long long int* address,
unsigned long long int val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes (old & val), and stores the result back to memory at the same
address. These three operations are performed in one atomic transaction. The function
returns old.
The 64-bit version of atomicAnd() is only supported by devices of computecapability
3.5 and higher.
B.12.2.2. atomicOr()
int atomicOr(int* address, int val);
unsigned int atomicOr(unsigned int* address,
unsigned int val);
unsigned long long int atomicOr(unsigned long long int* address,
unsigned long long int val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes (old | val), and stores the result back to memory at the same
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 120
C Language Extensions
address. These three operations are performed in one atomic transaction. The function
returns old.
The 64-bit version of atomicOr() is only supported by devices of computecapability
3.5 and higher.
B.12.2.3. atomicXor()
int atomicXor(int* address, int val);
unsigned int atomicXor(unsigned int* address,
unsigned int val);
unsigned long long int atomicXor(unsigned long long int* address,
unsigned long long int val);
reads the 32-bit or 64-bit word old located at the address address in global or shared
memory, computes (old ^ val), and stores the result back to memory at the same
address. These three operations are performed in one atomic transaction. The function
returns old.
The 64-bit version of atomicXor() is only supported by devices of computecapability
3.5 and higher.
The warp vote functions allow the threads of a given warp to perform a reduction-and-
broadcast operation. These functions take as input an integer predicate from each
thread in the warp and compare those values with zero. The results of the comparisons
are combined (reduced) across the active threads of the warp in one of the following
ways, broadcasting a single return value to each participating thread:
all(predicate):
Evaluate predicate for all active threads of the warp and return non-zero if and
only if predicate evaluates to non-zero for all of them.
any(predicate):
Evaluate predicate for all active threads of the warp and return non-zero if and
only if predicate evaluates to non-zero for any of them.
ballot(predicate):
Evaluate predicate for all active threads of the warp and return an integer whose
Nth bit is set if and only if predicate evaluates to non-zero for the Nth thread of the
warp and the Nth thread is active.
Notes
For each of these warp vote operations, the result excludes threads that are inactive (e.g.,
due to warp divergence). Inactive threads are represented by 0 bits in the value returned
by ballot() and are not considered in the reductions performed by all() and
any().
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 121
C Language Extensions
1. Synopsis
B.14.2. Description
The shfl() intrinsics permit exchanging of a variable between threads within a
warp without use of shared memory. The exchange occurs simultaneously for all
active threads within the warp, moving 4 bytes of data per thread. Exchange of 8-byte
quantities must be broken into two separate invocations of shfl().
Threads within a warp are referred to as lanes, and may have an index between 0 and
warpSize-1 (inclusive). Four source-lane addressing modes are supported:
shfl()
Direct copy from indexed lane
shfl_up()
Copy from a lane with lower ID relative to caller
shfl_down()
Copy from a lane with higher ID relative to caller
shfl_xor()
Copy from a lane based on bitwise XOR of own lane ID
Threads may only read data from another thread which is actively participating in the
shfl() command. If the target thread is inactive, the retrieved value is undefined.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 122
C Language Extensions
All of the shfl() intrinsics take an optional width parameter which alters the
behavior of the intrinsic. width must have a value which is a power of 2; resultsare
undefined if width is not a power of 2, or is a number greater than warpSize.
shfl() returns the value of var held by the thread whose ID is given by srcLane.
If width is less than warpSize then each subsection of the warp behaves as a separate
entity with a starting logical lane ID of 0. If srcLane is outside the range [0:width-1],
the value returned corresponds to the value of var held by the srcLane modulo width
(i.e. within the same subsection).
shfl_up() calculates a source lane ID by subtracting delta from the caller's lane ID.
The value of var held by the resulting lane ID is returned: in effect, var is shifted up the
warp by delta lanes. If width is less than warpSize then each subsection of the warp
behaves as a separate entity with a starting logical lane ID of 0. The source lane index
will not wrap around the value of width, so effectively the lower delta lanes will be
unchanged.
shfl_down() calculates a source lane ID by adding delta to the caller's lane ID. The
value of var held by the resulting lane ID is returned: this has the effect of shifting var
down the warp by delta lanes. If width is less than warpSize then each subsection
of the warp behaves as a separate entity with a starting logical lane ID of 0. As for
shfl_up(), the ID number of the source lane will not wrap around the value of
width and so the upper delta lanes will remain unchanged.
shfl_xor() calculates a source line ID by performing a bitwise XOR of thecaller's
lane ID with laneMask: the value of var held by the resulting lane ID is returned. If
width is less than warpSize then each group of width consecutive threads are able
to access elements from earlier groups of threads, however if they attempt to access
elements from later groups of threads their own value of var will be returned. This
mode implements a butterfly addressing pattern such as is used in tree reduction and
broadcast.
3. Return Value
All shfl() intrinsics return the 4-byte word referenced by var from the source lane
ID as an unsigned integer. If the source lane ID is out of range or the source thread has
exited, the calling thread's own var is returned.
4. Notes
All shfl() intrinsics share the same semantics with respect to code motion as the
vote intrinsics any() and all().
Threads may only read data from another thread which is actively participating in the
shfl() command. If the target thread is inactive, the retrieved value is undefined.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 123
C Language Extensions
width must be a power-of-2 (i.e., 2, 4, 8, 16 or 32). Results are unspecified for other
values.
Types other than int or float must first be cast in order to use the shfl() intrinsics.
5. Examples
1. Broadcast of a single value across a warp
#include <stdio.h>
int main() {
bcast<<< 1, 32 >>>(1234);
cudaDeviceSynchronize();
return 0;
}
int main() {
scan4<<< 1, 32 >>>();
cudaDeviceSynchronize();
return 0;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 124
C Language Extensions
int main() {
warpReduce<<< 1, 32 >>>();
cudaDeviceSynchronize();
return 0;
}
increments by one per warp the per-multiprocessor hardware counter of index counter.
Counters 8 to 15 are reserved and should not be used by applications.
The value of counters 0, 1, ..., 7 can be obtained via nvprof by nvprof --events
prof_trigger_0x where x is 0, 1, ..., 7. All counters are reset before each kernel launch
(note that when collecting counters, kernel launches are synchronous as mentioned in
Concurrent Execution between Host and Device).
B.16. Assertion
Assertion is only supported by devices of compute capability 2.x and higher. It is not
supported on MacOS, regardless of the device, and loading a module that references the
assert function on Mac OS will fail.
void assert(int expression);
stops the kernel execution if expression is equal to zero. If the program is run within a
debugger, this triggers a breakpoint and the debugger can be used to inspect the current
state of the device. Otherwise, each thread for which expression is equal to zero prints
a message to stderr after synchronization with the host via cudaDeviceSynchronize(),
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 125
C Language Extensions
will output:
test.cu:19: void testAssert(): block: [0,0,0], thread: [0,0,0] Assertion
`should_be_one` failed.
Assertions are for debugging purposes. They can affect performance and it is therefore
recommended to disable them in production code. They can be disabled at compile
time by defining the NDEBUG preprocessor macro before including assert.h. Note that
expression should not be an expression with side effects (something like (++i > 0),
for example), otherwise disabling the assertion will affect the functionality of the code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 126
C Language Extensions
The following fields are supported (see widely-available documentation for a complete
description of all behaviors):
B.17.2. Limitations
Final formatting of the printf() output takes place on the host system. This means
that the format string must be understood by the host-system's compiler and C library.
Every effort has been made to ensure that the format specifiers supported by CUDA's
printf function form a universal subset from the most common host compilers, but exact
behavior will be host-OS-dependent.
As described in Format Specifiers, printf() will accept all combinations of valid flags
and types. This is because it cannot determine what will and will not be valid on the
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 127
C Language Extensions
host system where the final output is formatted. The effect of this is that output may be
undefined if the program emits a format string which contains invalid combinations.
The printf() command can accept at most 32 arguments in addition to the format
string. Additional arguments beyond this will be ignored, and the format specifier
output as-is.
Owing to the differing size of the long type on 64-bit Windows platforms (four bytes
on 64-bit Windows platforms, eight bytes on other 64-bit platforms), a kernel which is
compiled on a non-Windows 64-bit machine but then run on a win64 machine will see
corrupted output for all format strings which include "%ld". It is recommended that the
compilation platform matches the execution platform to ensure safety.
The output buffer for printf() is set to a fixed size before kernel launch (see
Associated Host-Side API). It is circular and if more output is produced during kernel
execution than can fit in the buffer, older output is overwritten. It is flushed only when
one of these actions is performed:
‣ Kernel launch via <<<>>> or cuLaunchKernel() (at the start of the launch, and if
the CUDA_LAUNCH_BLOCKING environment variable is set to 1, at the end of the
launch as well),
‣ Synchronization via cudaDeviceSynchronize(), cuCtxSynchronize(),
cudaStreamSynchronize(), cuStreamSynchronize(),
cudaEventSynchronize(), or cuEventSynchronize(),
‣ Memory copies via any blocking version of cudaMemcpy*() or cuMemcpy*(),
‣ Module loading/unloading via cuModuleLoad() or cuModuleUnload(),
‣ Context destruction via cudaDeviceReset() or cuCtxDestroy().
‣ Prior to executing a stream callback added by cudaStreamAddCallback or
cuStreamAddCallback.
Note that the buffer is not flushed automatically when the program exits. The user must
call cudaDeviceReset() or cuCtxDestroy() explicitly, as shown in the examples
below.
Internally printf() uses a shared data structure and so it is possible that calling
printf() might change the order of execution of threads. In particular, a thread
which calls printf() might take a longer execution path than one which does notcall
printf(), and that path length is dependent upon the parameters of the printf().
Note, however, that CUDA makes no guarantees of thread execution order except at
explicit syncthreads() barriers, so it is impossible to tell whether execution order
has been modified by printf() or by other scheduling behaviour in the hardware.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 128
C Language Extensions
‣ cudaDeviceGetLimit(size_t* size,cudaLimitPrintfFifoSize)
‣ cudaDeviceSetLimit(cudaLimitPrintfFifoSize, size_t size)
B.17.4. Examples
The following code sample:
#include <stdio.h>
int main()
{
helloCUDA<<<1, 5>>>(1.2345f);
cudaDeviceSynchronize();
return 0;
}
will output:
Hello thread 2, f=1.2345
Hello thread 1, f=1.2345
Hello thread 4, f=1.2345
Hello thread 0, f=1.2345
Hello thread 3, f=1.2345
Notice how each thread encounters the printf() command, so there are as manylines
of output as there were threads launched in the grid. As expected, global values (i.e.,
float f) are common between all threads, and local values (i.e., threadIdx.x) are
distinct per-thread.
The following code sample:
#include <stdio.h>
int main()
{
helloCUDA<<<1, 5>>>(1.2345f);
cudaDeviceSynchronize();
return 0;
}
will output:
Self-evidently, the if() statement limits which threads will call printf, so that onlya
single line of output is seen.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 129
C Language Extensions
allocate and free memory dynamically from a fixed-size heap in global memory.
void* memcpy(void* dest, const void* src, size_t size);
copy size bytes from the memory location pointed by src to the memory location
pointed by dest.
void* memset(void* ptr, int value, size_t size);
set size bytes of memory block pointed by ptr to value (interpreted as an unsigned
char).
The CUDA in-kernel malloc() function allocates at least size bytes from the device
heap and returns a pointer to the allocated memory or NULL if insufficient memory
exists to fulfill the request. The returned pointer is guaranteed to be aligned to a 16-byte
boundary.
The CUDA in-kernel free() function deallocates the memory pointed to by ptr,which
must have been returned by a previous call to malloc(). If ptr is NULL, the call to
free() is ignored. Repeated calls to free() with the same ptr has undefined behavior.
The memory allocated by a given CUDA thread via malloc() remains allocated for the
lifetime of the CUDA context, or until it is explicitly released by a call to free(). It can
be used by any other CUDA threads even from subsequent kernel launches. Any CUDA
thread may free memory allocated by another thread, but care should be taken to ensure
that the same pointer is not freed more than once.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 130
C Language Extensions
The actual memory allocation for the heap occurs when a module is loaded into the
context, either explicitly via the CUDA driver API (see Module), or implicitly via the
CUDA runtime API (see CUDA C Runtime). If the memory allocation fails, the module
load will generate a CUDA_ERROR_SHARED_OBJECT_INIT_FAILED error.
Heap size cannot be changed once a module load has occurred and it does not resize
dynamically according to need.
Memory reserved for the device heap is in addition to memory allocated through host-
side CUDA API calls such as cudaMalloc().
3. Examples
int main()
{
// Set a heap size of 128 megabytes. Note that this must
// be done before any kernel is launched.
cudaDeviceSetLimit(cudaLimitMallocHeapSize, 128*1024*1024);
mallocTest<<<1, 5>>>();
cudaDeviceSynchronize();
return 0;
}
will output:
Thread 0 got pointer: 00057020
Thread 1 got pointer: 0005708c
Thread 2 got pointer: 000570f8
Thread 3 got pointer: 00057164
Thread 4 got pointer: 000571d0
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 131
C Language Extensions
Notice how each thread encounters the malloc() and memset() commands and so
receives and initializes its own allocation. (Exact pointer values will vary: these are
illustrative.)
// The first thread in the block does the allocation and initialization
// and then shares the pointer with all other threads through shared memory,
// so that access can easily be coalesced.
// 64 bytes per thread are allocated.
if (threadIdx.x == 0) {
size_t size = blockDim.x * 64;
data = (int*)malloc(size);
memset(data, 0, size);
}
syncthreads();
int main()
{
cudaDeviceSetLimit(cudaLimitMallocHeapSize, 128*1024*1024);
mallocTest<<<10, 128>>>();
cudaDeviceSynchronize();
return 0;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 132
C Language Extensions
#define NUM_BLOCKS 20
}
Simple example: store thread ID into each element
//global void usemem()
{
int* ptr = dataptr[blockIdx.x];
if (ptr != NULL)
ptr[threadIdx.x] += threadIdx.x;
}
int main()
{
cudaDeviceSetLimit(cudaLimitMallocHeapSize, 128*1024*1024);
// Allocate memory
allocmem<<< NUM_BLOCKS, 10 >>>();
// Use memory
usemem<<< NUM_BLOCKS, 10 >>>();
usemem<<< NUM_BLOCKS, 10 >>>();
usemem<<< NUM_BLOCKS, 10 >>>();
// Free memory
freemem<<< NUM_BLOCKS, 10 >>>();
cudaDeviceSynchronize();
return 0;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 133
C Language Extensions
‣ Dg is of type dim3 (see dim3) and specifies the dimension and size of the grid, such
that Dg.x * Dg.y * Dg.z equals the number of blocks being launched;
‣ Db is of type dim3 (see dim3) and specifies the dimension and size of each block,
such that Db.x * Db.y * Db.z equals the number of threads per block;
‣ Ns is of type size_t and specifies the number of bytes in shared memory that is
dynamically allocated per block for this call in addition to the statically allocated
memory; this dynamically allocated memory is used by any of the variables
declared as an external array as mentioned in shared ; Ns is an optional
argument which defaults to 0;
‣ S is of type cudaStream_t and specifies the associated stream; S is an optional
argument which defaults to 0.
The arguments to the execution configuration are evaluated before the actual function
arguments.
The function call will fail if Dg or Db are greater than the maximum sizes allowed for
the device as specified in Compute Capabilities, or if Ns is greater than the maximum
amount of shared memory available on the device, minus the amount of shared memory
required for static allocation.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 134
C Language Extensions
‣ If the initial register usage is higher than L, the compiler reduces it further until it
becomes less or equal to L, usually at the expense of more local memory usage and/
or higher number of instructions;
‣ If the initial register usage is lower than L
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 135
C Language Extensions
Optimal launch bounds for a given kernel will usually differ across major architecture
revisions. The sample code below shows how this is typically handled in device code
using the CUDA_ARCH macro introduced in Application Compatibility
#define THREADS_PER_BLOCK 256
#if CUDA_ARCH >= 200
#define MY_KERNEL_MAX_THREADS (2 * THREADS_PER_BLOCK)
#define MY_KERNEL_MIN_BLOCKS 3
#else
#define MY_KERNEL_MAX_THREADS THREADS_PER_BLOCK
#define MY_KERNEL_MIN_BLOCKS 2
#endif
// Device code
global void
launch_bounds (MY_KERNEL_MAX_THREADS, MY_KERNEL_MIN_BLOCKS)
MyKernel(...)
{
...
}
In the common case where MyKernel is invoked with the maximum number ofthreads
per block (specified as the first parameter of launch_bounds ()), it is tempting
to use MY_KERNEL_MAX_THREADS as the number of threads per block in the execution
configuration:
// Host code
MyKernel<<<blocksPerGrid, MY_KERNEL_MAX_THREADS>>>(...);
This will not work however since CUDA_ARCH is undefined in host code as
mentioned in Application Compatibility, so MyKernel will launch with 256 threads
per block even when CUDA_ARCH is greater or equal to 200. Instead the number of
threads per block should be determined:
‣ Either at compile time using a macro that does not depend on CUDA_ARCH , for
example
// Host code
MyKernel<<<blocksPerGrid, THREADS_PER_BLOCK>>>(...);
‣ Or at runtime based on the compute capability
// Host code
cudaGetDeviceProperties(&deviceProp, device);
int threadsPerBlock =
(deviceProp.major >= 2 ?
2 * THREADS_PER_BLOCK : THREADS_PER_BLOCK);
MyKernel<<<blocksPerGrid, threadsPerBlock>>>(...);
Register usage is reported by the --ptxas options=-v compiler option. Thenumber
of resident blocks can be derived from the occupancy reported by the CUDA profiler
(see Device Memory Accessesfor a definition of occupancy).
Register usage can also be controlled for all global functions in a file using the
maxrregcount compiler option. The value of maxrregcount is ignored for functions
with launch bounds.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 136
C Language Extensions
Examples:
struct S1_t { static const int value = 4; };
template <int X, typename T2>
device void foo(int *p1, int *p2) {
// no argument specified, loop will be completely unrolled
#pragma unroll
for (int i = 0; i < 12; ++i)
p1[i] += p2[i]*2;
// unroll value = 8
#pragma unroll (X+1)
for (int i = 0; i < 12; ++i)
p1[i] += p2[i]*4;
// unroll value = 4
#pragma unroll (T2::value)
for (int i = 0; i < 12; ++i)
p1[i] += p2[i]*16;
}
‣ vadd2, vadd4
‣ vsub2, vsub4
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 137
C Language Extensions
‣ vavrg2, vavrg4
‣ vabsdiff2, vabsdiff4
‣ vmin2, vmin4
‣ vmax2, vmax4
‣ vset2, vset4
PTX instructions, such as the SIMD video instructions, can be included in CUDA
programs by way of the assembler, asm(), statement.
The basic syntax of an asm() statement is:
asm("template-string" : "constraint"(output) : "constraint"(input)"));
This uses the vabsdiff4 instruction to compute an integer quad byte SIMD sum of
absolute differences. The absolute difference value is computed for each byte of the
unsigned integers A and B in SIMD fashion. The optional accumulate operation (.add)
is specified to sum these differences.
Refer to the document "Using Inline PTX Assembly in CUDA" for details on using
the assembly statement in your code. Refer to the PTX ISA documentation ("Parallel
Thread Execution ISA Version 3.0" for example) for details on the PTX instructions for
the version of PTX that you are using.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 138
Appendix C.
CUDA DYNAMIC PARALLELISM
1. Introduction
1. Overview
Dynamic Parallelism is an extension to the CUDA programming model enabling a CUDA
kernel to create and synchronize with new work directly on the GPU. The creation of
parallelism dynamically at whichever point in a program that it is needed offers exciting
new capabilities.
The ability to create work directly from the GPU can reduce the need to transfer
execution control and data between host and device, as launch configuration decisions
can now be made at runtime by threads executing on the device. Additionally,
data-dependent parallel work can be generated inline within a kernel at run-time,
taking advantage of the GPU's hardware schedulers and load balancers dynamically
and adapting in response to data-driven decisions or workloads. Algorithms and
programming patterns that had previously required modifications to eliminate
recursion, irregular loop structure, or other constructs that do not fit a flat, single-level of
parallelism may more transparently be expressed.
This document describes the extended capabilities of CUDA which enable Dynamic
Parallelism, including the modifications and additions to the CUDA programming
model necessary to take advantage of these, as well as guidelines and best practices for
exploiting this added capacity.
Dynamic Parallelism is only supported by devices of compute capability 3.5 and higher.
2. Glossary
Definitions for terms used in this guide.
Grid
A Grid is a collection of Threads. Threads in a Grid execute a Kernel Function and are
divided into Thread Blocks.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 139
CUDA Dynamic Parallelism
Thread Block
A Thread Block is a group of threads which execute on the same multiprocessor
(SMX). Threads within a Thread Block have access to shared memory and can be
explicitly synchronized.
Kernel Function
A Kernel Function is an implicitly parallel subroutine that executes under the CUDA
execution and memory model for every Thread in a Grid.
Host
The Host refers to the execution environment that initially invoked CUDA. Typically
the thread running on a system's CPU processor.
Parent
A Parent Thread, Thread Block, or Grid is one that has launched new grid(s), the Child
Grid(s). The Parent is not considered completed until all of its launched Child Grids
have also completed.
Child
A Child thread, block, or grid is one that has been launched by a Parent grid. A Child
grid must complete before the Parent Thread, Thread Block, or Grid are considered
complete.
Thread Block Scope
Objects with Thread Block Scope have the lifetime of a single Thread Block. They only
have defined behavior when operated on by Threads in the Thread Block that created
the object and are destroyed when the Thread Block that created them is complete.
Device Runtime
The Device Runtime refers to the runtime system and APIs available to enable Kernel
Functions to use Dynamic Parallelism.
1. Execution Environment
The CUDA execution model is based on primitives of threads, thread blocks, and grids,
with kernel functions defining the program executed by individual threads within a
thread block and grid. When a kernel function is invoked the grid's properties are
described by an execution configuration, which has a special syntax in CUDA. Support
for dynamic parallelism in CUDA extends the ability to configure, launch, and
synchronize upon new grids to threads that are running on the device.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 140
CUDA Dynamic Parallelism
launched, the runtime guarantees an implicit synchronization between the parent and
child.
Time
CPU Thread
Grid A Threads
Grid A - Parent
3. Synchronization
CUDA runtime operations from any thread, including kernel launches, are visible across
a thread block. This means that an invoking thread in the parent grid may perform
synchronization on the grids launched by that thread, by other threads in the thread
block, or on streams created within the same thread block. Execution of a thread block
is not considered complete until all launches by all threads in the block have completed.
If all threads in a block exit before all child launches have completed, a synchronization
operation will automatically be triggered.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 141
CUDA Dynamic Parallelism
Streams and events created within a grid exist within thread block scope but have
undefined behavior when used outside of the thread block where they were created. As
described above, all work launched by a thread block is implicitly synchronized when
the block exits; work launched into streams is included in this, with all dependencies
resolved appropriately. The behavior of operations on a stream that has been modified
outside of thread block scope is undefined.
Streams and events created on the host have undefined behavior when used within any
kernel, just as streams and events created by a parent grid have undefined behavior if
used within a child grid.
6. Device Management
There is no multi-GPU support from the device runtime; the device runtime is only
capable of operating on the device upon which it is currently executing. It is permitted,
however, to query properties for any CUDA capable device in the system.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 142
CUDA Dynamic Parallelism
syncthreads();
if (threadIdx.x == 0) {
child_launch<<< 1, 256 >>>(data);
cudaDeviceSynchronize();
}
syncthreads();
}
3. Constant Memory
Constants are immutable and may not be modified from the device, even between
parent and child launches. That is to say, the value of all constant variables must
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 143
CUDA Dynamic Parallelism
be set from the host prior to launch. Constant memory is inherited automatically by all
child kernels from their respective parents.
Taking the address of a constant memory object from within a kernel thread has the
same semantics as for all CUDA programs, and passing that pointer from parent to child
or from a child to parent is naturally supported.
5. Local Memory
Local memory is private storage for an executing thread, and is not visible outside of
that thread. It is illegal to pass a pointer to local memory as a launch argument when
launching a child kernel. The result of dereferencing such a local memory pointer from a
child will be undefined.
For example the following is illegal, with undefined behavior if x_array is accessedby
child_launch:
int x_array[10]; // Creates x_array in parent's local memory
child_launch<<< 1, 1 >>>(x_array);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 144
CUDA Dynamic Parallelism
of a child grid and when a child grid completes. This means that writes to memory prior
to a child kernel launch are reflected in texture memory accesses of the child. Similarly,
writes to memory by a child will be reflected in the texture memory accesses by a parent,
but only after the parent synchronizes on the child's completion. Concurrent accesses by
parent and child may result in inconsistent data.
3. Programming Interface
‣ Dg is of type dim3 and specifies the dimensions and size of the grid
‣ Db is of type dim3 and specifies the dimensions and size of each thread block
‣ Ns is of type size_t and specifies the number of bytes of shared memory that
is dynamically allocated per thread block for this call and addition to statically
allocated memory. Ns is an optional argument that defaults to 0.
‣ S is of type cudaStream_t and specifies the stream associated with this call. The
stream must have been allocated in the same thread block where the call is being
made. S is an optional argument that defaults to 0.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 145
CUDA Dynamic Parallelism
may begin execution at any time after launch, but is not guaranteed to begin execution
until the launching thread reaches an explicit launch-synchronization point.
2. Streams
Both named and unnamed (NULL) streams are available from the device runtime.
Named streams may be used by any thread within a thread-block, but stream handles
may not be passed to other blocks or child/parent kernels. In other words, a stream
should be treated as private to the block in which it is created. Stream handles are not
guaranteed to be unique between blocks, so using a stream handle within a block that
did not allocate it will result in undefined behavior.
Similar to host-side launch, work launched into separate streams may run concurrently,
but actual concurrency is not guaranteed. Programs that depend upon concurrency
between child kernels are not supported by the CUDA programming model and will
have undefined behavior.
The host-side NULL stream's cross-stream barrier semantic is not supported on the
device (see below for details). In order to retain semantic compatibility with the host
runtime, all device streams must be created using the cudaStreamCreateWithFlags()
API, passing the cudaStreamNonBlocking flag. The cudaStreamCreate() call is a
host-runtime- only API and will fail to compile for the device.
As cudaStreamSynchronize() and cudaStreamQuery() are unsupported by
the device runtime, cudaDeviceSynchronize() should be used instead when the
application needs to know that stream-launched child kernels have completed.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 146
CUDA Dynamic Parallelism
a block, but as all named streams must be created with the cudaStreamNonBlocking
flag, work launched into the NULL stream will not insert an implicit dependency on
pending work in any other streams.
3. Events
Only the inter-stream synchronization capabilities of CUDA events are
supported. This means that cudaStreamWaitEvent() is supported, but
cudaEventSynchronize(), cudaEventElapsedTime(), and cudaEventQuery() are
not. As cudaEventElapsedTime() is not supported, cudaEvents must be created via
cudaEventCreateWithFlags(), passing the cudaEventDisableTiming flag.
As for all device runtime objects, event objects may be shared between all threads
withinthe thread-block which created them but are local to that block and may not be
passed to other kernels, or between blocks within the same kernel. Event handles are not
guaranteed to be unique between blocks, so using an event handle within a block that
did not create it will result in undefined behavior.
4. Synchronization
The cudaDeviceSynchronize() function will synchronize on allwork launched by
any thread in the thread-block up to the point where cudaDeviceSynchronize() was
called. Note that cudaDeviceSynchronize() may be called from within divergent
code (see Block Wide Synchronization).
It is up to the program to perform sufficient additional inter-thread synchronization, for
example via a call to syncthreads(), if the calling thread is intended to synchronize
with child grids invoked from other threads.
5. Device Management
Only the device on which a kernel is running will be controllable from that kernel.
This means that device APIs such as cudaSetDevice() are not supported by
the device runtime. The active device as seen from the GPU (returned from
cudaGetDevice()) will have the same device number as seen from the host system.
The cudaGetDeviceProperty() call may request information about another device
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 147
CUDA Dynamic Parallelism
as this API allows specification of a device ID as a parameter of the call. Note that the
catch-all cudaGetDeviceProperties() API is not offered by the device runtime -
properties must be queried individually.
6. Memory Declarations
The device runtime does not support legacy module-scope (i.e., Fermi-style) textures
and surfaces within a kernel launched from the device. Module-scope (legacy)
textures may be created from the host and used in device code as for any kernel,
but may only be used by a top-level kernel (i.e., the one which is launched from the
host).
1 Dynamically created texture and surface objects are an addition to the CUDA memory model introduced with CUDA
5.0. Please see the CUDA Programming Guide for details.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 148
CUDA Dynamic Parallelism
smem[threadIdx.x] = data[threadIdx.x];
syncthreads();
permute_data(smem, n);
syncthreads();
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 149
CUDA Dynamic Parallelism
The APIs for these launch functions are different to those of the CUDA Runtime API,
and are defined as follows:
extern device cudaError_t cudaGetParameterBuffer(void **params);
extern device cudaError_t cudaLaunchDevice(void *kernel,
void *params, dim3 gridDim,
dim3 blockDim,
unsigned int sharedMemSize = 0,
cudaStream_t stream = 0);
cudaDeviceGetLimit
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 150
CUDA Dynamic Parallelism
cudaGetErrorString
cudaGetDeviceCount
cudaStreamWaitEvent
cudaEventDestroy
cudaFuncGetAttributes
cudaMemset2DAsync
cudaMemset3DAsync
cudaRuntimeGetVersion
cudaOccupancyMaxActiveBlocksPerMultiproce ssor
cudaOccupancyMaxPotentialBlockSize
cudaOccupancyMaxPotentialBlockSizeVariabl eSMem
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 151
CUDA Dynamic Parallelism
to the launched kernel. The parameter buffer can be NULL, i.e., no need to invoke
cudaGetParameterBuffer(), if the launched kernel does not take any parameters.
C.3.2.1.1. cudaLaunchDevice
At the PTX level, cudaLaunchDevice()needs to be declared in one of the two forms
shown below before it is used.
// PTX-level Declaration of cudaLaunchDevice() when .address_size is 64
.extern .func(.param .b32 func_retval0) cudaLaunchDevice
(
.param .b64 func,
.param .b64 parameterBuffer,
.param .align 4 .b8 gridDimension[12],
.param .align 4 .b8 blockDimension[12],
.param .b32 sharedMemSize,
.param .b64 stream
)
;
// PTX-level Declaration of cudaLaunchDevice() when .address_size is 32
.extern .func(.param .b32 func_retval0) cudaLaunchDevice
(
.param .b32 func,
.param .b32 parameterBuffer,
.param .align 4 .b8 gridDimension[12],
.param .align 4 .b8 blockDimension[12],
.param .b32 sharedMemSize,
.param .b32 stream
)
;
The CUDA-level declaration below is mapped to one of the aforementioned PTX-level
declarations and is found in the system header file cuda_device_runtime_api.h.
The function is defined in the cudadevrt system library, which must be linked witha
program in order to use device-side kernel launch functionality.
// CUDA-level declaration of cudaLaunchDevice()
extern "C" device
cudaError_t cudaLaunchDevice(void *func, void *parameterBuffer,
dim3 gridDimension, dim3 blockDimension,
unsigned int sharedMemSize,
cudaStream_t stream);
The first parameter is a pointer to the kernel to be is launched, and the second parameter
is the parameter buffer that holds the actual parameters to the launched kernel. The
layout of the parameter buffer is explained in Parameter Buffer Layout, below. Other
parameters specify the launch configuration, i.e., as grid dimension, block dimension,
shared memory size, and the stream associated with the launch (please refer to
Execution Configuration for the detailed description of launch configuration.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 152
CUDA Dynamic Parallelism
C.3.2.1.2. cudaGetParameterBuffer
cudaGetParameterBuffer() needs to be declared at the PTX level before it's used.
The PTX-level declaration must be in one of the two forms given below, depending on
address size:
// PTX-level Declaration of cudaGetParameterBuffer() when .address_size is 64
// When .address_size is 64
.extern .func(.param .b64 func_retval0) cudaGetParameterBuffer
(
.param .b64 alignment,
.param .b64 size
)
;
// PTX-level Declaration of cudaGetParameterBuffer() when .address_size is 32
.extern .func(.param .b32 func_retval0) cudaGetParameterBuffer
(
.param .b32 alignment,
.param .b32 size
)
;
The following CUDA-level declaration of cudaGetParameterBuffer() is mapped to
the aforementioned PTX-level declaration:
// CUDA-level Declaration of cudaGetParameterBuffer()
extern "C" device
void *cudaGetParameterBuffer(size_t alignment, size_t size);
The first parameter specifies the alignment requirement of the parameter buffer and
the second parameter the size requirement in bytes. In the current implementation, the
parameter buffer returned by cudaGetParameterBuffer() is always guaranteed to
be 64- byte aligned, and the alignment requirement parameter is ignored. However,
it is recommended to pass the correct alignment requirement value - which is
the largest alignment of any parameter to be placed in the parameter buffer - to
cudaGetParameterBuffer() to ensure portability in the future.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 153
CUDA Dynamic Parallelism
It is also possible to compile CUDA .cu source files first to object files, and then link
these together in a two-stage process:
$ nvcc -arch=sm_35 -dc hello_world.cu -o hello_world.o
$ nvcc -arch=sm_35 -rdc=true hello_world.o -o hello -lcudadevrt
Please see the Using Separate Compilation section of The CUDA Driver Compiler NVCC
guide for more details.
4. Programming Guidelines
1. Basics
The device runtime is a functional subset of the host runtime. API level device
management, kernel launching, device memcpy, stream management, and event
management are exposed from the device runtime.
Programming for the device runtime should be familiar to someone who already has
experience with CUDA. Device runtime syntax and semantics are largely the same as
that of the host API, with any exceptions detailed earlier in this document.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 154
CUDA Dynamic Parallelism
The following example shows a simple Hello World program incorporating dynamic
parallelism:
#include <stdio.h>
printf("World!\n");
}
return 0;
}
This program may be built in a single step from the command line as follows:
2. Performance
1. Synchronization
Synchronization by one thread may impact the performance of other threads in the same
Thread Block, even when those other threads do not call cudaDeviceSynchronize()
themselves. This impact will depend upon the underlying implementation.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 155
CUDA Dynamic Parallelism
calls when made from the device compared to from the host side. This overhead is, in
general, incurred for applications that link against the device runtime library.
1. Runtime
1. Memory Footprint
The device runtime system software reserves memory for various management
purposes, in particular one reservation which is used for saving parent-grid state
during synchronization, and a second reservation for tracking pending grid launches.
Configuration controls are available to reduce the size of these reservations in exchange
for certain launch limitations. See Configuration Options, below, for details.
The majority of reserved memory is allocated as backing-store for parent kernel state, for
use when synchronizing on a child launch. Conservatively, this memory must support
storing of state for the maximum number of live threads possible on the device. This
means that each parent generation at which cudaDeviceSynchronize() is callable
may require up to 150MB of device memory, depending on the device configuration,
which will be unavailable for program use even if it is not all consumed.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 156
CUDA Dynamic Parallelism
4. Configuration Options
Resource allocation for the device runtime system software is controlled via the
cudaDeviceSetLimit() API from the host program. Limits must be set before
any kernel is launched, and may not be changed while the GPU is actively running
programs.
The following named limits may be set:
Limit Behavior
cudaLimitDevRuntimeSyncDepth Sets the maximum depth at which
cudaDeviceSynchronize() may be
called. Launches may be performed deeper
than this, but explicit synchronization
deeper than this limit will return the
cudaErrorLaunchMaxDepthExceeded. The
default maximum sync depth is 2.
cudaLimitDevRuntimePendingLaunchCount Controls the amount of memory set aside for
buffering kernel launches which have not yet
begun to execute, due either to unresolved
dependencies or lack of execution resources.
When the buffer is full, the device runtime
system software will attempt to track new
pending launches in a lower performance
virtualized buffer. If the virtualized buffer
is also full, i.e. when all available heap
space is consumed, launches will not occur,
and the thread's last error will be set to
cudaErrorLaunchPendingCountExceeded.
The default pending launch count is 2048
launches.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 157
CUDA Dynamic Parallelism
6. SM Id and Warp Id
Note that in PTX %smid and %warpid are defined as volatile values. The device runtime
may reschedule thread blocks onto different SMs in order to more efficiently manage
resources. As such, it is unsafe to rely upon %smid or %warpid remaining unchanged
across the lifetime of a thread or thread block.
7. ECC Errors
No notification of ECC errors is available to code within a CUDA kernel. ECC errors
are reported at the host side once the entire launch tree has completed. Any ECC errors
which arise during execution of a nested program will either generate an exception or
continue execution (depending upon error and configuration).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 158
Appendix D.
MATHEMATICAL FUNCTIONS
The reference manual lists, along with their description, all the functions of the C/C++
standard library mathematical functions that are supported in device code, as well as all
intrinsic functions (that are only supported in device code).
This appendix provides accuracy information for some of these functions when
applicable.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 159
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 160
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 161
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 162
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 163
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 164
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 165
Mathematical Functions
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 166
Mathematical Functions
sinf(x) sinf(x)
cosf(x) cosf(x)
tanf(x) tanf(x)
sincosf(x,sptr,cptr) sincosf(x,sptr,cptr)
logf(x) logf(x)
log2f(x) log2f(x)
log10f(x) log10f(x)
expf(x) expf(x)
exp10f(x) exp10f(x)
powf(x,y) powf(x,y)
Functions suffixed with _rn operate using the round to nearest even rounding mode.
Functions suffixed with _rz operate using the round towards zero roundingmode.
Functions suffixed with _ru operate using the round up (to positive infinity)rounding
mode.
Functions suffixed with _rd operate using theround down (to negative infinity)
rounding mode.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 167
Mathematical Functions
additions and multiplications generated from the '*' and '+' operators will frequently be
combined into FMADs.
The accuracy of floating-point division varies depending on whether the code is
compiled with -prec-div=false or -prec-div=true. When the code is compiled
with -prec-div=false, both the regular division / operator and fdividef(x,y)
have the same accuracy, but for 2126 < y < 2128, fdividef(x,y) delivers a result of
zero, whereas the / operator delivers the correct result to within the accuracy stated
in Table 9. Also, for 2126 < y < 2128, if x is infinity, fdividef(x,y) delivers a NaN (as a
result of multiplying infinity by zero), while the / operator returns infinity. On the
other hand, the / operator is IEEE-compliant when the code is compiled with -prec-
div=true or without any -prec-div option at all since its default value is true.
fadd_[rn,rz,ru,rd](x,y) IEEE-compliant.
fsub_[rn,rz,ru,rd](x,y) IEEE-compliant.
fmul_[rn,rz,ru,rd](x,y) IEEE-compliant.
fmaf_[rn,rz,ru,rd](x,y,z) IEEE-compliant.
frcp_[rn,rz,ru,rd](x) IEEE-compliant.
fsqrt_[rn,rz,ru,rd](x) IEEE-compliant.
frsqrt_rn(x) IEEE-compliant.
fdiv_[rn,rz,ru,rd](x,y) IEEE-compliant.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 168
Mathematical Functions
dadd_[rn,rz,ru,rd](x,y) IEEE-compliant.
dsub_[rn,rz,ru,rd](x,y) IEEE-compliant.
dmul_[rn,rz,ru,rd](x,y) IEEE-compliant.
fma_[rn,rz,ru,rd](x,y,z) IEEE-compliant.
ddiv_[rn,rz,ru,rd](x,y)(x,y) IEEE-compliant.
drcp_[rn,rz,ru,rd](x) IEEE-compliant.
dsqrt_[rn,rz,ru,rd](x) IEEE-compliant.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 169
Appendix E.
C/C++ LANGUAGE SUPPORT
As described in Compilation with NVCC, CUDA source files compiled with nvcc can
include a mix of host code and device code. The CUDA frontend compiler aims to
emulate the host compiler behavior with respect to C++ input code. The input source
code is processed according to the C++ ISO/IEC 14882:2003 or C++ ISO/IEC 14882:2011
specifications, and the CUDA frontend compiler aims to emulate any host compiler
divergences from the ISO specification. In addition, the supported language is extended
with CUDA-specific constructs described in this document 7, and is subject to the
restrictions described below.
C++11 Language Features provides a support matrix for the C++11 features. Restrictions
lists the language restrictions. Polymorphic Function Wrappers and #unique_333
describe additional features. Code Samples gives code samples.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 170
C/C++ Language Support
Available
C++11 in nvcc
Language Feature (device
Proposal
code)
Variadic templates N2242 7.0
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 171
C/C++ Language Support
Available
C++11 in nvcc
Language Feature (device
Proposal
code)
Standard Layout Types N2342 7.0
N3206
Explicit virtual overrides 7.0
N3272
Minimal support for garbage collection and reachability-based leak N2670 N/A (see
detection Restrictions)
Allowing move constructors to throw [noexcept] N3050 7.0
Defining move special member functions N3053 7.0
Concurrency
Sequence points N2239
Atomic operations N2427
Strong Compare and Exchange N2748
Bidirectional Fences N2752
Memory model N2429
Data-dependency ordering: atomics and memory model N2664
Propagating exceptions N2179
Allow atomics use in signal handlers N2547
Thread-local storage N2659
Dynamic initialization and destruction with concurrency N2660
C99 Features in C++11
func predefined identifier N2340 7.0
C99 preprocessor N1653 7.0
long long N1811 7.0
Extended integral types N1988
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 172
C/C++ Language Support
2. Restrictions
1. Host Compiler Extensions
The use of host compiler specific language extensions in device code is unsupported.
2. Preprocessor Symbols
1. CUDA_ARCH
1. The type signature of the following entities shall not depend on whether
CUDA_ARCH is defined or not, or on a particular value of CUDA_ARCH :
‣ global functions and function templates
‣ device and constant variables
‣ textures and surfaces
Example:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 173
C/C++ Language Support
Example:
device int result;
template <typename T>
global void kern(T in)
{ result = in;
}
host device void foo(void)
{
#if !defined( CUDA_ARCH )
kern<<<1,1>>>(1); // error: "kern<int>" instantiation only
// when CUDA_ARCH is undefined!
#endif
}
int main(void)
{
foo();
cudaDeviceSynchronize();
return 0;
}
3. In separate compilation mode, the presence or absence of a definition of a function
or variable with external linkage shall not depend on whether CUDA_ARCH is
defined or on a particular value of CUDA_ARCH 8.
Example:
8 This does notapply to entities that may be defined in more than one translation unit, such as compiler generated
template instantiations.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 174
C/C++ Language Support
3. Qualifiers
1. Device Memory Qualifiers
The device , shared and constant qualifiers are not allowed on:
‣ class, struct, and union data members,
‣ formal parameters,
‣ local variables within a function that executes on the host.
shared and constant variables have implied static storage.
device and constant variable definitions are only allowed in namespace
scope (including global namespace scope).
device , constant and shared variables defined in namespace scope,
that are of class type, cannot have a non-empty constructor or a non-empty destructor. A
constructor for a class type is considered empty at a point in the translation unit, if it is
either a trivial constructor or it satisfies all of the following conditions:
‣ The constructor function has been defined.
‣ The constructor function has no parameters, the initializer list is empty and the
function body is an empty compound statement.
‣ Its class has no virtual functions and no virtual base classes.
‣ The default constructors of all base classes of its class can be considered empty.
‣ For all the nonstatic data members of its class that are of class type (or array thereof),
the default constructors can be considered empty.
A destructor for a class is considered empty at a point in the translation unit, if it is
either a trivial destructor or it satisfies all of the following conditions:
‣ The destructor function has been defined.
‣ The destructor function body is an empty compound statement.
‣ Its class has no virtual functions and no virtual base classes.
‣ The destructors of all base classes of its class can be considered empty.
‣ For all the nonstatic data members of its class that are of class type (or array thereof),
the destructor can be considered empty.
When compiling in the whole program compilation mode (see the nvcc user manual for
a description of this mode), device , shared , and constant variables
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 175
C/C++ Language Support
cannot be defined as external using the extern keyword. The only exception is for
dynamically allocated shared variables as described in shared .
When compiling in the separate compilation mode (see the nvcc user manual for a
description of this mode), device , shared , and constant variables can
be defined as external using the extern keyword. nvlink will generate an error when
it cannot find a definition for an external variable (unless it is a dynamically allocated
shared variable).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 176
C/C++ Language Support
device managed const int yyy = 10; // error: const qualified type
device managed int &zzz = xxx; // error: reference type
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 177
C/C++ Language Support
4. Pointers
Dereferencing a pointer either to global or shared memory in code that is executed
on the host, or to host memory in code that is executed on the device results in an
undefined behavior, most often in a segmentation fault and application termination.
The address obtained by taking the address of a device , shared or
constant variable can only be used in device code. The address of a device
or constant variable obtained through cudaGetSymbolAddress() as described
in Device Memory can only be used in host code.
As a consequence of the use of C++ syntax rules, void pointers (e.g., returnedby
malloc()) cannot be assigned to non-void pointers without atypecast.
5. Operators
1. Assignment Operator
constant variables can only be assigned from the host code through runtime
functions (Device Memory); they cannot be assigned from the device code.
shared variables cannot have an initialization as part of their declaration.
It is not allowed to assign values to any of the built-in variables defined in Built-in
Variables.
2. Address Operator
It is not allowed to take the address of any of the built-in variables defined in Built-in
Variables.
7. Exception Handling
Exception handling is only supported in host code, but not in device code.
8. Standard Library
Standard libraries are only supported in host code, but not in device code, unless
specified otherwise.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 178
C/C++ Language Support
9. Functions
1. External Linkage
A call within some device code of a function declared with the extern qualifier is only
allowed if the function is defined within the same compilation unit as the device code,
i.e., a single file or several files linked together with relocatable device code and nvlink.
Other D2;
}
host void bar(void)
{ Other D3;
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 179
C/C++ Language Support
int x = 33;
static int i6 = x; // error: dynamic initialization is not allowed
static S1_t i7 = {x}; // error: dynamic initialization is not allowed
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 180
C/C++ Language Support
10. Classes
1. Data Members
Static data members are not supported.
2. Function Members
Static member functions cannot be global functions.
3. Virtual Functions
When a function in a derived class overrides a virtual function in a base class, the
execution space qualifiers (i.e., host , device ) on the overridden and
overriding functions must match.
It is not allowed to pass as an argument to a global function an object of a class
with virtual functions.
The virtual function table is placed in global or constant memory by the compiler.
5. Anonymous Unions
Member variables of a namespace scope anonymous union cannot be referenced in a
global or device function.
6. Windows-Specific
The CUDA compiler follows the IA64 ABI for class layout, while the Microsoft host
compiler does not. This may cause the CUDA compiler to compute the class layout and
size differently than the Microsoft host compiler, for a class type 'T' that satisfies any of
the following conditions or for any class type that has T as a field type or as a base class
type (direct or indirect):
‣ T has virtual functions.
‣ T has a virtual base class.
‣ T has multiple inheritance with more than one direct or indirect empty base class.
‣ All direct and indirect base classes ('B') of T are empty and the type of the first field
of T ('F') uses B in its definition, such that B is laid out at offset 0 in the definition of
F.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 181
C/C++ Language Support
As long as affected class types are used exclusively in host or device code, the program
should work correctly; do not pass objects of such class types between between host and
device code (e.g., as arguments to global functions or through cudaMemcpy*()
calls) 10.
E.2.11. Templates
A type or template cannot be used in the type, non-type or template template argument
of a global function template instantiation if either:
‣ The type or template is defined within a host or host device .
‣ The type or template is a class member with private or protected access and its
parent class is not defined within a device or global function.
Example:
template <typename T>
global void myKernel(void) { }
class myClass {
private:
struct inner_t { };
public:
static void launch(void)
{
// error: inner_t is used in template argument
// but it is private
myKernel<inner_t><<<1,1>>>();
}
};
10 One way to debug suspected layout mismatch of a type C is to use printf to output the values of sizeof(C) and
offsetof(C, field) in host and device code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 182
C/C++ Language Support
Example:
const int xxx = 10; int yyy = 20; };
struct S1_t { static const
1. Lambda Expressions
The execution space qualifiers for all member functions12 of the closure class associated
with a lambda expression are derived by the compiler as follows. As described in
the C++11 standard, the compiler creates a closure type in the smallest block scope,
class scope or namespace scope that contains the lambda expression. The innermost
function scope enclosing the closure type is computed, and the corresponding function's
execution space qualifiers are assigned to the closure class member functions. If there is
no enclosing function scope, the execution space qualifier is host .
11 At present, the -std=c++11 flag is supported only for the following host compilers : gcc version >= 4.7, clang, icc >= 15
(without extended lambda), and xlc >= 13.1
12 including operator()
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 183
C/C++ Language Support
Examples of lambda expressions and computed execution space qualifiers are shown
below (in comments).
auto globalVar = [] { return 0; }; // host
void f1(void) {
auto l1 = [] { return 1; }; // host
}
The closure type of a lambda expression cannot be used in the type or non-type
argument of a global function template instantiation, unless the lambda is defined
within a device or global function.
Example:
void bar(void) {
auto temp1 = [] { };
foo<<<1,1>>>(temp1); // error: lambda closure type used in
// template type argument
foo<<<1,1>>>( S1_t<decltype(temp1)>()); // error: lambda closure type used in
// template type argument
}
E.2.14.2. std::initializer_list
By default, the CUDA compiler will implicitly consider the member functions of
std::initializer_list to have host device execution space qualifiers,
and therefore they can be invoked directly from device code. The nvcc flag --no-
host-device-initializer-list will disable this behavior; member functions of
std::initializer_list will then be considered as host functions and will not
be directly invokable from device code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 184
C/C++ Language Support
Example:
#include <initializer_list>
int i = 4;
foo({i,5,6}); // (b) initializer list with at least one
// non-constant element.
// This form may have better performance than (a).
}
3. Rvalue references
By default, the CUDA compiler will implicitly consider std::move and std::forward
function templates to have host device execution space qualifiers, and
therefore they can be invoked directly from device code. The nvcc flag --no-host-
device-move-forward will disable this behavior; std::move and std::forward
will then be considered as host functions and will not be directly invokable from
device code.
5. Constexpr variables
Let 'V' denote a namespace scope variable or a class static member variable that has
been marked constexpr and that does not have execution space annotations (e.g.,
device , constant , shared ). V is considered to be a host code variable.
If V is of scalar type 14 other than long double, the value of V can be directly used in
device code. In addition, the value of V can be used inside a constexpr device or
host device function, if the call to the function is a constant expression 15.
Device source code cannot contain a reference to V or take the address of V.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 185
C/C++ Language Support
Example:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 186
C/C++ Language Support
Example:
inline namespace N1 {
namespace N2 {
device int Gvar;
}
}
namespace N2 {
device int Gvar;
}
inline namespace {
namespace N2 {
template <typename T>
global void foo(void); // error
template <>
global void foo<int>(void) { } // error
7. thread_local
The thread_local storage specifier is not allowed in device code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 187
C/C++ Language Support
Example:
kernel<<<1,1>>>( [] device { } );
kernel<<<1,1>>>( [] host device { } );
kernel<<<1,1>>>( [] { } );
}
auto lam1 = [] { };
auto lam2 = [] host device { };
void foo_host(void)
{
// OK: instantiated with closure type of an extended device lambda
kernel<<<1,1>>>( [] device { } );
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 188
C/C++ Language Support
Example:
// ok
template <template <typename...> class Wrapper, typename... Pack>
global void foo1(Wrapper<Pack...>);
// error: pack parameter is not last in parameter list
template <typename... Pack, template <typename...> class Wrapper>
global void foo2(Wrapper<Pack...>);
// error: multiple parameter packs
template <typename... Pack1, int...Pack2, template<typename...> class Wrapper1,
template<int...> class Wrapper2>
global void foo3(Wrapper1<Pack1...>, Wrapper2<Pack2...>);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 189
C/C++ Language Support
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 190
C/C++ Language Support
namespace nvstd {
template <class _RetType, class ..._ArgTypes>
class function<_RetType(_ArgTypes...)>
{
public:
// constructors
device host function() noexcept;
device host function(nullptr_t) noexcept;
device host function(const function &);
device host function(function &&);
template<class _F>
device host function(_F);
// destructor
device host ~function();
// assignment operators
device host function& operator=(const function&);
device host function& operator=(function&&);
device host function& operator=(nullptr_t);
device host function& operator=(_F&&);
// swap
device host void swap(function&) noexcept;
// function capacity
device host explicit operator bool() const noexcept;
// function invocation
device _RetType operator()(_ArgTypes...) const;
};
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 191
C/C++ Language Support
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 192
C/C++ Language Support
Example:
void foo_host(void)
{
// not an extended lambda: no explicit execution space annotations
auto lam1 = [] { };
// not an extended lambda: explicitly annotated with only ' host '
auto lam4 = [] host { };
}
// not an extended lambda: explicitly annotated with only ' host '
auto lam4 = [] host { };
}
// lam1 and lam2 are not extended lambdas because they are not defined
// within a host or host device function.
auto lam1 = [] { };
auto lam2 = [] host device { };
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 193
C/C++ Language Support
These traits can be used in all compilation modes, irrespective of whether lambdas or
extended lambdas are enabled16.
Example:
void foo(void) {
auto lam1 = [] { };
auto lam2 = [] device { };
auto lam3 = [] host device { };
16 The traits will always return false if extended lambda mode is not active.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 194
C/C++ Language Support
Example:
void foo(void) {
// OK
auto lam1 = [] device { return 0; };
{
// OK
auto lam2 = [] device { return 0; };
// OK
auto lam3 = [] device host { return 0; };
}
auto outer = [] {
// Error: enclosing function (operator()) is a class member
// of an unnamed class
auto lam3 = [] device { return 0; };
};
}
auto outer = [] {
// Error: enclosing function (operator()) is a class member
// of an unnamed class
auto lam3 = [] device host { return 0; };
};
struct S1_t {
S1_t(void) {
// Error: cannot take address of enclosing function
auto lam4 = [] device { return 0; };
}
};
class C0_t {
void foo(void) {
// Error: enclosing function has private access in parent class
auto temp1 = [] device { return 10; };
}
struct S2_t {
void foo(void) {
// Error: enclosing class S2_t has private access in its
// parent class
auto temp1 = [] device { return 10; };
}
};
2. An extended
}; lambda cannot be defined in a class that is local to a function.
Example:
void foo(void) {
struct S1_t {
void bar(void) {
// Error: bar is member of a class that is local to a function.
auto lam4 = [] host device { return 0; };
}
};
}
3. If the enclosing function is an instantiation of a function template or a member
function template, and/or the function is a member of a class template, the
template(s) must satisfy the following constraints:
‣ The template must have at most one variadic parameter, and it must be listed
last in the template parameter list.
‣ The template parameters must be named.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 195
C/C++ Language Support
‣ The template instantiation argument types cannot involve types that are either
local to a function (except for closure types for extended lambdas), or are private
or protected class members.
Example:
int main()
{
foo<char, int, float> f1;
foo<char, int> f2;
bar1(f1, f2);
bar2(f1, 10);
bar3<int, 10>();
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 196
C/C++ Language Support
Example:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 197
C/C++ Language Support
Example
#include <type_traits>
void foo(void)
{
auto lam1 = [] device { return 10; };
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 198
C/C++ Language Support
Example
device int result;
template <typename T>
global void kernel(T in) { result = in(); }
void foo(void) {
int x1 = 1;
auto lam1 = [=] host device {
// Error: "x1" is only captured when CUDA_ARCH is defined.
#ifdef CUDA_ARCH
return x1 + 1;
#else
return 10;
#endif
};
kernel<<<1,1>>>(lam1);
}
9. As described previously, the CUDA compiler replaces an extended host
device lambda expression with an instance of a placeholder type in the
code sent to the host compiler. For an extended host device lambda
that did not specify a capture default, or which did not capture any variables, the
placeholder type provides a pointer-to-function conversion operator, but only when
g++ is the host compiler. Thus, attempting to convert such an extended host
device lambda into a function pointer in host code will cause host compilation
failure when the host compiler is not g++.
Example
void foo(void)
{
auto lam_d = [] device (double) { return 1; };
auto lam_hd = [] host device (double) { return 1; };
kern<<<1,1>>>(lam_d);
kern<<<1,1>>>(lam_hd);
The CUDA compiler will generate compiler diagnostics for a subset of cases described
in 1-5; no diagnostic will be generated for cases 6-9, but the host compiler may fail to
compile the generated code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 199
C/C++ Language Support
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 200
C/C++ Language Support
Example:
#include <cstdio>
struct S1_t {
int xxx;
host device S1_t(void) : xxx(10) { };
void doit(void) {
};
// Kernel launch fails at run time because 'this->xxx'
// is not accessible from the GPU
foo<<<1,1>>>(lam1);
cudaDeviceSynchronize();
}
};
int main(void) {
S1_t s1;
s1.doit();
}
C++17 solves this problem by adding a new "*this" capture mode. In this mode, the
compiler makes a copy of the object denoted by "*this" instead of capturing the pointer
this by value. The "*this" capture mode is described in more detail here: http://
www.open-std.org/jtc1/sc22/wg21/docs/papers/2016/p0018r3.html .
The CUDA compiler supports the "*this" capture mode for lambdas defined within
device and global functions and for extended device lambdas
defined in host code, when the --expt-extended-lambda nvcc flag is used.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 201
C/C++ Language Support
struct S1_t {
int xxx;
host device S1_t(void) : xxx(10) { };
void doit(void) {
};
// Kernel launch succeeds
foo<<<1,1>>>(lam1);
cudaDeviceSynchronize();
}
};
int main(void) {
S1_t s1;
s1.doit();
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 202
C/C++ Language Support
"*this" capture mode is not allowed for unannotated lambdas defined in host code, or for
extended host device lambdas. Examples of supported and unsupported
usage:
struct S1_t {
int xxx;
host device S1_t(void) : xxx(10) { };
void host_func(void) {
};
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 203
C/C++ Language Support
Example:
namespace N1 {
struct S1_t { };
template <typename T> void foo(T);
};
namespace N2 {
template <typename T> int foo(T);
In the example above, the CUDA compiler replaced the extended lambda with a
placeholder type that involves the N1 namespace. As a result, the namespace N1
participates in the ADL lookup for foo(in) in the body of N2::doit, and host
compilation fails because multiple overload candidates N1::foo and N2::foo are
found.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 204
C/C++ Language Support
5. Code Samples
1. Data Aggregation Class
class PixelRGBA {
public:
device PixelRGBA(): r_(0), g_(0), b_(0), a_(0) { }
private:
unsigned char r_, g_, b_, a_;
device
PixelRGBA operator+(const PixelRGBA& p1, const PixelRGBA& p2)
{
return PixelRGBA(p1.r_ + p2.r_, p1.g_ + p2.g_,
p1.b_ + p2.b_, p1.a_ + p2.a_);
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 205
C/C++ Language Support
int main()
{
...
useValues<int><<<blocks, threads>>>(buffer);
...
}
template <>
device bool func<int>(T x) // Specialization
{
return true;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 206
C/C++ Language Support
class Sub {
public:
device float operator() (float a, float b) const
{
return a - b;
}
};
// Device code
template<class O> global
void VectorOperation(const float * A, const float * B, float * C,
unsigned int N, O op)
{
unsigned int iElement = blockDim.x * blockIdx.x + threadIdx.x;
if (iElement < N)
C[iElement] = op(A[iElement], B[iElement]);
}
// Host code
int main()
{
...
VectorOperation<<<blocks, threads>>>(v1, v2, v3, N, Add());
...
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 207
Appendix F.
TEXTURE FETCHING
This appendix gives the formula used to compute the value returned by the texture
functions of Texture Functions depending on the various attributes of the texture
reference (see Texture and Surface Memory).
The texture bound to the texture reference is represented as an array T of
‣ N texels for a one-dimensional texture,
‣ N x M texels for a two-dimensional texture,
‣ N x M x L texels for a three-dimensional texture.
It is fetched using non-normalized texture coordinates x, y, and z, or the normalized
texture coordinates x/N, y/M, and z/L as described in Texture Memory. In this appendix,
the coordinates are assumed to be in the valid range. Texture Memory explained how
out-of-range coordinates are remapped to the valid range based on the addressing
mode.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 208
Texture Fetching
tex(x )
T[3]
T[0]
T[2]
T[1]
x
0 1 2 3 4 Non- Normalize
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 209
Texture Fetching
tex(x)
T[3]
T[0]
T[2]
T[1]
x
0 1 2 3 4 Non- Normalize
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 210
Texture Fetching
TL(x)
T[3]
T[0]
T[2]
T[1]
x
0 4/ 3 8/ 3 4
0 1/ 3 2/ 3 1
Figure 15 One-Dimensional Table Lookup Using Linear Filtering
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 211
Appendix G.
COMPUTE CAPABILITIES
The general specifications and features of a compute device depend on its compute
capability (see Compute Capability).
Table 12 gives the features and technical specifications associated to each compute
capability.
Floating-Point Standard reviews the compliance with the IEEE floating-point standard.
Sections Compute Capability 2.x, Compute Capability 3.x, Compute Capability 5.x, and
Compute Capability 6.x give more details on the architecture of devices of compute
capability 2.x, 3.x, 5.x, and 6.x respectively.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 212
Compute Capabilities
Maximum dimensionality of 3
grid of thread blocks
Maximum x-dimension of a grid 65535 231-1
of thread blocks
Maximum y- or z-dimension of 65535
a grid of thread blocks
Maximum dimensionality of 3
thread block
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 213
Compute Capabilities
Compute Capability
Technical Specifications 2.x 3.0 3.2 3.5 3.7 5.0 5.2 5.3 6.0 6.1 6.2
Maximum x- or y-dimension of 1024
a block
Maximum z-dimension of a 64
block
Maximum number of threads 1024
per block
Warp size 32
Maximum number of resident 8 16 32
blocks per multiprocessor
Maximum number of resident 48 64
warps per multiprocessor
Maximum number of resident 1536 2048
threads per multiprocessor
Number of 32-bit registers per 32 K 64 K 128 64 K
multiprocessor K
Maximum number of 32-bit 32 K 64 K 32 K 64 K 32 K 64 K 32 K
registers per thread block
Maximum number of 32-bit 63 255
registers per thread
Maximum amount of shared 48 KB 112 64 96 64 KB 96 64
memory per multiprocessor KB KB KB KB KB
Maximum amount of shared 48 KB
memory per thread block
Number of shared memory 32
banks
Amount of local memory per 512 KB
thread
Constant memory size 64 KB
Cache working set per
multiprocessor for constant 8 KB 10 KB
memory
Cache working set per 12 Between 24
multiprocessor for texture KB Between 12 KB and 48 KB KB and 48 KB
memory
Maximum width for a 1D
texture reference bound to a 65536
CUDA array
Maximum width for a 1D 227
texture reference bound to
linear memory
Maximum width and number of
layers for a 1D layered texture 16384 x 2048
reference
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 214
Compute Capabilities
Compute Capability
Technical Specifications 2.x 3.0 3.2 3.5 3.7 5.0 5.2 5.3 6.0 6.1 6.2
Maximum width and height for
a 2D texture reference bound 65536 x 65535
to a CUDA array
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 215
Compute Capabilities
Compute Capability
Technical Specifications 2.x 3.0 3.2 3.5 3.7 5.0 5.2 5.3 6.0 6.1 6.2
Maximum width (and height)
and number of layers for a
cubemap layered surface 32768 x 2046
reference
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 216
Compute Capabilities
The behavior of integer division by zero and integer overflow is left undefined by
IEEE-754. For compute devices, there is no mechanism for detecting that such integer
operation exceptions have occurred. Integer division by zero yields an unspecified,
machine-specific value.
http://developer.nvidia.com/content/precision-performance-floating-point-and-ieee-754-
compliance-nvidia-gpus includes more information on the floating point accuracy and
compliance of NVIDIA GPUs.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 217
Compute Capabilities
// Device code
global void MyKernel(int* foo, int* bar, int a)
{
...
}
// Host code
// Runtime API
// cudaFuncCachePreferShared: shared memory is 48 KB
// cudaFuncCachePreferL1: shared memory is 16 KB
// cudaFuncCachePreferNone: no preference
cudaFuncSetCacheConfig(MyKernel, cudaFuncCachePreferShared)
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 218
Compute Capabilities
If the size of the words accessed by each thread is more than 4 bytes, a memory
request by a warp is first split into separate 128-byte memory requests that are issued
independently:
‣ Two memory requests, one for each half-warp, if the size is 8 bytes,
‣ Four memory requests, one for each quarter-warp, if the size is 16 bytes.
Each memory request is then broken down into cache line requests that are issued
independently. A cache line request is serviced at the throughput of L1 or L2 cache in
case of a cache hit, or at the throughput of device memory, otherwise.
Note that threads can access any words in any order, including the same words.
If a non-atomic instruction executed by a warp writes to the same location in global
memory for more than one of the threads of the warp, only one thread performs a write
and which thread does it is undefined.
Figure 16 shows some examples of global memory accesses and corresponding memory
transactions.
Figure 18 shows some examples of memory read accesses that involve the broadcast
mechanism for devices of compute capability 3.x. The same examples apply for devices
of compute capability 2.x.
In this case, threads tid and tid+n access the same bank whenever s*n is a multiple of
the number of banks (i.e., 32) or, equivalently, whenever n is a multiple of 32/d where
d is the greatest common divisor of 32 and s. As a consequence, there will be no bank
conflict only if the warp size (i.e., 32) is less than or equal to 32/d, i.e., only if d is equal
to 1, i.e., s is odd.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 219
Compute Capabilities
Figure 17 shows some examples of strided access for devices of compute capability 3.x.
The same examples apply for devices of compute capability 2.x. However, the access
pattern for the example in the middle generates 2-way bank conflicts for devices of
compute capability 2.x.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 220
Compute Capabilities
‣ 192 CUDA cores for arithmetic operations (see Arithmetic Instructions for
throughputs of arithmetic operations),
‣ 32 special function units for single-precision floating-point transcendental functions,
‣ 4 warp schedulers.
When a multiprocessor is given warps to execute, it first distributes them among
the four schedulers. Then, at every instruction issue time, each scheduler issues two
independent instructions for one of its assigned warps that is ready to execute, if any.
A multiprocessor has a read-only constant cache that is shared by all functional units
and speeds up reads from the constant memory space, which resides in device memory.
There is an L1 cache for each multiprocessor and an L2 cache shared by all
multiprocessors. The L1 cache is used to cache accesses to local memory, including
temporary register spills. The L2 cache is used to cache accesses to local and global
memory. The cache behavior (e.g., whether reads are cached in both L1 and L2 or in L2
only) can be partially configured on a per-access basis using modifiers to the load or
store instruction. Some devices of compute capability 3.5 and devices of compute
capability 3.7 allow opt-in to caching of global memory in both L1 and L2 via compiler
options.
The same on-chip memory is used for both L1 and shared memory: It can be configured
as 48 KB of shared memory and 16 KB of L1 cache or as 16 KB of shared memory
and 48 KB of L1 cache or as 32 KB of shared memory and 32 KB of L1 cache, using
cudaFuncSetCacheConfig()/cuFuncSetCacheConfig():
// Device code
global void MyKernel()
{
...
}
// Host code
// Runtime API
// cudaFuncCachePreferShared: shared memory is 48 KB
// cudaFuncCachePreferEqual: shared memory is 32 KB
// cudaFuncCachePreferL1: shared memory is 16 KB
// cudaFuncCachePreferNone: no preference
cudaFuncSetCacheConfig(MyKernel, cudaFuncCachePreferShared)
The default cache configuration is "prefer none," meaning "no preference."
If a kernel is configured to have no preference, then it will default
to the preference of the current thread/context, which is set using
cudaDeviceSetCacheConfig()/cuCtxSetCacheConfig() (see the reference manual
for details). If the current thread/context also has no preference (which is again the
default setting), then whichever cache configuration was most recently used for any
kernel will be the one that is used, unless a different cache configuration is required to
launch the kernel (e.g., due to shared memory requirements). The initial configuration is
48 KB of shared memory and 16 KB of L1 cache.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 221
Compute Capabilities
Applications may query the L2 cache size by checking the l2CacheSize device property
(see Device Enumeration). The maximum L2 cache size is 1.5 MB.
Each multiprocessor has a read-only data cache of 48 KB to speed up reads from device
memory. It accesses this cache either directly (for devices of compute capability 3.5
or 3.7), or via a texture unit that implements the various addressing modes and data
filtering mentioned in Texture and Surface Memory. When accessed via the texture unit,
the read-only data cache is also referred to as texture cache.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 222
Compute Capabilities
64-Bit Mode
Successive 64-bit words map to successive banks.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 223
Compute Capabilities
A shared memory request for a warp does not generate a bank conflict between two
threads that access any sub-word within the same 64-bit word (even though the
addresses of the two sub-words fall in the same bank): In that case, for read accesses, the
64-bit word is broadcast to the requesting threads and for write accesses, each sub-word
is written by only one of the threads (which thread performs the write is undefined).
In this mode, the same access pattern generates fewer bank conflicts than on devices of
compute capability 2.x for 64-bit accesses and as many or fewer for 32-bit accesses.
32-Bit Mode
Successive 32-bit words map to successive banks.
A shared memory request for a warp does not generate a bank conflict between two
threads that access any sub-word within the same 32-bit word or within two 32-bit
words whose indices i and j are in the same 64-word aligned segment (i.e., a segment
whose first index is a multiple of 64) and such that j=i+32 (even though the addresses of
the two sub-words fall in the same bank): In that case, for read accesses, the 32-bit words
are broadcast to the requesting threads and for write accesses, each sub-word is written
by only one of the threads (which thread performs the write is undefined).
In this mode, the same access pattern generates as many or fewer bank conflicts than on
devices of compute capability 2.x.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 224
Compute Capabilities
There is also an L2 cache shared by all multiprocessors that is used to cache accesses to
local or global memory, including temporary register spills. Applications may query the
L2 cache size by checking the l2CacheSize device property (see Device Enumeration).
The cache behavior (e.g., whether reads are cached in both the unified L1/texture cache
and L2 or in L2 only) can be partially configured on a per-access basis using modifiers to
the load instruction.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 225
Compute Capabilities
the requesting threads and for write accesses, each address is written by only one of the
threads (which thread performs the write is undefined).
Figure 17 shows some examples of strided access.
Figure 18 shows some examples of memory read accesses that involve the broadcast
mechanism.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 226
Compute Capabilities
0 0 0 0 0 0
1 1 1 1 1 1
2 2 2 2 2 2
3 3 3 3 3 3
4 4 4 4 4 4
5 5 5 5 5 5
6 6 6 6 6 6
7 7 7 7 7 7
8 8 8 8 8 8
9 9 9 9 9 9
10 10 10 10 10 10
11 11 11 11 11 11
12 12 12 12 12 12
13 13 13 13 13 13
14 14 14 14 14 14
15 15 15 15 15 15
16 16 16 16 16 16
17 17 17 17 17 17
18 18 18 18 18 18
19 19 19 19 19 19
20 20 20 20 20 20
21 21 21 21 21 21
22 22 22 22 22 22
23 23 23 23 23 23
24 24 24 24 24 24
25 25 25 25 25 25
26 26 26 26 26 26
27 27 27 27 27 27
28 28 28 28 28 28
29 29 29 29 29 29
30 30 30 30 30 30
31 31 31 31 31 31
Left
Linear addressing with a stride of one 32-bit word (no bank conflict).
Middle
Linear addressing with a stride of two 32-bit words (two-way bank conflict).
Right
Linear addressing with a stride of three 32-bit words (no bank conflict).
Figure 17 Strided Shared Memory Accesses
Examples for devices of compute capability 3.x (in 32-bit mode) or compute capability 5.x
and 6.x
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 227
Compute Capabilities
0 0 0 0 0 0
1 1 1 1 1 1
2 2 2 2 2 2
3 3 3 3 3 3
4 4 4 4 4 4
5 5 5 5 5 5
6 6 6 6 6 6
7 7 7 7 7 7
8 8 8 8 8 8
9 9 9 9 9 9
10 10 10 10 10 10
11 11 11 11 11 11
12 12 12 12 12 12
13 13 13 13 13 13
14 14 14 14 14 14
15 15 15 15 15 15
16 16 16 16 16 16
17 17 17 17 17 17
18 18 18 18 18 18
19 19 19 19 19 19
20 20 20 20 20 20
21 21 21 21 21
21
22 22 22 22 22 22
23 23 23 23 23 23
24 24 24 24 24 24
25 25 25 25 25 25
26 26 26 26 26 26
27 27 27 27 27 27
28 28 28 28 28 28
29 29 29 29 29 29
30 30 30 30 30 30
31 31 31 31 31 31
Left
Conflict-free access via random permutation.
Middle
Conflict-free access since threads 3, 4, 6, 7, and 9 access the same word within bank 5.
Right
Conflict-free broadcast access (threads access the same word within a bank).
Figure 18 Irregular Shared Memory Accesses
Examples for devices of compute capability 3.x, 5.x, or 6.x.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 228
Compute Capabilities
2. Global Memory
Global memory behaves the same way as devices of compute capability 5.x (See Global
Memory).
3. Shared Memory
Shared memory behaves the same way as devices of compute capability 5.x (See Shared
Memory).
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 229
Appendix H.
DRIVER API
The driver API must be initialized with cuInit() before any function from the driver
API is called. A CUDA context must then be created that is attached to a specific device
and made current to the calling host thread as detailed in Context.
Within a CUDA context, kernels are explicitly loaded as PTX or binary objects by the
host code as described in Module. Kernels written in C must therefore be compiled
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 230
Driver API
separately into PTX or binary objects. Kernels are launched using API entry points as
described in Kernel Execution.
Any application that wants to run on future device architectures must load PTX,
not binary code. This is because binary code is architecture-specific and therefore
incompatible with future architectures, whereas PTX code is compiled to binary code at
load time by the device driver.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 231
Driver API
Here is the host code of the sample from Kernels written using the driver API:
int main()
{
int N = ...;
size_t size = N * sizeof(float);
// Initialize
cuInit(0);
// Create context
CUcontext cuContext;
cuCtxCreate(&cuContext, 0, cuDevice);
// Invoke kernel
int threadsPerBlock = 256;
int blocksPerGrid =
(N + threadsPerBlock - 1) / threadsPerBlock;
void* args[] = { &d_A, &d_B, &d_C, &N };
cuLaunchKernel(vecAdd,
blocksPerGrid, 1, 1, threadsPerBlock, 1, 1,
0, 0, args, 0);
...
}
Full code can be found in the vectorAddDrv CUDA sample.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 232
Driver API
H.1. Context
A CUDA context is analogous to a CPU process. All resources and actions performed
within the driver API are encapsulated inside a CUDA context, and the system
automatically cleans up these resources when the context is destroyed. Besides objects
such as modules and texture or surface references, each context has its own distinct
address space. As a result, CUdeviceptr values from different contexts reference
different memory locations.
A host thread may have only one device context current at a time. When a context
is created with cuCtxCreate(), it is made current to the calling host thread. CUDA
functions that operate in a context (most functions that do not involve device
enumeration or context management) will return CUDA_ERROR_INVALID_CONTEXT if a
valid context is not current to the thread.
Each host thread has a stack of current contexts. cuCtxCreate() pushes the new
context onto the top of the stack. cuCtxPopCurrent() may be called to detach the
context from the host thread. The context is then "floating" and may be pushed as the
current context for any host thread. cuCtxPopCurrent() also restores the previous
current context, if any.
A usage count is also maintained for each context. cuCtxCreate() creates a
context with a usage count of 1. cuCtxAttach() increments theusage count and
cuCtxDetach() decrements it. A context is destroyed when theusage count goes to 0
when calling cuCtxDetach() or cuCtxDestroy().
Usage count facilitates interoperability between third party authored code operating in
the same context. For example, if three libraries are loaded to use the same context, each
library would call cuCtxAttach() to increment the usage count and cuCtxDetach()
to decrement the usage count when the library is done using the context. For most
libraries, it is expected that the application will have created a context before loading
or initializing the library; that way, the application can create the context using its
own heuristics, and the library simply operates on the context handed to it. Libraries
that wish to create their own contexts - unbeknownst to their API clients who may or
may not have created contexts of their own - would use cuCtxPushCurrent() and
cuCtxPopCurrent() as illustrated in Figure 19.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 233
Driver API
Initialize
cuCtxCreate() context cuCtxPopCurrent()
Library Call
Use
cuCtxPushCurrent() context cuCtxPopCurrent()
H.2. Module
Modules are dynamically loadable packages of device code and data, akin to DLLs in
Windows, that are output by nvcc (see Compilation with NVCC). The names for all
symbols, including functions, global variables, and texture or surface references, are
maintained at module scope so that modules written by independent third parties may
interoperate in the same CUDA context.
This code sample loads a module and retrieves a handle to some kernel:
CUmodule cuModule;
cuModuleLoad(&cuModule, "myModule.ptx");
CUfunction myKernel;
cuModuleGetFunction(&myKernel, cuModule, "MyKernel");
This code sample compiles and loads a new module from PTX code and parses
compilation errors:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 234
Driver API
This code sample compiles, links, and loads a new module from multiple PTX codes and
parses link and compilation errors:
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 235
Driver API
Alignment requirements in device code for the built-in vector types are listed in Table
3. For all other basic types, the alignment requirement in device code matches the
alignment requirement in host code and can therefore be obtained using alignof().
The only exception is when the host compiler aligns double and long long (and
long on a 64-bit system) on a one-word boundary instead of a two-word boundary(for
example, using gcc's compilation flag -mno-align-double) since in device code these
types are always aligned on a two-word boundary.
CUdeviceptr is an integer, but represents a pointer, so its alignment requirement is
alignof(void*).
The following code sample uses a macro (ALIGN_UP()) to adjust the offset
of each parameter to meet its alignment requirement and another macro
(ADD_TO_PARAM_BUFFER()) to add each parameter to the parameter buffer passed to
the CU_LAUNCH_PARAM_BUFFER_POINTER option.
#define ALIGN_UP(offset, alignment) \
(offset) = ((offset) + (alignment) - 1) & ~((alignment) - 1)
char paramBuffer[1024];
size_t paramBufferSize = 0;
int i;
ADD_TO_PARAM_BUFFER(i, alignof(i));
float4 f4;
ADD_TO_PARAM_BUFFER(f4, 16); // float4's alignment is 16
char c;
ADD_TO_PARAM_BUFFER(c, alignof(c));
float f;
ADD_TO_PARAM_BUFFER(f, alignof(f));
CUdeviceptr devPtr;
ADD_TO_PARAM_BUFFER(devPtr, alignof(devPtr));
float2 f2;
ADD_TO_PARAM_BUFFER(f2, 8); // float2's alignment is 8
void* extra[] = {
CU_LAUNCH_PARAM_BUFFER_POINTER, paramBuffer,
CU_LAUNCH_PARAM_BUFFER_SIZE, ¶mBufferSize,
CU_LAUNCH_PARAM_END
};
cuLaunchKernel(cuFunction,
blockWidth, blockHeight, blockDepth,
gridWidth, gridHeight, gridDepth,
0, 0, 0, extra);
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 236
Driver API
is padded in device code with 12 bytes after field f since the alignment requirementfor
field f4 is 16.
typedef struct {
float f;
float4 f4;
} myStruct;
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 237
Appendix I.
CUDA ENVIRONMENT VARIABLES
Environment variables related to the Multi-Process Service are documented in the Multi-
Process Service section of the GPU Deployment and Management guide.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 238
CUDA Environment Variables
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 239
CUDA Environment Variables
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 240
Appendix J.
UNIFIED MEMORY PROGRAMMING
A processor refers to any independent execution unit with a dedicated MMU. This
includes both CPUs and GPUs of any type and architecture.
The underlying system manages data access and locality within a CUDA program
without need for explicit memory copy calls. This benefits GPU programming in two
primary ways:
‣ GPU programming is simplified by unifying memory spaces coherently across all
GPUs and CPUs in the system and by providing tighter and more straightforward
language integration for CUDA programmers.
‣ Data access speed is maximized by transparently migrating data towards the
processor using it.
In simple terms, Unified Memory eliminates the need for explicit data movement via the
cudaMemcpy*() routines without the performance penalty incurred by placing all data
into zero-copy memory. Data movement, of course, still takes place, so a program’s run
time typically does not decrease; Unified Memory instead enables the writing of simpler
and more maintainable code.
Unified Memory offers a “single-pointer-to-data” model that is conceptually similar to
CUDA’s zero-copy memory. One key difference between the two is that with zero-copy
allocations the physical location of memory is pinned in CPU system memory such that
a program may have fast or slow access to it depending on where it is being accessed
from. Unified Memory, on the other hand, decouples memory and execution spaces so
that all data accesses are fast.
The term Unified Memory describes a system that provides memory management
services to a wide range of programs, from those targeting the Runtime API down to
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 241
Unified Memory Programming
those using the Virtual ISA (PTX). Part of this system defines the managed memory
space that opts in to Unified Memory services.
Managed memory is interoperable and interchangeable with device-specific allocations,
such as those created using the cudaMalloc() routine. All CUDA operations that are
valid on device memory are also valid on managed memory; the primary difference is
that the host portion of a program is able to reference and access the memory as well.
On supporting platforms with devices of compute capability 6.x, Unified Memory will
enable applications to allocate and share data using the default system allocator. This
allows the GPU to access the entire system virtual memory without using a special
allocator.
The following code examples illustrate how the use of managed memory can change the
way in which host code is written. First, a simple program written without the benefit of
unified memory:
global void AplusB( int *ret, int a, int b) {
ret[threadIdx.x] = a + b + threadIdx.x;
}
int main() {
int *ret;
cudaMalloc(&ret, 1000 * sizeof(int));
AplusB<<< 1, 1000 >>>(ret, 10, 100);
int *host_ret = (int *)malloc(1000 * sizeof(int));
cudaMemcpy(host_ret, ret, 1000 * sizeof(int), cudaMemcpyDefault);
for(int i=0; i<1000; i++)
printf("%d: A+B = %d\n", i, host_ret[i]);
free(host_ret);
cudaFree(ret);
return 0;
}
This first example combines two numbers together on the GPU with a per-thread ID and
returns the values in an array. Without managed memory, both host- and device-side
storage for the return values is required (host_ret and ret in the example), as is an
explicit copy between the two using cudaMemcpy().
Compare this with the Unified Memory version of the program, which allows direct
access of GPU data from the host. Notice the cudaMallocManaged() routine, which
returns a pointer valid from both host and device code. This allows ret to be used
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 242
Unified Memory Programming
without a separate host_ret copy, greatly simplifying and reducing the size ofthe
program.
global void AplusB(int *ret, int a, int b) {
ret[threadIdx.x] = a + b + threadIdx.x;
}
int main() {
int *ret;
cudaMallocManaged(&ret, 1000 * sizeof(int));
AplusB<<< 1, 1000 >>>(ret, 10, 100);
cudaDeviceSynchronize();
for(int i=0; i<1000; i++)
printf("%d: A+B = %d\n", i, ret[i]);
cudaFree(ret);
return 0;
}
Finally, language integration allows direct reference of a GPU-declared managed
variable and simplifies a program further when global variables are used.
device managed int ret[1000];
global void AplusB(int a, int b) {
ret[threadIdx.x] = a + b + threadIdx.x;
}
int main() {
AplusB<<< 1, 1000 >>>(10, 100);
cudaDeviceSynchronize();
for(int i=0; i<1000; i++)
printf("%d: A+B = %d\n", i, ret[i]);
return 0;
}
Note the absence of explicit cudaMemcpy() commands and the fact that the return array
ret is visible on both CPU and GPU.
It is worth a comment on the synchronization between host and device. Notice how in
the non-managed example, the synchronous cudaMemcpy() routine is used both to
synchronize the kernel (that is, to wait for it to finish running), and to transfer the data
to the host. The Unified Memory examples do not call cudaMemcpy() and so requirean
explicit cudaDeviceSynchronize() before the host program can safely use theoutput
from the GPU.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 243
Unified Memory Programming
any processor regardless of locality. Note that maintaining coherence is the primary
requirement, ahead of performance; within the constraints of the host operating system,
the system is permitted to either fail accesses or move data in order to maintain global
coherence between processors.
GPU architectures of compute capability lower than 6.x do not support fine-grained
movement of the managed data to GPU on-demand. Whenever a GPU kernel is
launched all managed memory generally has to be transfered to GPU memory to avoid
faulting on memory access. With compute capability 6.x a new GPU page faulting
mechanism is introduced that provides more seamless Unified Memory functionality.
Combined with the system-wide virtual address space, page faulting provides several
benefits. First, page faulting means that the CUDA system software doesn’t need to
synchronize all managed memory allocations to the GPU before each kernel launch.
If a kernel running on the GPU accesses a page that is not resident in its memory, it
faults, allowing the page to be automatically migrated to the GPU memory on-demand.
Alternatively, the page may be mapped into the GPU address space for access over
the PCIe or NVLink interconnects (mapping on access can sometimes be faster than
migration). Note that Unified Memory is system-wide: GPUs (and CPUs) can fault on
and migrate memory pages either from CPU memory or from the memory of other
GPUs in the system.
4. Multi-GPU Support
For devices of compute capability lower than 6.x managed memory allocation behaves
identically to unmanaged memory allocated using cudaMalloc(): the current active
device is the home for the physical allocation, and all other GPUs receive peer mappings
to the memory. This means that other GPUs in the system will access the memory at
reduced bandwidth over the PCIe bus. Note that if peer mappings are not supported
between the GPUs in the system, then the managed memory pages are placed in CPU
system memory (“zero-copy” memory), and all GPUs will experience PCIe bandwidth
restrictions. See Managed Memory with Multi-GPU Programs on pre-6.x Architectures
for details.
Managed allocations on systems with devices of compute capability 6.x are visible to all
GPUs and can migrate to any processor on-demand.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 244
Unified Memory Programming
2. Programming Model
1. Managed Memory Opt In
Most platforms require a program to opt in to automatic data management by either
annotating a device variable with the managed keyword (see the
Language Integration section) or by using a new cudaMallocManaged() call to allocate
data.
Devices of compute capability lower than 6.x must always allocate managed memory on
the heap, either with an allocator or by declaring global storage. It is not possible either
to associate previously allocated memory with Unified Memory, or to have the Unified
Memory system manage a CPU or a GPU stack pointer.
Starting with CUDA 8.0 and on supporting systems with devices of compute capability
6.x, memory allocated with the default OS allocator (e.g. malloc or new) can beaccessed
from both GPU code and CPU code using the same pointer. On these systems, Unified
Memory is the default: there is no need to use a special allocator or the creation of a
specially managed memory pool.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 245
Unified Memory Programming
of the Unified Memory programming model (see Coherency and Concurrency). Below is
a simple example, showing the use of cudaMallocManaged():
global void printme(char *str) {
printf(str);
}
int main() {
// Allocate 100 bytes of memory, accessible to both Host and Device code
char *s;
cudaMallocManaged(&s, 100);
// Note direct Host-code use of "s"
strncpy(s, "Hello Unified Memory\n", 99);
// Here we pass "s" to a kernel without explicitly copying
printme<<< 1, 1 >>>(s);
cudaDeviceSynchronize();
// Free as for normal CUDA allocations
cudaFree(s);
return 0;
}
A program’s behavior is functionally unchanged when cudaMalloc() is replaced
with cudaMallocManaged(); however, the program should go on to eliminate explicit
memory copies and take advantage of automatic migration. Additionally, dual pointers
(one to host and one to device memory) can be eliminated.
Device code is not able to call cudaMallocManaged(). All managed memory must be
allocated from the host or at global scope (see the next section). Allocations on the device
heap using malloc() in a kernel will not be created in the managed memory space, and
so will not be accessible to CPU code.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 246
Unified Memory Programming
cudaDeviceSynchronize();
return 0;
}
In example above, the GPU program kernel is still active when the CPU touches y.
(Note how it occurs before cudaDeviceSynchronize().) The code runs successfully
on devices of compute capability 6.x due to the GPU page faulting capability which
lifts all restrictions on simultaneous access. However, such memory access is invalid on
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 247
Unified Memory Programming
pre-6.x architectures even though the CPU is accessing different data than the GPU. The
program must explicitly synchronize with the GPU before accessing y:
device managed int x, y=2;
global void kernel() {
x = 10;
}
int main() {
kernel<<< 1, 1 >>>();
cudaDeviceSynchronize();
y = 20; // Success on GPUs not supporing concurrent access
return 0;
}
As this example shows, on systems with pre-6.x GPU architectures, a CPU thread may
not access any managed data in between performing a kernel launch and a subsequent
synchronization call, regardless of whether the GPU kernel actually touches that same
data (or any managed data at all). The mere potential for concurrent CPU and GPU
access is sufficient for a process-level exception to be raised.
Note that if memory is dynamically allocated with cudaMallocManaged() or
cuMemAllocManaged() while the GPU is active, the behavior of the memory is
unspecified until additional work is launched or the GPU is synchronized. Attempting
to access the memory on the CPU during this time may or may not cause a segmentation
fault. This does not apply to memory allocated using the flag cudaMemAttachHost or
CU_MEM_ATTACH_HOST.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 248
Unified Memory Programming
‣ It is always permitted for the CPU to access non-managed zero-copy data while the
GPU is active.
‣ The GPU is considered active when it is running any kernel, even if that kernel does
not make use of managed data. If a kernel might use data, then access is forbidden,
unless device property concurrentManagedAccess is 1.
‣ There are no constraints on concurrent inter-GPU access of managed memory, other
than those that apply to multi-GPU access of non-managed memory.
‣ There are no constraints on concurrent GPU kernels accessing managed data.
Note how the last point allows for races between GPU kernels, as is currently the case
for non-managed GPU memory. As mentioned previously, managed memory functions
identically to non-managed memory from the perspective of the GPU. The following
code example illustrates these points:
int main() {
cudaStream_t stream1, stream2;
cudaStreamCreate(&stream1);
cudaStreamCreate(&stream2);
int *non_managed, *managed, *also_managed;
cudaMallocHost(&non_managed, 4); // Non-managed, CPU-accessible memory
cudaMallocManaged(&managed, 4);
cudaMallocManaged(&also_managed, 4);
// Point 1: CPU can access non-managed data.
kernel<<< 1, 1, 0, stream1 >>>(managed);
*non_managed = 1;
// Point 2: CPU cannot access any managed data while GPU is busy,
// unless concurrentManagedAccess = 1
// Note we have not yet synchronized, so "kernel" is still active.
*also_managed = 2; // Will issue segmentation fault
// Point 3: Concurrent GPU kernels can access the same data.
kernel<<< 1, 1, 0, stream2 >>>(managed);
// Point 4: Multi-GPU concurrent access is also permitted.
cudaSetDevice(1);
kernel<<< 1, 1 >>>(managed);
return 0;
}
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 249
Unified Memory Programming
into a specified stream or not. This enables opportunities for concurrency based on
program-specific data access patterns. The function to control this behaviour is:
cudaError_t cudaStreamAttachMemAsync(cudaStream_t stream,
void *ptr,
size_t length=0,
unsigned int flags=0);
The cudaStreamAttachMemAsync() function associates length bytes of memory
starting from ptr with the specified stream. (Currently, length must always be
0 to indicate that the entire region should be attached.) Because of this association,
the Unified Memory system allows CPU access to this memory region so long as all
operations in stream have completed, regardless of whether other streams are active.In
effect, this constrains exclusive ownership of the managed memory region by an active
GPU to per-stream activity instead of whole-GPU activity.
Most importantly, if an allocation is not associated with a specific stream, it is visible
to all running kernels regardless of their stream. This is the default visibility for a
cudaMallocManaged() allocation or a managed variable; hence, the simple-case
rule that the CPU may not touch the data while any kernel is running.
By associating an allocation with a specific stream, the program makes a guarantee
that only kernels launched into that stream will touch that data. No error checking is
performed by the Unified Memory system: it is the programmer’s responsibility to
ensure that guarantee is honored.
In addition to allowing greater concurrency, the use of cudaStreamAttachMemAsync()
can (and typically does) enable data transfer optimizations within the Unified Memory
system that may affect latencies and other overhead.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 250
Unified Memory Programming
Note that associating a variable with a stream does not change the associating of any
other variable. E.g. associating x with stream1 does not ensure that only x is accessed
by kernels launched in stream1, thus an error is caused by this code:
device managed int x, y=2;
global void kernel() {
x = 10;
}
int main() {
cudaStream_t stream1;
cudaStreamCreate(&stream1);
cudaStreamAttachMemAsync(stream1, &x);// Associate “x” with stream1.
cudaDeviceSynchronize(); // Wait for “x” attachment to occur.
kernel<<< 1, 1, 0, stream1 >>>(); // Note: Launches into stream1.
y = 20; // ERROR: “y” is still associated
globally
// with all streams by default
return 0;
}
Note how the access to y will cause an error because, even though x has been associated
with a stream, we have told the system nothing about who can see y. The system
therefore conservatively assumes that kernel might access it and prevents the CPU
from doing so.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 251
Unified Memory Programming
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 252
Unified Memory Programming
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 253
Unified Memory Programming
CUDA’s language integration: variables annotated with the managed keyword can
be referenced directly from both host and device code.
The following example, seen earlier in Simplifying GPU Programming, illustrates a
simple use of managed global declarations:
// Managed variable declaration is an extra annotation with device
device managed int x;
global void kernel() {
// Reference "x" directly - it's a normal variable on the GPU.
printf( "GPU sees: x = %d\n" , x);
}
int main() {
// Set "x" from Host code. Note it's just a normal variable on the CPU.
x = 1234;
// Launch a kernel which uses "x" from the GPU.
kernel<<< 1, 1 >>>();
cudaDeviceSynchronize();
return 0;
}
The capability available with managed variables is that the symbol is available in
both device code and in host code without the need to dereference a pointer, and the
data is shared by all. This makes it particularly easy to exchange data between host and
device programs without the need for explicit allocations or copying.
Semantically, the behavior of managed variables is identical to that of
storage allocated via cudaMallocManaged(). Data is hosted in physical GPU
storage and is visible to all GPUs in the system as well as the CPU. Stream
visibility defaults to cudaMemAttachGlobal, but may be constrained using
cudaStreamAttachMemAsync().
A valid CUDA context is necessary for the correct operation of managed variables.
Accessing managed variables can trigger CUDA context creation if a context for
the current device hasn’t already been created. In the example above, accessing x before
the kernel launch triggers context creation on device 0. In the absence of that access, the
kernel launch would have triggered context creation.
C++ objects declared as managed are subject to certain specific constraints,
particularly where static initializers are concerned. Please refer to C/C++ Language
Support in the CUDA C Programming Guide for a list of these constraints.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 254
Unified Memory Programming
2. Pointer Attributes
To determine if a given pointer refers to managed memory, a program can call
cudaPointerGetAttributes() and check the value of the isManaged attribute. This
attribute is set to 1 if the pointer refers to managed memory and to 0 if not.
5. Advanced Topics
1.Managed Memory with Multi-GPU Programs on pre-6.x
Architectures
On systems with devices of compute capabilities lower than 6.x managed allocations
are automatically visible to all GPUs in a system via the peer-to-peer capabilities of the
GPUs.
On Linux the managed memory is allocated in GPU memory as long as all GPUs that
are actively being used by a program have the peer-to-peer support. If at any time the
application starts using a GPU that doesn’t have peer-to-peer support with any of the
other GPUs that have managed allocations on them, then the driver will migrate all
managed allocations to system memory.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 255
Unified Memory Programming
On Windows if peer mappings are not available (for example, between GPUs of
different architectures), then the system will automatically fall back to using zero-copy
memory, regardless of whether both GPUs are actually used by a program. If only one
GPU is actually going to be used, it is necessary to set the CUDA_VISIBLE_DEVICES
environment variable before launching the program. This constrains which GPUs are
visible and allows managed memory to be allocated in GPU memory.
Alternatively, on Windows users can also set CUDA_MANAGED_FORCE_DEVICE_ALLOC to
a non-zero value to force the driver to always use device memory for physical storage.
When this environment variable is set to a non-zero value, all devices used in that
process that support managed memory have to be peer-to-peer compatible with each
other. The error ::cudaErrorInvalidDevice will be returned if a device that supports
managed memory is used and it is not peer-to-peer compatible with any of the other
managed memory supporting devices that were previously used in that process, even
if ::cudaDeviceReset has been called on those devices. These environment variables are
described in Appendix CUDA Environment Variables. Note that starting from CUDA
8.0 CUDA_MANAGED_FORCE_DEVICE_ALLOC has no effect on Linux operating systems.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 256
Unified Memory Programming
To achieve the same level of performance as what's possible without using Unified
Memory, the application has to guide the Unified Memory driver subsystem into
avoiding the aforementioned pitfalls. It is worthy to note that the Unified Memory
driver subsystem can detect common data access patterns and achieve some of these
objectives automatically without application participation. But when the data access
patterns are non-obvious, explicit guidance from the application is crucial. CUDA
8.0 introduces useful APIs for providing the runtime with memory usage hints
(cudaMemAdvise()) and for explicit prefetching (cudaMemPrefetchAsync()). These
tools allow the same capabilities as explicit memory copy and pinning APIs without
reverting to the limitations of explicit GPU memory allocation.
where the memory region specified by devPtr pointer and count number of bytes, with
ptr rounded down to the nearest page boundary and count rounded up to the nearest
page boundary, is migrated to the dstDevice by enqueueing a migration operation in
stream. Passing in cudaCpuDeviceId for dstDevice will cause data to be migrated to
CPU memory.
Consider a simple code example below:
void foo(cudaStream_t s) {
char *data;
cudaMallocManaged(&data, N);
init_data(data, N); // execute on CPU
cudaMemPrefetchAsync(data, N, myGpuId, s); // prefetch to GPU
mykernel<<<..., s>>>(data, N, 1, compare); // execute on GPU
cudaMemPrefetchAsync(data, N, cudaCpuDeviceId, s); // prefetch to CPU
cudaStreamSynchronize(s);
use_data(data, N);
cudaFree(data);
}
Without performance hints the kernel mykernel will fault on first access to data
which creates additional overhead of the fault processing and generally slows down
the application. By prefetching data in advance it is possible to avoid page faultsand
achieve better performance.
This API follows stream ordering semantics, i.e. the migration does not begin until all
prior operations in the stream have completed, and any subsequent operation in the
stream does not begin until the migration has completed.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 257
Unified Memory Programming
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 258
Unified Memory Programming
avoided. On the other hand, if the data is not in its preferred location or if a direct
mapping cannot be established, then it will be migrated to the processor accessing
it. It is important to note that setting the preferred location does not prevent data
prefetching done using cudaMemPrefetchAsync.
‣ cudaMemAdviseSetAccessedBy: This advice implies that the data will be accessed
by device. This does not cause data migration and has no impact on the location of
the data per se. Instead, it causes the data to always be mapped in the specified
processor’s page tables, as long as the location of the data permits a mapping to
be established. If the data gets migrated for any reason, the mappings are updated
accordingly. This advice is useful in scenarios where data locality is not important,
but avoiding faults is. Consider for example a system containing multiple GPUs
with peer-to-peer access enabled, where the data located on one GPU is occasionally
accessed by other GPUs. In such scenarios, migrating data over to the other GPUs is
not as important because the accesses are infrequent and the overhead of migration
may be too high. But preventing faults can still help improve performance, and
so having a mapping set up in advance is useful. Note that on CPU access of this
data, the data may be migrated to CPU memory because the CPU cannot access
GPU memory directly. Any GPU that had the cudaMemAdviceSetAccessedBy flag
set for this data will now have its mapping updated to point to the page in CPU
memory.
Each advice can be also unset by using one of the following values:
cudaMemAdviseUnsetReadMostly, cudaMemAdviseUnsetPreferredLocation and
cudaMemAdviseUnsetAccessedBy.
cudaMemRangeGetAttribute(void *data,
size_t dataSize,
enum cudaMemRangeAttribute attribute,
const void *devPtr,
size_t count);
This function queries an attribute of the memory range starting at devPtr with a
size of count bytes. The memory range must refer to managed memory allocatedvia
cudaMallocManaged or declared via managed variables. It is possible to query the
following attributes:
‣ cudaMemRangeAttributeReadMostly: the result returned will be 1 if all pages in
the given memory range have read-duplication enabled, or 0 otherwise.
‣ cudaMemRangeAttributePreferredLocation: the result returned will
be a GPU device id or cudaCpuDeviceId if all pages in the memory range
have the corresponding processor as their preferred location, otherwise
cudaInvalidDeviceId will be returned. An application can use this query API to
make decision about staging data through CPU or GPU depending on the preferred
location attribute of the managed pointer. Note that the actual location of the pages
in the memory range at the time of the query may be different from the preferred
location.
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 259
Unified Memory Programming
www.nvidia.com
CUDA C Programming Guide PG-02829-001_v8.0 | 260
Notice
ALL NVIDIA DESIGN SPECIFICATIONS, REFERENCE BOARDS, FILES, DRAWINGS,
DIAGNOSTICS, LISTS, AND OTHER DOCUMENTS (TOGETHER AND SEPARATELY,
"MATERIALS") ARE BEING PROVIDED "AS IS." NVIDIA MAKES NO WARRANTIES,
EXPRESSED, IMPLIED, STATUTORY, OR OTHERWISE WITH RESPECT TO THE
MATERIALS, AND EXPRESSLY DISCLAIMS ALL IMPLIED WARRANTIES OF
NONINFRINGEMENT, MERCHANTABILITY, AND FITNESS FOR A PARTICULAR
PURPOSE.
Information furnished is believed to be accurate and reliable. However, NVIDIA
Corporation assumes no responsibility for the consequences of use of such
information or for any infringement of patents or other rights of third parties
that may result from its use. No license is granted by implication of otherwise
under any patent rights of NVIDIA Corporation. Specifications mentioned in this
publication are subject to change without notice. This publication supersedes and
replaces all other information previously supplied. NVIDIA Corporation products
are not authorized as critical components in life support devices or systems
without express written approval of NVIDIA Corporation.
Trademarks
NVIDIA and the NVIDIA logo are trademarks or registered trademarks of NVIDIA
Corporation in the U.S. and other countries. Other company and product names
may be trademarks of the respective companies with which they are associated.
Copyright
© 2007-2017 NVIDIA Corporation. All rights reserved.
www.nvidia.com