Galaxy simulation

This version
2005-01-03


ood2 is now pre-multiplied into dirx and diry

This version
2005-01-03

ood2 is now pre-multiplied into dirx and diry

took out wrap-around from calculations

uses fft fluid solver with gravity forcefield

todo:

use forward projection instead of sime-lagrangian lookups

In this version:

  use forward projection instead of sime-lagrangian lookups
  combine ood2 lookup with dirx and diry

This version:

  Made gravity travel at the speed of light instead of infitiely fast

  Method:

    Store previous mass frames
    Instead of using current mass frame for each pixel's calculation use a past mass
    frame depending on the distance of the pixel from the one it is affecting.

  I'll call the technique Pass Mass Lookup

pprizm (prizm ep) - i miss u already (benbecular label, scotland, techno)

This version:

  woohoo - looks like i can do gravity using fft

  I've had to negate the gravity strength and swap the x and the y lookups for some unknown reason

Operation
-

The universe is DIM x DIM, and DIM must be a power of two.
The universe is tiling (wraps around at the edges)

(because it uses &(DIM-1) as a quick way of wrapping around,
and because the FFT expects powers of two (i think) )

The universe is represented by a mass and an x and y direction (velocity)
for each pixel in a DIM x DIM square.

The gravity calculation involves calculating the force on each pixel
from the rest of the universe:

  for each pixel A
  {
    force_x = 0
    force_y = 0

    for each pixel B except A
    {
      force_x += horizontal force on pixel A of pixel B
      force_y += vertical force on pixel A of pixel B
    }

    //  (force_x, force_y) is the total force on this pixel
    //  from the rest of the universe
  }

After the force is calculated it is used to modify the velocity of the pixel.

Because the universe wraps around, each of the other copies of the universe
has an influence.

  for each pixel A
  {
    force_x = 0
    force_y = 0

    for each tile
    {
      for each pixel B except A  //  i think you need to include A in all galaxies apart from real one?
      {
        force_x += horizontal force on pixel A of pixel B
        force_y += vertical force on pixel A of pixel B
      }
    }

    //  (force_x, force_y) is the total force on this pixel
    //  from the rest of the universe
  }

The force is a convolution of the mass array and the gravity force array.

To move the mass based on the force, start with an empty array (zero's)

  NEW_MASS = zero array

  for each pixel
  {
    M,DX,DY = pixel's mass and velocity

    NEW_X = pixel's x position + DX
    NEW_Y = pixel's y position + DY

    Place mass and velocity into NEW_MASS with sub-pixel accuracy
  }

To do a fast convolution of A and B (both power of two
The full sequence is:

  Create gravity lookup tables LOOKUPX and LOOKUPY

  each frame
  {
    FORCE_X = convolution of LOOKUPX and MASS
    FORCE_Y = convolution of LOOKUPY and MASS
    Apply FORCE_X,FORCE_Y to MASS
  }






  Create gravity lookup tables LOOKUPX and LOOKUPY
  Calculate fft of LOOKUPX as V0
  Calculate fft of LOOKUPY as V1

  each frame
  {
  }


The gravity force array can be precalculated and saved off.






Notes
-

. Horizontal and vertical bias because of the grid nature of the universe
  . Horiz and vertical streams coming out of tiny heavy galaxies when matter should be fired out in all directions
  . Long thin galaxies form which suspiciously travel exactly in the x or y dir

003:

  use forward projection instead of sime-lagrangian lookups
  combine ood2 lookup with dirx and diry

004:

  use fft image convolution to calculate gravity (avoids n^4 calc and reduces it to n^2)

007

trying FFT for gravity convolution

just put init_fft in, need to calculate kernel etc.