Esta pergunta é uma extensão de duas discussões que surgiram recentemente nas respostas a " C ++ vs Fortran for HPC ". E é um pouco mais um desafio do que uma pergunta ... Um dos argumentos mais ouvidos em favor do Fortran é que os compiladores são apenas melhores. Como a maioria dos compiladores...
Eu tenho procurado em bibliotecas de álgebra linear C ++ para um projeto que eu tenho trabalhado. Algo que ainda não entendi é a conexão do BLAS e LAPACK a outras bibliotecas de álgebra linear. Examinando este artigo sobre bibliotecas de álgebra linear, achei interessante que: algumas...
Dadas duas implementações BLAS diferentes, podemos esperar que elas façam exatamente os mesmos cálculos de ponto flutuante e retornem os mesmos resultados? Ou pode acontecer, por exemplo, que alguém calcule um produto escalar como e um como ( x 1 y 1 + x 2 y 2 ) + ( x 3 y 3 + x 4( ( x1 1y1 1+...
Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com
Desculpe pelo longo post, mas eu queria incluir tudo o que achava relevante desde o início. O que eu quero Estou implementando uma versão paralela dos métodos de subespaço de Krylov para matrizes densas. Principalmente GMRES, QMR e CG. Percebi (após a criação de perfil) que minha rotina DGEMV era...
Onde quer que eu tenha visto, os tutoriais / documentos do PETSc etc. dizem que é útil para álgebra linear e geralmente especifica que sistemas esparsos serão beneficiados. E matrizes densas? Estou preocupado sobre como resolver para densa .AA x = bUMAx=bAx=bUMAUMAA Eu escrevi meu próprio código...
A nVidia, por exemplo, possui CUBLAS, que promete aceleração de 7-14x. Ingenuamente, isso não chega nem perto do rendimento teórico de qualquer placa de GPU da nVidia. Quais são os desafios para acelerar a álgebra linear em GPUs, e já existem roteiros mais rápidos de álgebra linear...
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...
O produto complexo interior tem duas definições diferentes decididas pelo convenções: ˉ u T v ou u T ˉ v . No BLAS, encontrei as rotinas cdotu, zdotu e cdotc, zdotc. Os dois primeiros rotinas realmente computação u T V (um produto interno falso!) E os dois últimos rotinas conjugar o primeiro vector...
A compilação do PETSc com uma biblioteca externa BLAS / LAPACK afeta significativamente o desempenho em matrizes esparsas ou apenas usa essas bibliotecas para matemática densa de
Eu implementei o CG no FORTRAN vinculando-o ao Intel MKL. Quando houver declarações como: ( Consulte Wikipedia ) p=r; x=x+alpha*p r=r-alpha*Ap; ou similares no QMR (em quantidade muito maior) v_tld = r; y = v_tld; rho = norm( y ); w_tld = r; z = w_tld; xi = norm( z ); (and more) Faz...
Eu tenho lido sobre a implementação de multiplicação densa de matrizes quando a matriz não se encaixa no cache. Um dos gráficos que eu vi (slide 9 desses slides ) mostra quedas repentinas no desempenho usando o algoritmo ingênuo. Essas quedas atingem cerca de 50% da velocidade, ocorrem enquanto a...