Memory access pattern

From HandWiki

In computing, a memory access pattern or IO access pattern is the pattern with which a system or program reads and writes memory on secondary storage. These patterns differ in the level of locality of reference and drastically affect cache performance,[1] and also have implications for the approach to parallelism[2][3] and distribution of workload in shared memory systems.[4] Further, cache coherency issues can affect multiprocessor performance,[5] which means that certain memory access patterns place a ceiling on parallelism (which manycore approaches seek to break).[6] Computer memory is usually described as "random access", but traversals by software will still exhibit patterns that can be exploited for efficiency. Various tools exist to help system designers[7] and programmers understand, analyse and improve the memory access pattern, including VTune and Vectorization Advisor,[8][9][10][11][12] including tools to address GPU memory access patterns[13]

Memory access patterns also have implications for security,[14][15] which motivates some to try and disguise a program's activity for privacy reasons.[16][17]

Examples

Sequential and Linear patterns are incorrectly drawn as counterparts to each other by some publications; while real-world workloads contain almost innumerable patterns.[18]

Sequential

Main page: Sequential access

The simplest extreme is the sequential access pattern, where data is read, processed, and written out with straightforward incremented/decremented addressing. These access patterns are highly amenable to prefetching.

Strided

Strided or simple 2D, 3D access patterns (e.g., stepping through multi-dimensional arrays) are similarly easy to predict, and are found in implementations of linear algebra algorithms and image processing. Loop tiling is an effective approach.[19] Some systems with DMA provided a strided mode for transferring data between subtile of larger 2D arrays and scratchpad memory.[20]

Linear

A linear access pattern is closely related to "strided", where a memory address may be computed from a linear combination of some index. Stepping through indices sequentially with a linear pattern yields strided access. A linear access pattern for writes (with any access pattern for non-overlapping reads) may guarantee that an algorithm can be parallelised, which is exploited in systems supporting compute kernels.

Nearest neighbor

Nearest neighbor memory access patterns appear in simulation, and are related to sequential or strided patterns. An algorithm may traverse a data structure using information from the nearest neighbors of a data element (in one or more dimensions) to perform a calculation. These are common in physics simulations operating on grids.[21] Nearest neighbor can also refer to inter-node communication in a cluster; physics simulations which rely on such local access patterns can be parallelized with the data partitioned into cluster nodes, with purely nearest-neighbor communication between them, which may have advantages for latency and communication bandwidth. This use case maps well onto torus network topology.[22]

2D spatially coherent

In 3D rendering, access patterns for texture mapping and rasterization of small primitives (with arbitrary distortions of complex surfaces) are far from linear, but can still exhibit spatial locality (e.g., in screen space or texture space). This can be turned into good memory locality via some combination of morton order[23] and tiling for texture maps and frame buffer data (mapping spatial regions onto cache lines), or by sorting primitives via tile based deferred rendering.[24] It can also be advantageous to store matrices in morton order in linear algebra libraries.[25]

Scatter

A scatter memory access pattern combines sequential reads with indexed/random addressing for writes.[26] Compared to gather, It may place less load on a cache hierarchy since a processing element may dispatch writes in a "fire and forget" manner (bypassing a cache altogether), whilst using predictable prefetching (or even DMA) for its source data.

However, it may be harder to parallelise since there is no guarantee the writes do not interact,[27] and many systems are still designed assuming that a hardware cache will coalesce many small writes into larger ones.

In the past, forward texture mapping attempted to handle the randomness with "writes", whilst sequentially reading source texture information.

The PlayStation 2 console used conventional inverse texture mapping, but handled any scatter/gather processing "on-chip" using EDRAM, whilst 3D model (and a lot of texture data) from main memory was fed sequentially by DMA. This is why it lacked support for indexed primitives, and sometimes needed to manage textures "up front" in the display list.

Gather

In a gather memory access pattern, reads are randomly addressed or indexed, whilst the writes are sequential (or linear).[26] An example is found in inverse texture mapping, where data can be written out linearly across scan lines, whilst random access texture addresses are calculated per pixel.

Compared to scatter, the disadvantage is that caching (and bypassing latencies) is now essential for efficient reads of small elements, however it is easier to parallelise since the writes are guaranteed to not overlap. As such the gather approach is more common for gpgpu programming,[27] where the massive threading (enabled by parallelism) is used to hide read latencies.[27]

