Ciência computacional

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 CSPRNG, que não é …

8 performance  c  random-number-generation 

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 …

8 linear-algebra  optimization 

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 função. Tentei procurar por …

8 special-functions  numerics 

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

8 matlab  eigensystem  sparse-matrix  eigenvalues 

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 - 10 00 0⋮0 0bn - 2uman …

8 matrix  linear-solver 

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) = ∇ ⋅ …

8 pde  finite-element  reference-request 

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 …

8 reference-request  nonlinear-equations 

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 quero escrever o código …

8 parallel-computing  gpu 

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 geral, mas …

8 pde  matlab  finite-difference 

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 / reais.

8 simulation  numerical-analysis 

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 …

8 ode  c++  time-integration  cuda  runge-kutta 

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 qualquer outro pacote python para calcular esse caminho? …

8 optimization  python  scipy 

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 …

8 optimization  convex-optimization 

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 suaves de . Como os valores próprios não são retornados …

8 linear-algebra  matrix  eigensystem  numpy 

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 correspondente isendainda não foi chamado). Para ficar claro, essa …

8 parallel-computing  software  mpi