Ciência computacional

Não vejo menção a isso nas listas de empregos. Já vi programação inteira mencionada, MIP, programação não-linear com números inteiros mistos, LP, programação dinâmica etc., mas sem SDP. É muito mais moderno na academia do que na indústria? Da minha exposição limitada a acadêmicos e participantes da indústria em sistemas …

10 convex-optimization  semidefinite-programming 

Na minha resposta a uma pergunta no MSE referente a uma simulação física hamiltoniana 2D, sugeri o uso de um integrador simplético de ordem superior . Então achei que seria uma boa ideia demonstrar os efeitos de diferentes etapas no tempo na precisão global de métodos com ordens diferentes, e …

10 python  time-integration 

O campo da dinâmica dos fluidos computacional (CFD) é dedicado à resolução das equações de Navier-Stokes (ou alguma simplificação delas). Um subconjunto de CFD, modelos oceânicos e atmosféricos resolve numericamente as mesmas equações para aplicações realistas. Quais são as diferenças e trade-offs entre as abordagens gerais de CFD e os …

10 fluid-dynamics  simulation  computational-physics  numerical-modelling 

Há algum tempo eu escrevi um código que tentava calcular o sem usar as funções da biblioteca. Ontem, eu estava revisando o código antigo e tentei torná-lo o mais rápido possível (e correto). Aqui está minha tentativa até agora:l o g( X )log(x)log(x) const double ee = exp(1); double series_ln_taylor(double …

10 numerical-analysis 

Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico do Computational Science Stack Exchange. Fechado há 5 anos . Existe um grau de complexidade maior que e menor que O ( n …

10 algorithms  complexity  efficiency 

Eu gostaria de usar o método de 8ª ordem Runge-Kutta (89) em um aplicativo de mecânica / astrodinâmica celeste, escrito em C ++, usando uma máquina Windows. Portanto, pergunto-me se alguém conhece uma boa biblioteca / implementação documentada e livre para usar? Tudo bem se estiver escrito em C, desde …

10 c++  reference-request  libraries  c  runge-kutta 

Por exemplo, as bibliotecas de matriz esparsa C ++ que eu usei - Eigen e SuiteSparse, elas parecem não ter nenhuma funcionalidade SVD para matriz esparsa. Então, apenas curioso, SVD é mais difícil que QR / LU para matriz esparsa?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Eu uso o FEA há alguns anos, mas usá-lo e usá-lo corretamente são duas coisas diferentes: o fator de segurança não é a solução para tudo. Tenho a sensação de que não vou usá-lo corretamente, a menos que tenha uma resposta clara para essa pergunta: Estou ciente de que os …

10 finite-element  simulation 

Costumo encontrar o ditado geral de que métodos de pontos interiores são difíceis de aquecer. Existe uma explicação intuitiva por trás desse conselho? Existem situações em que se pode esperar benefícios com o início quente em um método de ponto interior? Alguém pode recomendar algumas referências úteis sobre o assunto?

10 optimization  constrained-optimization 

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 )f(x)f(x) …

10 monte-carlo  nonlinear-equations  stochastic  numerics 

Estou tentando entender a diferença entre essas duas placas gráficas para computação acadêmica, especificamente para o componente DGEMM. Se olharmos para as estatísticas brutas, ambas têm o mesmo chip GK110, estatísticas comparáveis ​​em praticamente todas as categorias e, acredito, têm a mesma arquitetura principal. Antes de qualquer desconto, o K20X …

10 performance  gpu  efficiency 

Sou iniciante na FE. Minha aplicação é a precificação de derivativos financeiros em que o espaço é tridimensional. Então, adicionando tempo, o problema tem seis dimensões. Eu tentei olhar em volta (Fenics, escript, deal.II, ...), mas meu entendimento é que esses softwares são limitados a 3 + 1 (espaço 3d …

10 pde  finite-element 

Gostaria de saber como as condições de Dirichlet são normalmente aplicadas ao usar o método de volume finito em uma grade não uniforme centralizada em célula, Minha implementação atual simplesmente impõe a condição de contorno ao fixar o valor da primeira célula, ϕ1 1= gD( xeu)ϕ1=gD(xL) \phi_1 = g_D(x_L) onde …

10 boundary-conditions  finite-volume  poisson 

Estou trabalhando em um sistema de equações linear esparso grande e mal condicionado. Quero usar aritmética duplo-duplo ou aritmético quad-duplo para resolvê-los. Eu sei que existe um pacote chamado MPACK desenvolvido por Nakata, Maho, que pode executar cálculos algébricos lineares numéricos sob aritmética quádrupla. No entanto, ele é projetado para …

10 sparse-matrix  matrix  condition-number 

Estou resolvendo um sistema não linear de equações acopladas e calculei o jacobiano do sistema discretizado. O resultado é realmente complicado, abaixo estão (apenas!) As 3 primeiras colunas de uma matriz ,3 × 93×93\times 9 (A complexidade surge, em parte, porque o esquema numérico requer ajuste exponencial para estabilidade.) Eu …

10 symbolic-computation  jacobian  design-patterns  sympy