Essa pergunta foi feita um bilhão de vezes no Stackoverflow, no entanto, o foco sempre foi a codificação não numérica. Estou procurando um projeto para contribuir dentro dos limites da computação numérica e da computação de alto desempenho. Idealmente, eu preferiria um projeto pequeno, mas isso não...
Para aplicativos que exigem recursos computacionais significativos, o alto desempenho pode ser um fator crítico quando se trata de fornecer resultados científicos ou obter "avanços" em tempo razoável. Quanto tempo e esforço os desenvolvedores de software devem investir na otimização de um...
Vamos começar com um problema do formulário (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 com um conjunto de determinadas condições de contorno ( Dirichlet , Neumann , Robin , Periodic , Bloch-Periodic ). Isso corresponde à localização dos valores próprios e dos vetores próprios para algum operador ,...
Estou tentando resolver a equação Ideal MHD usando métodos semi-discretos, reconstruções espaciais ENO e escalonamento de tempo TVD RK. Estou recebendo diferentes soluções de estado estacionário com ordem temporal diferente. Está
Estou interessado em implementar uma malha móvel para um problema de difusão de advecção. Os métodos de malha móvel adaptativa fornecem um bom exemplo de como fazer isso para a equação de Burger em 1D usando diferença finita. Alguém seria capaz de oferecer um exemplo elaborado para resolver a...
Vamos considerar uma condição inicial suave e a equação do calor em uma dimensão: no intervalo aberto , e vamos assumir que queremos resolvê-lo numericamente com diferenças finitas.∂tu=∂xxu∂tu=∂xxu \partial_t u = \partial_{xx} u]0,1[]0,1[]0,1[ Eu sei que para o meu problema a ser bem colocada,...
Os estudos de caso da minha tese de doutorado exigem que eu tenha extensões de intervalo das sub - rotinas Fortran no CHEMKIN-II (desculpas pelo link; é o melhor que eu poderia encontrar para um pacote que não é mais distribuído pelo Sandia National Laboratory). CHEMKIN-II é um pacote químico usado...
Estou trabalhando na resolução das equações de poroelasticidade unidimensional acopladas (modelo de biot), dadas como: −(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} =
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,...
Eu tenho essa confusão sobre a regra Armijo usada na pesquisa de linha. Eu estava lendo a pesquisa de linha de rastreamento anterior, mas não entendi o que é essa regra Armijo. Alguém pode elaborar qual é a regra de Armijo? A Wikipédia não parece explicar bem.
Normalmente, escrevo código serial e, quando o faço, escrevo testes de unidade com alguma estrutura de teste no estilo xUnit (MATLAB xUnit, PyUnit / nose ou estrutura de teste C ++ do Google). Com base em uma pesquisa superficial do Google, não tenho visto muito sobre como os profissionais...
Estou desenvolvendo uma simulação FEM. Para os primeiros testes, usarei o mesher auto-escrito simples e a visualização do gráfico de malha. Mas quero preparar meu programa para usar os dados gerados por uma malha existente e enviá-los para as ferramentas de visualização existentes. Existe um...
Gostaria de modelar o fluxo laminar de água das raízes ao caule de uma planta. No final das raízes, os tubos variam de milímetro a centímetro de diâmetro e comprimento. À medida que nos aproximamos do caule, as raízes aumentam em comprimento e diâmetro. Quero criar domínios aleatórios em 3D que...
Eu tenho usado métodos diferentes para calcular a classificação de uma matriz e a solução de um sistema matricial de equações. Me deparei com a função linalg.svd. Comparando isso com o meu próprio esforço de resolver o sistema com a Eliminação Gaussiana, parece ser mais rápido e preciso. Estou...
Não estou familiarizado com o formato de arquivo usado no HDF5, mas gostaria de saber se os arquivos HDF5 são adequados para o controle de revisão com o git (ou, por exemplo, Mercurial ou Subversion)? Acho que o que quero dizer é: os arquivos HDF5 são adequados para diferenças baseadas em linhas ou...
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
Suponha que eu tenha a seguinte função interessante: Tem algumas propriedades desagradáveis, como sua derivada não sendo contínua em múltiplos racionais de . Eu suspeito que um formulário fechado não existe.f(x)=∑k≥1coskxk2(2−coskx).f(x)=∑k≥1coskxk2(2−coskx). f(x) = \sum_{k\geq1} \frac{\cos k...
TL, DR Qual é a melhor prática aceita nos círculos científicos da computação para armazenar grandes quantidades de dados hierarquicamente estruturados? Por exemplo, o SQL não funciona bem com grandes matrizes esparsas. Existe uma boa ferramenta disponível para estruturar, armazenar e analisar esse...
Segundo a Wikipedia, a taxa de convergência é expressa como uma proporção específica de normas de vetores. Estou tentando entender a diferença entre taxas "lineares" e "quadráticas", em diferentes pontos do tempo (basicamente, "no início" da iteração e "no final"). Pode-se afirmar que: ek +...
Eu estou querendo saber se o algoritmo Thomas é a maneira mais rápida (provàvelmente?) De resolver um sistema tridiagonal esparso com dominação diagonal simétrica em termos de complexidade algorítmica (sem procurar pacotes de implementação como LAPACK etc.). Eu sei que o algoritmo Thomas e o...