Suponha que a seguinte matriz seja dada com sua transposta . O produto gera ,AAA⎡⎣⎢0.500−0.500−0.500−0.3330.667−0.333−0.167−0.1670.833⎤⎦⎥[0.500−0.333−0.167−0.5000.667−0.167−0.500−0.3330.833] \left[\begin{array}{ccc} 0.500 & -0.333 & -0.167\\ -0.500 & 0.667 & -0.167\\ -0.500 & -0.333 &...
(σ2(x)u′′(x))′′=f(x),0⩽x⩽1(σ2(x)u″(x))″=f(x),0⩽x⩽1 \left( \sigma^{2}(x) u ''(x) \right)'' = f(x), \;\;\; 0 \leqslant x \leqslant 1 u(0)=u(1)=0u(0)=u(1)=0u(0) = u(1) = 0u′′(0)=u′′(1)=0u″(0)=u″(1)=0u''(0) = u''(1) = 0σ(x)⩾σ0>0σ(x)⩾σ0>0\sigma(x) \geqslant \sigma_{0} > 0Au=fAu=fAu =...
Por favor, liste o pacote Python (petsc4py, etc ...) e os solucionadores diretos esparsos que ele suporta. Uma resposta (wiki da comunidade) por pacote, por
(Espero que esta pergunta se encaixe neste site; caso contrário, aceite minhas desculpas). Fiz uma certa simulação e obtive uma série temporal y (t), t = 0, 1, ... 20. Depois de tentar algumas funções, descobri que: y(t) =~ 1 / (A t + B) Onde A e B são coeficientes calculados usando regressão...
Suponha que A seja uma matriz esparsa geral e eu quero calcular os autovalores. Não sei como detectar a multiplicidade dos valores próprios. Até onde eu sei, para um caso especial, encontrando as raízes polinomiais pelo método da matriz companheira, podemos aplicar o RRQR para detectar a...
Dado o sistema onde , li que, caso a iteração Jacobi seja usada como solucionador, o método não convergirá se tiver um valor diferente de zero componente no nulo-espaço de . Então, como alguém poderia formalmente afirmar que, desde que tenha um componente diferente de zero, abrangendo o espaço nulo...
Quanto mais lento o petsc4py vs c / c ++ / fortran? Sei que isso dependerá significativamente do código que está sendo executado, mas e algo simples como um produto de vetor de
Estou implementando o artigo " Transporte de massa ideal para registro e entortamento ", meu objetivo é colocá-lo online, pois não consigo encontrar nenhum código de transporte de massa euleriano on-line e isso seria interessante pelo menos para a comunidade de pesquisa em processamento de...
Este é um problema simples de enumeração de simetria. Eu dou o histórico completo aqui, mas nenhum conhecimento de química quântica é necessário. O integrante dois partícula é: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) ψ k ( x ) ψ l ( x ' )⟨ I j | k l ⟩⟨ij|kl⟩\langle ij|kl\rangle E tem os...
Eu tenho conjuntos de dados de campo multiponto, cada conjunto de dados de pontos se refere a uma única célula de uma malha não estruturada. O objetivo é interpolar os dados para o centro celular, direta ou indiretamente, da maneira mais precisa. Se eu usar a interpolação com Distância inversa...
Estou tentando resolver uma equação do tipo: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Onde f(x)f(x)f(x) possui um pólo simples em 000 , para os menores NNN autovalores e autovetores. As condições...
Existem algoritmos para detecção da comunidade para gráficos bipartidos (redes de dois modos) implementados no igraph, networkX, R ou Python etc.? Em particular, existe uma implementação em que alguém seria capaz de restringir a detecção de comunidades apenas em um dos dois...
A seguinte equação da matriz em Σ - para determinadas matrizes B e C - aparece no meu trabalho como uma caracterização de uma matriz de covariância. Aprendi que essa equação é conhecida, particularmente na teoria de controle de tempo contínuo, como a equação de Lyapunov , e que existem vários...
Tenho experiência em codificar OpenMP para máquinas de memória compartilhada (em C e FORTRAN) para realizar tarefas simples como adição de matriz, multiplicação etc. (Apenas para ver como ele compete com o LAPACK). Conheço o OpenMP o suficiente para realizar tarefas simples sem a necessidade de...
A revista Association for Computing Machinery Transactions on Mathematics Software (ACM TOMS) publica muitos artigos sobre algoritmos numéricos que incluem implementações de software. De acordo com sua política editorial , a submissão de um artigo sobre algoritmos inclui o código-fonte para uma...
Existe alguma esperança em resolver o seguinte sistema linear de maneira eficiente com um método iterativo? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n >
Suponha que é uma matriz simétrica real e sua decomposição de autovalor V Λ V T é dada. É fácil ver o que acontece com os autovalores da soma A + c I, em que c é uma constante escalar (consulte esta pergunta ). Podemos tirar alguma conclusão no caso geral A + D, onde D é uma matriz diagonal...
Atualmente, estou produzindo imagens 1d e 2d usando o Matlab para comparação de precisão em relação a um determinado modelo. Preciso comparar meus métodos com o modelo gaussiano .wfn padrão e farei isso analisando a densidade de uma molécula e também o laplaciano. Estou interessado em produzir...
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...
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...