Ciência computacional

Estou trabalhando com a biblioteca OpenFOAM C ++ Computational Continuum Mechanics (ela pode lidar com interação fluido-sólido, fluxos MHD ...) que usa malhas não estruturadas arbitrárias. Isso foi motivado pela idéia de usar a vantagem da geração rápida (geralmente automática) de malhas não estruturadas para simular problemas em geometrias complexas. …

12 mesh-generation  fluid-dynamics 

Nos métodos de decomposição de domínio (DD) e multigrid (MG), pode-se compor a aplicação das atualizações de bloco ou correções grosseiras como aditivas ou multiplicativas . Para solucionadores pontuais, essa é a diferença entre as iterações de Jacobi e Gauss-Seidel. O mais suave multiplicativo para atuando como é aplicado comoS …

12 linear-algebra  performance  multigrid  domain-decomposition 

Em uma classe específica de detectores, nossos dados são apresentados como pares de pontos em duas dimensões, e queremos agrupar esses pontos em linhas. Os dados são barulhentos e são armazenados em uma direção, mas não na outra. Não podemos garantir um acerto em cada compartimento, mesmo quando cada elemento …

12 regression  computer-vision  pattern-recognition 

Estou resolvendo um problema físico usando esquema numérico implícito. Isso me leva a resolver uma equação linear com matriz tridiagonal. Eu codifiquei esse algoritmo da Wikipedia. Gostaria de saber se existe uma biblioteca eficiente que permita resolver esse tipo de equação de maneira otimizada. Uma observação importante é que a …

12 linear-algebra  libraries  c++ 

A diferenciação automática nos permite avaliar numericamente a derivada de um programa em uma entrada específica. Existe um teorema de que esse cálculo pode ser feito a um custo inferior a cinco vezes o custo para executar o programa original. Esse fator de cinco é um limite superior. Em que …

12 automatic-differentiation  efficiency 

Eu tenho usado o GSL como base de muitas de minhas simulações, mas é um pouco exagerado para meus propósitos e define seu próprio tipo complexo por motivos legados. Em vez de codificar meu próprio solucionador Runge-Kutta ODE, o que provavelmente não seria muito eficiente, há algum solucionador ODE de …

12 c  ode 

Os altos fluxos de número de Reynolds produzem camadas de limite muito finas. Se a resolução de parede for usada na simulação de redemoinho grande, a proporção da imagem pode ser da ordem de 10610610^6 . Muitos métodos se tornam instáveis ​​nesse regime porque a constante inf-sup se degrada como …

12 pde  finite-element  fluid-dynamics  stability  incompressible 

Eu tenho uma pergunta sobre o ajuste quadrático de um conjunto de pontos e normais correspondentes (ou equivalentes, tangentes). Ajustar superfícies quadráticas para apontar dados é bem explorado. Alguns trabalhos são os seguintes: Montagem direta com restrição de tipo de superfícies quadráticas , James Andrews, Carlo H. Sequin Design e …

12 computational-geometry  regression  geometry  curve-fitting 

Estou interessado em maximizar globalmente uma função de muitos ( ) parâmetros reais (resultado de uma simulação complexa). No entanto, a função em questão é relativamente cara de avaliar, exigindo cerca de 2 dias para cada conjunto de parâmetros. Estou comparando opções diferentes e queria saber se alguém tinha sugestões.≈ …

12 optimization 

Gostaria de conhecer um bom método para interpolar dados entre duas grades não estruturadas, onde uma grade é uma versão mais grossa da outra. A eficiência é muito importante para mim, pois estou resolvendo um problema temporário de PDE em que preciso transferir dados entre as grades a cada etapa …

12 interpolation 

Eu atualmente tentando calcular de forma barata uma estimativa boa classificação para uma matriz . Portanto, eu calculo uma decomposição QR dinâmica de colunas usandoUMAAA [Q,R,E]=qr(A) no Matlab. Estimo a classificação de usandoUMAAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Isso funciona bem e um gráfico sobre todas as entradas diagonais …

12 matlab  lapack 

Dadas duas matrizes e , eu gostaria de encontrar os vetores e , de modo que, Em forma de matriz, estou tentando minimizar a norma Frobenius de A - \ mbox {diag} (x) \ cdot B \ cdot \ mbox {diag} (y) = A - B \ circ (xy ^ …

12 linear-algebra  matrix  reference-request 

As proteínas são tipicamente representadas em forma de desenho animado, com folhas β como flechas e hélices α como bobinas: Gostaria de saber, existe em algum lugar uma referência que descreva a construção dessa representação? Ou seja, que objetos matemáticos são usados ​​para construir esses gráficos e em quais átomos …

12 computational-chemistry  visualization 

Preciso resolver o mesmo sistema linear esparso (300x300 a 1000x1000) com muitos lados do lado direito (300 a 1000). Além deste primeiro problema, eu também gostaria de resolver sistemas diferentes, mas com os mesmos elementos diferentes de zero (apenas valores diferentes), ou seja, muitos sistemas esparsos com padrão de esparsidade …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

Existem muitas filosofias em diferentes disciplinas de engenharia de software sobre como as bibliotecas devem lidar com erros ou outras condições excepcionais. Alguns dos que eu já vi: Retorne um código de erro com o resultado retornado por um argumento de ponteiro. É isso que o PETSc faz. Retornar erros …

11 software