Otimização da simulação de corpo N, procurando nome ou trabalho existente


8

durante o desenvolvimento da minha simulação de corpo N com visualização no WebGL, desenvolvi uma otimização e fico imaginando se ele tem um nome. Acho improvável que isso nunca tenha sido feito antes.

Funciona assim: Durante o primeiro passo, faça uma iteração de todos os pares. Durante essa iteração, para cada partícula:

  1. Armazene todas as interações próximas em uma lista - representando todas as partículas próximas a essa. Essas interações serão avaliadas a partir de então a cada etapa do processo. Essa lista normalmente contém algumas entradas.
  2. Itere em todas as outras partículas e calcule uma força líquida que você armazena com a partícula. Essa força resultante é assim lembrada entre os intervalos de tempo e aplicada continuamente à partícula.

Então, à medida que a simulação continua após seu primeiro passo no tempo, de maneira rotativa, cada passo atualiza um pequeno número de listas de partículas de interações próximas e forças distantes . Para que, durante um certo número de etapas, digamos 1000, todas as interações íntimas das partículas e forças líquidas distantes tenham sido atualizadas. Os que você não atualizar apenas verificará suas interações estreitas e aplicará a força distante líquida. Neste exemplo, a complexidade computacional de cada etapa do tempo é algo como vez de .N2/1000N2

Um truque para também tornar isso razoavelmente preciso é identificar melhor "interações íntimas". Às vezes, a proximidade não é o melhor indicador - você também pode considerar a massa e a velocidade relativa e assim por diante. "Interações mais significativas" pode ser uma palavra melhor. Ou "interações mais prováveis ​​de mudar em breve".

Essa otimização permite muito mais partículas interagindo do que o método dos pares, mas não sei como descrevê-lo em termos de O (), pois ele não cria uma solução completa a cada intervalo de tempo, mas reutiliza (ligeiramente incorreto) o antigo informações e espalha o esforço computacional ao longo do tempo.

( Isenção de responsabilidade: minha simulação webgl também possui partículas "vfx" que são afetadas apenas pela gravidade e não retribuem o efeito, por isso não é tão incrivelmente rápido quanto parece )

Então, essa técnica de otimização tem um nome?


Bem-vindo ao SciComp.SE. Sua visualização está ótima, BTW. O que você está descrevendo não parece ser uma técnica de otimização, mas sim um solucionador de força não tão bruta para o problema do corpo n. Talvez relacionado ao método multipolo rápido .
nicoguaro

Respostas:


6

O método que você descreve é ​​um pouco semelhante ao algoritmo de Barnes-Hut . A principal diferença é que você tem um único nível de interações próximas, enquanto a Barnes-Hut tem .registroN

O(NregistroN)

EDITAR:

(N/3.8)14


Muito legal! Parece correto .. embora eu ache que o ganho em velocidade computacional possa ser significativamente maior se você aceitar um alto nível de imprecisão, o que pode ser bom se a simulação, como no meu caso, é principalmente para efeito visual e não para fins científicos.
Magnus Wolffelt

O(N)
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.