using System.Linq; using System.Numerics.Tensors; using System.Runtime.CompilerServices; using Content.Server.Atmos.Reactions; using Content.Shared.Atmos; using Content.Shared.Atmos.Reactions; using JetBrains.Annotations; namespace Content.Server.Atmos.EntitySystems { public sealed partial class AtmosphereSystem { private GasReactionPrototype[] _gasReactions = []; /// /// List of gas reactions ordered by priority. /// public IEnumerable GasReactions => _gasReactions; public override void InitializeGases() { base.InitializeGases(); } /// /// Caches all gas reactions into an array ordered by priority. /// public void CacheGases() { _gasReactions = ProtoMan.EnumeratePrototypes().ToArray(); Array.Sort(_gasReactions, (a, b) => b.Priority.CompareTo(a.Priority)); } public override float GetMass(GasMixture mix) { return GetMass(mix.Moles); } public override float GetMass(float[] moles) { Span tmp = stackalloc float[moles.Length]; TensorPrimitives.Multiply(moles, GasMolarMasses, tmp); // Conversion of grams to kilograms. return TensorPrimitives.Sum(tmp) * Atmospherics.gToKg; } [MethodImpl(MethodImplOptions.AggressiveInlining)] protected override float GetHeatCapacityCalculation(float[] moles, bool space) { // Little hack to make space gas mixtures have heat capacity, therefore allowing them to cool down rooms. if (space && MathHelper.CloseTo(TensorPrimitives.Sum(moles), 0f)) { return Atmospherics.SpaceHeatCapacity; } Span tmp = stackalloc float[moles.Length]; TensorPrimitives.Multiply(moles, GasMolarHeatCapacities, tmp); // Adjust heat capacity by speedup, because this is primarily what // determines how quickly gases heat up/cool. return MathF.Max(TensorPrimitives.Sum(tmp), Atmospherics.MinimumHeatCapacity); } public override bool IsMixtureFuel(GasMixture mixture, float epsilon = Atmospherics.Epsilon) { Span tmp = stackalloc float[Atmospherics.AdjustedNumberOfGases]; TensorPrimitives.Multiply(mixture.Moles, GasFuelMask, tmp); return TensorPrimitives.Sum(tmp) > epsilon; } public override bool IsMixtureOxidizer(GasMixture mixture, float epsilon = Atmospherics.Epsilon) { Span tmp = stackalloc float[Atmospherics.AdjustedNumberOfGases]; TensorPrimitives.Multiply(mixture.Moles, GasOxidizerMask, tmp); return TensorPrimitives.Sum(tmp) > epsilon; } /// /// Return speedup factor for pumped or flow-based devices that depend on MaxTransferRate. /// public float PumpSpeedup() { return Speedup; } /// /// Add 'dQ' Joules of energy into 'mixture'. /// public void AddHeat(GasMixture mixture, float dQ) { var c = GetHeatCapacity(mixture); float dT = dQ / c; mixture.Temperature += dT; } /// /// Divides a source gas mixture into several recipient mixtures, scaled by their relative volumes. Does not /// modify the source gas mixture. Used for pipe network splitting. Note that the total destination volume /// may be larger or smaller than the source mixture. /// public void DivideInto(GasMixture source, List receivers) { var totalVolume = 0f; foreach (var receiver in receivers) { if (!receiver.Immutable) totalVolume += receiver.Volume; } float? sourceHeatCapacity = null; var buffer = new float[Atmospherics.AdjustedNumberOfGases]; foreach (var receiver in receivers) { if (receiver.Immutable) continue; var fraction = receiver.Volume / totalVolume; // Set temperature, if necessary. if (MathF.Abs(receiver.Temperature - source.Temperature) > Atmospherics.MinimumTemperatureDeltaToConsider) { // Often this divides a pipe net into new and completely empty pipe nets if (receiver.TotalMoles == 0) receiver.Temperature = source.Temperature; else { sourceHeatCapacity ??= GetHeatCapacity(source); var receiverHeatCapacity = GetHeatCapacity(receiver); var combinedHeatCapacity = receiverHeatCapacity + sourceHeatCapacity.Value * fraction; if (combinedHeatCapacity > Atmospherics.MinimumHeatCapacity) receiver.Temperature = (GetThermalEnergy(source, sourceHeatCapacity.Value * fraction) + GetThermalEnergy(receiver, receiverHeatCapacity)) / combinedHeatCapacity; } } // transfer moles TensorPrimitives.Multiply(source.Moles, fraction, buffer); TensorPrimitives.Add(receiver.Moles, buffer, receiver.Moles); } } /// /// Releases gas from this mixture to the output mixture. /// If the output mixture is null, then this is being released into space. /// It can't transfer air to a mixture with higher pressure. /// public bool ReleaseGasTo(GasMixture mixture, GasMixture? output, float targetPressure) { var outputStartingPressure = output?.Pressure ?? 0; var inputStartingPressure = mixture.Pressure; if (outputStartingPressure >= MathF.Min(targetPressure, inputStartingPressure - 10)) // No need to pump gas if the target is already reached or input pressure is too low. // Need at least 10 kPa difference to overcome friction in the mechanism. return false; if (!(mixture.TotalMoles > 0) || !(mixture.Temperature > 0)) return false; // We calculate the necessary moles to transfer with the ideal gas law. var pressureDelta = MathF.Min(targetPressure - outputStartingPressure, (inputStartingPressure - outputStartingPressure) / 2f); var transferMoles = pressureDelta * (output?.Volume ?? Atmospherics.CellVolume) / (mixture.Temperature * Atmospherics.R); // And now we transfer the gas. var removed = mixture.Remove(transferMoles); if(output != null) Merge(output, removed); return true; } /// /// Pump gas from this mixture to the output mixture. /// Amount depends on target pressure. /// /// The mixture to pump the gas from /// The mixture to pump the gas to /// The target pressure to reach /// Whether we could pump air to the output or not public bool PumpGasTo(GasMixture mixture, GasMixture output, float targetPressure) { var outputStartingPressure = output.Pressure; var pressureDelta = targetPressure - outputStartingPressure; if (pressureDelta < 0.01) // No need to pump gas, we've reached the target. return false; if (!(mixture.TotalMoles > 0) || !(mixture.Temperature > 0)) return false; // We calculate the necessary moles to transfer with the ideal gas law. var transferMoles = pressureDelta * output.Volume / (mixture.Temperature * Atmospherics.R); // And now we transfer the gas. var removed = mixture.Remove(transferMoles); Merge(output, removed); return true; } /// /// Scrubs specified gases from a gas mixture into a gas mixture. /// public void ScrubInto(GasMixture mixture, GasMixture destination, IReadOnlyCollection filterGases) { var buffer = new GasMixture(mixture.Volume){Temperature = mixture.Temperature}; foreach (var gas in filterGases) { buffer.AdjustMoles(gas, mixture.GetMoles(gas)); mixture.SetMoles(gas, 0f); } Merge(destination, buffer); } /// /// Checks whether a gas mixture is probably safe. /// This only checks temperature and pressure, not gas composition. /// /// Mixture to be checked. /// Whether the mixture is probably safe. public bool IsMixtureProbablySafe(GasMixture? air) { // Note that oxygen mix isn't checked, but survival boxes make that not necessary. if (air == null) return false; switch (air.Pressure) { case <= Atmospherics.WarningLowPressure: case >= Atmospherics.WarningHighPressure: return false; } switch (air.Temperature) { case <= 260: case >= 360: return false; } return true; } /// /// Compares two TileAtmospheres to see if they are within acceptable ranges for group processing to be enabled. /// public GasCompareResult CompareExchange(TileAtmosphere sample, TileAtmosphere otherSample) { if (sample.AirArchived == null || otherSample.AirArchived == null) return GasCompareResult.NoExchange; return CompareExchange(sample.AirArchived, otherSample.AirArchived); } /// /// Compares two gas mixtures to see if they are within acceptable ranges for group processing to be enabled. /// public GasCompareResult CompareExchange(GasMixture sample, GasMixture otherSample) { var moles = 0f; for(var i = 0; i < Atmospherics.TotalNumberOfGases; i++) { var gasMoles = sample.Moles[i]; var delta = MathF.Abs(gasMoles - otherSample.Moles[i]); if (delta > Atmospherics.MinimumMolesDeltaToMove && (delta > gasMoles * Atmospherics.MinimumAirRatioToMove)) return (GasCompareResult)i; // We can move gases! moles += gasMoles; } if (moles > Atmospherics.MinimumMolesDeltaToMove) { var tempDelta = MathF.Abs(sample.Temperature - otherSample.Temperature); if (tempDelta > Atmospherics.MinimumTemperatureDeltaToSuspend) return GasCompareResult.TemperatureExchange; // There can be temperature exchange. } // No exchange at all! return GasCompareResult.NoExchange; } [PublicAPI] public override ReactionResult React(GasMixture mixture, IGasMixtureHolder? holder) { var reaction = ReactionResult.NoReaction; var temperature = mixture.Temperature; var energy = GetThermalEnergy(mixture); foreach (var prototype in GasReactions) { if (energy < prototype.MinimumEnergyRequirement || temperature < prototype.MinimumTemperatureRequirement || temperature > prototype.MaximumTemperatureRequirement) continue; var doReaction = true; for (var i = 0; i < Atmospherics.TotalNumberOfGases; i++) { var req = prototype.MinimumRequirements[i]; if (!(mixture.GetMoles(i) < req)) continue; doReaction = false; break; } if (!doReaction) continue; reaction = prototype.React(mixture, holder, this, HeatScale); if(reaction.HasFlag(ReactionResult.StopReactions)) break; } return reaction; } /// /// Adds an array of moles to a . /// Guards against negative moles by clamping to zero. /// /// The to add moles to. /// The of moles to add. /// Thrown when the length of the /// is not the same as the length of the gas array. [PublicAPI] public static void AddMolsToMixture(GasMixture mixture, ReadOnlySpan molsToAdd) { TensorPrimitives.Add(mixture.Moles, molsToAdd, mixture.Moles); TensorPrimitives.Max(mixture.Moles, 0f, mixture.Moles); } public enum GasCompareResult { NoExchange = -2, TemperatureExchange = -1, } } }