Feature/maneuvers (#79)

* Implement maneuvers * Implement attitudes * Upgrade API and Scenario * Improve test coverage * Fix Cosmographia export * Fix orientation computation * Fix warnings * Upgrade native astrodynamics framework
IO-Aerospace-software-engineering · Mar 5, 2024 · 2622e1d · 2622e1d
1 parent e3a1097
commit 2622e1d
Show file tree

Hide file tree

Showing 153 changed files with 2,278 additions and 2,294 deletions.
diff --git a/IO.Astrodynamics.CLI.Tests/BodyInformationTests.cs b/IO.Astrodynamics.CLI.Tests/BodyInformationTests.cs
@@ -1,6 +1,5 @@
 using System.Text;
 using IO.Astrodynamics.CLI.Commands;
-using IO.Astrodynamics.CLI.Commands.Parameters;
 
 namespace IO.Astrodynamics.CLI.Tests;
 

diff --git a/IO.Astrodynamics.CLI.Tests/FrameTests.cs b/IO.Astrodynamics.CLI.Tests/FrameTests.cs
@@ -1,7 +1,6 @@
 using System.Text;
 using IO.Astrodynamics.CLI.Commands;
 using IO.Astrodynamics.CLI.Commands.Parameters;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Tests;
 

diff --git a/IO.Astrodynamics.CLI.Tests/OrbitalParametersTests.cs b/IO.Astrodynamics.CLI.Tests/OrbitalParametersTests.cs
@@ -1,7 +1,6 @@
 using System.Text;
 using IO.Astrodynamics.CLI.Commands;
 using IO.Astrodynamics.CLI.Commands.Parameters;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Tests;
 

diff --git a/IO.Astrodynamics.CLI.Tests/PropagateTests.cs b/IO.Astrodynamics.CLI.Tests/PropagateTests.cs
@@ -1,6 +1,5 @@
 // Copyright 2024. Sylvain Guillet (sylvain.guillet@tutamail.com)
 
-using System.Globalization;
 using System.Text;
 using IO.Astrodynamics.CLI.Commands;
 using IO.Astrodynamics.CLI.Commands.Parameters;

diff --git a/IO.Astrodynamics.CLI/Commands/BodyInformationCommand.cs b/IO.Astrodynamics.CLI/Commands/BodyInformationCommand.cs
@@ -1,14 +1,8 @@
 using System;
-using System.Globalization;
 using System.IO;
 using System.Threading.Tasks;
 using Cocona;
 using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
-using IO.Astrodynamics.CLI.Commands.Parameters;
-using IO.Astrodynamics.Frames;
-using IO.Astrodynamics.OrbitalParameters;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;
 

diff --git a/IO.Astrodynamics.CLI/Commands/EphemerisCommand.cs b/IO.Astrodynamics.CLI/Commands/EphemerisCommand.cs
@@ -4,11 +4,8 @@
 using System.Threading.Tasks;
 using Cocona;
 using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
 using IO.Astrodynamics.Frames;
-using IO.Astrodynamics.Surface;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;
 

diff --git a/IO.Astrodynamics.CLI/Commands/FrameConverterCommand.cs b/IO.Astrodynamics.CLI/Commands/FrameConverterCommand.cs
@@ -1,14 +1,9 @@
 using System;
-using System.Globalization;
 using System.IO;
 using System.Threading.Tasks;
 using Cocona;
-using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
 using IO.Astrodynamics.Frames;
-using IO.Astrodynamics.OrbitalParameters;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;
 

diff --git a/IO.Astrodynamics.CLI/Commands/GeometryFinderCommand.cs b/IO.Astrodynamics.CLI/Commands/GeometryFinderCommand.cs
@@ -2,11 +2,8 @@
 using System.IO;
 using System.Threading.Tasks;
 using Cocona;
-using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
 using IO.Astrodynamics.Frames;
-using IO.Astrodynamics.OrbitalParameters;
 using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;

diff --git a/IO.Astrodynamics.CLI/Commands/OrbitalParametersConverterCommand.cs b/IO.Astrodynamics.CLI/Commands/OrbitalParametersConverterCommand.cs
@@ -1,15 +1,9 @@
 using System;
-using System.Globalization;
 using System.IO;
 using System.Threading.Tasks;
 using Cocona;
-using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
 using IO.Astrodynamics.Frames;
-using IO.Astrodynamics.Math;
-using IO.Astrodynamics.OrbitalParameters;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;
 

diff --git a/IO.Astrodynamics.CLI/Commands/OrientationCommand.cs b/IO.Astrodynamics.CLI/Commands/OrientationCommand.cs
@@ -3,8 +3,6 @@
 using System.IO;
 using System.Threading.Tasks;
 using Cocona;
-using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
 using IO.Astrodynamics.Frames;
 using IO.Astrodynamics.OrbitalParameters;

diff --git a/IO.Astrodynamics.CLI/Commands/Parameters/WindowParameters.cs b/IO.Astrodynamics.CLI/Commands/Parameters/WindowParameters.cs
@@ -1,7 +1,5 @@
 // Copyright 2024. Sylvain Guillet (sylvain.guillet@tutamail.com)
 
-using System;
-using System.Text;
 using Cocona;
 
 namespace IO.Astrodynamics.CLI.Commands.Parameters;

diff --git a/IO.Astrodynamics.CLI/Commands/PropagateCommand.cs b/IO.Astrodynamics.CLI/Commands/PropagateCommand.cs
@@ -1,19 +1,14 @@
 using System;
 using System.Collections.Generic;
-using System.Globalization;
 using System.IO;
 using System.Linq;
-using System.Threading;
 using System.Threading.Tasks;
 using Cocona;
 using IO.Astrodynamics.Body;
 using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
 using IO.Astrodynamics.Cosmographia;
-using IO.Astrodynamics.Frames;
 using IO.Astrodynamics.Mission;
-using IO.Astrodynamics.OrbitalParameters;
-using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;
 
@@ -64,7 +59,7 @@ public async Task Propagate(
         //Simulate
         var qualifiedOutputDirectory = new DirectoryInfo(outputDirectory);
         Console.WriteLine("Orbit propagation in progress. This operation could take a long time. Duration expected for one orbit 7000km x 7000km on mid range computer with geopotential N=30 = 600ms");
-        await scenario.SimulateWithoutManeuverAsync(qualifiedOutputDirectory, useAtmosphericDrag, useSolarRadiationPressure);
+        await scenario.SimulateAsync(qualifiedOutputDirectory, useAtmosphericDrag, useSolarRadiationPressure,TimeSpan.FromSeconds(1.0));
 
         //Export to cosmographia
         CosmographiaExporter cosmographiaExporter = new CosmographiaExporter();

diff --git a/IO.Astrodynamics.CLI/Commands/TimeConverterCommand.cs b/IO.Astrodynamics.CLI/Commands/TimeConverterCommand.cs
@@ -2,11 +2,7 @@
 using System.Globalization;
 using System.Threading.Tasks;
 using Cocona;
-using IO.Astrodynamics.Body;
-using IO.Astrodynamics.Body.Spacecraft;
 using IO.Astrodynamics.CLI.Commands.Parameters;
-using IO.Astrodynamics.Frames;
-using IO.Astrodynamics.OrbitalParameters;
 using IO.Astrodynamics.Time;
 
 namespace IO.Astrodynamics.CLI.Commands;

diff --git a/IO.Astrodynamics.CLI/Helpers.cs b/IO.Astrodynamics.CLI/Helpers.cs
@@ -45,7 +45,7 @@ internal static CelestialItem CreateCelestialItem(int objectId)
         }
         else
         {
-            localizableObject = new Spacecraft(objectId, $"spc{objectId}", 1.0, 1.0, new Clock($"clk{objectId}", 1 / 65536.0), null);
+            localizableObject = new Spacecraft(objectId, $"spc{objectId}", 1.0, 1.0, new Clock($"clk{objectId}", 65536), null);
         }
 
         return localizableObject;
@@ -60,7 +60,7 @@ internal static IOrientable CreateOrientable(int objectId)
         }
         else
         {
-            celestialItem = new Spacecraft(objectId, "spc", 1.0, 1.0, new Clock("clk", 1 / 65536.0), null);
+            celestialItem = new Spacecraft(objectId, "spc", 1.0, 1.0, new Clock("clk", 65536), null);
         }
 
         return celestialItem;
