a981a46170
Sector renderer is now also rendering asteroid belts. All gravity calculations have been moved out of monobehaviors.
409 lines
No EOL
16 KiB
C#
409 lines
No EOL
16 KiB
C#
using System;
|
|
using System.Collections;
|
|
using System.Collections.Generic;
|
|
using System.Globalization;
|
|
using System.Linq;
|
|
using MessagePack;
|
|
using Unity.Mathematics;
|
|
using static Unity.Mathematics.math;
|
|
using Random = Unity.Mathematics.Random;
|
|
|
|
public class ZoneGenerator
|
|
{
|
|
public static ZonePack GenerateZone(
|
|
GameContext context,
|
|
string name,
|
|
float2 position,
|
|
IEnumerable<GalaxyMapLayerData> mapLayers,
|
|
IEnumerable<SimpleCommodityData> resources,
|
|
float mass = 0,
|
|
float radius = 0)
|
|
{
|
|
var zone = new ZoneData
|
|
{
|
|
Name = name,
|
|
Position = position
|
|
};
|
|
|
|
var pack = new ZonePack
|
|
{
|
|
Data = zone
|
|
};
|
|
|
|
var zoneRadius = radius != 0 ? radius : context.Random.NextPowerDistribution(
|
|
context.GlobalData.MinimumZoneRadius,
|
|
context.GlobalData.MaximumZoneRadius,
|
|
context.GlobalData.ZoneRadiusExponent,
|
|
1 / lerp(
|
|
context.GlobalData.ResourceDensityMinimum,
|
|
context.GlobalData.ResourceDensityMaximum,
|
|
mapLayers
|
|
.First(m => m.Name == "Radius")
|
|
.Evaluate(zone.Position, context.GlobalData)));
|
|
zone.Radius = zoneRadius;
|
|
|
|
var zoneMass = mass != 0 ? mass : context.Random.NextPowerDistribution(
|
|
context.GlobalData.MinimumZoneMass,
|
|
context.GlobalData.MaximumZoneMass,
|
|
context.GlobalData.ZoneMassExponent,
|
|
1 / lerp(
|
|
context.GlobalData.ResourceDensityMinimum,
|
|
context.GlobalData.ResourceDensityMaximum,
|
|
mapLayers
|
|
.First(m => m.Name == "Mass")
|
|
.Evaluate(zone.Position, context.GlobalData)));
|
|
zone.Mass = zoneMass;
|
|
|
|
//Debug.Log($"Generating zone at position {zone.Position} with radius {zoneRadius} and mass {zoneMass}");
|
|
|
|
var planets = GenerateEntities(context, zone, zoneMass, zoneRadius);
|
|
|
|
// Create collections to map between zone generator output and database entries
|
|
var orbitMap = new Dictionary<GeneratorPlanet, OrbitData>();
|
|
var orbitInverseMap = new Dictionary<OrbitData, GeneratorPlanet>();
|
|
|
|
// Create orbit database entries
|
|
pack.Orbits = planets.Select(planet =>
|
|
{
|
|
var data = new OrbitData
|
|
{
|
|
Distance = planet.Distance,
|
|
Period = planet.Period,
|
|
Phase = planet.Phase
|
|
};
|
|
orbitMap[planet] = data;
|
|
orbitInverseMap[data] = planet;
|
|
return data;
|
|
}).ToList();
|
|
|
|
// Link OrbitData parents to database GUIDs
|
|
foreach (var data in pack.Orbits)
|
|
data.Parent = orbitInverseMap[data].Parent != null
|
|
? orbitMap[orbitInverseMap[data].Parent].ID
|
|
: Guid.Empty;
|
|
|
|
// Cache resource densities
|
|
var resourceMaps = mapLayers.ToDictionary(m => m.ID, m => m.Evaluate(zone.Position, context.GlobalData));
|
|
|
|
pack.Planets = planets.Where(p=>!p.Empty).Select(planet =>
|
|
{
|
|
Dictionary<Guid, float> planetResources = new Dictionary<Guid, float>();
|
|
BodyType bodyType = planet.Belt ? BodyType.Asteroid :
|
|
planet.Mass > context.GlobalData.SunMass ? BodyType.Sun :
|
|
planet.Mass > context.GlobalData.GasGiantMass ? BodyType.GasGiant :
|
|
planet.Mass > context.GlobalData.PlanetMass ? BodyType.Planet : BodyType.Planetoid;
|
|
|
|
foreach (var r in resources)
|
|
{
|
|
if ((bodyType & r.ResourceBodyType) != 0)
|
|
{
|
|
float quantity = ResourceValue(r, context, r.ResourceDensity.Aggregate(1f, (m, rdm) => m * resourceMaps[rdm]));
|
|
if (r.Floor < quantity) planetResources.Add(r.ID, quantity);
|
|
}
|
|
}
|
|
|
|
BodyData planetData;
|
|
switch (bodyType)
|
|
{
|
|
case BodyType.Asteroid:
|
|
planetData = new AsteroidBeltData();
|
|
break;
|
|
case BodyType.Planetoid:
|
|
case BodyType.Planet:
|
|
planetData = new PlanetData();
|
|
break;
|
|
case BodyType.GasGiant:
|
|
planetData = new GasGiantData();
|
|
break;
|
|
case BodyType.Sun:
|
|
planetData = new SunData();
|
|
break;
|
|
default:
|
|
throw new ArgumentOutOfRangeException();
|
|
}
|
|
|
|
planetData.Mass.Value = planet.Mass;
|
|
planetData.Orbit = orbitMap[planet].ID;
|
|
planetData.Resources = planetResources;
|
|
planetData.Name.Value = planetData.ID.ToString().Substring(0, 8);
|
|
if (planetData is AsteroidBeltData beltData)
|
|
{
|
|
beltData.Asteroids =
|
|
Enumerable.Range(0, (int) (pow(planetData.Mass.Value, context.GlobalData.BeltMassExponent) / context.GlobalData.BeltMassRatio * orbitMap[planet].Distance))
|
|
.Select(_ => new Asteroid
|
|
{
|
|
Distance = context.Random.NextFloat(orbitMap[planet].Distance * .5f, orbitMap[planet].Distance * 1.5f),
|
|
Phase = context.Random.NextFloat(),
|
|
Size = pow(context.Random.NextFloat(), context.GlobalData.AsteroidSizeExponent) *
|
|
(context.GlobalData.AsteroidSizeMax - context.GlobalData.AsteroidSizeMin) +
|
|
context.GlobalData.AsteroidSizeMin,
|
|
RotationSpeed = context.Random.NextFloat() * context.GlobalData.AsteroidRotationSpeed
|
|
}).ToArray();
|
|
}
|
|
else if (planetData is GasGiantData gas)
|
|
{
|
|
var times = Enumerable.Range(0, 5).Select(i => (float) i / 4).ToArray();
|
|
gas.Colors.Value = times.Select(time => float4(context.Random.NextFloat3(),time)).ToArray();
|
|
gas.AlbedoRotationSpeed.Value = -3;
|
|
gas.FirstOffsetRotationSpeed.Value = 5;
|
|
gas.FirstOffsetDomainRotationSpeed.Value = 0;
|
|
gas.SecondOffsetRotationSpeed.Value = 10;
|
|
gas.SecondOffsetDomainRotationSpeed.Value = -25;
|
|
if (planetData is SunData sun)
|
|
{
|
|
sun.FogTintColor.Value = sun.LightColor.Value = context.Random.NextFloat3();
|
|
gas.FirstOffsetDomainRotationSpeed.Value = 5;
|
|
}
|
|
}
|
|
return planetData;
|
|
}).ToList();
|
|
|
|
return pack;
|
|
}
|
|
|
|
static float ResourceValue(SimpleCommodityData resource, GameContext context, float density)
|
|
{
|
|
return context.Random.NextPowerDistribution(resource.Minimum, resource.Maximum, resource.Exponent,
|
|
1 / lerp(context.GlobalData.ResourceDensityMinimum, context.GlobalData.ResourceDensityMaximum, density));
|
|
}
|
|
|
|
public static GeneratorPlanet[] GenerateEntities(GameContext context, ZoneData data, float mass, float radius)
|
|
{
|
|
var random = new Random();
|
|
random.InitState(unchecked((uint)data.Name.GetHashCode()));
|
|
|
|
var root = new GeneratorPlanet
|
|
{
|
|
Context = context,
|
|
Mass = mass,
|
|
ChildDistanceMaximum = radius * .75f,
|
|
ChildDistanceMinimum = context.GlobalData.PlanetRadius(mass)
|
|
};
|
|
|
|
// There is some chance of generating a rosette or binary system
|
|
// Probabilities which are fixed for the entire galaxy are in GlobalData, contained in the GameContext
|
|
var rosette = random.NextFloat() < context.GlobalData.RosetteProbability;
|
|
|
|
if (rosette)
|
|
{
|
|
// Create a rosette with a number of vertices between 2 and 9 inclusive
|
|
root.ExpandRosette((int)(random.NextFloat(1, 5) + random.NextFloat(1, 5)));
|
|
|
|
// Create a small number of less massive "captured" planets orbiting past the rosette
|
|
root.ExpandSolar(
|
|
count: (int)(random.NextFloat(1, 3) * random.NextFloat(1, 2)),
|
|
massMulMin: .6f,
|
|
massMulMax: .8f,
|
|
distMulMin: 1.25f,
|
|
distMulMax: 1.75f,
|
|
jupiterJump: 1,
|
|
massFraction: .1f);
|
|
|
|
var averageChildMass = root.Children.Sum(p => p.Mass) / root.Children.Count;
|
|
foreach(var p in root.Children.Where(c=>c.Mass>context.GlobalData.GasGiantMass))
|
|
{
|
|
var m = p.Mass / averageChildMass;
|
|
// Give each child in the rosette its own mini solar system
|
|
p.ExpandSolar(
|
|
count: (int) (random.NextFloat(1, 3 * m) + random.NextFloat(1, 3 * m)),
|
|
massMulMin: 0.75f,
|
|
massMulMax: 2.5f,
|
|
distMulMin: 1 + m * .25f,
|
|
distMulMax: 1.05f + m * .5f,
|
|
jupiterJump: random.NextFloat() * random.NextFloat() * 10 + 1,
|
|
massFraction: .5f
|
|
);
|
|
}
|
|
}
|
|
else
|
|
{
|
|
// Create a regular old boring solar system
|
|
root.ExpandSolar(
|
|
count: random.NextInt(5, 15),
|
|
massMulMin: 0.5f,
|
|
massMulMax: 2.0f,
|
|
distMulMin: 1.1f,
|
|
distMulMax: 1.25f,
|
|
jupiterJump: random.NextFloat() * random.NextFloat() * 10 + 1,
|
|
massFraction: .25f
|
|
);
|
|
}
|
|
|
|
// Get all children that are above the satellite creation mass floor and not rosette members
|
|
var satelliteCandidates = rosette
|
|
? root.AllPlanets().Where(p => p != root && p.Parent != root && p.Mass > context.GlobalData.SatelliteCreationMassFloor)
|
|
: root.AllPlanets().Where(p => p != root && p.Mass > context.GlobalData.SatelliteCreationMassFloor);
|
|
|
|
var binaries = new List<GeneratorPlanet>();
|
|
foreach (var planet in satelliteCandidates)
|
|
{
|
|
// There's a chance of generating satellites for each qualified planet
|
|
if (random.NextFloat() < context.GlobalData.SatelliteCreationProbability)
|
|
{
|
|
// Sometimes the satellite is so massive that it forms a binary system (like Earth!)
|
|
if (random.NextFloat() < context.GlobalData.BinaryCreationProbability)
|
|
{
|
|
planet.ExpandRosette(2);
|
|
binaries.AddRange(planet.Children);
|
|
}
|
|
// Otherwise, terrestrial planets get a couple satellites while gas giants get many
|
|
else planet.ExpandSolar(
|
|
count: planet.Mass < context.GlobalData.GasGiantMass ? random.NextInt(1,3) : random.NextInt(4,10),
|
|
massMulMin: .75f,
|
|
massMulMax: 1.5f,
|
|
distMulMin: 1.05f,
|
|
distMulMax: 1.25f,
|
|
jupiterJump: 1,
|
|
massFraction: .15f); // Planetary satellites are not nearly as massive as planets themselves
|
|
}
|
|
}
|
|
|
|
// Get all children that are below the belt creation mass floor and not rosette members, also exclude binaries
|
|
var beltCandidates = rosette
|
|
? root.AllPlanets().Where(p => p != root && p.Parent != root && p.Mass < context.GlobalData.BeltMassCeiling && !binaries.Contains(p))
|
|
: root.AllPlanets().Where(p => p != root && p.Mass < context.GlobalData.BeltMassCeiling && !binaries.Contains(p));
|
|
foreach(var planet in beltCandidates.Reverse())
|
|
if (random.NextFloat() < context.GlobalData.BeltProbability && !planet.Parent.Children.Any(p=>p.Belt))
|
|
planet.Belt = true;
|
|
|
|
return root.AllPlanets().ToArray();
|
|
}
|
|
}
|
|
|
|
// All orbiting bodies are referred to here as "Planets"
|
|
// This goes back to the original etymology for the word planet, which is traveler,
|
|
// Because they were objects we observed to move in the sky instead of remaining fixed like the stars
|
|
public class GeneratorPlanet
|
|
{
|
|
public GameContext Context;
|
|
public float Distance;
|
|
public float Phase;
|
|
public float Mass;
|
|
public float Period;
|
|
public float ChildDistanceMinimum;
|
|
public float ChildDistanceMaximum;
|
|
public bool Empty = false;
|
|
public bool Belt = false;
|
|
public List<GeneratorPlanet> Children = new List<GeneratorPlanet>(); // Planets orbiting this one are referred to as children
|
|
public GeneratorPlanet Parent;
|
|
|
|
// Recursively gather all planets in the hierarchy
|
|
public IEnumerable<GeneratorPlanet> AllPlanets()
|
|
{
|
|
return new[]{this}.Concat(Children.SelectMany(c=>c.AllPlanets()));
|
|
}
|
|
|
|
// Create children that fill a single orbit, equally spaced
|
|
// https://en.wikipedia.org/wiki/Klemperer_rosette
|
|
public void ExpandRosette(int vertices)
|
|
{
|
|
//Debug.Log("Expanding Rosette");
|
|
|
|
// Rosette children replace the parent, parent orbital node is left empty
|
|
Empty = true;
|
|
|
|
// Masses in a rosette alternate, so every sequential pair has the same shared mass
|
|
var sharedMass = Mass / vertices * 2;
|
|
|
|
// Proportion of the masses of sequential pairs is random, but only for even vertex counts
|
|
var proportion = vertices % 2 == 0 ? Context.Random.NextFloat(.5f,.95f) : .5f;
|
|
|
|
// Place children at a fixed distance in the center of the range
|
|
var dist = (ChildDistanceMinimum + ChildDistanceMaximum) / 2;
|
|
|
|
// Position of first child
|
|
var p0 = new float2(0, dist);
|
|
|
|
// Position of second child
|
|
var p1 = OrbitData.Evaluate(1.0f / vertices) * dist;
|
|
|
|
// Maximum child distance is half the distance to the neighbor minus the neighbor's radius
|
|
var p0ChildDist = (distance(p0, p1) * proportion - Context.GlobalData.PlanetRadius(sharedMass * (1 - proportion))) * .75f;
|
|
var p1ChildDist = (distance(p0, p1) * (1 - proportion) - Context.GlobalData.PlanetRadius(sharedMass * proportion)) * .75f;
|
|
|
|
for (int i = 0; i < vertices; i++)
|
|
{
|
|
var child = new GeneratorPlanet
|
|
{
|
|
Context = Context,
|
|
Parent = this,
|
|
Mass = sharedMass * (i % 2 == 0 ? proportion : 1 - proportion), // Masses alternate
|
|
Distance = dist,
|
|
Phase = (float) i / vertices,
|
|
ChildDistanceMaximum = (i % 2 == 0 ? p0ChildDist : p1ChildDist)
|
|
};
|
|
child.ChildDistanceMinimum = Context.GlobalData.PlanetRadius(child.Mass) * 2;
|
|
child.Period = Context.GlobalData.OrbitalPeriod(child.Distance);
|
|
Children.Add(child);
|
|
}
|
|
|
|
ChildDistanceMinimum = dist + p0ChildDist;
|
|
}
|
|
|
|
// Create children that mimic the distribution of planetary masses in the solar system
|
|
public void ExpandSolar(int count, float massMulMin, float massMulMax, float distMulMin, float distMulMax, float jupiterJump, float massFraction)
|
|
{
|
|
//Debug.Log("Expanding Solar");
|
|
// Expansion is impossible when child space is full
|
|
if (count == 0 || ChildDistanceMaximum < ChildDistanceMinimum)
|
|
return;
|
|
|
|
var masses = new float[count];
|
|
var distances = new float[count];
|
|
|
|
// Accumulate total mass
|
|
// Initialize first mass and distance, actual number doesn't matter since total mass will be divided proportionally
|
|
float massTotal = distances[0] = masses[0] = 1;
|
|
|
|
// Masses and distances multiply from one planet to the next
|
|
// Mass typically increases exponentially as you go further out
|
|
// jupiterJump is an additional mass multiplier applied after half of the planets
|
|
for (var i = 1; i < count; i++)
|
|
{
|
|
massTotal += masses[i] = masses[i - 1] * Context.Random.NextFloat(massMulMin, massMulMax) * (count/2==i ? jupiterJump : 1);
|
|
distances[i] = distances[i - 1] * Context.Random.NextFloat(distMulMin, distMulMax);
|
|
}
|
|
|
|
// Add some randomness to child masses
|
|
for (var i = 0; i < masses.Length; i++)
|
|
masses[i] *= Context.Random.NextFloat(.1f, 1f);
|
|
|
|
// Normalize the masses of the children
|
|
for (var i = 0; i < count; i++)
|
|
masses[i] = masses[i] / massTotal * Mass * massFraction;
|
|
|
|
// Map child distances to range between minimum and maximum
|
|
if (count > 1)
|
|
{
|
|
var oldDistances = (float[]) distances.Clone();
|
|
for (var i = 0; i < count; i++)
|
|
distances[i] = lerp(ChildDistanceMinimum, ChildDistanceMaximum,
|
|
(oldDistances[i] - oldDistances[0]) / (oldDistances[count - 1] - oldDistances[0])) + Context.GlobalData.PlanetRadius(masses[i]);
|
|
}
|
|
|
|
for (var i = 0; i < count; i++)
|
|
{
|
|
// Only instantiate children above the mass floor
|
|
if (masses[i] > Context.GlobalData.MassFloor)
|
|
{
|
|
var child = new GeneratorPlanet
|
|
{
|
|
Context = Context,
|
|
Parent = this,
|
|
Mass = masses[i],
|
|
Distance = distances[i],
|
|
Phase = Context.Random.NextFloat()
|
|
};
|
|
child.Period = Context.GlobalData.OrbitalPeriod(child.Distance);
|
|
child.ChildDistanceMinimum = Context.GlobalData.PlanetRadius(child.Mass) * 2;
|
|
// Maximum child distance of child is the smallest distance to either of its neighbors
|
|
child.ChildDistanceMaximum = min(i == 0 ? child.Distance - ChildDistanceMinimum : child.Distance - distances[i - 1],
|
|
i < count - 1 ? distances[i + 1] - child.Distance : float.PositiveInfinity);
|
|
if (float.IsNaN(child.Distance))
|
|
throw new NotFiniteNumberException($"Planet created with NaN distance, something went very wrong!");
|
|
if (child.ChildDistanceMaximum > child.ChildDistanceMinimum)
|
|
Children.Add(child);
|
|
}
|
|
}
|
|
}
|
|
} |