Perguntas com a marcação «performance»

Trabalho em ciência da computação e, como resultado, passo uma quantidade não trivial do meu tempo tentando aumentar a taxa de transferência científica de muitos códigos, além de entender a eficiência desses códigos. Vamos supor que eu avaliei a troca desempenho versus legibilidade / reutilização...

Nos meus cursos de análise numérica, aprendi a analisar a eficiência dos algoritmos contando o número de operações de ponto flutuante (flops) necessárias, em relação ao tamanho do problema. Por exemplo, no texto de Trefethen & Bau sobre Álgebra Linear Numérica, existem até imagens em 3D das...

Fortran tem um lugar especial na programação numérica. Você certamente pode criar software bom e rápido em outros idiomas, mas o Fortran continua com um desempenho muito bom, apesar da idade. Além disso, é mais fácil criar programas rápidos no Fortran. Criei programas rápidos em C ++, mas você...

Eu estava comparando alguns dos meus códigos com os códigos "estoque" do MATLAB. Estou surpreso com os resultados. Corri um código de exemplo (Sparse Matrix) n = 5000; a = diag(rand(n,1)); b = rand(n,1); disp('For a\b'); tic;a\b;toc; disp('For LU'); tic;LULU;toc; disp('For Conj...

O ATLAS é um substituto gratuito do BLAS / LAPACK que se ajusta à máquina quando compilado. MKL é a biblioteca comercial fornecida pela Intel. Essas duas bibliotecas são comparáveis ​​no que diz respeito ao desempenho ou o MKL tem a vantagem de algumas tarefas? Se sim, quais...

Antecedentes: Eu acho que posso querer portar algum código que calcule produtos vetoriais exponenciais da matriz usando um método de subespaço Krylov do MATLAB ao Python. (Especificamente, a função expmvp de Jitse Niesen , que usa um algoritmo descrito neste documento .) No entanto, eu sei que, a...

Suponha que eu esteja trabalhando em um código científico em C ++. Em uma discussão recente com um colega, foi argumentado que os modelos de expressão poderiam ser uma coisa muito ruim, potencialmente tornando o software compilável apenas em determinadas versões do gcc. Supostamente, esse problema...

Qual é a abordagem preferida e eficiente para interpolar dados multidimensionais? Coisas com as quais me preocupo: desempenho e memória para construção, avaliação de lote único dimensões de manipulação de 1 a 6 ordem linear ou superior capacidade de obter gradientes (se não lineares) grade...

Começarei com minha experiência pessoal em nosso laboratório. Nos 9 e 10 dias da ifort, costumávamos ser bastante agressivos com as otimizações, compilando com -O3 e sinalizadores específicos do processador (-xW -xSSE4.2 por exemplo). Mas a partir do ifort 11, começamos a notar: 1. algumas...

Existe um método para resolver sistemas lineares da forma onde é uma matriz SPD fixa e são matrizes diagonais positivas?O(n3+n2k)O(n3+n2k)O(n^3+n^2 k)kkk(Di+A)xi=bi(Di+A)xi=bi(D_i + A) x_i = b_iAAADiDiD_i Por exemplo, se cada é escalar, basta calcular o SVD Uma . No entanto, isso se divide em D...

Como as matrizes Python / Numpy escalam com o aumento das dimensões da matriz? Isso se baseia em algum comportamento que notei ao comparar o código Python para esta pergunta: Como expressar essa expressão complicada usando fatias numpy O problema envolveu principalmente a indexação para preencher...

A linguagem C ++ fornece programação genérica e metaprogramação por meio de modelos. Essas técnicas chegaram a muitos pacotes de computação científica em larga escala (por exemplo, MPQC , LAMMPS , CGAL , Trilinos ). Mas o que eles realmente contribuíram para a computação científica em valor que vai...

Quando cálculos limitados de largura de banda de memória são executados em ambientes de memória compartilhada (por exemplo, encadeados via OpenMP, Pthreads ou TBB), há um dilema de como garantir que a memória seja distribuída corretamente na memória física , de modo que cada encadeamento acesse a...

Estou procurando alguma implementação de código aberto (qualquer um de Python, C, C ++, Fortran é bom) de aproximação racional a uma função. Algo ao longo do artigo [1]. Dou-lhe uma função e ele me devolve dois polinômios, cuja razão é a aproximação no intervalo determinado e o erro está oscilando...

Nossas simulações de dinâmica de plasma geralmente produzem muita informação. Durante as simulações, registramos várias propriedades físicas em uma grade (x, y, z, t) que é tão grande quanto (8192x1024x1024x1500), por pelo menos 10 propriedades. Esta informação é processada após a simulação ser...

Estou comparando dois métodos iterativos para inverter matrizes quadradas aleatórias. Como as matrizes são aleatórias, todo caso de teste leva quantidades diferentes de iterações e diferentes tempos decorridos. Minha pergunta é, além do tempo médio da CPU, o valor médio das iterações realizadas...

Eu sei que o ATLAS é capaz de otimizar-se para a máquina na qual é compilado e, portanto, o máximo de benefícios é encontrado ao compilar a partir do código-fonte. Existe algum benefício em compilar o LAPACK a partir da fonte? Seria muito mais fácil instalar o pacote

Pergunta: Quais métodos estão disponíveis para calcular com precisão e eficiência a estrutura de esparsidade de uma matriz de elementos finitos? Info: Estou trabalhando em um solucionador de equações de pressão de Poisson, usando o método de Galerkin com base em Lagrange quadrático, escrito em C,...

Ao avaliar o número de FLOPs em uma função simples, geralmente é possível simplesmente descer a expressão que calcula os operadores aritméticos básicos. No entanto, no caso de declarações matemáticas envolvendo divisão uniforme, não se pode fazer isso e espera poder comparar com contagens de FLOP...

Neste programa de exemplo, estou fazendo a mesma coisa (pelo menos acho) de duas maneiras diferentes. Estou executando isso no meu pc Linux e monitorando o uso de memória com o top. Usando o gfortran, acho que na primeira maneira (entre "1" e "2") a memória usada é de 8,2 GB, enquanto na segunda...