Combined gather and scatter

An algorithm may gather data from one source, perform some computation in local or on chip memory, and scatter results elsewhere. This is essentially the full operation of a GPU pipeline when performing 3D rendering- gathering indexed vertices and textures, and scattering shaded pixels in screen space. Rasterization of opaque primitives using a depth buffer is "commutative", allowing reordering, which facilitates parallel execution. In the general case synchronisation primitives would be needed.

Random

Main page: Random access

At the opposite extreme is a truly random memory access pattern. A few multiprocessor systems are specialised to deal with these.[28] The PGAS approach may help by sorting operations by data on the fly (useful when the problem *is* figuring out the locality of unsorted data).[21] Data structures which rely heavily on pointer chasing can often produce poor locality of reference, although sorting can sometimes help. Given a truly random memory access pattern, it may be possible to break it down (including scatter or gather stages, or other intermediate sorting) which may improve the locality overall; this is often a prerequisite for parallelizing.

Approaches

Data-oriented design

Data-oriented design is an approach intended to maximise the locality of reference, by organising data according to how it is traversed in various stages of a program, contrasting with the more common object oriented approach (i.e., organising such that data layout explicitly mirrors the access pattern).[29]

Contrast with locality of reference

Locality of reference refers to a property exhibited by memory access patterns. A programmer will change the memory access pattern (by reworking algorithms) to improve the locality of reference,[30] and/or to increase potential for parallelism.[26] A programmer or system designer may create frameworks or abstractions (e.g., C++ templates or higher-order functions) that encapsulate a specific memory access pattern.[31][32]

Different considerations for memory access patterns appear in parallelism beyond locality of reference, namely the separation of reads and writes. E.g.: even if the reads and writes are "perfectly" local, it can be impossible to parallelise due to dependencies; separating the reads and writes into separate areas yields a different memory access pattern, maybe initially appear worse in pure locality terms, but desirable to leverage modern parallel hardware.[26]

Locality of reference may also refer to individual variables (e.g., the ability of a compiler to cache them in registers), whilst the term memory access pattern only refers to data held in an indexable memory (especially main memory).

See also

