don't extrapolate points in powermap that have to neighbours with higher values

VectoCore/VectoCore/InputData/Reader/ComponentData/IEPCMapReader.cs

@@ -170,20 +170,52 @@ namespace TUGraz.VectoCore.InputData.Reader.ComponentData
 			}
 
 
-			
+			//Don't extrapolate speedbuckets which have a left AND right neighbour with significantly higher torque values
+			//
+
+			var ignoredSpeedBucketsRecuperation = new HashSet<PerSecond>();
+			var ignoredSpeedBucketsDrive = new HashSet<PerSecond>();
+			var ignoreThreshold = 0.8;
+			var orderedBuckets = speedBuckets.OrderBy(x => x.Key).ToList();
+			for (var i = 1; i < speedBuckets.Count - 1; i++) {
+				var current = orderedBuckets[i];
+				var prev = orderedBuckets[i - 1];
+				var next = orderedBuckets[i + 1];
+
+				//Drive
+				if (current.Value.MinBy(x => x.Torque).Torque / prev.Value.MinBy(x => x.Torque).Torque < ignoreThreshold
+					&& 
+					current.Value.MinBy(x => x.Torque).Torque / next.Value.MinBy(x => x.Torque).Torque < ignoreThreshold) {
+					ignoredSpeedBucketsDrive.Add(current.Key);
+				}
+
+                //Recuperation
+				if (current.Value.MaxBy(x => x.Torque).Torque / prev.Value.MaxBy(x => x.Torque).Torque < ignoreThreshold
+					&&
+					current.Value.MaxBy(x => x.Torque).Torque / next.Value.MaxBy(x => x.Torque).Torque < ignoreThreshold)
+				{
+					ignoredSpeedBucketsRecuperation.Add(current.Key);
+				}
+            }
+
+
+
 			var maxTargetTorque = fullLoadCurve.MaxGenerationTorque * extrapolationfactor;
 			var minTargetTorque = fullLoadCurve.MaxDriveTorque * extrapolationfactor;
 			var ratedSpeed = ElectricMotorRatedSpeedHelper.GetRatedSpeed(fullLoadCurve.FullLoadEntries,
 				e => e.MotorSpeed, e => e.FullDriveTorque);
 			PerSecond prevSpeed = null;
