Ciência computacional

Eu tenho um vetor que pode ser decomposto no espaço eigens do operador esparso hermitiano :MVVVMMM V=∑ivim^iV=∑ivim^iV = \sum_i v_i \hat{m}_i Existe uma maneira de encontrar o (o próprio vetor próprio) que corresponde ao maior (em magnitude)?vim^im^i\hat{m}_iviviv_i Eu quero essencialmente os maiores termos da soma, incluindo os autovetores de , …

8 sparse-matrix  eigensystem 

Estou procurando um método ou algoritmo para produzir um valor que descreva como duas malhas diferentes são geometricamente, mas que têm topologias diferentes. Um exemplo seria alguns dados CAD que tiveram uma malha FE aplicada e a peça depois que ela foi fabricada e, em seguida, digitalizada em uma malha …

8 optimization  finite-element 

dada uma matriz simétrica , e uma matriz arbitrária X ∈ R n × n , e um vetor v ∈ R n × 1 , é possível calcular a seguinte expressão no tempo O ( n 2 ) ?Y∈Rn×nY∈Rn×nY \in \mathbb{R}^{n \times n}X∈Rn×nX∈Rn×nX \in \mathbb{R}^{n \times n}v∈Rn×1v∈Rn×1v \in \mathbb{R}^{n …

8 optimization  algorithms  performance  matrix  complexity 

Para um protetor de tela bobo que estou tentando desenvolver, gostaria de gerar aleatoriamente uma matriz 2D livre de divergência de vetores 2D e usá-la para gerar um gráfico de convolução integral de linha. Ouvi que uma maneira de fazer isso é gerar ruído aleatório e depois projetar o componente …

8 fourier-analysis 

Agora que o Formulize / Eureqa começou a cobrar US $ 2500 por ano por usá-lo e ter aleijado a versão de teste, alguém sabe de alguma substituição que possa fazer coisas semelhantes, como encontrar uma equação com base em pontos de dados? Era gratuito, mas agora está fora da …

8 machine-learning  regression  symbolic-computation 

Ao escrever código MPI 3.0, que deve passar argce argvà MPI_Initchamada ou não, e por quê? EDIT: Como a resposta de Bill Barth levantou algumas perguntas, quero fazer algumas observações: A aprovação argc/ argvnão é necessária desde o MPI 1.1. A pergunta é especificamente sobre por que deveria / não …

8 c++  c  mpi 

Durante minha peregrinação no Mathematica.se , notei gradualmente que um certo tipo de problema de solução de equações diferenciais está nos "incomodando" o tempo todo, ou seja, o problema do valor limite (BVP) das equações diferenciais ordinárias não lineares (ODEs). O método de disparo , usado pela função MathematicaNDSolve , …

8 ode  nonlinear-equations 

Estou entrando no meu terceiro ano do meu programa de doutorado e até agora meu código de pesquisa (problemas inversos numéricos / imagens médicas / processamento de imagens / etc) consiste principalmente de scripts e funções desorganizados do MATLAB, com alguns C ++ e CUDA MEX'd rotinas lançadas para obter …

8 open-source 

Quero publicar o pacote de software que escrevi para meu trabalho de pós-graduação, juntamente com um artigo descrevendo o pacote. Até agora, estive considerando licenças liberais como o BSD. No entanto, agora estou ciente de que a publicação do seu código pode envolver um acordo com os termos de licenciamento …

8 publications 

Eu codifiquei um problema do Stokes Flow usando elementos finitos e estou no processo de verificar se ele funciona. Só não tenho certeza de qual taxa de convergência devo esperar enquanto refino a malha globalmente. Conheço problemas escalares usando funções de base linear que eu esperaria convergência de ordem ( …

8 pde  finite-element  convergence 

Estou perguntando sobre a solução da equação de Helmholtz em um domínio periódico com velocidade de onda constante por partes em diferentes regimes de frequência. Uma abordagem possível é resolver esse problema: escrever equações integrais nas superfícies de contorno em termos da função de Green do sistema. Como o domínio …

8 integral-equations  helmholtz-equation 

Fiz uma pergunta semelhante em physics.stackexchange , ignorando este site. Estou basicamente procurando uma maneira eficiente de implementar campos gravitacionais. Eu tenho um espaço 2D enorme, com milhares de objetos. Preciso simular como esses objetos são afetados pela gravidade um do outro. Eu pensei que seria possível classificar os objetos …

8 optimization  simulation  approximation 

Fundo: Atualmente, estou executando uma grande quantidade de experimentos de variação de parâmetros. Eles estão sendo executados no Python 2.6+, usando numpy. Esses experimentos levarão cerca de 2 semanas para serem executados. Aproximadamente, estou variando três parâmetros (variáveis ​​independentes) em um intervalo de valores. Estou corrigindo mais 6 variáveis ​​independentes …

8 software  visualization  data-analysis 

Preciso calcular a integral Eu= ∫R0 0f( R ) Jn( zn mrR) rdrI=∫0Rf(r)Jn(znmrR)rdrI = \int_0^R f(r)J_n\left(\frac{z_{nm}r}{R}\right)rdr onde é a função de Bessel da ordem do primeiro tipo, é seu m t h zero ef ( r ) é uma função real que é um pouco semelhante a J n (mas …

8 quadrature  integration  special-functions  numerics 

Quando aprendi sobre o SOR, ele foi dado principalmente como um dos primeiros exemplos de métodos iterativos e, posteriormente, os métodos iterativos que eu acabaria usando seriam os métodos do subespaço de Krylov. Algum dos métodos iterativos como Gauss-Seidel e SOR já foi usado na prática? Você conhece algum pacote …

8 linear-algebra