Ciência computacional

Até onde eu sei, as duas grandes estruturas genéricas de software de ciências computacionais do Departamento de Energia dos EUA são PETSc e Trilinos . Eles parecem semelhantes à primeira vista, além das diferenças de linguagem (C versus C ++). Quais são as principais diferenças entre as duas estruturas e …

24 petsc  software  trilinos 

Ao passar da forma forte de uma PDE para a forma FEM, parece que sempre se deve fazer isso declarando primeiro a forma variacional. Para fazer isso, você multiplica a forma forte por um elemento em algum espaço (Sobolev) e integra-se na sua região. Isso eu posso aceitar. O que …

24 pde  finite-element  elliptic-pde 

Um problema muito comum na cadeia de Markov Monte Carlo envolve probabilidades de computação que são soma de grandes termos exponenciais, ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... onde os componentes de lata podem variar de muito pequeno a muito grande. Minha abordagem foi fatorar o maior termo exponencial para que:K …

24 monte-carlo  floating-point  statistics 

Não estou executando nenhum código paralelo agora, mas prevejo a execução de código paralelo no futuro usando um híbrido de OpenMP e MPI. Depuradores foram ferramentas inestimáveis ​​para mim ao executar projetos em série. Alguém pode recomendar um depurador paralelo (ou vários depuradores) para usar na depuração de software paralelo? …

24 software  parallel-computing 

Descrição da experiência: Na interpolação de Lagrange, a equação exata é amostrada em pontos (ordem polinomial ) e é interpolada em 101 pontos. Aqui varia de 2 a 64. Cada vez que , e gráficos de erro são preparados. Vê-se que, quando a função é amostrado em pontos equidistantes, o …

24 finite-difference  interpolation  error-estimation  polynomials 

Pelo que entendi, existem duas categorias principais de métodos iterativos para resolver sistemas lineares de equações: Métodos estacionários (Jacobi, Gauss-Seidel, SOR, Multigrid) Métodos de Krylov Subespaço (Gradiente Conjugado, GMRES, etc.) Entendo que a maioria dos métodos estacionários funciona relaxando (suavizando) iterativamente os modos de Fourier do erro. Pelo que entendi, …

24 linear-algebra  convergence  krylov-method  conjugate-gradient 

A reprodutibilidade se tornou cada vez mais importante na pesquisa em ciências computacionais. (Por exemplo, consulte este artigo de Roger Peng na Science ; também conheço outros artigos e sites da Web.) No entanto, não está claro para mim quanta informação devo incluir em um artigo de revista (ou online) …

23 software  publications  documentation  reproducibility 

Um cluster de computação do Amazon EC2 custa entre US $ 800 e US $ 1000 (dependendo do ciclo de trabalho) por núcleo de CPU física ao longo de três anos. Em nossa última rodada de aquisição de hardware, meu laboratório pegou 48 núcleos no valor de hardware muito semelhantes …

23 cloud-computing  hpc 

Muitas abordagens numéricas do CFD podem ser estendidas para ordem arbitrariamente alta (por exemplo, métodos descontínuos de Galerkin, métodos WENO, diferenciação espectral etc.). Como devo escolher uma ordem de precisão apropriada para um determinado problema?

23 fluid-dynamics  accuracy 

Qual é o estado da arte na aproximação de integrais altamente oscilatórias em uma dimensão e dimensões superiores à precisão arbitrária?

23 quadrature  precision  numerical-analysis 

Eu tenho uma série de pontos de dados que espero (aproximadamente) seguir uma função que assintota uma linha em geral . Essencialmente, aproxima de zero como , e provavelmente o mesmo pode ser dito de todas as derivadas , , etc. Mas não sei qual é a forma funcional de …

23 algorithms 

Esta pergunta foi migrada da Cross Validated porque pode ser respondida no Computational Science Stack Exchange. Migrou há 7 anos . Quero aprender a programação do MATLAB para que eu possa realizar algumas pesquisas / análises por conta própria e também, para que eu possa estudar / modificar alguns scripts …

23 software  matlab  languages  octave 

Atualmente, essa questão não se encaixa no nosso formato de perguntas e respostas. Esperamos que as respostas sejam apoiadas por fatos, referências ou conhecimentos, mas essa pergunta provavelmente solicitará debates, argumentos, pesquisas ou discussões prolongadas. Se você acha que essa pergunta pode ser melhorada e possivelmente reaberta, visite o centro …

22 languages 

O que seria uma boa discretização de diferenças finitas para a seguinte equação: ∂ρ∂t+ ∇ ⋅ ( ρ u ) = 0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? Podemos pegar o caso 1D: ∂ρ∂t+ ddx( ρ u ) = 0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho u\right)=0 Por alguma razão, …

22 finite-difference  fluid-dynamics 

Estou usando o SNES do pacote de resolução não-linear do PETSc para resolver um sistema de equações não-lineares obtido discretizando uma equação diferencial parcial. Como posso determinar por que o solucionador não está convergindo e o que posso fazer para resolver com êxito minhas equações?

22 petsc  pde  implicit-methods