-			foreach (var speedBucket in speedBuckets.OrderBy(x => x.Key))
+			foreach (var speedBucket in orderedBuckets)
 			{
+				
+
 				var maxRecuperationEntry = speedBucket.Value.MaxBy(x => x.Torque);
 				var maxDriveEntry = speedBucket.Value.MinBy(x => x.Torque); //drive torque < 0
 				var recuperationFactor = maxTargetTorque / maxRecuperationEntry.Torque;
 				var driveFactor = minTargetTorque / maxDriveEntry.Torque;
 
-				if (!recuperationFactor.IsSmallerOrEqual(1)) {
+				//Recuperation
+				if (!recuperationFactor.IsSmallerOrEqual(1) && !ignoredSpeedBucketsRecuperation.Contains(speedBucket.Key)) {
 					var nrExtrapolationPointsRecuperation = (uint)Math.Ceiling(speedBucket.Value.Count(x => x.Torque.IsGreater(0)) * (recuperationFactor - 1));
 					for (var i = 1; i <= nrExtrapolationPointsRecuperation; i++)
 					{
@@ -193,7 +225,9 @@ namespace TUGraz.VectoCore.InputData.Reader.ComponentData
 					}
 				}
 
-				if (!driveFactor.IsSmallerOrEqual(1)) {
+				//Drive
+				if (!driveFactor.IsSmallerOrEqual(1) && !ignoredSpeedBucketsDrive.Contains(speedBucket.Key)) {
+
 					var nrExtrapolationPointsDrive = (uint)Math.Ceiling(speedBucket.Value.Count(x => x.Torque.IsSmaller(0)) * (driveFactor - 1));
 
 					for (var i = 1; i <= nrExtrapolationPointsDrive; i++)

VectoCore/VectoCoreTest/InputData/IEPCMapReaderTest.cs

@@ -28,27 +28,31 @@ namespace TUGraz.VectoCore.Tests.InputData
 	[TestFixture]
 	public class IEPCMapReaderTest
 	{
-		public const string basePath = "../../../../../"; //from cwd to repo root
+		public const string repoRoot = "../../../../../"; //from cwd to repo root
+		public const string iepcMapTestFiles = @"TestData\Components\IEPC\IEPCMapReader\";
 
 
 
 
-
-		[OneTimeSetUp]
+        [OneTimeSetUp]
 		public void OneTimeSetup()
 		{
 			
 		}
 
 
-		[TestCase(@$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_2.viepco", 4.65f, 2.74f)]
-		[TestCase(@$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_1.viepco", 8.31f, 2.74f)]
-		[TestCase(@$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_3.viepco", 2.74f, 2.74f)]
+		[TestCase(@$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_2.viepco", 4.65f, 2.74f)]
+		[TestCase(@$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_1.viepco", 8.31f, 2.74f)]
+		[TestCase(@$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_3.viepco", 2.74f, 2.74f)]
+
+
+		[TestCase(@$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{repoRoot}VectoCore/VectoCoreTest/TestData/Components/IEPC/IEPCMapReader/IEPC_Gbx3_1_different_number_of_points.viepco", 8.31f, 2.74f)]
+		[TestCase(@$"{repoRoot}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{repoRoot}VectoCore/VectoCoreTest/TestData/Components/IEPC/IEPCMapReader/IEPC_Gbx3_1_following_curve.viepco", 8.31f, 2.74f)]
 
 
-		[TestCase(@$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{basePath}VectoCore/VectoCoreTest/TestData/Components/IEPC/IEPCMapReader/IEPC_Gbx3_1_different_number_of_points.viepco", 8.31f, 2.74f)]
-		[TestCase(@$"{basePath}Generic Vehicles/Engineering Mode/GenericIEPC-S/IEPC-S_Gbx3Speed/IEPC_Gbx3_FLD_max.viepcp", @$"{basePath}VectoCore/VectoCoreTest/TestData/Components/IEPC/IEPCMapReader/IEPC_Gbx3_1_following_curve.viepco", 8.31f, 2.74f)]
-		public void ReadIEPCMap(string fullLoadCurvePath, string powerMapPath, double ratio, double fldMeasuredRatio)
+		[TestCase($"{iepcMapTestFiles}IEPC_Gbx1_FLD_max.viepcp",$"{iepcMapTestFiles}IEPC_Gbx1.viepco", 4.65f, 4.65f)]
+		[TestCase($"{iepcMapTestFiles}IEPC_Gbx1_FLD_max.viepcp", $"{iepcMapTestFiles}IEPC_Gbx1_below_fld_outlierpoints.viepco", 4.65f, 4.65f)]
+        public void ReadIEPCMap(string fullLoadCurvePath, string powerMapPath, double ratio, double fldMeasuredRatio)
 		{
 			Assert.That(File.Exists(fullLoadCurvePath), Path.GetFullPath(fullLoadCurvePath));
 			Assert.That(File.Exists(powerMapPath), Path.GetFullPath(fullLoadCurvePath));
@@ -63,7 +67,6 @@ namespace TUGraz.VectoCore.Tests.InputData
 
 			var fld = IEPCFullLoadCurveReader.Create(fullLoadCurveData, 1, fldMeasuredRatio);
 			var powerMapInput = IEPCMapReader.GetEntries(powerMapData, ratio);
-			//PrintMaps("powerMapInput", powerMapInput, "powerMapInput", x => x.MotorSpeed.AsRPM, y => y.Torque.Value());
 
 
 			var effMap = IEPCMapReader.Create(powerMapData, 1, ratio, fld);
@@ -84,10 +87,6 @@ namespace TUGraz.VectoCore.Tests.InputData
 
 
 
-
-			//PrintMaps("powerMap", effMap.Entries, "powerMap", x => x.MotorSpeed.AsRPM, y => y.Torque.Value());
-
-
 			PrintMaps("FLD_input_powermap.png",
 				new SeriesProperties<ElectricMotorFullLoadCurve.FullLoadEntry>()
 				{
@@ -229,7 +228,7 @@ namespace TUGraz.VectoCore.Tests.InputData
 				
 
 
-				var dirInfo = Directory.CreateDirectory($@"{nameof(IEPCMapReaderTest)}//{String.Join("", TestContext.CurrentContext.Test.Name.Replace(basePath, "").Split(Path.GetInvalidFileNameChars()))}");
+				var dirInfo = Directory.CreateDirectory($@"{nameof(IEPCMapReaderTest)}//{String.Join("", TestContext.CurrentContext.Test.Name.Replace(repoRoot, "").Split(Path.GetInvalidFileNameChars()))}");
 				chart.SaveImage($"{dirInfo.FullName}\\{(fileName.EndsWith(".png") ? fileName : fileName + ".png")}",
 					ChartImageFormat.Png);
 				TestContext.WriteLine($"{dirInfo.FullName}");

VectoCore/VectoCoreTest/TestData/Components/IEPC/IEPCMapReader/IEPC_Gbx1_FLD_max.viepcp

+n_out   , T_drive_out , T_recuperation_out
+0.00    , 6765.75     , -7190.72
+8.82    , 6765.75     , -7190.72
+89.03   , 6765.75     , -7190.72
+178.06  , 6765.75     , -7190.72
+266.88  , 6765.75     , -7190.72
+355.91  , 6765.75     , -7190.72
+444.95  , 6765.75     , -7190.72
+533.98  , 5638.13     , -5992.27
+623.01  , 4832.64     , -5136.19
+711.83  , 4228.59     , -4494.20
+800.86  , 3758.78     , -3994.88
+889.89  , 3382.88     , -3595.36
+978.92  , 3075.35     , -3268.52
+1067.96 , 2819.06     , -2996.13
+1156.77 , 2602.20     , -2765.65
+1245.81 , 2416.37     , -2568.14
+1334.84 , 2255.25     , -2396.91
+1423.87 , 2114.34     , -2247.15
+1512.90 , 1989.94     , -2114.94
+1601.72 , 1879.34     , -1997.39
+1690.75 , 1780.47     , -1892.31
+1779.78 , 1691.44     , -1797.68

VectoCore/VectoCoreTest/TestData/Components/IEPC/IEPCMapReader/IEPC_Gbx1_below_fld_outlierpoints.viepco

+n_out   , T_out    , P_el
+0,-7190.72,0
+0,-6831.19,0
+0,-6471.65,0
+0,-6112.11,0
+0,-5752.58,0
+0,-5393.04,0
+0,-5033.51,-0.13
+0,-4673.97,-0.66
+0,-4314.43,-0.98
+0,-3954.9,-1.19
+0,-3595.36,-1.32
+0,-3235.82,-1.38
+0,-2876.29,-1.39
+0,-2516.75,-1.35
+0,-2157.22,-1.26
+0,-1797.68,-1.13
+0,-1438.14,-0.97
+0,-1078.61,-0.77
+0,-719.07,-0.53
+0,-359.54,-0.27
+0,-71.91,-0.02
+0,67.66,0.13
+0,338.29,0.56
+0,676.58,1.12
+0,1014.86,1.7
+0,1353.15,2.31
+0,1691.44,2.93
+0,2029.73,3.58
+0,2368.01,4.25
+0,2706.3,4.93
+0,3044.59,5.64
+0,3382.88,6.37
+0,3721.16,7.12
+0,4059.45,7.88
+0,4397.74,8.67
+0,4736.03,9.47
+0,5074.31,10.29
+0,5412.6,11.13
+0,5750.89,11.99
+0,6089.18,12.87
+0,6427.46,13.76
+0,6765.75,14.67
+8.9,-7190.72,0
+8.9,-6831.19,0
+8.9,-6471.65,0
+8.9,-6112.11,0
+8.9,-5752.58,0
+8.9,-5393.04,0
+8.9,-5033.51,-0.5
+8.9,-4673.97,-0.84
+8.9,-4314.43,-1.04
+8.9,-3954.9,-1.16
+8.9,-3595.36,-1.22
+8.9,-3235.82,-1.24
+8.9,-2876.29,-1.21
+8.9,-2516.75,-1.15
+8.9,-2157.22,-1.06
+8.9,-1797.68,-0.94
+50,-1438.14,-0.8
+8.9,-1078.61,-0.62
+8.9,-719.07,-0.43
+8.9,-359.54,-0.21
+8.9,-71.91,-0.02
+8.9,67.66,0.1
+8.9,338.29,0.43
+8.9,676.58,0.85
+8.9,1014.86,1.3
+8.9,1353.15,1.75
+8.9,1691.44,2.23
+8.9,2029.73,2.71
+8.9,2368.01,3.21
+8.9,2706.3,3.73
+8.9,3044.59,4.26
+8.9,3382.88,4.8
+8.9,3721.16,5.36
+8.9,4059.45,5.93
+8.9,4397.74,6.52
+8.9,4736.03,7.12
+8.9,5074.31,7.73
+8.9,5412.6,8.35
+8.9,5750.89,8.99
+8.9,6089.18,9.64
+8.9,6427.46,10.3
+8.9,6765.75,10.98
+88.99,-7190.72,-58.1
+88.99,-6831.19,-55.37
+88.99,-6471.65,-52.62
+88.99,-6112.11,-49.85
+88.99,-5752.58,-47.06
+88.99,-5393.04,-44.25
+88.99,-5033.51,-41.42
+88.99,-4673.97,-38.57
+88.99,-4314.43,-35.7
+88.99,-3954.9,-32.81
+88.99,-3595.36,-29.9
+88.99,-3235.82,-26.98
+88.99,-2876.29,-24.03
+88.99,-2516.75,-21.07
+88.99,-2157.22,-18.08
+88.99,-1797.68,-15.08
+88.99,-1438.14,-12.06
+88.99,-1078.61,-9.02
+88.99,-719.07,-5.96
+88.99,-359.54,-2.88
+88.99,-71.91,-0.37
+88.99,67.66,0.9
+88.99,338.29,3.68
+88.99,676.58,7.16
+88.99,1014.86,10.66
+88.99,1353.15,14.17
+88.99,1691.44,17.71
+88.99,2029.73,21.26
+88.99,2368.01,24.84
+88.99,2706.3,28.43
+88.99,3044.59,32.04
+88.99,3382.88,35.67
+88.99,3721.16,39.32
+88.99,4059.45,42.99
+88.99,4397.74,46.68
+88.99,4736.03,50.38
+88.99,5074.31,54.11
+88.99,5412.6,57.85
+88.99,5750.89,61.61
+88.99,6089.18,65.39
+88.99,6427.46,69.19
+88.99,6765.75,73.01
+177.98,-7190.72,-120.57
+177.98,-6831.19,-114.73
+177.98,-6471.65,-108.87
+177.98,-6112.11,-102.99
+177.98,-5752.58,-97.08
+177.98,-5393.04,-91.15
+177.98,-5033.51,-85.2
+177.98,-4673.97,-79.23
+177.98,-4314.43,-73.23
+177.98,-3954.9,-67.21
+177.98,-3595.36,-61.17
+177.98,-3235.82,-55.1
+177.98,-2876.29,-49.01
+177.98,-2516.75,-42.9
+177.98,-2157.22,-36.77
+177.98,-1797.68,-30.61
+177.98,-1438.14,-24.43
+177.98,-1078.61,-18.23
+177.98,-719.07,-12.01
+177.98,-359.54,-5.76
+177.98,-71.91,-0.66
+177.98,67.66,1.84
+177.98,338.29,7.35
+177.98,676.58,14.22
+177.98,1014.86,21.11
+177.98,1353.15,28.02
+177.98,1691.44,34.96
+177.98,2029.73,41.92
+177.98,2368.01,48.9
+177.98,2706.3,55.9
+177.98,3044.59,62.93
+177.98,3382.88,69.98
+177.98,3721.16,77.05
+177.98,4059.45,84.15
+177.98,4397.74,91.27
+177.98,4736.03,98.41
+177.98,5074.31,105.58
+177.98,5412.6,112.76
+177.98,5750.89,119.97
+177.98,6089.18,127.21
+177.98,6427.46,134.46
+177.98,6765.75,141.74
+266.97,-7190.72,-182.44
+266.97,-6831.19,-173.55
+266.97,-6471.65,-164.62
+266.97,-6112.11,-155.67
+266.97,-5752.58,-146.69
+266.97,-5393.04,-137.68
+266.97,-5033.51,-128.65
+266.97,-4673.97,-119.58
+266.97,-4314.43,-110.49
+266.97,-3954.9,-101.36
+266.97,-3595.36,-92.21
+266.97,-3235.82,-83.03
+266.97,-2876.29,-73.82
+266.97,-2516.75,-64.58
+266.97,-2157.22,-55.31
+266.97,-1797.68,-46.02
+266.97,-1438.14,-36.7
+266.97,-1078.61,-27.34
+266.97,-719.07,-17.96
+266.97,-359.54,-8.54
+266.97,-71.91,-0.79
+266.97,67.66,2.86
+266.97,338.29,11.12
+266.97,676.58,21.39
+266.97,1014.86,31.68
+266.97,1353.15,41.99
+266.97,1691.44,52.34
+266.97,2029.73,62.71
+266.97,2368.01,73.12
+266.97,2706.3,83.55
+266.97,3044.59,94.01
+266.97,3382.88,104.5
+266.97,3721.16,115.02
+266.97,4059.45,125.57
+266.97,4397.74,136.15
+266.97,4736.03,146.76
+266.97,5074.31,157.4
+266.97,5412.6,168.06
+266.97,5750.89,178.76
+266.97,6089.18,189.48
+266.97,6427.46,200.23
+266.97,6765.75,211.02
+355.96,-7190.72,-243.71
+355.96,-6831.19,-231.81
+355.96,-6471.65,-219.87
+355.96,-6112.11,-207.89
+355.96,-5752.58,-195.88
+355.96,-5393.04,-183.83
+355.96,-5033.51,-171.74
+355.96,-4673.97,-159.61
+355.96,-4314.43,-147.45
+355.96,-3954.9,-135.25
+355.96,-3595.36,-123.01
+355.96,-3235.82,-110.74
+355.96,-2876.29,-98.43
+355.96,-2516.75,-86.08
+355.96,-2157.22,-73.7
+355.96,-1797.68,-61.27
+355.96,-1438.14,-48.82
+355.96,-1078.61,-36.32
+355.96,-719.07,-23.78
+355.96,-359.54,-11.19
+355.96,-71.91,-0.63
+355.96,67.66,3.98
+355.96,338.29,15.01
+355.96,676.58,28.68
+355.96,1014.86,42.37
+355.96,1353.15,56.1
+355.96,1691.44,69.87
+355.96,2029.73,83.67
+355.96,2368.01,97.52
+355.96,2706.3,111.4
+355.96,3044.59,125.31
+355.96,3382.88,139.27
+355.96,3721.16,153.26
+355.96,4059.45,167.29
+355.96,4397.74,181.36
+355.96,4736.03,195.46
+355.96,5074.31,209.61
+355.96,5412.6,223.79
+355.96,5750.89,238.01
+355.96,6089.18,252.26
+355.96,6427.46,266.56
+355.96,6765.75,280.89
+444.95,-7190.72,-304.37
+444.95,-6831.19,-289.5
+444.95,-6471.65,-274.59
+444.95,-6112.11,-259.63
+444.95,-5752.58,-244.61
+444.95,-5393.04,-229.56
+444.95,-5033.51,-214.45
+444.95,-4673.97,-199.3
+444.95,-4314.43,-184.1
+444.95,-3954.9,-168.85
+444.95,-3595.36,-153.55
+444.95,-3235.82,-138.2
+444.95,-2876.29,-122.81
+444.95,-2516.75,-107.37
+444.95,-2157.22,-91.88
+444.95,-1797.68,-76.35
+444.95,-1438.14,-60.76
+444.95,-1078.61,-45.13
+444.95,-719.07,-29.44
+444.95,-359.54,-13.67
+444.95,-71.91,0
+444.95,67.66,5.21
+444.95,338.29,19.04
+444.95,676.58,36.13
+444.95,1014.86,53.24
+444.95,1353.15,70.39
+444.95,1691.44,87.59
+444.95,2029.73,104.83
+444.95,2368.01,122.13
+444.95,2706.3,139.47
+444.95,3044.59,156.87
+444.95,3382.88,174.31
+444.95,3721.16,191.8
+444.95,4059.45,209.33
+444.95,4397.74,226.92
+444.95,4736.03,244.56
+444.95,5074.31,262.24
+444.95,5412.6,279.97
+444.95,5750.89,297.76
+444.95,6089.18,315.58
+444.95,6427.46,333.46
+444.95,6765.75,351.39
+533.94,-7190.72,-364.37
+533.94,-6831.19,-346.59
+533.94,-6471.65,-328.74
+533.94,-6112.11,-310.84
+533.94,-5752.58,-292.87
+533.94,-5393.04,-274.84
+533.94,-5033.51,-256.75
+533.94,-4673.97,-238.6
+533.94,-4314.43,-220.39
+533.94,-3954.9,-202.12
+533.94,-3595.36,-183.79
+533.94,-3235.82,-165.4
+533.94,-2876.29,-146.94
+533.94,-2516.75,-128.43
+533.94,-2157.22,-109.85
+533.94,-1797.68,-91.21
+533.94,-1438.14,-72.51
+533.94,-1078.61,-53.75
+533.94,-719.07,-34.91
+533.94,-359.54,-15.95
+533.94,-71.91,0
+533.94,67.66,6.58
+533.94,338.29,23.25
+533.94,676.58,43.76
+533.94,1014.86,64.29
+533.94,1353.15,84.87
+533.94,1691.44,105.52
+533.94,2029.73,126.22
+533.94,2368.01,146.98
+533.94,2706.3,167.81
+533.94,3044.59,188.7
+533.94,3382.88,209.65
+533.94,3721.16,230.66
+533.94,4059.45,251.73
+533.94,4397.74,272.87
+533.94,4736.03,294.07
+533.94,5074.31,315.33
+533.94,5412.6,336.65
+533.94,5750.89,358.03
+533.94,6089.18,379.47
+533.94,6427.46,400.98
+533.94,6765.75,422.54
+622.93,-7190.72,-423.69
+622.93,-6831.19,-403.04
+622.93,-6471.65,-382.31
+622.93,-6112.11,-361.5
+622.93,-5752.58,-340.62
+622.93,-5393.04,-319.66
+622.93,-5033.51,-298.62
+622.93,-4673.97,-277.51
+622.93,-4314.43,-256.31
+622.93,-3954.9,-235.05
+622.93,-3595.36,-213.7
+622.93,-3235.82,-192.28
+622.93,-2876.29,-170.79
+622.93,-2516.75,-149.21
+622.93,-2157.22,-127.57
+622.93,-1797.68,-105.84
+622.93,-1438.14,-84.03
+622.93,-1078.61,-62.14
+622.93,-719.07,-40.16
+622.93,-359.54,-18
+622.93,-71.91,0
+622.93,67.66,8.09
+622.93,338.29,27.64
+622.93,676.58,51.59
+622.93,1014.86,75.56
+622.93,1353.15,99.58
+622.93,1691.44,123.68
+622.93,2029.73,147.86
+622.93,2368.01,172.1
+622.93,2706.3,196.43
+622.93,3044.59,220.84
+622.93,3382.88,245.32
+622.93,3721.16,269.88
+622.93,4059.45,294.51
+622.93,4397.74,319.23
+622.93,4736.03,344.02
+622.93,5074.31,368.89
+622.93,5412.6,393.83
+622.93,5750.89,418.85
+622.93,6089.18,443.95
+622.93,6427.46,469.13
+622.93,6765.75,494.38
+711.92,-7190.72,-482.31
+711.92,-6831.19,-458.83
+711.92,-6471.65,-435.26
+711.92,-6112.11,-411.59
+711.92,-5752.58,-387.83
+711.92,-5393.04,-363.97
+711.92,-5033.51,-340.02
+711.92,-4673.97,-315.97
+711.92,-4314.43,-291.83
+711.92,-3954.9,-267.59
+711.92,-3595.36,-243.26
+711.92,-3235.82,-218.84
+711.92,-2876.29,-194.32
+711.92,-2516.75,-169.71
+711.92,-2157.22,-145
+711.92,-1797.68,-120.19
+711.92,-1438.14,-95.29
+711.92,-1078.61,-70.28
+711.92,-719.07,-45.14
+711.92,-359.54,-19.76
+711.92,-71.91,0
+711.92,67.66,9.76
+711.92,338.29,32.25
+711.92,676.58,59.66
+711.92,1014.86,87.07
+711.92,1353.15,114.55
+711.92,1691.44,142.11
+711.92,2029.73,169.77
+711.92,2368.01,197.52
+711.92,2706.3,225.37
+711.92,3044.59,253.31
+711.92,3382.88,281.34
+711.92,3721.16,309.47
+711.92,4059.45,337.7
+711.92,4397.74,366.02
+711.92,4736.03,394.44
+711.92,5074.31,422.95
+711.92,5412.6,451.56
+711.92,5750.89,480.26
+711.92,6089.18,509.05
+711.92,6427.46,537.94
+711.92,6765.75,566.93
+800.91,-7190.72,-540.17
+800.91,-6831.19,-513.92
+800.91,-6471.65,-487.55
+800.91,-6112.11,-461.06
+800.91,-5752.58,-434.46
+800.91,-5393.04,-407.75
+800.91,-5033.51,-380.92
+800.91,-4673.97,-353.97
+800.91,-4314.43,-326.91
+800.91,-3954.9,-299.73
+800.91,-3595.36,-272.44
+800.91,-3235.82,-245.03
+800.91,-2876.29,-217.51
+800.91,-2516.75,-189.87
+800.91,-2157.22,-162.12
+800.91,-1797.68,-134.25
+800.91,-1438.14,-106.26
+800.91,-1078.61,-78.13
+800.91,-719.07,-49.84
+800.91,-359.54,-21.19
+800.91,-71.91,0
+800.91,67.66,11.59
+800.91,338.29,37.1
+800.91,676.58,67.98
+800.91,1014.86,98.85
+800.91,1353.15,129.79
+800.91,1691.44,160.84
+800.91,2029.73,191.99
+800.91,2368.01,223.26
+800.91,2706.3,254.64
+800.91,3044.59,286.14
+800.91,3382.88,317.75
+800.91,3721.16,349.48
+800.91,4059.45,381.32
+800.91,4397.74,413.28
+800.91,4736.03,445.36
+800.91,5074.31,477.54
+800.91,5412.6,509.85
+800.91,5750.89,542.27
+800.91,6089.18,574.8
+800.91,6427.46,607.45
+800.91,6765.75,640.21
+889.9,-7190.72,-597.26
+889.9,-6831.19,-568.28
+889.9,-6471.65,-539.16
+889.9,-6112.11,-509.9
+889.9,-5752.58,-480.5
+889.9,-5393.04,-450.96
+889.9,-5033.51,-421.28
+889.9,-4673.97,-391.47
+889.9,-4314.43,-361.52
+889.9,-3954.9,-331.43
+889.9,-3595.36,-301.2
+889.9,-3235.82,-270.83
+889.9,-2876.29,-240.33
+889.9,-2516.75,-209.68
+889.9,-2157.22,-178.9
+889.9,-1797.68,-147.97
+889.9,-1438.14,-116.9
+889.9,-1078.61,-85.66
+889.9,-719.07,-54.21
+889.9,-359.54,-22.23
+889.9,-71.91,0
+889.9,67.66,13.6
+889.9,338.29,42.2
+889.9,676.58,76.58
+889.9,1014.86,110.92
+889.9,1353.15,145.34
+889.9,1691.44,179.88
+889.9,2029.73,214.55
+889.9,2368.01,249.35
+889.9,2706.3,284.29
+889.9,3044.59,319.36
+889.9,3382.88,354.57
+889.9,3721.16,389.92
+889.9,4059.45,425.41
+889.9,4397.74,461.03
+889.9,4736.03,496.79
+889.9,5074.31,532.69
+889.9,5412.6,568.73
+889.9,5750.89,604.9
+889.9,6089.18,641.21
+889.9,6427.46,677.66
+889.9,6765.75,714.24
+978.89,-7190.72,-653.54
+978.89,-6831.19,-621.88
+978.89,-6471.65,-590.05
+978.89,-6112.11,-558.06
+978.89,-5752.58,-525.9
+978.89,-5393.04,-493.58
+978.89,-5033.51,-461.09
+978.89,-4673.97,-428.44
+978.89,-4314.43,-395.63
+978.89,-3954.9,-362.65
+978.89,-3595.36,-329.51
+978.89,-3235.82,-296.21
+978.89,-2876.29,-262.74
+978.89,-2516.75,-229.11
+978.89,-2157.22,-195.31
+978.89,-1797.68,-161.34
+978.89,-1438.14,-127.18
+978.89,-1078.61,-92.83
+978.89,-719.07,-58.2
+978.89,-359.54,-22.82
+978.89,-71.91,0
+978.89,67.66,15.79
+978.89,338.29,47.58
+978.89,676.58,85.48
+978.89,1014.86,123.3
+978.89,1353.15,161.22
+978.89,1691.44,199.26
+978.89,2029.73,237.46
+978.89,2368.01,275.81
+978.89,2706.3,314.32
+978.89,3044.59,352.99
+978.89,3382.88,391.83
+978.89,3721.16,430.82
+978.89,4059.45,469.98
+978.89,4397.74,509.3
+978.89,4736.03,548.78
+978.89,5074.31,588.42
+978.89,5412.6,628.23
+978.89,5750.89,668.2
+978.89,6089.18,708.32
+978.89,6427.46,748.61
+978.89,6765.75,789.06
+1067.88,-7190.72,-708.97
+1067.88,-6831.19,-674.68
+1067.88,-6471.65,-640.19
+1067.88,-6112.11,-605.51
+1067.88,-5752.58,-570.63
+1067.88,-5393.04,-535.56
+1067.88,-5033.51,-500.3
+1067.88,-4673.97,-464.85
+1067.88,-4314.43,-429.21
+1067.88,-3954.9,-393.37
+1067.88,-3595.36,-357.35
+1067.88,-3235.82,-321.13
+1067.88,-2876.29,-284.71
+1067.88,-2516.75,-248.11
+1067.88,-2157.22,-211.31
+1067.88,-1797.68,-174.3
+1067.88,-1438.14,-137.08
+1067.88,-1078.61,-99.61
+1067.88,-719.07,-61.78
+1067.88,-359.54,-22.86
+1067.88,-71.91,0
+1067.88,67.66,18.17
+1067.88,338.29,53.24
+1067.88,676.58,94.71
+1067.88,1014.86,136.03
+1067.88,1353.15,177.45
+1067.88,1691.44,219.01
+1067.88,2029.73,260.75
+1067.88,2368.01,302.67
+1067.88,2706.3,344.78
+1067.88,3044.59,387.07
+1067.88,3382.88,429.54
+1067.88,3721.16,472.21
+1067.88,4059.45,515.06
+1067.88,4397.74,558.11
+1067.88,4736.03,601.34
+1067.88,5074.31,644.76
+1067.88,5412.6,688.37
+1067.88,5750.89,732.17
+1067.88,6089.18,776.15
+1067.88,6427.46,820.33
+1067.88,6765.75,864.69
+1156.87,-7190.72,-763.53
+1156.87,-6831.19,-726.65
+1156.87,-6471.65,-689.55
+1156.87,-6112.11,-652.22
+1156.87,-5752.58,-614.67
+1156.87,-5393.04,-576.89
+1156.87,-5033.51,-538.89
+1156.87,-4673.97,-500.67
+1156.87,-4314.43,-462.22
+1156.87,-3954.9,-423.56
+1156.87,-3595.36,-384.67
+1156.87,-3235.82,-345.56
+1156.87,-2876.29,-306.22
+1156.87,-2516.75,-266.66
+1156.87,-2157.22,-226.86
+1156.87,-1797.68,-186.83
+1156.87,-1438.14,-146.54
+1156.87,-1078.61,-105.95
+1156.87,-719.07,-64.9
+1156.87,-359.54,-22.24
+1156.87,-71.91,0
+1156.87,67.66,20.74
+1156.87,338.29,59.2
+1156.87,676.58,104.27
+1156.87,1014.86,149.12
+1156.87,1353.15,194.05
+1156.87,1691.44,239.16
+1156.87,2029.73,284.45
+1156.87,2368.01,329.96
+1156.87,2706.3,375.67
+1156.87,3044.59,421.6
+1156.87,3382.88,467.75
+1156.87,3721.16,514.11
+1156.87,4059.45,560.69
+1156.87,4397.74,607.48
+1156.87,4736.03,654.5
+1156.87,5074.31,701.73
+1156.87,5412.6,749.18
+1156.87,5750.89,796.84
+1156.87,6089.18,844.73
+1156.87,6427.46,892.83
+1156.87,6765.75,941.14
+1245.86,-7190.72,-817.17
+1245.86,-6831.19,-777.76
+1245.86,-6471.65,-738.08
+1245.86,-6112.11,-698.15
+1245.86,-5752.58,-657.96
+1245.86,-5393.04,-617.52
+1245.86,-5033.51,-576.82
+1245.86,-4673.97,-535.86
+1245.86,-4314.43,-494.65
+1245.86,-3954.9,-453.18
+1245.86,-3595.36,-411.45
+1245.86,-3235.82,-369.47
+1245.86,-2876.29,-327.22
+1245.86,-2516.75,-284.72
+1245.86,-2157.22,-241.95
+1245.86,-1797.68,-198.9
+1245.86,-1438.14,-155.54
+1245.86,-1078.61,-111.82
+1245.86,-719.07,-67.51
+1245.86,-359.54,-20.8
+1245.86,-71.91,0
+1245.86,67.66,23.52
+1245.86,338.29,65.48
+1245.86,676.58,114.2
+1245.86,1014.86,162.59
+1245.86,1353.15,211.06
+1245.86,1691.44,259.72
+1245.86,2029.73,308.59
+1245.86,2368.01,357.69
+1245.86,2706.3,407.04
+1245.86,3044.59,456.62
+1245.86,3382.88,506.46
+1245.86,3721.16,556.54
+1245.86,4059.45,606.87
+1245.86,4397.74,657.45
+1245.86,4736.03,708.28
+1245.86,5074.31,759.35
+1245.86,5412.6,810.68
+1245.86,5750.89,862.25
+1245.86,6089.18,914.07
+1245.86,6427.46,966.13
+1245.86,6765.75,1018.44
+1334.85,-7190.72,-869.86
+1334.85,-6831.19,-827.96
+1334.85,-6471.65,-785.77
+1334.85,-6112.11,-743.28
+1334.85,-5752.58,-700.5
+1334.85,-5393.04,-657.42
+1334.85,-5033.51,-614.05
+1334.85,-4673.97,-570.39
+1334.85,-4314.43,-526.44
+1334.85,-3954.9,-482.19
+1334.85,-3595.36,-437.65
+1334.85,-3235.82,-392.82
+1334.85,-2876.29,-347.69
+1334.85,-2516.75,-302.26
+1334.85,-2157.22,-256.52
+1334.85,-1797.68,-210.46
+1334.85,-1438.14,-164.04
+1334.85,-1078.61,-117.17
+1334.85,-719.07,-69.54
+1334.85,-359.54,-18.27
+1334.85,-71.91,0
+1334.85,67.66,26.49
+1334.85,338.29,72.1
+1334.85,676.58,124.5
+1334.85,1014.86,176.46
+1334.85,1353.15,228.49
+1334.85,1691.44,280.71
+1334.85,2029.73,333.18
+1334.85,2368.01,385.9
+1334.85,2706.3,438.89
+1334.85,3044.59,492.16
+1334.85,3382.88,545.71
+1334.85,3721.16,599.54
+1334.85,4059.45,653.64
+1334.85,4397.74,708.03
+1334.85,4736.03,762.71
+1334.85,5074.31,817.66
+1334.85,5412.6,872.89
+1334.85,5750.89,928.41
+1334.85,6089.18,984.2
+1334.85,6427.46,1040.27
+1334.85,6765.75,1096.63
+1423.84,-7190.72,-921.56
+1423.84,-6831.19,-877.23
+1423.84,-6471.65,-832.56
+1423.84,-6112.11,-787.56
+1423.84,-5752.58,-742.22
+1423.84,-5393.04,-696.56
+1423.84,-5033.51,-650.56
+1423.84,-4673.97,-604.23
+1423.84,-4314.43,-557.56
+1423.84,-3954.9,-510.57
+1423.84,-3595.36,-463.24
+1423.84,-3235.82,-415.58
+1423.84,-2876.29,-367.58
+1423.84,-2516.75,-319.24
+1423.84,-2157.22,-270.55
+1423.84,-1797.68,-221.48
+1423.84,-1438.14,-171.99
+1423.84,-1078.61,-121.95
+1423.84,-719.07,-70.93
+1423.84,-359.54,-14.1
+1423.84,-71.91,0
+1423.84,67.66,29.68
+1423.84,338.29,79.05
+1423.84,676.58,135.2
+1423.84,1014.86,190.76
+1423.84,1353.15,246.36
+1423.84,1691.44,302.17
+1423.84,2029.73,358.24
+1423.84,2368.01,414.6
+1423.84,2706.3,471.27
+1423.84,3044.59,528.24
+1423.84,3382.88,585.52
+1423.84,3721.16,643.12
+1423.84,4059.45,701.03
+1423.84,4397.74,759.26
+1423.84,4736.03,817.8
+1423.84,5074.31,876.66
+1423.84,5412.6,935.84
+1423.84,5750.89,995.34
+1423.84,6089.18,1055.15
+1423.84,6427.46,1115.27
+1423.84,6765.75,1175.71
+1512.83,-7190.72,-972.24
+1512.83,-6831.19,-925.52
+1512.83,-6471.65,-878.43
+1512.83,-6112.11,-830.96
+1512.83,-5752.58,-783.11
+1512.83,-5393.04,-734.89
+1512.83,-5033.51,-686.3
+1512.83,-4673.97,-637.33
+1512.83,-4314.43,-587.99
+1512.83,-3954.9,-538.28
+1512.83,-3595.36,-488.19
+1512.83,-3235.82,-437.72
+1512.83,-2876.29,-386.87
+1512.83,-2516.75,-335.63
+1512.83,-2157.22,-283.99
+1512.83,-1797.68,-231.92
+1512.83,-1438.14,-179.35
+1512.83,-1078.61,-126.12
+1512.83,-719.07,-71.6
+1512.83,-359.54,-6.45
+1512.83,-71.91,0
+1512.83,67.66,33.08
+1512.83,338.29,86.37
+1512.83,676.58,146.31
+1512.83,1014.86,205.5
+1512.83,1353.15,264.7
+1512.83,1691.44,324.11
+1512.83,2029.73,383.81
+1512.83,2368.01,443.83
+1512.83,2706.3,504.18
+1512.83,3044.59,564.87
+1512.83,3382.88,625.91
+1512.83,3721.16,687.3
+1512.83,4059.45,749.05
+1512.83,4397.74,811.14
+1512.83,4736.03,873.59
+1512.83,5074.31,936.4
+1512.83,5412.6,999.56
+1512.83,5750.89,1063.07
+1512.83,6089.18,1126.93
+1512.83,6427.46,1191.14
+1512.83,6765.75,1255.71
+1601.82,-3235.82,-459.19
+1601.82,-2876.29,-405.51
+1601.82,-2516.75,-351.39
+1601.82,-2157.22,-296.81
+1601.82,-1797.68,-241.73
+1601.82,-1438.14,-186.07
+1601.82,-1078.61,-129.61
+1601.82,-719.07,-71.46
+1601.82,-359.54,0
+1601.82,-71.91,0
+1601.82,67.66,36.7
+1601.82,338.29,94.04
+1601.82,676.58,157.85
+1601.82,1014.86,220.7
+1601.82,1353.15,283.52
+1601.82,1691.44,346.56
+1601.82,2029.73,409.9
+1601.82,2368.01,473.59
+1601.82,2706.3,537.65
+1601.82,3044.59,602.09
+1601.82,3382.88,666.91
+1601.82,3721.16,732.12
+1690.81,-3954.9,-591.53
+1690.81,-3595.36,-535.99
+1690.81,-3235.82,-479.97
+1690.81,-2876.29,-423.47
+1690.81,-2516.75,-366.48
+1690.81,-2157.22,-308.96
+1690.81,-1797.68,-250.88
+1690.81,-1438.14,-192.11
+1690.81,-1078.61,-132.36
+1690.81,-719.07,-70.38
+1690.81,-359.54,0
+1690.81,-71.91,0
+1690.81,67.66,40.53
+1690.81,338.29,102.1
+1690.81,676.58,169.83
+1690.81,1014.86,236.37
+1690.81,1353.15,302.84
+1690.81,1691.44,369.53
+1690.81,2029.73,436.54
+1690.81,2368.01,503.92
+1690.81,2706.3,571.71
+1690.81,3044.59,639.91
+1690.81,3382.88,708.54
+1690.81,3721.16,777.59
+1690.81,4059.45,847.08
+1690.81,4397.74,917
+1690.81,4736.03,987.35
+1690.81,5074.31,1058.14
+1690.81,5412.6,1129.36
+1690.81,5750.89,1201.01
+1690.81,6089.18,1273.09
+1690.81,6427.46,1345.61
+1690.81,6765.75,1418.56
+1779.8,-2516.75,-380.85
+1779.8,-2157.22,-320.4
+1779.8,-1797.68,-259.31
+1779.8,-1438.14,-197.41
+1779.8,-1078.61,-134.3
+1779.8,-719.07,-68.24
+1779.8,-359.54,0
+1779.8,-71.91,0
+1779.8,67.66,44.59
+1779.8,338.29,110.53
+1779.8,676.58,182.27
+1779.8,1014.86,252.55
+1779.8,1353.15,322.69
+1779.8,1691.44,393.04
+1779.8,2029.73,463.74