Tenho experiência em codificar OpenMP para máquinas de memória compartilhada (em C e FORTRAN) para realizar tarefas simples como adição de matriz, multiplicação etc. (Apenas para ver como ele compete com o LAPACK). Conheço o OpenMP o suficiente para realizar tarefas simples sem a necessidade de...
Olá comunidade scicomp, Trabalhei na área de algoritmos de grafos usando frameworks como NetworkX (Python), JUNG e YFiles (Java). Agora estou entrando na área da computação paralela e de alto desempenho. Para um novo projeto, estou procurando uma biblioteca de gráficos C ++ com os seguintes...
Simultaneamente à minha pesquisa sobre os métodos do subespaço de Krylov, tenho a opção de explorar a matemática por trás do HPC um passo à frente ou a teoria da computação (hardware, SO, compiladores etc.). Atualmente, eu sei que tanto o suficiente para apenas sobreviver. Por exemplo, eu sei como...
Pelo que entendi, o método multigrid resolve um sistema linear, resolvendo uma versão mais grossa do mesmo problema (eliminando o erro de baixa frequência) e projetando-se de volta à grade fina para suavizar os erros de alta frequência. Para sistemas grandes, posso ver como um método iterativo pode...
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...
Estou tentando diagonalizar algumas matrizes densas e mal condicionadas. Na precisão da máquina, os resultados são imprecisos (retornando valores próprios negativos, os vetores próprios não possuem as simetrias esperadas). Eu mudei para a função Eigensystem [] do Mathematica para tirar proveito da...
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,...
Recentemente, me deparei com um comentário afirmando que quase todos os pesquisadores que usam os métodos de Monte Carlo estão fazendo errado. Em seguida, foi elaborado que a simples escolha de sementes diferentes para diferentes instâncias de um PRNG, como o Mersenne Twister, não é suficiente para...
Atualmente, estou trabalhando na solução de sistemas simétricos muito grandes (mas não positivos definidos), gerados por alguns algoritmos. Essas matrizes têm uma boa esparsidade de bloco que pode ser usada para solução paralela. Mas não consigo decidir se devo usar uma abordagem direta (como...
Estou procurando portar algum código que resolva um conjunto de equações diferenciais parciais (PDE) pelo método de volume finito no formato IMPLICIT (para discretização do tempo). Como resultado, existe um sistema tridiagonal de equações nas direções x, y, z, que é tratado pelo esquema ADI /...
Eu já desenvolvi uma solução funcional do Método dos Elementos Finitos para resolver problemas de transferência de calor usando GPU e OpenCL usando o método Gradiente Conjugado. A principal desvantagem desse método é a alta demanda por memória. Além disso, no caso de placas gráficas, a memória...
Os algoritmos Parareal, PITA e PFASST são técnicas de todo o domínio para paralelizar a solução de problemas dependentes do tempo no tempo. Quais são os princípios orientadores por trás desses métodos? Quais são as principais diferenças entre eles? Posso dizer que um é baseado no outro? Como? E...
Estou resolvendo um problema de várias escalas usando o HMM (Heterogeneous Multiscale Method) . Essencialmente, meu procedimento específico usa o seguinte processo iterativo: Resolva muitos sistemas matriciais locais. Calcule um valor de interesse das soluções dos sistemas locais. Montar um...
Sempre ouvi dizer que a fácil paralelização era uma das vantagens dos métodos DG, mas não vejo realmente por que nenhuma dessas razões também se aplica ao Galerkin
Gostaria de saber se existe uma lista de discussão muito boa ou grupos do Google para a Computational Science, onde podemos discutir perguntas em vez de apenas fazer e responder perguntas. De fato, estou mais interessado em computação paralela e soluções numéricas de PDEs. Mas não sei o que e como...
Precisamos calcular matrizes de covariância com tamanhos que variam de a . Temos acesso a GPUs e clusters, nos perguntamos qual é a melhor abordagem paralela para acelerar esses cálculos.100000 × 10000010000 × 1000010000×1000010000\times10000100000 ×
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...
Normalmente, quando eu proponho um gerador seqüencial de números aleatórios em C, uso a chamada srand(time(NULL)) então use rand() mod N para obter um número inteiro aleatório entre 0 e N-1. No entanto, quando faço isso em paralelo, as chamadas ao tempo (NULL) são tão próximas uma da outra...
Escrevo códigos em Fortran e C para vários algoritmos matriciais. No entanto, quando eu perfil meus códigos usando o VTune, geralmente encontro alguma terminologia que não consigo apreciar completamente. Existe um bom recurso para aprender a criar perfis do zero para um nível bastante...
Alguém conhece referências que resumem a história e as ideias por trás da supercomputação, incluindo menção de desenvolvimentos em linguagens de programação paralela, aplicativos, empresas iniciantes (algumas foram estabelecidas, mas não sobreviveram), necessidades da indústria...