Ciência computacional

Eu tenho um problema físico governado pela equação de Poisson em duas dimensões Tenho medidas dos dois componentes do gradiente ∂ u / ∂ x e ∂ u / ∂ y ao longo de alguma parte do limite, Γ m , então gostaria de impor ∂ u−∇2u=f(x,y),inΩ−∇2u=f(x,y),inΩ -\nabla^2 u = …

11 finite-element  python  fenics 

Digamos que eu tenho uma função que desejo integrar em um tetraedro . Se fosse arbitrário, a quadratura de Gauss seria uma boa solução, mas por acaso sei que é harmônico. Quanto a quadratura de Gauss pode ser acelerada usando essas informações?f:R3→Rf:R3→Rf : \mathbf{R}^3 \to \mathbf{R}T⊂R3T⊂R3T \subset \mathbf{R}^3ffffff Por exemplo, …

11 quadrature 

Eu tenho um sistema linear não homogêneo Ax=bAx=b Ax=b onde é uma matriz real com . O espaço nulo de é garantido como de dimensão zero, portanto a equação possui um inverso exclusivo . Como o resultado entra no lado direito de uma ODE, que pretendo resolver usando um método …

11 linear-algebra  ode 

Teoricamente, como o tempo para fazer um cálculo da teoria funcional da densidade (DFT) com o número de elétrons? Estou interessado em implementações "típicas" de DFT, como VASP, ABINIT, etc., não códigos O (N).

11 density-functional-theory 

A computação paralela é freqüentemente modelada usando uma taxa local determinística de computação, sobrecarga de latência e largura de banda da rede. Na realidade, são espacialmente variáveis ​​e não determinísticas. Estudos como Skinner e Kramer (2005) observam distribuições multimodais, mas a análise de desempenho parece sempre usar distribuições determinísticas ou …

11 performance  hpc  statistics 

Deseja melhorar este post? Forneça respostas detalhadas para esta pergunta, incluindo citações e uma explicação de por que sua resposta está correta. Respostas sem detalhes suficientes podem ser editadas ou excluídas. Sou totalmente novo na noção de ciência da computação e estou procurando um bom ponto de partida. Entendo que …

11 languages 

Assim, os estados de Cholesky decomposição teorema que que qualquer verdadeira simétrica definida positiva matriz tem uma decomposição de Cholesky M = G G ⊤ onde L é uma matriz triangular inferior.MMMM=LL⊤M=LL⊤M= LL^\topLLL Dado , já sabemos que existem algoritmos rápidos para calcular seu fator L de Cholesky .MMMLLL Agora, …

11 linear-algebra  algorithms 

Estou resolvendo para uma enorme matriz definida positiva esparsa A usando o método do gradiente conjugado (CG). É possível calcular o determinante de A usando as informações produzidas durante a resolução?Ax=bAx=bAx=bAAAAAA

11 linear-algebra  optimization  krylov-method  conjugate-gradient 

Olá comunidade scicomp, Trabalhei na área de algoritmos de grafos usando frameworks como NetworkX (Python), JUNG e YFiles (Java). Agora estou entrando na área da computação paralela e de alto desempenho. Para um novo projeto, estou procurando uma biblioteca de gráficos C ++ com os seguintes recursos: possui uma interface …

11 algorithms  parallel-computing  libraries  graph-theory 

A maioria das bibliotecas de matemática possui várias versões de funções de logaritmo. Na maioria das vezes, assumimos que eles são perfeitos, mas na verdade muitos deles oferecem apenas um certo número de dígitos de precisão. Para algumas funções, existem variantes numericamente mais estáveis. Por exemplo, Fortran, R, Java e …

11 floating-point  accuracy 

CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome Qual a relação entre eles? Quais são as vantagens / desvantagens deles, respectivamente? BTW, existe alguma outra biblioteca de otimização convexa de alta qualidade para Python / C ++ de propósito geral?

11 python  convex-optimization 

Existe uma lista publicada de melhores práticas para garantir a longevidade do código, visando resultados científicos reproduzíveis? (por exemplo, código aberto, práticas de documentação, seleção de dependências, seleção de idioma, máquinas virtuais etc.). Conhece estudos (ou faltam exemplos / anedotas) que tentaram estimar a meia-vida de códigos científicos típicos ou …

11 software  publications  reproducibility 

Eu gostaria de executar um PCA em um conjunto de dados composto por aproximadamente 40.000 amostras, cada uma exibindo cerca de 10.000 recursos. O uso da função princomp do Matlab leva consistentemente mais de meia hora; nesse ponto, eu mato o processo. Gostaria de encontrar uma implementação / algoritmo que …

11 high-dimensional  data-analysis 

A nVidia, por exemplo, possui CUBLAS, que promete aceleração de 7-14x. Ingenuamente, isso não chega nem perto do rendimento teórico de qualquer placa de GPU da nVidia. Quais são os desafios para acelerar a álgebra linear em GPUs, e já existem roteiros mais rápidos de álgebra linear disponíveis?

11 linear-algebra  lapack  blas  gpu 

Uma pergunta bastante simples: fazer uma integral multidimensional, considerando que alguém decidiu que algum tipo de método de Monte Carlo é apropriado, existe alguma vantagem de que uma integração regular de MC usando números pseudo-aleatórios tenha uma integração quase-Monte Carlo usando uma sequência quase aleatória ? Se sim, como eu …

11 monte-carlo