1/57
Práce s pamětí
Vláknové programování
část IX
Lukáš Hejmánek, Petr Holub
{xhejtman,hopet}@ics.muni.cz
Laboratoř pokročilých síťových technologií
PV192
2014–04–29
2/57
Práce s pamětí
Přehled přednášky
Práce s pamětí
Paměť RAM
Cache
Rychlosti čtení paměti
Zápis do paměti
Optimalizace
3/57
Práce s pamětí
Hierarchie pamětí
• Oproti původnímu von Neumannovu konceptu existuje více druhů
pamětí
• Různé druhy paměti se liší rychlostí, kapacitou, cenou
• Nelze mít veškerou paměť pouze rychlou a s dostatečnou kapacitou
• Čím má paměť vetší kapacitu, tím je obvykle pomalejší
• Ať už z ekonomických důvodů
• Nebo z fyzikálních omezení
• Pevný disk ⇒ Flash paměť ⇒ Paměť DRAM ⇒ Cache ⇒ Registry
4/57
Práce s pamětí
Hierarchie pamětí
Odkud Cyklů
Registr ≤ 1
L1d ∼ 3
L2 ∼ 14
DRAM ∼ 240
Údaje pro Pentium M
5/57
Práce s pamětí
Architektura paměti
• Tradiční schéma pamětí RAM a CPU
Northbridge
Southbridge
CPUCPU
RAM
PCI-E SATA
USB
FSB
1 2
• Problémy takového přístupu
• Všechny procesory sdílí jednu sběrnici (quadpumped 800–1333MHz)
• Komunikace s pamětí všech procesorů probíhá přes jeden Northbridge
• Paměť má jen jeden přístupový port
• Komunikace mezi procesory a zařízeními probíhá přes jeden
Northbridge
6/57
Práce s pamětí
Architektura paměti
• NUMA schéma pamětí RAM a CPU
Southbridge
CPUCPU
RAM
PCI-E SATA
USB
1 2
CPU3 CPU4
RAM
RAM
RAM
• Výhody takového přístupu
• Komunikace s pamětí nepřetěžuje sběrnici s Northbridge
• Více připojitelné paměti
• Nevýhody
• Neuniformní přístup do paměti
• Vysoký NUMA faktor je problém
• Ne každý procesor může přímo komunikovat se zařízením
7/57
Práce s pamětí
Principy DRAM
• Buňka paměti DRAM
DL
AL
M C
• Kondenzátor C svým nabitím uchovává bitovou informaci 1 nebo 0
• Čtení/zápis bitu je řízen AL (adresní linka)
• Čtení/zápis bitu je provedeno na DL (datová linka)
8/57
Práce s pamětí
Principy DRAM
• Čtení a zápis do/z kondenzátoru není instantní
0
10
20
30
40
50
60
70
80
90
100
1RC 2RC 3RC 4RC 5RC 6RC 7RC 8RC 9RC
Charge DischargePercentageCharge
9/57
Práce s pamětí
Principy DRAM
• AL a DL samostatně ke každému bitu by vyžadovalo extrémní
množství propojů (Kolik je 32 Gb?)
• Buňky jsou v matici, adresujeme sloupce, řádky
a0
a1
a2
a3
Data
Column Address Selection
RowAddressSelection
• Pinů pro zvlášť pro sloupce a zvlášť pro řádky je pořád příliš
• Probíhá adresace řádků a následně (v dalším cyklu) adresace sloupců
10/57
Práce s pamětí
Čtení z paměti
CLK
RAS
CAS
Address Row
Addr
Col
Addr
Data
Out
Data
Out
Data
Out
Data
Out
DQ
RCD CL
11/57
Práce s pamětí
Čtení z paměti
CLK
RAS
CAS
Address Row
Addr
Col
Addr
Data
Out
Data
Out
DQ
RCD CL
Row
Addr
Col
Addr
Data
Out
Data
Out
RCD CL
RP
WE
Act.to.Prech. (RAS)
• CAS Latence (CL)
• RAS to CAS delay (RCD)
• RAS Precharge (RP)
• Active to Precharge delay (RAS)
• Command rate
• 2-3-2-8-T1
12/57
Práce s pamětí
SDR, DDR, DDR2, DDR3
• Princip SDR paměti
DRAM
Cell
Array
f f
• Princip DDR paměti
DRAM
Cell
Array
f f
I/O
Buffer
f
• Princip DDR2 paměti
DRAM
Cell
Array
f 2f
I/O
Buffer
2f
• Princip DDR3 paměti
DRAM
Cell
Array
f 4f
I/O
Buffer
4f
13/57
Práce s pamětí
Problémy DDR pamětí
• Zvyšující se frekvence klade nároky na paralelní sběrnice
• DDR2 modul má 240 pinů, každý by měl být zhruba stejné délky
• Nejvýše dva DDR2 moduly lze osadit na jeden kanál
• DDR3 na rychlosti 1600MB/s podporovala pouze 1 modul na kanále
• Možná řešení
• Paměťový řadič přímo v CPU (AMD Opteron, Intel Nehalem)
• Výměna DDR3 za FB-DRAM – sériové moduly, 69 pinů
14/57
Práce s pamětí
Cache paměť
• Zavádí hierarchii podobně jako paměti obecně
• Zde opět platí, že nelze mít velkou a zároveň rychlou paměť
Main Memory
L3 Cache
L2 Cache L1i Cache
L1d Cache CPU Core
Bus
15/57
Práce s pamětí
Cache paměť
• Ve světě multijaderných, hyperthreadovaných CPU jsou cache
paměti různým způsobem sdíleny
Main Memory
Bus
L3 Cache
L2 Cache
L1i CacheL1d Cache
CPU CoreCPU Core
L1i CacheL1d Cache
CPU CoreCPU Core
L3 Cache
L2 Cache
L1i CacheL1d Cache
CPU CoreCPU Core
L1i CacheL1d Cache
CPU CoreCPU Core
16/57
Práce s pamětí
Organizace cache
• Cache je organizována po řádcích (cacheline)
• Velikost řádku cache je obvykle 64bytů (někdy 32 nebo 128)
• Při přístupu do paměti je do cache načten zpravidla celý řádek
• Hierarchii cache rozlišujeme exkluzivní a inkluzivní
• Exkluzivní
• Oblast paměti pokrývá vždy pouze jedna úroveň cache (např. L2)
• Případ AMD, VIA
• Inkluzivní
• Vyšší úroveň cache pokrývá celou nižší úroveň cache
• Případ Intelu
17/57
Práce s pamětí
Organizace cache
• Paměť RAM je zpravidla rozdělena do:
• Cachovatelná oblast
• Necachovatelná oblast
• Write combining oblast
• Cachovatelná oblast
• Pro čtení dat použita cache (postupně všechny úrovně)
• Data jsou do cache přednačítána
• Data jsou zapisována do cache a rovnou do paměti RAM (write
through režim) nebo do paměti RAM až později (write back režim)
18/57
Práce s pamětí
Organizace cache
• Přednačítání do plné cache
• Nějaký řádek je nutné uvolnit (eviction)
• Exkluzivní cache – řádek lze uvolnit zrušením nemodifikovaného
řádku nebo zápisem modifikovaného řádku do cache vyšší úrovně
případně do RAM
• Inkluzivní cache – řádek lze uvolnit zrušením nemodifikovaného
řádku nebo zrušením modifikovaného řádku (pokud je řádek
cachován ve vyšší úrovni)
19/57
Práce s pamětí
Mapování cache na RAM
• Statické v blocích
• Např. řádka 1 pokrývá RAM 0-999B, řádka 2 pokrývá RAM
1000-1999B, atd.
• Cache moc nepomůže u kódu a dat blízko sebe
• Asociativní
• Podobně jako hash tabulka
• Každý řádek cache má u sebe adresu, kterou cachuje
• Vyrobit velkou asociativní paměť skoro nelze
• Řešení
• Kombinace bloku (cache set) a asociativního mapování (tag)
• Pojem 8-way cache
Tag Cache Set Offset
064/32
T S O
20/57
Práce s pamětí
Cache paměť
• V případě Intel procesorů, L1 cache neobsahuje x86(64) kód, ale
mikrokód, AMD cachuje x86 kód
• L2 cache obsahuje již x86 kód
• L1 cache používá virtuální adresy (překlad virtuální adresy na
fyzickou by příliš trval)
• L2 cache používá fyzické adresy
21/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
• Mějme program:
1 struct l {
2 struct l *next;
3 long int pad[NPAD];
4 } l;
5
6 void
7 iterate_list(struct l *head)
8 {
9 struct l *e = head;
10 long i;
11 for(i = 0; i < nelem*200; i++) {
12 e = e->next;
13 asm volatile("" : : "m" (*e));
14 }
15 }
22/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
1
2
3
4
5
6
7
8
9
10
CPUCycles/ListElement
Sequential Reading
NPAD=0
L1d L2
23/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
20
40
60
80
100
120
140
CPUCycles/ListElement
NPAD 0
NPAD 7
NPAD 15
NPAD 31
Sequential Reading
24/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
• Virtuální paměť přináší drobnou komplikaci
• Překlad virtuální adresy na fyzickou znamená nutnost projít tabulky
stránek
• Procházení tabulek stránek = náhodný přístup do paměti
• Cache na překlad adres – TLB
• C2D – 256 TLB položek pro data, 128 TLB položek pro instrukce
• Opteron 10. generace, 1024 TLB položek
25/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
50
100
150
200
250
300
CPUCycles/ListElement
NPAD 7
NPAD 512
Sequential Reading
26/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
• TLB jsou vázány na konkrétní stránkovací tabulky
• Přepnutí adresního prostoru = zahození všech položek TLB
• Důsledek:
TLB má málo položek, protože se adresní prostor přepíná velmi často
• Použití větších stránek (2 MB místo 4 kB) redukuje nutný počet
TLB položek
27/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
• Přidáme více práce při procházení seznamu
• Funkce inc přičte do prvku jedničku
• Funkce add přičte do prvku následující prvek
1 void
2 iterate_list_inc(struct l *head)
3 {
4 struct l *e = head;
5 long i;
6 for(i = 0; i < nelem*200; i++) {
7 e->pad[0]++;
8 e = e->next;
9 asm volatile("" : : "m" (*e));
10 }
11 }
12
13 void
14 iterate_list_add(struct l *head)
15 {
16 struct l *e = head;
17 long i;
18 for(i = 0; i < nelem*200; i++) {
19 e->pad[0]+=e->next->pad[0];
20 e = e->next;
21 asm volatile("" : : "m" (*e));
22 }
23 }
28/57
Práce s pamětí
Jednovláknový sekvenční přístup do paměti
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
10
20
30
40
50
60
70
CPUCycles/ListElement
NPAD 7
NPAD 7 Inc
NPAD 7 Add
29/57
Práce s pamětí
Jednovláknový náhodný přístup do paměti
• Pro sekvenční přístup je schopno CPU sledovat, jaká data budou
potřeba v blízké budoucnosti.
• Pro NPAD=0, paměť v burst režimu pošle celou cacheline (64B) =
8 rvků seznamu
• Pro NPAD=7 a další už burst režim pošle pouze jeden prvek
• U sekvenčního přístupu je velká šance, že se nebude
přeprogramovávat řádek paměti
• U náhodného přístupu je predikce problematická, často je potřeba
kompletně naprogramovat přístup k paměti
30/57
Práce s pamětí
Jednovláknový náhodný přístup do paměti
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
50
100
150
200
250
300
CPUCycles/ListElement
Sequential
Random
Sequential Reading
31/57
Práce s pamětí
Jednovláknový náhodný přístup do paměti
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
50
100
150
200
250
300
CPUCycles/ListElement
NPAD 7 Seq
NPAD 7 Rand
Sequential vs. Random reading
32/57
Práce s pamětí
Critical Word Load
• Cacheline není naplněna instantně ale postupně
• Sběrnicí od paměti nebo cache vyšší úrovně nelze poslat 64 bytů
v jednom cyklu
• Procesor je schopen pracovat i s částečně naplněnou cacheline
• Důsledek:
• Záleží na pozici slova v cache, často používané položky je lepší
umísťovat na začátek
33/57
Práce s pamětí
Critical Word Load
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
2
4
6
8
10
12
14
16
18
20
22
24
26
28
30
CPUCycles/ListElement
NPAD 7
NPAD7 Critical word last
Critical Word Load
34/57
Práce s pamětí
Zápisy do paměti
• L1 cache je prakticky vždy nesdílená mezi jádry multiprocesorového
systému
• Cachováním položek paměti může dojít k nekonzistencím
• Existují protokoly koherence cache
• Je nepraktické umožnit procesorům přímý přístup do cache jiného
procesoru
• Procesory odposlouchávají na paměťové sběrnici a mají představu o
tom, co ostatní procesory cachují
35/57
Práce s pamětí
MESI protokol
• Každý řádek cache je v jednom ze čtyř stavů:
• Modified: lokální procesor modifikoval řádku cache, která se
nenachází v žádné cache jiného procesoru
• Exclusive: řádka není modifikovaná a není uložena v cache jiného
procesoru
• Shared: řádka není modifikovaná a může existovat v cache jiného
procesoru
• Invalid: řádka není platná
36/57
Práce s pamětí
Zápisy do paměti
M
S
E
I
local read
remote read
local write
remote write
37/57
Práce s pamětí
Přechody stavů
• Request for ownership (RFO) zpráva
• Modified ⇒ Invalid
• Shared ⇒ Modified
38/57
Práce s pamětí
Zápisy do paměti
• Mějme seznam, který leží v paměti následovně
Cache line
• Každý seznam prochází jiné vlákno
• Použijeme inc variantu
tj. položky cache se mění střídavě různými CPU
• Velké množství RFO zpráv (Modified ⇒ Invalid)
39/57
Práce s pamětí
Zápisy do paměti
QuadCore Opteron, 64kB L1, 512kB L2, 2MB L3 Cache
NPAD 1, Inc
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
100
200
300
400
CPUCycles/ListElement
1 Thread*
2 Threads*
4 Threads*
8 Threads*
1 Thread
2 Threads
4 Threads
8 Threads
RFO Costs
bez * je původní seznam
40/57
Práce s pamětí
Práce ve vláknech
QuadCore Opteron, 64kB L1, 512kB L2, 2MB L3 Cache
NPAD 7
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
50
100
150
200
CPUCycles/ListElement
1 Thread
2 Threads
4 Threads
8 Threads
Sequential read acces, Multiple threads
NPAD 7
41/57
Práce s pamětí
Práce ve vláknech
QuadCore Opteron, 64kB L1, 512kB L2, 2MB L3 Cache
NPAD 7, Increment
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
100
200
300
400
CPUCycles/ListElement
1 Thread
2 Threads
4 Threads
8 Threads
Sequential Inc access, Multipe threads
NPAD 7
42/57
Práce s pamětí
Práce ve vláknech
QuadCore Opteron, 64kB L1, 512kB L2, 2MB L3 Cache
NPAD 7, Add next
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
50
100
150
200
250
300
350
400
CPUCycles/ListElement
1 Thread
2 Threads
4 Threads
8 Threads
Sequential Add next, Multiple threads
NPAD 7
43/57
Práce s pamětí
Násobení matic
• Mějme 2 velké matice 1000×1000 typu double
• Potřebujeme je vynásobit
• Naivní algoritmus
1 for(i=0; i < N; ++i)
2 for(j=0; j < N; ++j)
3 for(k=0; k < N; ++k)
4 res[i][j] += mul1[i][k] * mul2[k][j];
• Lze vidět, že matici mul2 procházíme po sloupcích = nesekvenční
přístup
44/57
Práce s pamětí
Násobení matic
= *
45/57
Práce s pamětí
Násobení matic
• Matici mul2 můžeme transponovat
1 for(i=0; i < N; ++i)
2 for(j=0; j < N; ++j)
3 tmp[i][j] = mul2[j][i];
4
5 for(i=0; i < N; ++i)
6 for(j=0; j < N; ++j)
7 for(k=0; k < N; ++k)
8 res[i][j] += mul1[i][k] * tmp[j][k];
• Čas nutný na transpozici započítáme do běhu programu
Original Transposed
Opteron 2356
Cycles 54,415,069,438 11,303,294,945
Relative 100% 20%
C2D E6850
Cycles 10,543,642,668 3,822,534,279
Relative 100% 36%
46/57
Práce s pamětí
Násobení matic
• Ne vždy je možné udělat kopii matice, není to ani optimální řešení
• Lze rozdělit matici do bloků SM×SM a násobit bloky
1 for(i=0; i < N; i+= SM)
2 for(j=0; j < N; j+=SM)
3 for(k=0; k < N; k+= SM)
4 for(i2 = 0, rres = &res[i][j], rmul1 = &mul1[i][k];
5 i2 < SM; ++i2, rres += N, rmul1 += N)
6 for(k2 = 0, rmul2 = &mul2[k][j]; k2 < SM; ++k2,
7 rmul2 += N)
8 for(j2 = 0; j2 < SM; ++j2)
9 rres[j2] += rmul1[k2] * rmul2[j2];
47/57
Práce s pamětí
Násobení matic
= *
= *
= *
48/57
Práce s pamětí
Násobení matic
• Dále je možno využít vektorových instrukcí, které zpracují několik
operací najednou
1 for(i=0; i < N; i+= SM)
2 for(j=0; j < N; j+=SM)
3 for(k=0; k < N; k+= SM)
4 for(i2 = 0, rres = &res[i][j], rmul1 = &mul1[i][k];
5 i2 < SM; ++i2, rres += N, rmul1 += N) {
6 _mm_prefetch(&rmul1[8], _MM_HINT_NTA);
7 for(k2 = 0, rmul2 = &mul2[k][j]; k2 < SM; ++k2,
8 rmul2 += N) {
9 __m128d m1d = _mm_load_sd(&rmul1[k2]);
10 m1d = _mm_unpacklo_pd(m1d, m1d);
11 for(j2 = 0; j2 < SM; j2+=2) {
12 __m128d m2 = _mm_load_pd(&rmul2[j2]);
13 __m128d r2 = _mm_load_pd(&rres[j2]);
14 _mm_store_pd(&rres[j2],
15 _mm_add_pd(_mm_mul_pd(m2, m1d), r2));
16 }
17 }
18 }
49/57
Práce s pamětí
Násobení matic
Original Transposed Sub-matrix Vectorized
Opteron 2356
Cycles 54,415,069,438 11,303,294,945 5,237,161,135 5,178,595,240
Relative 100% 20% 9% 9%
C2D E6850
Cycles 10,543,642,668 3,822,534,279 3,056,154,579 1,442,368,584
Relative 100% 36% 28% 13%
50/57
Práce s pamětí
Blízkost dat
• Rozdíly mezi malou strukturou zaujímající 1 cacheline a mezi větší
přes více cachelines je znát
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
100
200
300
400
500
600
CPUCycles/ListElement
1 Thread
2 Threads
4 Threads
8 Threads
1 Thread*
2 Threads*
4 Threads*
8 Threads*
Inc Sequential NPAD 15, Multiple threads
1 cacheline vs. 2 cachelines*
51/57
Práce s pamětí
Blízkost dat
• Blízkost dat se vztahuje i na NUMA systémy, přístup do nelokální
paměti je dražší
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
10
20
30
40
50
60
CPUCycles/ListElement
1 Thread
2 Threads
4 Threads
1 Thread*
2 Threads
4 Threads
Sequential Inc, Multiple threads
NUMA factor
52/57
Práce s pamětí
Blízkost dat
• Blízkost malých dynamických kusů paměti narušuje malloc()
• Každý alokovaný kus paměti má 16 bytů hlavičku a případně patičku
pro zarovnání na 16 bytů
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
50
100
150
CPUCycles/ListElement
NPAD 0
NPAD 7
NPAD 15
NPAD 31
NPAD 0(*)
NPAD 7(*)
NPAD 15(*)
NPAD 31(*)
Sequential Reading
bez * je ruční alokace v malloc() bloku
53/57
Práce s pamětí
Blízkost dat
• U blízkosti dat je nutné respektovat jejich vhodné zarovnání,
zejména zarovnání na cacheline a slova na hranici slov
• Tj. 8bytový long na hranici 8 bytů
• Struktury na 64 bytů
• Kódy funkcí na 64 bytů
• Pro zarovnání dat existuje volání int posix_memalign(void
**memptr, size_t align, size_t size)
• Pro zarovnávání golbálních dat gcc nabízí __attribute__
((aligned(64)))
54/57
Práce s pamětí
Blízkost dat
2
10
2
12
2
14
2
16
2
18
2
20
2
22
2
24
2
26
2
28
Set Size
0
10
20
30
40
50
CPUCycles/ListElement
1 Thread
2 Threads
4 Threads
1 Thread
2 Threads
4 Threads
Sequential Inc, Multiple threads
Aligned/Unaligned access
55/57
Práce s pamětí
Blízkost instrukcí
• Inlining funkcí nevede vždy k nárůstu rychlosti
• Penalta za cache miss je proti volání funkce obrovská
• Optimalizace větvení
• Překladač může optimalizovat podmínky podle toho, zda se většinou
provede či nikoliv
• Konstrukty v jádře if(likely(a>1)) { },
if(unlikely(a>1)) { }
• #define likely(expr) __builtin_expect(!!(expr),1)
• #define unlikely(expr)
__builtin_expect(!!(expr),0)
56/57
Práce s pamětí
Nástroje pro analýzu
• valgrind – cachegrind – analýza využití cache, cache miss
apod.
• valgrind – massif – analýza využití paměti v čase
• podobně memusage
• gcc je schopno analyzovat úspěšnost branchprediction
-fprofile-generate, -fprofile-use
• Informace o cache hits/miss, TLB hit/miss lze získat přímo z CPU
http://user.it.uu.se/~mikpe/linux/perfctr/
57/57
Práce s pamětí
Výhled do budoucnosti
• TLB
• Tagovaná TLB cache, není nutné vylití TLB při přepnutí kontextu
• Transakční paměť
• Podpora LL/SC instrukcí (Nemají ABA problém)
• Transakční L1t cache
• Rozšíření MESI protokolu – není technologicky náročné