Recentemente, vi esta pergunta: como medir o erro de um método de diferenças finitas Eu sou estudante de ciências da simulação e, infelizmente, para mim, não está totalmente claro qual norma usar em que contexto. Frequentemente, usamos a norma euclidiana ou a norma L2, mas por que escolher normas...
Eu tenho parte de um problema descrito pela equação de conservação do momento: ∂ρ∂t+ 1pecadoθ∂∂θ( ρ u sinθ ) = 0∂ρ∂t+1sinθ∂∂θ(ρusinθ)=0\frac{\partial \rho}{\partial t} + \frac{1}{\sin\theta} \frac{\partial}{\partial \theta}(\rho u \sin \theta) =0 Onde e ρ = f ( θ , t ) (velocidade constante).u...
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...
Eu tenho um monte de processos MPI trocando mensagens de solicitação para frente e para trás. Os processos não sabem quais outros processos enviarão a eles mensagens ou quantos. Diante dessa situação, quero uma maneira eficiente de saber se todos os outros processos se consideram concluídos no...
A equação do calor é discretizada no espaço com FV (ou FEM), e uma equação semi-discreta é obtida (sistema de EDOs). Essa abordagem, conhecida como método de linhas , permite alternar facilmente de uma discretização temporal para outra, sem duplicação de código. Em particular, você pode reutilizar...
Peço desculpas antecipadamente se esta pergunta é boba. Eu preciso calcular a raiz de u - f( u ) = 0você-f(você)=0 0\begin{equation} u -f(u) =0 \end{equation} Onde é um vetor real ef ( u ) é uma função com valor vetorial real. Comecei com o método de Newton (que funcionava), mas depois percebi...
Eu tenho lido sobre a implementação de multiplicação densa de matrizes quando a matriz não se encaixa no cache. Um dos gráficos que eu vi (slide 9 desses slides ) mostra quedas repentinas no desempenho usando o algoritmo ingênuo. Essas quedas atingem cerca de 50% da velocidade, ocorrem enquanto a...
Recentemente, herdei um grande corpo de código legado que resolve um problema muito rígido e transitório. Eu gostaria de demonstrar que os tamanhos de etapas espaciais e temporais são pequenos o suficiente para que a natureza qualitativa da solução computada não seja alterada se forem diminuídas....
Quero acelerar meu programa python usando cython, mas meu loop interno ainda está fazendo chamadas lentas de python para o gerador de números aleatórios! Vários anos atrás, esse mesmo problema foi levantado por alguém que apoiava os sábios e parecia não haver uma boa solução naquele momento. Não é...
Supondo que uma simulação numérica direta seja realizada, qual é um bom método para estimar o impacto de pequena escala em larga escala na dinâmica de fluidos? Por exemplo, é pertinente comparar duas séries com tamanho de grade diferente ou duas séries com viscosidade diferente? Existem algumas...
Eu preciso resolver numericamente o sistema de duas equações diferenciais parciais acopladas. ∂x1∂t∂x2∂t=c1∇2x1+f1(x1,x2)=c2∇2x2+K∂x1∂t∂x1∂t=c1∇2x1+f1(x1,x2)∂x2∂t=c2∇2x2+K∂x1∂t\begin{align} \frac{\partial x_1}{\partial t} &= c_1\nabla ^2 x_1 + f_1(x_1,x_2)\\ \frac{\partial x_2}{\partial t} &=...
Estou tentando aprender mais sobre a malha de volume e decidi tentar implementar uma simples malha de volume. A estratégia que escolhi é subdividir meu espaço usando uma octree, refinada com base em alguns critérios. Verifiquei se meu octree está equilibrado, de modo que a diferença entre duas...
Existe uma implementação de código aberto e segura de thread de solucionadores de otimização convexos em C / C ++? Algumas bibliotecas como NLopt , Ipopt , OPT ++ não atendem aos meus requisitos. O OPT ++ e o Ipopt não são seguros para threads, e o NLopt não parece ter uma maneira específica /...
Eu tenho duas variáveis ke tcomo funções de duas outras variáveis p1e p2. Eu também sei seus valores máximos. Não tenho expressão analítica para isso. Quero encontrar os valores de ke tquais são os mais próximos de seus valores máximos. Existe uma maneira de otimizar o k = f1(p1, p2)e t = f2(p1,...
Estou procurando exemplos de equações de Helmholtz e Biharmonic em coordenadas cartesianas com soluções exatas, a fim de comparar minhas soluções numéricas com ela. Consegui encontrar alguns exemplos na internet, onde o problema com as condições de contorno foi definido com precisão. Infelizmente,...
Estou tentando reimplementar redes neurais em Python. Eu implementei a função de custo e o algoritmo de retropropagação corretamente. Eu os verifiquei executando seu código equivalente do Octave. Mas quando estou tentando usar a scipy.optimize.fmin_cgfunção, as iterações levam muito tempo para...
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...
Implementei um Navier Stokes baseado na discretização FEM e PETSc para resolver o sistema linear de equações. Para criar um procedimento de solução eficiente, sigo o artigo "Pré-condicionamento eficiente das equações de Navier-Stokes linearizadas para fluxo incompressível" (Silvester et al.), Que...
Se eu quiser usar uma matriz grande, como mpz_t A [100000], recebi "Falha na segmentação (core despejado)" durante minha compilação. Existe alguma maneira mais fácil de resolver
Estou tentando entender o algoritmo de rebinning da integração de Monte Carlo do VEGAS ( publicação original ( pré - impressão de LKlevin) e notas de implementação ). Vou tentar explicar primeiro o que acho que entendi e depois fazer minhas perguntas. Por simplicidade, vamos assumir que temos uma...