@@ -113,6 +113,7 @@ internal static OrbitalParameters.OrbitalParameters ConvertToOrbitalParameters(s
         {
             originalFrame = "j2000";
         }
+
         var inputFrame = new Frame(originalFrame);
         var inputEpoch = Helpers.ConvertDateTimeInput(epoch);
 
@@ -173,11 +174,11 @@ internal static Window ConvertWindowInput(string begin, string end)
         return new Window(ConvertDateTimeInput(begin), ConvertDateTimeInput(end));
     }
 
-    internal static Window ConvertWindowInput(EpochParameters begin,EpochParameters end)
+    internal static Window ConvertWindowInput(EpochParameters begin, EpochParameters end)
     {
         return ConvertWindowInput(begin.Epoch, end.Epoch);
     }
-    
+
     internal static Window ConvertWindowInput(WindowParameters windowParameters)
     {
         return ConvertWindowInput(windowParameters.Begin, windowParameters.End);

diff --git a/IO.Astrodynamics.CLI/IO.Astrodynamics.CLI.csproj b/IO.Astrodynamics.CLI/IO.Astrodynamics.CLI.csproj
@@ -9,7 +9,7 @@
         <FileVersion>0.0.1</FileVersion>
         <PackAsTool>true</PackAsTool>
         <ToolCommandName>astro</ToolCommandName>
-        <Version>0.0.1</Version>
+        <Version>0.1.0</Version>
         <Title>Astrodynamics command line interface</Title>
         <Authors>Sylvain Guillet</Authors>
         <Description>This CLI allows end user to exploit IO.Astrodynamics framework </Description>

diff --git a/IO.Astrodynamics.Performance/Scenario.cs b/IO.Astrodynamics.Performance/Scenario.cs
@@ -19,7 +19,7 @@ namespace IO.Astrodynamics.Performance;
 [SkewnessColumn]
 [KurtosisColumn]
 [StatisticalTestColumn]
-[ShortRunJob]
+// [ShortRunJob]
 public class Scenario
 {
     private readonly GeopotentialGravitationalField _geopotential;
@@ -29,7 +29,7 @@ public class Scenario
     private readonly CelestialBody _sun;
     private readonly CelestialBody _moon;
     private readonly VVIntegrator _integrator;
-    private readonly Propagator.Propagator _propagator;
+    private readonly Propagator.SpacecraftPropagator _spacecraftPropagator;
     // IO.Astrodynamics.Tests.Mission.ScenarioTests _scenario = new IO.Astrodynamics.Tests.Mission.ScenarioTests();
 
     public Scenario()
@@ -39,7 +39,7 @@ public Scenario()
         _moon = new CelestialBody(301);
         _sun = new CelestialBody(10);
         _geopotential = new GeopotentialGravitationalField(new StreamReader("Data/SolarSystem/EGM2008_to70_TideFree"));
-        Clock clk = new Clock("My clock", 1.0 / 256.0);
+        Clock clk = new Clock("My clock", 256);
         Spacecraft spc = new Spacecraft(-1001, "MySpacecraft", 100.0, 10000.0, clk,
             new StateVector(new Vector3(6800000.0, 0.0, 0.0), new Vector3(0.0, 7656.2204182967143, 0.0), _earth, DateTimeExtension.J2000, Frames.Frame.ICRF));
         _srp = new SolarRadiationPressure(spc);
@@ -53,11 +53,11 @@ public Scenario()
         _integrator = new VVIntegrator(forces, TimeSpan.FromSeconds(1.0), new StateVector(new Vector3(6800000.0 - Random.Shared.NextDouble(), 0.0, 0.0),
             new Vector3(0.0, 8000.0 - Random.Shared.NextDouble(), 0.0), _earth,
             DateTimeExtension.J2000, Frame.ICRF));
-        _propagator = new Propagator.Propagator(new Window(DateTimeExtension.J2000, DateTimeExtension.J2000 + spc.InitialOrbitalParameters.Period()), spc,
+        _spacecraftPropagator = new Propagator.SpacecraftPropagator(new Window(DateTimeExtension.J2000, DateTimeExtension.J2000 + spc.InitialOrbitalParameters.Period()), spc,
             new[] { _moon, _earth, _sun }, true, true, TimeSpan.FromSeconds(1.0));
     }
 
-    // [Benchmark(Description = "Spacecraft propagator C++")]
+    // [Benchmark(Description = "Spacecraft spacecraftPropagator C++")]
     public void Propagate()
     {
         // _scenario.PropagateWithoutManeuver();
@@ -96,13 +96,13 @@ public void VVIntegration()
         // var res = _integrator.Integrate(sv);
     }
 
-    // [Benchmark(Description = "Propagator per orbit (GeoPotentials // Moon and sun perturbation // Atmospheric drag // Solar radiation) ")]
+    [Benchmark(Description = "SpacecraftPropagator per orbit (GeoPotentials // Moon and sun perturbation // Atmospheric drag // Solar radiation) ")]
     public void Propagator()
     {
-        var res = _propagator.Propagate();
+        var res = _spacecraftPropagator.Propagate();
     }
 
-    [Benchmark(Description = "IO Vector", Baseline = true)]
+    // [Benchmark(Description = "IO Vector", Baseline = true)]
     public void Vector()
     {
         Vector3 vectorA = new Vector3(1.0, 2.0, 3.0);
@@ -169,7 +169,7 @@ public void Vector()
         mag = tot.Magnitude();
     }
 
-    [Benchmark(Description = "Numerics Vector")]
+    // [Benchmark(Description = "Numerics Vector")]
     public void VectorNumerics()
     {
         var vectorA = new Vector<double>([0.0, 1.0, 2.0, 3.0]);