/// <summary>
/// Generates nodes and edges endlessly for demonstration.
/// </summary>
public void StartGeneration()
{
lock (_nodeLock) {
Color colour = Color.FromArgb(120, Color.White);
// Add basis nodes.
for (int i = 0; i < 100; i++)
{
Add(new Node(PseudoRandom.UInt64().ToString(), colour));
}
// Connect some basis nodes.
for (int i = 0; i < 80; i++)
{
Node a, b;
do
{
a = _nodes[PseudoRandom.Int32(_nodes.Count - 1)];
b = _nodes[PseudoRandom.Int32(_nodes.Count - 1)];
} while (a == b || a.IsConnectedTo(b));
Connect(a, b);
}
// Add group nodes.
for (int i = 0; i < 200; i++)
{
Node node = new Node(PseudoRandom.UInt64().ToString(), colour);
Connect(node, _nodes[PseudoRandom.Int32(10)]);
Add(node);
}
// Add outlier nodes.
for (int i = 0; i < 200; i++)
{
Node node = new Node(PseudoRandom.UInt64().ToString(), colour);
Connect(node, _nodes[PseudoRandom.Int32(_nodes.Count - 1)]);
Add(node);
}
// Connect more nodes.
for (int i = 0; i < 50; i++)
{
Node a, b;
do
{
a = _nodes[PseudoRandom.Int32(_nodes.Count - 1)];
b = _nodes[PseudoRandom.Int32(_nodes.Count - 1)];
} while (a == b || a.IsConnectedTo(b));
Connect(a, b);
}
}
}
/// <summary>
/// Generates the specified gravitational system.
/// </summary>
/// <param name="type">The system type to generate.</param>
public void Generate(SystemType type)
{
// Reset frames elapsed.
Frames = 0;
lock (_bodyLock) {
switch (type)
{
// Clear bodies.
case SystemType.None:
Array.Clear(_bodies, 0, _bodies.Length);
break;
// Generate slow particles.
case SystemType.SlowParticles: {
for (int i = 0; i < _bodies.Length; i++)
{
double distance = PseudoRandom.Double(1e6);
double angle = PseudoRandom.Double(Math.PI * 2);
Vector location = new Vector(Math.Cos(angle) * distance, PseudoRandom.Double(-2e5, 2e5), Math.Sin(angle) * distance);
double mass = PseudoRandom.Double(1e6) + 3e4;
Vector velocity = PseudoRandom.Vector(5);
_bodies[i] = new Body(location, mass, velocity);
}
}
break;
// Generate fast particles.
case SystemType.FastParticles: {
for (int i = 0; i < _bodies.Length; i++)
{
double distance = PseudoRandom.Double(1e6);
double angle = PseudoRandom.Double(Math.PI * 2);
Vector location = new Vector(Math.Cos(angle) * distance, PseudoRandom.Double(-2e5, 2e5), Math.Sin(angle) * distance);
double mass = PseudoRandom.Double(1e6) + 3e4;
Vector velocity = PseudoRandom.Vector(5e3);
_bodies[i] = new Body(location, mass, velocity);
}
}
break;
// Generate massive body demonstration.
case SystemType.MassiveBody: {
_bodies[0] = new Body(Vector.Zero, 1e10);
Vector location1 = PseudoRandom.Vector(8e3) + new Vector(-3e5, 1e5 + _bodies[0].Radius, 0);
double mass1 = 1e6;
Vector velocity1 = new Vector(2e3, 0, 0);
_bodies[1] = new Body(location1, mass1, velocity1);
for (int i = 2; i < _bodies.Length; i++)
{
double distance = PseudoRandom.Double(2e5) + _bodies[1].Radius;
double angle = PseudoRandom.Double(Math.PI * 2);
double vertical = Math.Min(2e8 / distance, 2e4);
Vector location = (new Vector(Math.Cos(angle) * distance, PseudoRandom.Double(-vertical, vertical), Math.Sin(angle) * distance) + _bodies[1].Location);
double mass = PseudoRandom.Double(5e5) + 1e5;
double speed = Math.Sqrt(_bodies[1].Mass * _bodies[1].Mass * G / ((_bodies[1].Mass + mass) * distance));
Vector velocity = Vector.Cross(location, Vector.YAxis).Unit() * speed + velocity1;
location = location.Rotate(0, 0, 0, 1, 1, 1, Math.PI * 0.1);
velocity = velocity.Rotate(0, 0, 0, 1, 1, 1, Math.PI * 0.1);
_bodies[i] = new Body(location, mass, velocity);
}
}
break;
// Generate orbital system.
case SystemType.OrbitalSystem: {
_bodies[0] = new Body(1e10);
for (int i = 1; i < _bodies.Length; i++)
{
double distance = PseudoRandom.Double(1e6) + _bodies[0].Radius;
double angle = PseudoRandom.Double(Math.PI * 2);
Vector location = new Vector(Math.Cos(angle) * distance, PseudoRandom.Double(-2e4, 2e4), Math.Sin(angle) * distance);
double mass = PseudoRandom.Double(1e6) + 3e4;
double speed = Math.Sqrt(_bodies[0].Mass * _bodies[0].Mass * G / ((_bodies[0].Mass + mass) * distance));
Vector velocity = Vector.Cross(location, Vector.YAxis).Unit() * speed;
_bodies[i] = new Body(location, mass, velocity);
}
}
break;
// Generate binary system.
case SystemType.BinarySystem: {
double mass1 = PseudoRandom.Double(9e9) + 1e9;
double mass2 = PseudoRandom.Double(9e9) + 1e9;
double angle0 = PseudoRandom.Double(Math.PI * 2);
double distance0 = PseudoRandom.Double(1e5) + 3e4;
double distance1 = distance0 / 2;
double distance2 = distance0 / 2;
Vector location1 = new Vector(Math.Cos(angle0) * distance1, 0, Math.Sin(angle0) * distance1);
Vector location2 = new Vector(-Math.Cos(angle0) * distance2, 0, -Math.Sin(angle0) * distance2);
double speed1 = Math.Sqrt(mass2 * mass2 * G / ((mass1 + mass2) * distance0));
double speed2 = Math.Sqrt(mass1 * mass1 * G / ((mass1 + mass2) * distance0));
Vector velocity1 = Vector.Cross(location1, Vector.YAxis).Unit() * speed1;
Vector velocity2 = Vector.Cross(location2, Vector.YAxis).Unit() * speed2;
_bodies[0] = new Body(location1, mass1, velocity1);
_bodies[1] = new Body(location2, mass2, velocity2);
for (int i = 2; i < _bodies.Length; i++)
{
double distance = PseudoRandom.Double(1e6);
double angle = PseudoRandom.Double(Math.PI * 2);
Vector location = new Vector(Math.Cos(angle) * distance, PseudoRandom.Double(-2e4, 2e4), Math.Sin(angle) * distance);
double mass = PseudoRandom.Double(1e6) + 3e4;
double speed = Math.Sqrt((mass1 + mass2) * (mass1 + mass2) * G / ((mass1 + mass2 + mass) * distance));
speed /= distance >= distance0 / 2 ? 1 : (distance0 / 2 / distance);
Vector velocity = Vector.Cross(location, Vector.YAxis).Unit() * speed;
_bodies[i] = new Body(location, mass, velocity);
}
}
break;
// Generate planetary system.
case SystemType.PlanetarySystem: {
_bodies[0] = new Body(1e10);
int planets = PseudoRandom.Int32(10) + 5;
int planetsWithRings = PseudoRandom.Int32(1) + 1;
int k = 1;
for (int i = 1; i < planets + 1 && k < _bodies.Length; i++)
{
int planetK = k;
double distance = PseudoRandom.Double(2e6) + 1e5 + _bodies[0].Radius;
double angle = PseudoRandom.Double(Math.PI * 2);
Vector location = new Vector(Math.Cos(angle) * distance, PseudoRandom.Double(-2e4, 2e4), Math.Sin(angle) * distance);
double mass = PseudoRandom.Double(1e8) + 1e7;
double speed = Math.Sqrt(_bodies[0].Mass * _bodies[0].Mass * G / ((_bodies[0].Mass + mass) * distance));
Vector velocity = Vector.Cross(location, Vector.YAxis).Unit() * speed;
_bodies[k++] = new Body(location, mass, velocity);
// Generate rings.
const int RingParticles = 100;
if (--planetsWithRings >= 0 && k < _bodies.Length - RingParticles)
{
for (int j = 0; j < RingParticles; j++)
{
double ringDistance = PseudoRandom.Double(1e1) + 1e4 + _bodies[planetK].Radius;
double ringAngle = PseudoRandom.Double(Math.PI * 2);
Vector ringLocation = location + new Vector(Math.Cos(ringAngle) * ringDistance, 0, Math.Sin(ringAngle) * ringDistance);
double ringMass = PseudoRandom.Double(1e3) + 1e3;
double ringSpeed = Math.Sqrt(_bodies[planetK].Mass * _bodies[planetK].Mass * G / ((_bodies[planetK].Mass + ringMass) * ringDistance));
Vector ringVelocity = Vector.Cross(location - ringLocation, Vector.YAxis).Unit() * ringSpeed + velocity;
_bodies[k++] = new Body(ringLocation, ringMass, ringVelocity);
}
continue;
}
// Generate moons.
int moons = PseudoRandom.Int32(4);
while (moons-- > 0 && k < _bodies.Length)
{
double moonDistance = PseudoRandom.Double(1e4) + 5e3 + _bodies[planetK].Radius;
double moonAngle = PseudoRandom.Double(Math.PI * 2);
Vector moonLocation = location + new Vector(Math.Cos(moonAngle) * moonDistance, PseudoRandom.Double(-2e3, 2e3), Math.Sin(moonAngle) * moonDistance);
double moonMass = PseudoRandom.Double(1e6) + 1e5;
double moonSpeed = Math.Sqrt(_bodies[planetK].Mass * _bodies[planetK].Mass * G / ((_bodies[planetK].Mass + moonMass) * moonDistance));
Vector moonVelocity = Vector.Cross(moonLocation - location, Vector.YAxis).Unit() * moonSpeed + velocity;
_bodies[k++] = new Body(moonLocation, moonMass, moonVelocity);
}
}
// Generate asteroid belt.
while (k < _bodies.Length)
{
double asteroidDistance = PseudoRandom.Double(4e5) + 1e6;
double asteroidAngle = PseudoRandom.Double(Math.PI * 2);
Vector asteroidLocation = new Vector(Math.Cos(asteroidAngle) * asteroidDistance, PseudoRandom.Double(-1e3, 1e3), Math.Sin(asteroidAngle) * asteroidDistance);
double asteroidMass = PseudoRandom.Double(1e6) + 3e4;
double asteroidSpeed = Math.Sqrt(_bodies[0].Mass * _bodies[0].Mass * G / ((_bodies[0].Mass + asteroidMass) * asteroidDistance));
Vector asteroidVelocity = Vector.Cross(asteroidLocation, Vector.YAxis).Unit() * asteroidSpeed;
_bodies[k++] = new Body(asteroidLocation, asteroidMass, asteroidVelocity);
}
}
break;
// Generate distribution test.
case SystemType.DistributionTest: {
Array.Clear(_bodies, 0, _bodies.Length);
double distance = 4e4;
double mass = 5e6;
int side = (int)Math.Pow(_bodies.Length, 1.0 / 3);
int k = 0;
for (int a = 0; a < side; a++)
{
for (int b = 0; b < side; b++)
{
for (int c = 0; c < side; c++)
{
_bodies[k++] = new Body(distance * (new Vector(a - side / 2, b - side / 2, c - side / 2)), mass);
}
}
}
}
break;
}
}
}
