Qual é a política de várias transferências assíncronas sobrepostas no MPI? Eu tenho um programa com várias irecvoperações assíncronas abertas . Acho que as transferências que podem ocorrer (o correspondente isendfoi chamado) aguardam outras transferências que ainda não estão prontas (o...
Estou procurando o caminho mínimo entre os mínimos de uma superfície potencial que já é conhecida em uma grade. (fonte: http://www.math.nus.edu.sg/~matrw/string/ ) Qualquer ponto no caminho é o mínimo potencial em todas as direções perpendicular ao caminho. Existe algum método SciPy ou...
Quero resolver a tarefa de otimização (convexa): m a xr , zrmumaxr,zrmax_{r,z}\quad r sujeito às duas restrições a seguir \ | z \ | \ leq 1 r \ geq0 ‖ z ‖ ≤ 1 r ≥ 0r ∥ xEu∥ - xTEuz≤ 0∀ i = 1 , … , Nr__xEu__-xEuTz≤0 0∀Eu=1,…,Nr\|x_i\| - x_i^Tz \leq 0 \qquad \forall i=1,\dots, N ∥ z∥ ≤...
Eu tenho matrizes dimensionais dependendo do parâmetro do vetor .H ( → k ) → knnnH^( k⃗ )H^(k→)\mathrm{\hat{H}}(\vec{k})k⃗ k→\vec{k} Agora, as rotinas de autovalores retornam autovalores em nenhuma ordem específica (geralmente são classificadas), mas quero rastrear os autovalores como funções...
Suponha que você tenha um sistema de PDEs para resolver. Pelo menos por simplicidade, vamos supor que seja independente do tempo, quase linear (linear em suas derivadas) resolvido em uma grade retangular no espaço (x, y) e com condições de contorno especificadas ao redor. Minha pergunta é mais...
Estou tentando modelar a condução de calor dentro de um cilindro de madeira usando métodos implícitos de diferenças finitas. A equação geral de calor que estou usando para formas cilíndricas e esféricas é: Onde p é o fator de forma, p = 1 para cilindro ep = 2 para esfera. As condições de...
Estou apenas me perguntando se o Cálculo Externo Discreto, como um novo método numérico, é bom para resolver problemas numéricos - em elasticidade, fluidos ou outras áreas físicas /
Estou procurando bibliotecas de GPU para acelerar um código que tenho, cujas regiões primárias de "alto desempenho" possuem decomposição de valor singular, fatoração QR e valores de Eigen, computação de Eigenvectors. Eu procurei no site da Nvidia em suas ferramentas como CuBLAS, Magma etc., mas não...
Em muitas áreas de aplicação, é necessário resolver um sistema de equações não lineares Por vezes, a formulação ‖ F ( x ) ‖ 2 → min é usada. Claramente, cada solução x de F ( x ) = 0 é também uma solução do segundo problema; o inverso também é verdadeiro (se houver uma...
Meu projeto atual é uma reprogramação de um modelo de dobramento de proteínas envolvendo a solução de milhares de ODEs em C ++. Tenho parado e começado a progredir enquanto escrevia o solucionador para rodar totalmente na GPU. Finalmente consegui integrar, mas quando tento resolver dC / dt = -C...
Eu tenho uma matriz esparsa e simétrica A (nxn). O método Lanczos transforma a matriz A em matriz tridiagonal e simétrica T e os vetores Lanczos na matriz V. A partir daí, como computo k autovalores menores ou maiores e autovetores
Eu tenho esse problema no meu livro: Sugira um algoritmo eficiente para resolver sistemas de equações lineares com matriz tridimensional cíclica, que é da forma: sem trocar nenhuma linha e coluna. Estimar a complexidade.⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢uma1c20 0⋮0 0d2b1uma2⋱⋮⋯0 00 0b2⋱cn - 2⋯⋯⋯0 0⋱uman - 2cn -...
Gostaria de calcular, como parte da solução da equação de Laplace, usando o Método Multipolar Rápido, a segunda derivada das funções legendre associadas do primeiro tipo . Especificamente, estou procurando implementações em C ou apenas a relação de recorrência correta a ser usada para escrever a...
Estou pensando em acelerar os produtos vetoriais de matriz, mas tudo o que leio é sobre como fazer isso para matrizes muito grandes. No meu caso, as matrizes são pequenas, mas o número de vezes que isso deve ser feito é muito grande. Quais métodos, se houver, existem para otimizar isso? Seria mais...
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ê...
Eu usei MT19937em um equipamento de teste para gerar valores [0, - 1] de 32 bits uniformemente (não assinados) , com base na implementação mt19937.c do Authors original , para gerar um suprimento (essencialmente inesgotável) de fluxo estatisticamente aleatório de octetos de bits. Isso substitui um...
Eu fiz uma pergunta um pouco semelhante anteriormente, mas talvez possa ter sido muito específico para alguém realmente responder. Aqui está um pouco mais geral de uma pergunta com a qual estou lutando. Considere o seguinte sistema: ∂ u 2- ∇ ⋅ ( D1( u2) ∇ u1) = ∇ ⋅ f1(...
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 não do mesmo modo, é bastante...
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...
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...