A pesquisa reproduzível em computação visa disponibilizar a outros pesquisadores o código necessário para gerar os resultados em um artigo computacional, para que eles possam executar esse código para reproduzir os resultados nesse artigo. Gostaria de tornar toda a minha pesquisa reproduzível, mas...
é umamatriz esparsa n × n simétrica positiva positiva definida (SPD). G é uma matriz diagonal esparsa. n é grande (UMAUMAAn × nn×nn \times nGGGnnn > 10000) e o número de não-zeros no G é geralmente 100 ~ 1000.nnnGGG foi fatorado na forma cholesky comoUMAUMAA .L D LTeuDeuTLDL^T Como atualizar...
Estou com um pouco de dificuldade em tentar entender um artigo. O artigo usa o método espectral para resolver um autovalor proveniente de um sistema de ODEs acoplados. Escreverei apenas uma equação agora, porque é suficiente para chegar ao cerne da (s) minha (s) pergunta (s). A equação é V[ r ]...
Para a solução de grandes sistemas lineares usando métodos iterativos, é frequentemente interessante introduzir o pré-condicionamento, por exemplo, resolva M - 1 ( A x = b ) , onde M é usado aqui para o pré-condicionamento à esquerda do sistema. Normalmente, devemos ter esse M - 1 ≈ A - 1 e...
Suponha que eu tenha uma equação de matriz definida recursivamente como A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Então a equação para A [1] se parece com uma fração contínua, para a qual existem métodos altamente eficientes que evitam o recálculo tedioso (consulte "Receitas numéricas" para...
Quantidades físicas como pressão, densidade, energia, temperatura e concentração devem sempre ser positivas, mas os métodos numéricos às vezes calculam valores negativos durante o processo de solução. Isso não está certo porque as equações calcularão valores complexos ou infinitos (normalmente...
Na maioria das (todas?) Implementações do Fast Multipole Method (FMM), os octrees são usados para decompor o domínio relevante. Teoricamente, os octrees fornecem um limite volumétrico simples, útil para provar o tempo de execução O (n) de um FMM. Além dessa lógica teórica, há benefícios em usar...
Eu gostaria de usar f2pycom o Fortran moderno. Em particular, estou tentando fazer o exemplo básico a seguir funcionar. Este é o menor exemplo útil que eu poderia gerar. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out)...
Estou tentando implementar a seguinte função no ponto flutuante de precisão dupla com baixo erro relativo : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log(exp(x)+exp(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Isso é usado extensivamente em aplicativos estatísticos para adicionar...
Gostaria de aprender como os métodos wavelet podem ser aplicados ao PDE, mas infelizmente não conheço um bom recurso para aprender sobre este tópico. Parece que muitas introduções às wavelets se concentram na teoria da interpolação, por exemplo, na montagem de um sinal por uma superposição de...
Eu quero desenvolver um software de computação científica paralelo do zero. Eu quero alguns pensamentos sobre qual idioma começar. O programa envolve ler / gravar dados em arquivos txt e realizar cálculos pesados em paralelo, com muitas fatorações de LU e o uso de solucionadores lineares...
Quais são os melhores algoritmos (escalabilidade e eficiência) para gerar malhas quadradas não estruturadas em 2D? Onde posso encontrar um bom gerador de malha quadrada não estruturado? (de código aberto
Estou tentando resolver alguns problemas de otimização não-linear sem restrições na GPU (CUDA). A função objetivo é uma função não-linear suave e seu gradiente é relativamente barato para calcular analiticamente, portanto, não preciso me preocupar com a aproximação numérica. Quero resolver esse...
Enquanto lia alguma literatura sobre os solucionadores de PDE, me deparei com o termo pseudo passo a passo hoje. Parece ser um termo comum, no entanto, não consegui encontrar uma boa definição ou um artigo introdutório para ele. Portanto: o que é pseudo-time-stepping, e como é geralmente...
Deseja melhorar este post? Forneça respostas detalhadas para esta pergunta, incluindo citações e uma explicação de por que sua resposta está correta. Respostas sem detalhes suficientes podem ser editadas ou excluídas. Refinamento de malha adaptável (AMR) é uma técnica...
Estou interessado nos casos em que o gradiente conjugado funciona muito melhor que o método GMRES. Em geral, o CG é a escolha preferível em muitos casos de SPD (simétrico-positivo-definido) porque requer menos armazenamento e o limite teórico da taxa de convergência para o CG é o dobro do GMRES....
Eu tenho uma função 2D não-convexa limitada que eu gostaria de encontrar o mínimo de. A função é bastante suave. Avaliando é caro. Um erro aceitável é de cerca de 3% do domínio da função em cada eixo. Tentei executar a implementação do algoritmo DIRECT na biblioteca NLOPT, mas não proporcionou uma...
Estou à procura de um software de geração de malha que é gratuito e de código aberto, fornece uma interface de script sã para especificação de domínio, trabalha para geometrias complexas, pode gerar malhas 2D e 3D, Que opções eu
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...
Suponha que eu tenho alguma função e quero encontrar tal que . Eu poderia usar o método Newton-Raphson. Mas isso requer que eu conheça a função derivada . Uma expressão analítica para pode estar indisponível. Por exemplo, pode ser definido por um código complicado de computador que consulta um...