Meu sistema é um problema simétrico de FE com multiplicadores de lagrange (por exemplo, fluxo incompressível de Stokes): ( ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} onde é o caso típico (eu até me certifiquei de que as equações sejam numeradas para que os multiplicadores de...
Quanto mais lento o petsc4py vs c / c ++ / fortran? Sei que isso dependerá significativamente do código que está sendo executado, mas e algo simples como um produto de vetor de
12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Atualmente, estou fazendo isso com a biblioteca do tatu , mas ela acaba sendo muito lenta. O problema é que eu preciso fazer isso por um trilhão de matrizes e acontece que calcular os dois...
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...
Aqui está um uso típico de ponteiros de função em C. Eu gostaria de fazer algo semelhante no Fortran. Tenho algumas idéias, mas gostaria de saber se existe alguma maneira canônica de fazer isso. Os ponteiros de função e os contextos transmitidos pelo usuário são armazenados e depois chamados...
Eu gostaria de usar o método de 8ª ordem Runge-Kutta (89) em um aplicativo de mecânica / astrodinâmica celeste, escrito em C ++, usando uma máquina Windows. Portanto, pergunto-me se alguém conhece uma boa biblioteca / implementação documentada e livre para usar? Tudo bem se estiver escrito em C,...
Qual é calculado mais rapidamente, ou log a c ou b √umababa^bregistroumaclogac\log_a c ? a,becsão reais positivos comb>1.c√bcb\sqrt[b]{c}umaaabbbcccb > 1b>1b>1 Que tipos de algoritmos você usará na comparação? Quais são as suas complexidades? Por exemplo, quando ou c ≈ a bc ≡ abc≡abc...
Eu sou novo no pacote PETSc. Eu tenho uma matriz A de 4000 x 4000 no formato de mercado de matriz e quero que o PETSc resolva isso usando vários processadores. Sei como resolver o sistema em um único processador, mas não sei como distribuir a matriz e os vetores entre diferentes processadores....
Eu tenho um modelo considerável (~ 5000 linhas) escrito em C. É um programa serial, sem geração de números aleatórios em nenhum lugar. Utiliza a biblioteca FFTW para funções que utilizam a FFT - não conheço os detalhes da implementação da FFTW, mas presumo que as funções nela também sejam...
Estou procurando referências credíveis informando quanto recursos os supercomputadores gastam em coordenar versus realizar um trabalho real relacionado a tarefas. Os recursos podem estar disponíveis para processamento, mas até Watts parece uma unidade válida. Acredito que um dos meus professores...
Nos últimos anos, surgiram várias bibliotecas / projetos de software que oferecem uma forma ou outra de paralelismo de memória compartilhada de uso geral orientado a dados. A idéia principal é que, em vez de escrever um código explicitamente encadeado, os programadores implementam seus algoritmos...
Eu sou muito novo na programação GPGPU, então, desculpe-me se a pergunta não for particularmente apropriada. Pelo que entendi, a programação da GPU é um trabalho muito complicado de engenharia, quando comparado à programação usual da CPU. É preciso ter muito cuidado com questões de divergência,...
Existem muitos métodos numéricos bem conhecidos para resolver equações do tipo por exemplo, método de bissecção, método de Newton, etc.f( x ) = 0 ,x ∈ Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, Na minha aplicação, é calculado com um método estocástico (o resultado é uma média).f( X...
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
Como exemplificado pela resposta de Jed Brown a Custos de pesquisas versus cálculos , o uso de operações de ponto flutuante vetorizado versus não vetorizado resulta em código muito mais rápido. Muitos compiladores modernos afirmam que podem executar a vetorização automática. Como faço para que...
Preciso resolver Ax = b, mas percebo que, mesmo que seja escasso, o armazenamento dos coeficientes da matriz do meu problema precisará de muita memória. Então, agora estou pensando em usar um método sem matriz, porque os mesmos coeficientes aparecem muito tempo na matriz, para que eu possa usar meu...
A maior parte da minha programação são códigos de pesquisa pontuais em C para uso próprio. Nunca distribuí nenhum código para outros colaboradores que não sejam colaboradores próximos. Eu desenvolvi um algoritmo que estou publicando em uma revista científica. Desejo fornecer o código fonte e talvez...
Basicamente, o FEM parece ser um problema praticamente "resolvido". Existem inúmeras estruturas poderosas, como Trilinos, PETSc, FEniCS, Libmesh ou MOOSE. Uma coisa que eles têm em comum: são extremamente "pesados". Primeiro, a instalação normalmente é super dolorosa. Segundo, a interface / API...
Estou mergulhando no mundo fascinante da análise de elementos finitos e gostaria de resolver um grande problema termomecânico (apenas térmico mecânico, sem retorno).→→\rightarrow Para o problema mecânico, eu já entendi da resposta de Geoff , que precisarei usar o solucionador iterativo devido ao...
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...