Esiste un'implementazione a livello di riferimento di BLAS in C / C ++?

L'implementazione BLAS di netlib è un riferimento eccellente, essendo per lo più non ottimizzata e ben documentata (ad es. Zgemm ). Tuttavia, è in Fortran 77, rendendolo in qualche modo inaccessibile a quelli con una formazione di programmazione più moderna. Esiste un'implementazione a livello di riferimento di BLAS, come netlib, in C / C ++?

Hai esaminato l' implementazione della GNU Scientific Library ? Trovo che il codice sorgente sia sufficientemente leggibile e le routine sono ben documentate.

Una notevole implementazione in linguaggio C di BLAS è ATLAS. Tra le funzioni utili:

1. Le routine di algebra implementate sia in versione C semplice che in versioni assistite da assemblatore altamente ottimizzate per più architetture e varianti.
2. Il sistema di compilazione presenta un "sintonizzatore automatico" che compila più varianti della libreria ATLAS per stabilire quale sarà la più veloce su una determinata macchina.

http://math-atlas.sourceforge.net/

Netlib produce anche CLAPACK , che include BLAS, ma è solo il codice fortran eseguito attraverso f2c ed è quindi un po 'goffo (ad esempio zgemm ).

Per un'implementazione ad alte prestazioni che non è solo tra le più performanti (migliore dell'85% del picco su 60 core di Intel Xeon Phi), ma è anche la più magnificamente scritta, dai un'occhiata a BLIS:

https://github.com/flame/blis

Ho implementazioni di alcuni bit di BLAS / LAPACK in RNP e RNP2 .

Attualmente stiamo lavorando a un corso online aperto di massa, "LAFF-On High Performance Computing" che utilizza dgemm come esempio che guida attraverso diversi livelli di parallelizzazione: livello di istruzione, OpenMP, MPI.

Questa non è un'implementazione di riferimento per il BLAS, ma è un riferimento per come codificare il BLAS (per le prestazioni). Per essere informato, visitare www.ulaff.net