References

  1. "data oriented design". http://www.dice.se/wp-content/uploads/2014/12/Introduction_to_Data-Oriented_Design.pdf. 
  2. Jang, Byunghyun; Schaa, Dana; Mistry, Perhaad; Kaeli, David (2010-05-27). "Exploiting Memory Access Patterns to Improve Memory Performance in Data-Parallel Architectures". IEEE Transactions on Parallel and Distributed Systems (New York: IEEE) 22 (1): 105–118. doi:10.1109/TPDS.2010.107. NLM unique id 101212014. ISSN 1045-9219. 
  3. Jeffers, James; Reinders, James; Sodani, Avinash (2016-05-31). xeon phi optimization. ISBN 9780128091951. https://books.google.com/books?id=DDpUCwAAQBAJ&q=scatter+memory+access+pattern&pg=PA231. 
  4. "Analysis of Energy and Performance of Code Transformations for PGAS-based Data Access Patterns". http://nic.uoregon.edu/pgas14/papers/pgas14_submission_17.pdf. 
  5. "enhancing cache coherent architectures with memory access patterns for embedded many-core systems". http://www.cc.gatech.edu/~bader/papers/EnhancingCache-SoC12.pdf. 
  6. "intel terascale". https://cseweb.ucsd.edu/classes/fa12/cse291-c/talks/SCC-80-core-cern.pdf. 
  7. analysis of memory access patterns. WWC '98. 29 November 1998. p. 105. ISBN 9780769504506. http://dl.acm.org/citation.cfm?id=838115. 
  8. "QUAD a memory access pattern analyser". http://ce-publications.et.tudelft.nl/publications/207_quad__a_memory_access_pattern_analyser.pdf. 
  9. "Dymaxion: Optimizing Memory Access Patterns for Heterogeneous Systems". http://www.cs.virginia.edu/~skadron/Papers/sc11_dymaxion_dist.pdf. 
  10. evaluation of a memory access classification scheme for pointer intensive and numeric programs. 1996. 
  11. Matsubara, Yuki; Sato, Yukinori (2014). "Online Memory Access Pattern Analysis on an Application Profiling Tool". 2014 Second International Symposium on Computing and Networking. pp. 602–604. doi:10.1109/CANDAR.2014.86. ISBN 978-1-4799-4152-0. 
  12. "Putting Your Data and Code in Order: Data and layout". https://software.intel.com/en-us/articles/putting-your-data-and-code-in-order-data-and-layout-part-2. 
  13. CuMAPz: a tool to analyze memory access patterns in CUDA. Dac '11. 5 June 2011. pp. 128–133. doi:10.1145/2024724.2024754. ISBN 9781450306362. http://dl.acm.org/citation.cfm?id=2024754. 
  14. "Memory Access Pattern Protection for Resource-constrained Devices". https://www.cardis.org/proceedings/cardis_2012/CARDIS2012_14.pdf. 
  15. "understanding cache attacks". https://www.rocq.inria.fr/secret/Anne.Canteaut/Publications/RR-5881.pdf. 
  16. "protecting data in the cloud". https://news.mit.edu/2013/protecting-data-in-the-cloud-0702. 
  17. "boosting-cloud-security-with----oblivious-ram". 24 September 2013. http://www.information-age.com/technology/security/123457364/boosting-cloud-security-with----oblivious-ram---. proposed RAM design avoiding memory-access-pattern vulnerabilities
  18. "Storage Performance Benchmarking Guidelines - Part I: Workload Design". http://www.snia.org/sites/default/files/PerformanceBenchmarking.Nov2010.pdf. "In practice, IO access patterns are as numerous as the stars" 
  19. "optimizing for tiling and data locality". http://www.cs.utexas.edu/users/mckinley/papers/par-mem-ics-92.pdf. paper covers loop tiling and implication for parallel code
  20. "Optimal 2D Data Partitioning for DMA Transfers on MPSoCs". http://www-verimag.imag.fr/~maler/Papers/dma2dim.pdf. 
  21. 21.0 21.1 "partitioned global address space programming". https://www.youtube.com/watch?v=NU4VfjISk2M. covers cases where PGAS is a win, where data may not be already sorted, e.g., dealing with complex graphs - see "science across the irregularity spectrum".
  22. "Quantifying Locality In The Memory Access Patterns of HPC Applications". http://www.sdsc.edu/~allans/sc05_locality.pdf. mentions nearest neighbor access patterns in clusters
  23. "The Design and Analysis of a Cache Architecture for Texture Mapping". https://www.cs.cmu.edu/afs/cs/academic/class/15869-f11/www/readings/hakura97_texcaching.pdf. see morton order,texture access pattern
  24. "morton order to accelerate texturing". http://john.cs.olemiss.edu/~rhodes/papers/Nocentino10.pdf. 
  25. "Morton-order Matrices Deserve Compilers' Support Technical Report 533". http://www.cs.indiana.edu/pub/techreports/TR533.pdf. discusses the importance of morton order for matrices
  26. 26.0 26.1 26.2 26.3 "gpgpu scatter vs gather". http://http.developer.nvidia.com/GPUGems2/gpugems2_chapter31.html. 
  27. 27.0 27.1 27.2 GPU gems. 2011-01-13. ISBN 9780123849892. https://books.google.com/books?id=lGMzmbUhpiAC&q=scatter+memory+access+pattern&pg=PA51. deals with "scatter memory access patterns" and "gather memory access patterns" in the text
  28. "Cray and HPCC: Benchmark Developments and Results from the Past Year". https://cug.org/5-publications/proceedings_attendee_lists/2005CD/S05_Proceedings/pages/Authors/Wichmann/Wichmann_paper.pdf. see global random access results for Cray X1. vector architecture for hiding latencies, not so sensitive to cache coherency
  29. "data oriented design". http://www.dice.se/wp-content/uploads/2014/12/Introduction_to_Data-Oriented_Design.pdf. 
  30. "optimize-data-structures-and-memory-access-patterns-to-improve-data-locality". https://software.intel.com/en-us/articles/optimize-data-structures-and-memory-access-patterns-to-improve-data-locality. 
  31. "Template-based Memory Access Engine for Accelerators in SoCs". http://www.cs.utah.edu/~zfang/asp_dac11.pdf. 
  32. "Multi-Target Vectorization With MTPS C++ Generic Library". http://www.metz.supelec.fr/metz/recherche/publis_pdf/Supelec753.pdf. a C++ template library for producing optimised memory access patterns