Ciência computacional

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 imagens bidimensionais de …

11 matlab  computational-chemistry  publications 

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 …

11 linear-algebra 

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 > 10^6 Ax=bAx=bAx=b com A=(Δ−K)A=(Δ−K) A=(\Delta - K) , ondeΔΔ\Delta é uma matriz muito esparsa …

11 linear-solver  dense-matrix  linear-system 

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 implementação do algoritmo descrito no referido artigo. Este …

11 software  documentation 

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 consultar a documentação. Recentemente, …

11 parallel-computing  python 

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 …

11 matrix  matrix-equations 

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 modos?

11 python  graph-theory 

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 de contorno são: ψ(0)=0ψ(0)=0\psi(0) = 0 e ψ(R)=0ψ(R)=0\psi(R)=0 , e só estou olhando para …

11 finite-difference  ode  accuracy 

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 …

11 interpolation  unstructured-mesh 

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 matriz?

11 petsc  python 

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 …

11 linear-algebra  optimization  matrix  iterative-method  linear-solver 

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 multiplicidade das raízes.

11 linear-algebra 

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 ' ) …

11 algorithms  performance 

(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 …

11 regression  publications 

(σ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 = fAAA Seguindo o esquema do FEM, reduzo meu problema a um problema de otimização J(u)=(Au,u)−2(f,u)→minuJ(u)=(Au,u)−2(f,u)→minu J(u) = (Au,u) - 2(f,u) \to \min_{u} …

11 finite-element