Ciência computacional

Por favor, desculpe a pergunta demorada, ela só precisa de explicações para chegar ao problema real. Aqueles familiarizados com os algoritmos mencionados provavelmente poderiam pular diretamente para o primeiro tablau simplex. Para resolver problemas de desvio mínimo absoluto (também conhecido como otimização ), o algoritmo Barrodale-Roberts é um método simplex …

9 optimization  linear-programming 

Atualmente, estou tentando resolver um problema de minimização com restrições não-lineares, conforme implementado na função matlab "fmincon". Minhas expectativas são: minimizar (fun1, x0, uB, lB, fun2) onde x0 é o estado inicial, fun1 é uma função que precisa ser minimizada, uB são limites superiores, lB são limites inferiores e lB …

9 c++  nonlinear-programming 

Então, eu estou debatendo se devo ou não aprender o Python. Ao falar com meus professores, o Matlab parece ser a linguagem comum usada na matemática aplicada / ciência da computação no que diz respeito à academia ; enquanto na indústria , meus professores (especialmente aqueles que trabalharam na indústria) …

9 matlab  python  c++ 

Quando é preferível usar os polinômios de Bernstein para aproximar uma função contínua, em vez de usar os únicos métodos preliminares de Análise Numérica: "Polinômios de Lagrange", "Operadores de diferenças finitas simples". A questão é sobre comparar esses métodos.

9 finite-element  finite-difference  interpolation 

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 maneira mais …

9 linear-algebra  matrix  condition-number 

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 aumenta. As seguintes idéias vieram à …

9 interpolation  fourier-analysis  fftw 

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 …

9 pde  matlab 

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 deve . Olhando para a Wikipedia, vejo …

9 optimization  algorithms 

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 hospedagem …

9 matlab  libraries  c++ 

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 às …

9 finite-element  iterative-method  conjugate-gradient 

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: http://www.math.ualberta.ca/~bowman/m655/lab3d.pdf Usando as chamadas regras para negociação, que consiste em …

9 fluid-dynamics  simulation 

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 …

9 performance  eigensystem  gpu 

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 o código executável no …

9 linear-solver  libraries  c++  c 

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 estado da …

9 reference-request  precision 

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 geralmente dependem de aplicativos? …

9 optimization  reference-request