Administrace komentářů

1) První a základní optimalizace openMP programu je odstranění kritických sekcí. Dělat porovnávání různejch alokátorů, když tam mám odstranitelnou kritickou sekci je IMHO blbina.

2) Optimalizovat multithread aplikace bez toho, aniž by ses chtěl zabejvat jednotlivejma threada je jak jezdit s autem, ale nechtít šahat na volant. Jo, taky někam dojedeš... U multithread aplikací je dosti klíčové, jak se chovají thready, jak si navzájem zneplatňujou cache, jak na sebe musej čekat atd...

3) Použití TLS není žádná kontrola nad jednotlivými thready. Konkrétně v novym C++ 11 to znamená, že deklaruješ proměnou jako thread_local a to je vše. V gcc lze použít __thread. Je to na kód snad úspornější, než označování kritické sekce a vyžaduje to stejnou míru "starosti" o konkurenci threadů (tady v tý proměnný může dojít ke konfliktu, jak ji ochráním?) jen ji řešíš jiným (a daleko vhodnějším) způsobem.

4) Nevim jaks testoval, ale muj rychlej test ukazuje, že vlastní alokátor je škálovatelnější a násobně rychlejší než malloc

Muj alokátor, openmp (4thready, 2core + multithreading) 0.580s Muj alokátor singlethread 0.609s Malloc, openmp (4thready, 2core + multithreading) 10,867s Malloc singlethread 6.647s

na dvanactijadru jsou vysledky podobný, malloc neškáluje a i v singlethread je násobně pomalejší, custom allokátor škáluje krásně

kompilovano s

gcc -fopenmp -O1 test.c -std=c99

Kód

#include "stdlib.h"

void* alloc1(int bytes) {
     return malloc(bytes);
}

const int CHUNK_SIZE = 1024*1024;


void* alloc2(int bytes) {
   static __thread void* alloc_start = 0;
   static __thread void* alloc_end = 0;
   static __thread void* alloc_ptr = 0;


   if(alloc_start + bytes > alloc_end) {
     void* new_buffer = malloc(CHUNK_SIZE);
     *(void**)new_buffer = alloc_start;           //previous allocated chunk
     *(int *)(new_buffer + sizeof(void*)) = 0;    //number of allocated substrings
     alloc_start = new_buffer;
     alloc_end = new_buffer + CHUNK_SIZE;
     alloc_ptr = alloc_start + sizeof(void*) + sizeof(int);
   }
   void *ret = alloc_ptr;
   alloc_ptr+=bytes;
   return ret;
}

int main(void) {
#pragma omp parallel for
  for(int r=0;r<1024*1024*200;r++)
    alloc2(r % 21 * r % 13);
}