Eu implementei um solucionador multigrid do ciclo V usando uma correção linear de defeitos (LDC) e um esquema de aproximação total (FAS). Meu problema é o seguinte: Usando o LDC, o resíduo é reduzido em um fator de ~ 0,03 por ciclo. A implementação do FAS também converge com um fator linear, mas o...
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...
Eu queria saber qual é a melhor maneira de aprender a programar no Fortran? Estive procurando livros, os exemplos que me pareciam bem foram "Fortran 95/2003 Explained" e "Fortran For Scientists And Engineers" . Não sei o que vocês acham disso? Outra coisa que eu estou pensando é sobre qual...
Estou escrevendo um pequeno pacote MATLAB que resolverá certos problemas numéricos de classe. Existem 3 estágios do algoritmo e o usuário tem 5 opções para cada estágio. Eu implementei todo o problema usando funções e 3 casos de switch (um para cada estágio do algoritmo). Funciona bem, mas estou...
Eu tenho tentado debugar este erro nos últimos dias que eu quis saber se alguém tem o conselho em como proceder. Estou resolvendo a equação de Poisson para uma distribuição de carga escalonada (um problema comum em física eletrostática / semicondutora) em uma malha de volume finito não uniforme,...
O:=minx∈Rnf(x).O:=minx∈Rnf(x).\mathcal{O} := \min_{x \in \mathbb{R}^n} f(x).xoptxoptx_\text{opt}x0x0x_0xopt.xopt.x_\text{opt}.xxxϵ−ϵ−\epsilon-OO\mathcal{O}||x−xopt||2||x0−xopt||2≤ϵ.||x−xopt||2||x0−xopt||2≤ϵ.\begin{equation} \frac{||x - x_{\text{opt}}||_2}{||x_0 - x_\text{opt}||_2} \leq...
Essas perguntas podem ser um pouco fora de tópico no comp-sci. se necessário, sugira onde ele se encaixa. A questão é sobre como gerenciar todas as simulações executadas com eficiência. digamos, por exemplo, uma simulação requer a fixação de 2 parâmetros que devem ser definidos em um determinado...
Eu tenho um grande problema de autovalor cúbico: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Eu poderia resolver isso convertendo para um problema de autovalor linear, mas resultaria em...
Quais são as principais diferenças entre o FEM e o XFEM? Quando (não) devemos usar o XFEM em vez do FEM e vice-versa? Em outras palavras, quando encontro um novo problema, como posso saber qual
Existe alguma regra simples para dizer se vale a pena fazer SOR em vez de Gauss-Seidel? (e possível maneira de estimar o parâmetro de realxação )ωω\omega Quero dizer, apenas olhando para a matriz , ou conhecimento de um problema específico que a matriz representa? Eu estava lendo a resposta para...
Muitas bibliotecas adaptáveis do FEM usam estruturas de dados de malha mais avançadas para lidar com a adição / remoção de nós, arestas, triângulos, tetraedros, etc. Por exemplo, a biblioteca p4est usa estruturas de dados octree para refinamento de malha adaptável; você normalmente não encontra...
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
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } O código...
Qual a diferença entre esses dois métodos? Um problema pode ser resolvido por um método capaz de ser resolvido pelo outro? Ambos / ou um deles podem ser paralelos ao OpenMP e / ou
Eu tenho um projeto em que preciso usar campos quadráticos Especificamente, números do formulário com .a+b−3−−−√a+b−3a + b \sqrt{-3}a , b ∈Qa,b∈Qa,b \in \mathbb{Q} Por exemplo, aqui estão os números primos nos números inteiros de Eisenstein : Eu não quero usar sálvia. Eu gostaria de escrever...
No relatório técnico sobre Galahad [1], os autores declaram, no contexto de problemas gerais de programação não linear, Em nossa opinião, nunca houve muita dúvida de que os métodos SQP [programação quadrática sequencial] seriam mais bem-sucedidos [do que os métodos Lagrangeanos Aumentados] a...
Estou realizando uma pesquisa de linha como parte de um algoritmo quase-Newton BFGS. Em uma etapa da pesquisa de linha, uso uma interpolação cúbica para me aproximar do minimizador local. Seja a função de interesse. Eu quero encontrar um tal que .x ∗ f ′ ( x ∗ ) ≈ 0f: R → R , f∈ C1 1f:R→R,f∈C1f :...
Desejo modelar uma vara de pescar (ou uma corda) juntando segmentos curtos. (Os segmentos podem ter comprimento (curto) igual, mas cada segmento deve ter sua própria massa individual.) Um segmento influenciará o próximo pelo torque entre os segmentos. Por enquanto, as juntas podem ser consideradas...
Estou interessado em calcular a solução de um sistema Lage de ODEs usando um método de krylov como em [1]. Esse método envolve funções relacionadas ao exponencial (as chamadas funções ). Consiste essencialmente em calcular a ação da função de matriz, construindo um subespaço de Krylov usando a...
Alguns livros antigos que eu vi dizem que o número mínimo de estágios de um método Runge-Kutta explícito de uma ordem especificada é desconhecido para ordens . Isso ainda é verdade?≥9≥9\geq 9 Quais bibliotecas existem para trabalhar automaticamente com os métodos Runge-Kutta de alta...