Suponha que eu tenha duas matrizes Nx2, Mx2 representando N, M2 vetores, respectivamente. Existe uma maneira simples e boa de calcular distâncias entre cada par de vetores (n, m)? A maneira mais fácil, porém ineficiente, é: d = zeros(N, M); for i = 1:N, for j = 1:M, d(i,j) = norm(n(i,:) -...
Estou tentando descobrir se existe uma maneira mais rápida de calcular todos os autovalores e autovetores de uma matriz de adjacência muito grande e esparsa do que usar scipy.sparse.linalg.eigsh Até onde eu sei, esse método usa apenas a escassez e atributos de simetria da matriz. Uma matriz de...
Nos métodos de decomposição de domínio (DD) e multigrid (MG), pode-se compor a aplicação das atualizações de bloco ou correções grosseiras como aditivas ou multiplicativas . Para solucionadores pontuais, essa é a diferença entre as iterações de Jacobi e Gauss-Seidel. O mais suave multiplicativo...
Considere dois computadores com configurações diferentes de hardware e software. Ao executar exatamente o mesmo código serial de Navier-Stokes em cada plataforma, leva um tempo x e y para executar uma iteração para o computador 1 e 2, respectivamente. Nesse caso, , é a diferença de tempo de...
A computação paralela é freqüentemente modelada usando uma taxa local determinística de computação, sobrecarga de latência e largura de banda da rede. Na realidade, são espacialmente variáveis e não determinísticas. Estudos como Skinner e Kramer (2005) observam distribuições multimodais, mas a...
Este é um problema simples de enumeração de simetria. Eu dou o histórico completo aqui, mas nenhum conhecimento de química quântica é necessário. O integrante dois partícula é: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) ψ k ( x ) ψ l ( x ' )⟨ I j | k l ⟩⟨ij|kl⟩\langle ij|kl\rangle E tem os...
Estou pesquisando as referências mais importantes da literatura e de slides para modelar os custos atuais e futuros de energia de operações de ponto flutuante e transferências de dados na CPU, memória, rede e armazenamento. Marquei esta pergunta como wiki da comunidade e prefiro que você limite...
A multiplicação da matriz (Mat * Mat e Mat * Vec) é dimensionada com o número de não-zeros ou com o tamanho da matriz? Ou alguma combinação dos dois. Que tal com forma. Por exemplo, eu tenho uma matriz 100 x 100 com 100 valores ou uma matriz 1000 x 1000 com 100 valores. Ao esquadrar essas...
Eu tenho que fazer uma transformação de coordenadas entre dois sistemas de referência (eixos). Para isso, três matrizes ( ) devem ser multiplicadas devido ao uso de alguns eixos intermediários. Pensei em duas abordagens para resolver isso:3 × 33×33\times3 Método # 1 : Fazendo a multiplicação...
O C ++ 11 introduz a semântica de movimentação que pode, por exemplo, melhorar o desempenho do código em situações em que o C ++ 03 precisaria executar uma construção de cópia ou atribuição de cópia. Este artigo relata que o código a seguir experimenta uma velocidade de 5x quando compilado com C +...
Estou tentando entender a diferença entre essas duas placas gráficas para computação acadêmica, especificamente para o componente DGEMM. Se olharmos para as estatísticas brutas, ambas têm o mesmo chip GK110, estatísticas comparáveis em praticamente todas as categorias e, acredito, têm a mesma...
Eu tenho um problema em que preciso encontrar todos os pares de autovalores positivos (como o valor próprio é positivo) de uma matriz não simétrica pequena (geralmente menor que 60x60). Posso parar de calcular quando o valor próprio é menor que um determinado limite. Eu sei que os autovalores são...
O título diz a maior parte. Estou procurando uma biblioteca leve e fácil de usar que possa ser usada em projetos Android (NDK). Para coisas densas, gosto de usar o Eigen, mas não encontrei muitas bibliotecas abrangentes (e documentadas!) Para coisas esparsas que "simplesmente funcionam" em um...
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } O código...
Qual é a maneira mais eficiente de obter uma matriz esparsa complexa do meu código Fortran para o PETSc? Entendo que isso depende do problema, então tentei fornecer o maior número possível de detalhes relevantes abaixo. Eu estou brincando com o solucionador de autovalor FEAST [1] para problemas do...
Eu já vi muitas publicações em periódicos de Física Computacional usarem métricas diferentes para o desempenho de seu código. Especialmente para o código GPGPU, parece haver uma grande variedade de resultados de tempo publicados pelas pessoas. Em particular, eu já vi: Comparações de...
É comum dizer que o Linux é mais rápido e por boas razões. Mas, como declarado no título, quanto a escolha do SO importa para o desempenho do código de computação científica? Para algo em que o sistema operacional está mais diretamente envolvido (como jogar jogos ou renderização de vídeo), com...
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...
Estou analisando o desempenho dos solucionadores lineares do PETSc. Como eu entendo, speedup = Tempo sequencialTempo Paralelo.acelerar=Tempo SeqüencialTempo Paralelo.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Eu sei que a execução do código paralelo em um processador...
Estou pensando em usar redes neurais artificiais (RNA) para prever as taxas de reação no meu fluido, em vez de resolver o sistema completo de ODEs rígidas. Algumas pessoas do meu laboratório já fizeram algum trabalho para não começar do zero, mas estou tendo problemas com meus aplicativos. Acho que...