Perguntas com a marcação «blas»

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. …

74 fortran  c  blas  benchmarking 

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 bibliotecas são independentes …

26 linear-algebra  software  lapack  blas 

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 ) + ( …

17 floating-point  blas 

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 sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

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 …

15 parallel-computing  matrix  lapack  blas 

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 para CG e …

14 linear-algebra  petsc  blas 

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 …

11 linear-algebra  software  blas  complex 

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, …

11 blas  education 

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 disponíveis?

11 linear-algebra  lapack  blas  gpu 

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 matriz?

10 petsc  sparse-matrix  lapack  blas 

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 …

8 linear-algebra  matrix  blas 

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 …

8 performance  fortran  blas