Existe uma implementação de nível de referência do BLAS em C / C ++?

11

A implementação do netlib BLAS é uma excelente referência, sendo principalmente não otimizada e bem documentada (por exemplo, zgemm ). No entanto, está no Fortran 77, tornando-o um tanto inacessível para aqueles com uma educação em programação mais moderna. Existe uma implementação de nível de referência do BLAS, como netlib, em C / C ++?

blas education

— Max Hutchinson
fonte

10

Você já viu a implementação da Biblioteca Científica GNU ? Acho que o código fonte é suficientemente legível e as rotinas estão bem documentadas.

— Juan M. Bello-Rivas
fonte

Parece bom para mim. A documentação está um pouco ausente, mas os nomes das variáveis são escolhidos o suficiente para que eu fique claro. Provavelmente vou anexar o comentário de abertura das rotinas netlib BLAS. O que exatamente você discorda? Você tem uma alternativa?

— Max Hutchinson

6

Uma implementação notável da linguagem C do BLAS é o ATLAS. Entre recursos úteis:

As rotinas de álgebra implementadas tanto em C simples quanto em versões altamente otimizadas para montadores, com assistência otimizada, para múltiplas arquiteturas e variantes.
O sistema de compilação apresenta um "auto-tuner" que compila várias variantes da biblioteca ATLAS para estabelecer qual será a mais rápida na máquina especificada.

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

— oakad
fonte

Eu olhei para o ATLAS, mas perdi isso. O caminho para a implementação de referência é "src / blas / reference", com "ref" inserido entre o caractere de tipo e o nome da rotina e com os argumentos de caractere anexados.

— Max Hutchinson

3

O Netlib também produz CLAPACK , que inclui BLAS, mas é apenas o código fortran executado no f2c e, portanto, é um pouco desajeitado (por exemplo, zgemm ).

— Max Hutchinson
fonte

1

Para uma implementação de alto desempenho que não esteja apenas entre os de melhor desempenho (melhor que 85% do pico em 60 núcleos do Intel Xeon Phi), mas também seja a mais bem escrita, dê uma olhada no BLIS:

https://github.com/flame/blis

— Robert van de Geijn
fonte

Este não parece ser o que afirma ser uma implementação de nível de referência

— mabraham

Gostamos de reduzi-lo como uma implementação de referência de alto desempenho.

— Robert van de Geijn 01/04/19

0

Tenho implementações de alguns bits do BLAS / LAPACK no RNP e RNP2 .

— Victor Liu
fonte

Isso é interessante, mas é definitivamente pós-BLAS e menos direto que o GSL e a implementação de referência no ATLAS.

— precisa

-1

Atualmente, estamos trabalhando em um curso on-line aberto em grande escala, "LAFF-On High-Performance Computing", que usa dgemm como exemplo que leva a diferentes níveis de paralelização: nível de instrução, OpenMP, MPI.

Esta não é uma implementação de referência para o BLAS, mas é uma referência de como codificar o BLAS (para desempenho). Para se manter informado, visite www.ulaff.net

— Robert van de Geijn
fonte