No meu campo de pesquisa, a especificação de erros experimentais é comumente aceita e as publicações que não os fornecem são altamente criticadas. Ao mesmo tempo, muitas vezes acho que os resultados de cálculos numéricos são fornecidos sem nenhuma explicação de erros numéricos, mesmo que (ou talvez...
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...
Estou paralelizando código para resolver numericamente um modelo de equilíbrio populacional em 5 dimensões. Atualmente, tenho um código paralelo MPICH2 muito bom no FORTRAN, mas à medida que aumentamos os valores dos parâmetros, as matrizes ficam muito grandes para serem executadas no modo de...
Quando temos um modelo numérico que representa um sistema físico real e que exibe caos (por exemplo, modelos de dinâmica de fluidos, modelos climáticos), como podemos saber que o modelo está funcionando como deveria? Não podemos comparar diretamente dois conjuntos de resultados do modelo, porque...
Eu vejo muito poucas bibliotecas / pacotes de computação de ponto não flutuante. Dadas as várias imprecisões da representação de ponto flutuante, surge a questão de por que não existem pelo menos alguns campos em que essa precisão aumentada pode valer os meandros de trabalhar com ponto...
Existem muitos métodos numéricos bem conhecidos para resolver equações do tipo por exemplo, método de bissecção, método de Newton, etc.f( x ) = 0 ,x ∈ Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, Na minha aplicação, é calculado com um método estocástico (o resultado é uma média).f( X...
Estou trabalhando com funções que, em geral, são muito mais suaves e se comportam melhor no espaço de log-log - então é aí que eu executo interpolação / extrapolação, etc., e isso funciona muito bem. Existe uma maneira de integrar essas funções numéricas no espaço de log-log? ou seja, espero usar...
Cada matriz real pode ser reduzir a verdadeira forma Schur T = U T A L usando um ortogonal similiary transformar U . Aqui a matriz T é uma forma quase triangular com 1 por 1 ou 2 por 2 blocos na diagonal principal. Cada um por um do bloco corresponde a um valor próprio real Um e cada dois por dois...
Os algoritmos existentes para resolver EDOs manipulam funções , em que . Mas em muitos sistemas físicos, a equação diferencial é autônoma, então , , com deixado de fora. Com essa suposição simplificadora, que melhorias podem ser vistas nos métodos numéricos existentes? Por exemplo, se , o problema...
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
Antecedentes: Eu criei apenas uma solução numérica para o 2d Navier-Stokes, para um curso. Era uma solução para o fluxo da cavidade acionado pela tampa. O curso, no entanto, discutiu vários esquemas para discretizações espaciais e discretizações de tempo. Também fiz mais cursos de manipulação de...
Sou formado em ciências da computação e matemática. Eu amo os dois assuntos. Estou pensando em fazer uma carreira de pós-graduação, talvez em computação científica. Qual é a diferença real entre computação científica e análise numérica? Eles são estudados como
De acordo com Wolfram Alpha e o sistema de álgebra computacional Sage , a seguinte identidade é válida: porque( arctan( l1- l2d) ) = 11 + ( l1- l2)2d2--------√porque(arctan(eu1-eu2d))=11+(eu1-eu2)2d2 \cos\left(\arctan\left(\frac{l_1-l_2}{d}\right)\right) = \frac{1}{\sqrt{1 +...
Gostaria de calcular, como parte da solução da equação de Laplace, usando o Método Multipolar Rápido, a segunda derivada das funções legendre associadas do primeiro tipo . Especificamente, estou procurando implementações em C ou apenas a relação de recorrência correta a ser usada para escrever a...
Preciso calcular a integral Eu= ∫R0 0f( R ) Jn( zn mrR) rdrI=∫0Rf(r)Jn(znmrR)rdrI = \int_0^R f(r)J_n\left(\frac{z_{nm}r}{R}\right)rdr onde é a função de Bessel da ordem do primeiro tipo, é seu m t h zero ef ( r ) é uma função real que é um pouco semelhante a J n (mas não do mesmo modo, é bastante...