Estou mergulhando no mundo fascinante da análise de elementos finitos e gostaria de resolver um grande problema termomecânico (apenas térmico mecânico, sem retorno).→→\rightarrow Para o problema mecânico, eu já entendi da resposta de Geoff , que precisarei usar o solucionador iterativo devido ao...
Estou tentando otimizar um distribuidor de fluxo em um tanque, para que a distribuição de velocidade e temperatura em qualquer seção transversal seja relativamente uniforme. Existem muitos parâmetros que posso ajustar para a máxima uniformidade da seção transversal, como o número de tubos de...
Suponha que eu tenha o sistema linear grande e esparso original: . Agora, eu não tenho pois A é muito grande para fatorar ou qualquer tipo de decomposição de , mas suponha que eu tenha a solução encontrada com uma solução iterativa.A - 1 A x 0A x0 0= b0 0Ax0=b0A\textbf{x}_0=\textbf{b}_0UMA-...
Gostaria de prever tempos de execução para operações de álgebra linear densa em uma arquitetura específica usando uma biblioteca específica. Eu gostaria de aprender um modelo que se aproxima da função Fo p: :Fop::F_{op} \;::\; tamanhos de entrada tempo de execução→→ \rightarrow para operações...
Como exemplificado pela resposta de Jed Brown a Custos de pesquisas versus cálculos , o uso de operações de ponto flutuante vetorizado versus não vetorizado resulta em código muito mais rápido. Muitos compiladores modernos afirmam que podem executar a vetorização automática. Como faço para que...
O Answer fornece um software para calcular o expoente Lyapunov condicional (CLE) para osciladores acoplados na sincronização do caos. No entanto, é difícil de seguir e não há saída gráfica do gráfico (e é em C, mais complexo). Alguém sabe como modificar a LET Toolbox, que é muito boa para o sistema...
Sabe-se que o método de Newton para resolver equações não lineares converge quadraticamente quando o palpite inicial é "suficientemente próximo" da solução. O que é "suficientemente próximo"? Existe literatura sobre a estrutura dessa bacia de
Gostaria de saber se existe um algoritmo que, dado um conjunto de pontos e um ângulo, calcule o casco convexo se o ângulo for e se calcule um envelope que segue mais de perto o "perímetro "α = 0α=0 0\alpha = 0α > 0α>0 0\alpha > 0 E se houver uma definição de um perímetro sem interseção de...
Folheando alguns livros didáticos, notei que o problema de agrupar inicialmente um mínimo durante uma pesquisa de linha tende a ser uma reflexão tardia (pelo menos nos meus textos de graduação). Existem técnicas ou práticas recomendadas bem estabelecidas para esse tipo de problema ou as soluções...
Quero minimizar uma função objetiva complicada e não tenho certeza se é convexa. Existe um bom algoritmo que tenta provar que não é convexo? É claro que o algoritmo poderia falhar em provar isso; nesse caso, eu não saberia se é convexo ou não, e isso está correto; Eu só quero tentar descartar a...
Tenho a impressão, de recursos muito diferentes e conversas com pesquisas, de que há uma demanda crescente por cálculos de alta precisão em equações diferenciais parciais numéricas. Aqui, alta precisão significa mais precisão do que apenas a precisão dupla de 64 bits padrão. Eu me pergunto sobre o...
Eu tenho um problema em que preciso encontrar todos os pares de autovalores positivos (como o valor próprio é positivo) de uma matriz não simétrica pequena (geralmente menor que 60x60). Posso parar de calcular quando o valor próprio é menor que um determinado limite. Eu sei que os autovalores são...
Estou tentando implementar o método Broyden-Fletcher-Goldfarb-Shanno para encontrar o mínimo de uma função. Preciso de duas suposições iniciais & e uma aproximação inicial da Matriz Hessiana . Os únicos requisitos que encontro para são que, se o Hessian é definido simétrico positivo, também...
Existe uma maneira preferida de como implementar uma avaliação rápida (aproximada) do polinômio de interpolação de Chebyshev em grade uniforme (considerando os valores de função nos nós de Chebyshev)? Meu problema é que a interpolação se torna lenta quando o grau do polinômio interpolador...
A maior parte da minha programação são códigos de pesquisa pontuais em C para uso próprio. Nunca distribuí nenhum código para outros colaboradores que não sejam colaboradores próximos. Eu desenvolvi um algoritmo que estou publicando em uma revista científica. Desejo fornecer o código fonte e talvez...
Realizando uma simulação numérica direta de turbulência isotrópica com o método pseudo-espectral de Fourier (Orzag & Patterson, PRL, 1972) usando FFT. Para obter um plano de fundo do método, amplamente utilizado na comunidade de turbulências, você pode ver este curso:...
Quero resolver onde K é minha matriz de rigidez. No entanto, algumas restrições podem estar ausentes e, portanto, algum movimento corporal rígido ainda pode estar presente no sistema (devido ao valor próprio zero). Como estou usando o CG para resolver o sistema linear, isso não é aceitável, já que...
A partir da definição do número da condição, parece que é necessária uma inversão da matriz para calculá-la. Será que para uma matriz quadrada genérica (ou melhor, se a definição positiva simétrica positiva) é possível explorar alguma decomposição da matriz para calcular o número da condição em um...
Estou usando a caixa de ferramentas Matlab pde para resolver uma certa equação elíptica em 2D. A solução está correta, embora eu precise plotá-la ao longo de uma determinada linha, ou seja, cortar uma fatia plana da malha 3D que representa a solução. Não consigo descobrir uma maneira que envolva...
Para o meu trabalho de pesquisa, desenvolvi uma pequena biblioteca C ++ destinada a facilitar a comunicação entre os códigos computacionais C ++ e o Octave / Matlab (quando este último é usado para fins de pós-processamento). Gostaria de liberar essa biblioteca sob a GPL em alguns sites de...