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 /...
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...
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 /...
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....
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...
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...
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
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...
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...
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...
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
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 adaptativo, é...
Dada uma função desconhecida , podemos avaliar seu valor em qualquer ponto de seu domínio, mas não temos sua expressão. Em outras palavras, é como uma caixa preta para nós.f:Rd→Rf:Rd→Rf:\mathbb R^d \to \mathbb Rfff Qual é o nome do problema de encontrar o minimizador de ? Quais são alguns métodos...
Como o MATLAB, por exemplo, calcula o SVD de uma determinada matriz? Suponho que a resposta provavelmente envolva calcular os autovetores e autovalores de A*A'. Se for esse o caso, eu também gostaria de saber como isso os
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 = f(x,y), \; in \;...
No trabalho, eu essencialmente trabalho como consultor independente. Para a gerência e os clientes, preciso estimar quanto tempo levará para desenvolver o software como parte de minha pesquisa em ciências computacionais. No entanto, minhas estimativas de tempo geralmente estão desativadas. Eu sei...
Eu tenho um conjunto de dados de pontos tridimensionais para os quais gostaria de construir uma malha, usando python. Todo o software que eu vi exige que você forneça as bordas. Existe um programa em python que toma como entrada um conjunto de pontos em 3D e produz uma malha? Se possível, eu...
Começando com a equação de advecção em forma de conservação. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x onde é uma velocidade que depende do espaço, e é uma concentração de uma espécie que é conservada.ua(x)a(x)a(x)uuu Discretizar o fluxo (onde o fluxo , é definido nas bordas das células entre os...
Atualmente, estou desenvolvendo um método de decomposição de domínio para a solução do problema de espalhamento. Basicamente, estou resolvendo um sistema de Helmholtz BVPs iterativamente. Discreto as equações usando o método dos elementos finitos sobre malhas triangulares ou tetraédricas. Estou...
12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Atualmente, estou fazendo isso com a biblioteca do tatu , mas ela acaba sendo muito lenta. O problema é que eu preciso fazer isso por um trilhão de matrizes e acontece que calcular os dois...