Examples with TimeScalarFunction org.orekit.time.TimeScalarFunction used on opensource projects

Search in sources :

Example 1 with TimeScalarFunction

use of org.orekit.time.TimeScalarFunction in project Orekit by CS-SI.

the class IERSConventionsTest method testGMST82Sofa.

@Test
public void testGMST82Sofa() throws OrekitException {
    // the reference value has been computed using the March 2012 version of the SOFA library
    // http://www.iausofa.org/2012_0301_C.html, with the following code
    // 
    // double utc1, utc2, ut11, ut12, gmst;
    // 
    // // 2004-02-14:00:00:00Z, MJD = 53049, UT1-UTC = -0.4093509
    // utc1  = DJM0 + 53049.0;
    // utc2  = 0.0;
    // iauUtcut1(utc1, utc2, -0.4093509, &ut11, &ut12);
    // gmst = iauGmst82(ut11, ut12);
    // printf("iaugmst82(%.20g, %.20g)\n  --> %.20g\n", ut11, ut12, gmst);
    // 
    // // 2004-02-29:00:00:00Z, MJD = 53064, UT1-UTC = -0.4175723
    // utc1 = DJM0 + 53064.0;
    // utc2 = 0.0;
    // iauUtcut1(utc1, utc2, -0.4175723, &ut11, &ut12);
    // gmst = iauGmst82(ut11, ut12);
    // printf("iaugmst82(%.20g, %.20g)\n  --> %.20g\n", ut11, ut12, gmst);
    // 
    // the output of this test reads:
    // iaugmst82(2453049.5, -4.7378576388888813016e-06)
    // --> 2.5021977627453466653
    // iaugmst82(2453064.5, -4.8330127314815448519e-06)
    // --> 2.7602390405411441066
    final TimeScale ut1 = TimeScalesFactory.getUT1(IERSConventions.IERS_1996, true);
    final TimeScalarFunction gmst82 = IERSConventions.IERS_1996.getGMSTFunction(ut1);
    AbsoluteDate date = new AbsoluteDate(2004, 2, 14, TimeScalesFactory.getUTC());
    double gmst = MathUtils.normalizeAngle(gmst82.value(date), 0.0);
    Assert.assertEquals(2.5021977627453466653, gmst, 2.0e-13);
    date = new AbsoluteDate(2004, 2, 29, TimeScalesFactory.getUTC());
    gmst = MathUtils.normalizeAngle(gmst82.value(date), 0.0);
    Assert.assertEquals(2.7602390405411441066, gmst, 4.0e-13);
}
Also used : TimeScale(org.orekit.time.TimeScale) TimeScalarFunction(org.orekit.time.TimeScalarFunction) FieldAbsoluteDate(org.orekit.time.FieldAbsoluteDate) AbsoluteDate(org.orekit.time.AbsoluteDate) Test(org.junit.Test)

Example 2 with TimeScalarFunction

use of org.orekit.time.TimeScalarFunction in project Orekit by CS-SI.

the class IERSConventionsTest method testERASofa.

@Test
public void testERASofa() throws OrekitException {
    // the reference value has been computed using the March 2012 version of the SOFA library
    // http://www.iausofa.org/2012_0301_C.html, with the following code
    // 
    // double tai1, tai2, ut11, ut12, era, taiutc;
    // taiutc = 32.0;
    // 
    // // 2004-02-14:00:00:00Z, MJD = 53049, UT1-UTC = -0.4093509
    // tai1 = DJM0 + 53049.0;
    // tai2 = taiutc / DAYSEC;
    // iauTaiut1(tai1, tai2, -0.4093509 - taiutc, &ut11, &ut12);
    // era = iauEra00(ut11, ut12);
    // printf("iauera00(%.20g, %.20g)\n  --> %.20g\n", ut11, ut12, era);
    // 
    // // 2004-02-29:00:00:00Z, MJD = 53064, UT1-UTC = -0.4175723
    // tai1 = DJM0 + 53064.0;
    // tai2 = taiutc / DAYSEC;
    // iauTaiut1(tai1, tai2, -0.4175723 - taiutc, &ut11, &ut12);
    // era = iauEra00(ut11, ut12);
    // printf("iauera00(%.20g, %.20g)\n  --> %.20g\n", ut11, ut12, era);
    // 
    // the output of this test reads:
    // iauera00(2453049.5, -4.7378576388888813016e-06)
    // --> 2.5012766511308228701
    // iauera00(2453064.5, -4.8330127314815448519e-06)
    // --> 2.7593087452455264952
    final TimeScale ut1 = TimeScalesFactory.getUT1(IERSConventions.IERS_2003, true);
    final TimeScalarFunction era00 = IERSConventions.IERS_2003.getEarthOrientationAngleFunction(ut1);
    AbsoluteDate date = new AbsoluteDate(2004, 2, 14, TimeScalesFactory.getUTC());
    double era = MathUtils.normalizeAngle(era00.value(date), 0.0);
    Assert.assertEquals(2.5012766511308228701, era, 1.0e-15);
    date = new AbsoluteDate(2004, 2, 29, TimeScalesFactory.getUTC());
    era = MathUtils.normalizeAngle(era00.value(date), 0.0);
    Assert.assertEquals(2.7593087452455264952, era, 1.0e-15);
}
Also used : TimeScale(org.orekit.time.TimeScale) TimeScalarFunction(org.orekit.time.TimeScalarFunction) FieldAbsoluteDate(org.orekit.time.FieldAbsoluteDate) AbsoluteDate(org.orekit.time.AbsoluteDate) Test(org.junit.Test)

Example 3 with TimeScalarFunction

use of org.orekit.time.TimeScalarFunction in project Orekit by CS-SI.

the class IERSConventionsTest method testGMST06Sofa.

@Test
public void testGMST06Sofa() throws OrekitException {
    // the reference value has been computed using the March 2012 version of the SOFA library
    // http://www.iausofa.org/2012_0301_C.html, with the following code
    // 
    // double utc1, utc2, tai1, tai2, tt1, tt2, ut11, ut12, gmst;
    // 
    // // 2004-02-14:00:00:00Z, MJD = 53049, UT1-UTC = -0.4093509
    // utc1  = DJM0 + 53049.0;
    // utc2  = 0.0;
    // iauUtctai(utc1, utc2, &tai1, &tai2);
    // iauTaitt(tai1, tai2, &tt1, &tt2);
    // iauUtcut1(utc1, utc2, -0.4093509, &ut11, &ut12);
    // gmst = iauGmst06(ut11, ut12, tt1, tt2);
    // printf("iaugmst06(%.20g, %.20g, %.20g, %.20g)\n  --> %.20g\n",
    // ut11, ut12, tt1, tt2, gmst);
    // 
    // // 2004-02-29:00:00:00Z, MJD = 53064, UT1-UTC = -0.4175723
    // utc1 = DJM0 + 53064.0;
    // utc2 = 0.0;
    // iauUtctai(utc1, utc2, &tai1, &tai2);
    // iauTaitt(tai1, tai2, &tt1, &tt2);
    // iauUtcut1(utc1, utc2, -0.4175723, &ut11, &ut12);
    // gmst = iauGmst06(ut11, ut12, tt1, tt2);
    // printf("iaugmst06(%.20g, %.20g, %.20g, %.20g)\n  --> %.20g\n",
    // ut11, ut12, tt1, tt2, gmst);
    // 
    // the output of this test reads:
    // iaugmst06(2453049.5, -4.7378576388888813016e-06, 2453049.5, 0.00074287037037037029902)
    // --> 2.5021977784096232078
    // iaugmst06(2453064.5, -4.8330127314815448519e-06, 2453064.5, 0.00074287037037037029902)
    // --> 2.7602390556555129741
    final TimeScale ut1 = TimeScalesFactory.getUT1(IERSConventions.IERS_2010, true);
    final TimeScalarFunction gmst06 = IERSConventions.IERS_2010.getGMSTFunction(ut1);
    AbsoluteDate date = new AbsoluteDate(2004, 2, 14, TimeScalesFactory.getUTC());
    double gmst = MathUtils.normalizeAngle(gmst06.value(date), 0.0);
    Assert.assertEquals(2.5021977784096232078, gmst, 1.0e-15);
    date = new AbsoluteDate(2004, 2, 29, TimeScalesFactory.getUTC());
    gmst = MathUtils.normalizeAngle(gmst06.value(date), 0.0);
    Assert.assertEquals(2.7602390556555129741, gmst, 3.0e-15);
}
Also used : TimeScale(org.orekit.time.TimeScale) TimeScalarFunction(org.orekit.time.TimeScalarFunction) FieldAbsoluteDate(org.orekit.time.FieldAbsoluteDate) AbsoluteDate(org.orekit.time.AbsoluteDate) Test(org.junit.Test)

Example 4 with TimeScalarFunction

use of org.orekit.time.TimeScalarFunction in project Orekit by CS-SI.

the class IERSConventionsTest method testGST00ASofa.

@Test
public void testGST00ASofa() throws OrekitException {
    // the reference value has been computed using the March 2012 version of the SOFA library
    // http://www.iausofa.org/2012_0301_C.html, with the following code
    // 
    // double utc1, utc2, tai1, tai2, tt1, tt2, ut11, ut12, gst;
    // 
    // // 2004-02-14:00:00:00Z, MJD = 53049, UT1-UTC = -0.4093509
    // utc1  = DJM0 + 53049.0;
    // utc2  = 0.0;
    // iauUtctai(utc1, utc2, &tai1, &tai2);
    // iauTaitt(tai1, tai2, &tt1, &tt2);
    // iauUtcut1(utc1, utc2, -0.4093509, &ut11, &ut12);
    // gst = iauGst00a(ut11, ut12, tt1, tt2);
    // printf("iaugst00a(%.20g, %.20g, %.20g, %.20g)\n  --> %.20g\n",
    // ut11, ut12, tt1, tt2, gst);
    // 
    // // 2004-02-29:00:00:00Z, MJD = 53064, UT1-UTC = -0.4175723
    // utc1 = DJM0 + 53064.0;
    // utc2 = 0.0;
    // iauUtctai(utc1, utc2, &tai1, &tai2);
    // iauTaitt(tai1, tai2, &tt1, &tt2);
    // iauUtcut1(utc1, utc2, -0.4175723, &ut11, &ut12);
    // gst = iauGst00a(ut11, ut12, tt1, tt2);
    // printf("iaugst00a(%.20g, %.20g, %.20g, %.20g)\n  --> %.20g\n",
    // ut11, ut12, tt1, tt2, gst);
    // 
    // the output of this test reads:
    // iaugst00a(2453049.5, -4.7378576388888813016e-06, 2453049.5, 0.00074287037037037029902)
    // --> 2.5021491024360624778
    // iaugst00a(2453064.5, -4.8330127314815448519e-06, 2453064.5, 0.00074287037037037029902)
    // --> 2.7601901615614221619
    Utils.setLoaders(IERSConventions.IERS_2003, Utils.buildEOPList(IERSConventions.IERS_2003, ITRFVersion.ITRF_2008, new double[][] { { 53047, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53048, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53049, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53050, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53051, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53052, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53053, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53054, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, // ...
    { 53059, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53060, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53061, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53062, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53063, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53064, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53065, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53066, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 } }));
    EOPHistory eopHistory = FramesFactory.getEOPHistory(IERSConventions.IERS_2003, true);
    final TimeScalarFunction gst00a = IERSConventions.IERS_2003.getGASTFunction(TimeScalesFactory.getUT1(eopHistory), eopHistory);
    AbsoluteDate date = new AbsoluteDate(2004, 2, 14, TimeScalesFactory.getUTC());
    double gst = MathUtils.normalizeAngle(gst00a.value(date), 0.0);
    Assert.assertEquals(2.5021491024360624778, gst, 2.0e-13);
    date = new AbsoluteDate(2004, 2, 29, TimeScalesFactory.getUTC());
    gst = MathUtils.normalizeAngle(gst00a.value(date), 0.0);
    Assert.assertEquals(2.7601901615614221619, gst, 3.0e-13);
}
Also used : TimeScalarFunction(org.orekit.time.TimeScalarFunction) EOPHistory(org.orekit.frames.EOPHistory) FieldAbsoluteDate(org.orekit.time.FieldAbsoluteDate) AbsoluteDate(org.orekit.time.AbsoluteDate) Test(org.junit.Test)

Example 5 with TimeScalarFunction

use of org.orekit.time.TimeScalarFunction in project Orekit by CS-SI.

the class IERSConventionsTest method testGST94Sofa.

@Test
public void testGST94Sofa() throws OrekitException {
    // the reference value has been computed using the March 2012 version of the SOFA library
    // http://www.iausofa.org/2012_0301_C.html, with the following code
    // 
    // double utc1, utc2, ut11, ut12, tai1, tai2, tt1, tt2, gmst82, eqeq94;
    // 
    // // 2004-02-14:00:00:00Z, MJD = 53049, UT1-UTC = -0.4093509
    // utc1  = DJM0 + 53049.0;
    // utc2  = 0.0;
    // iauUtcut1(utc1, utc2, -0.4093509, &ut11, &ut12);
    // iauUtctai(utc1, utc2, &tai1, &tai2);
    // iauTaitt(tai1, tai2, &tt1, &tt2);
    // gmst82 = iauGmst82(ut11, ut12);
    // eqeq94 = iauEqeq94(tt1, tt2);
    // printf("iauGmst82(%.20g, %.20g)\n  --> %.20g\n", ut11, ut12, gmst82);
    // printf("iauEqeq94(%.20g, %.20g)\n  --> %.20g\n", tt1, tt2, eqeq94);
    // printf(" gmst82 + eqeq94  --> %.20g\n", gmst82 + eqeq94);
    // 
    // // 2004-02-29:00:00:00Z, MJD = 53064, UT1-UTC = -0.4175723
    // utc1 = DJM0 + 53064.0;
    // utc2 = 0.0;
    // iauUtcut1(utc1, utc2, -0.4175723, &ut11, &ut12);
    // iauUtctai(utc1, utc2, &tai1, &tai2);
    // iauTaitt(tai1, tai2, &tt1, &tt2);
    // gmst82 = iauGmst82(ut11, ut12);
    // eqeq94 = iauEqeq94(tt1, tt2);
    // printf("iauGmst82(%.20g, %.20g)\n  --> %.20g\n", ut11, ut12, gmst82);
    // printf("iauEqeq94(%.20g, %.20g)\n  --> %.20g\n", tt1, tt2, eqeq94);
    // printf(" gmst82 + eqeq94  --> %.20g\n", gmst82 + eqeq94);
    // 
    // the output of this test reads:
    // iauGmst82(2453049.5, -4.7378576388888813016e-06)
    // --> 2.5021977627453466653
    // iauEqeq94(2453049.5, 0.00074287037037037029902)
    // --> -4.8671604682267536886e-05
    // gmst82 + eqeq94  --> 2.5021490911406645274
    // iauGmst82(2453064.5, -4.8330127314815448519e-06)
    // --> 2.7602390405411441066
    // iauEqeq94(2453064.5, 0.00074287037037037029902)
    // --> -4.8893054690771762302e-05
    // gmst82 + eqeq94  --> 2.7601901474864534158
    // As can be seen in the code above, we didn't call iauGst94, because as
    // stated in the function header comment in SOFA source files:
    // "accuracy has been compromised for the sake of
    // convenience in that UT is used instead of TDB (or TT) to compute
    // the equation of the equinoxes."
    // So we rather performed the date conversion and then called ourselves iauGmst82
    // with a date in UTC and eqe94 with a date in TT, restoring full accuracy in the
    // SOFA computation for this test.
    // The thresholds below have been set up large enough to have the test pass with
    // the default Orekit setting and the reference SOFA setting. There are implementation
    // differences between the two libraries, as the Delaunay parameters are not the
    // same (in Orekit, we are compliant with page 23 in chapter 5 of IERS conventions 1996.
    // If the content of the IERS-conventions/1996/nutation-arguments.txt file by Orekit is
    // changed to match SOFA setting instead of IERS conventions as follows:
    // 
    // # Mean Anomaly of the Moon
    // F1 ≡ l = 485866.733″ + 1717915922.633″t + 31.310″t² + 0.064″t³
    // 
    // # Mean Anomaly of the Sun
    // F2 ≡ l' = 1287099.804″ + 129596581.224″t − 0.577″t² − 0.012t³
    // 
    // # L − Ω
    // F3 ≡ F = 335778.877″ + 1739527263.137″t − 13.257″t² + 0.011″t³
    // 
    // # Mean Elongation of the Moon from the Sun
    // F4 ≡ D = 1072261.307″ + 1602961601.328″t − 6.891″t² + 0.019″t³
    // 
    // # Mean Longitude of the Ascending Node of the Moon
    // F5 ≡ Ω = 450160.280″ − 6962890.539″t + 7.455″t² + 0.008″t³
    // 
    // then the thresholds for the test can be reduced to 4.2e-13 for the first test,
    // and 4.7e-13 for the second test.
    // We decided to stick with IERS published reference values for the default Orekit setting.
    // The differences are nevertheless quite small (9e-11 radians is sub-millimeter level
    // in low Earth orbit).
    Utils.setLoaders(IERSConventions.IERS_1996, Utils.buildEOPList(IERSConventions.IERS_1996, ITRFVersion.ITRF_2008, new double[][] { { 53047, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53048, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53049, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53050, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53051, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53052, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53053, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53054, -0.4093509, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, // ...
    { 53059, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53060, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53061, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53062, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53063, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53064, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53065, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 }, { 53066, -0.4175723, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 } }));
    EOPHistory eopHistory = FramesFactory.getEOPHistory(IERSConventions.IERS_1996, true);
    final TimeScalarFunction gst94 = IERSConventions.IERS_1996.getGASTFunction(TimeScalesFactory.getUT1(eopHistory), eopHistory);
    AbsoluteDate date = new AbsoluteDate(2004, 2, 14, TimeScalesFactory.getUTC());
    double gst = MathUtils.normalizeAngle(gst94.value(date), 0.0);
    // 4.2e-13 with SOFA values
    Assert.assertEquals(2.5021490911406645274, gst, 4.4e-11);
    date = new AbsoluteDate(2004, 2, 29, TimeScalesFactory.getUTC());
    gst = MathUtils.normalizeAngle(gst94.value(date), 0.0);
    // 4.7e-13 with SOFA values
    Assert.assertEquals(2.7601901474864534158, gst, 9.0e-11);
}
Also used : TimeScalarFunction(org.orekit.time.TimeScalarFunction) EOPHistory(org.orekit.frames.EOPHistory) FieldAbsoluteDate(org.orekit.time.FieldAbsoluteDate) AbsoluteDate(org.orekit.time.AbsoluteDate) Test(org.junit.Test)

Aggregations

TimeScalarFunction (org.orekit.time.TimeScalarFunction)11 Test (org.junit.Test)10 AbsoluteDate (org.orekit.time.AbsoluteDate)9 FieldAbsoluteDate (org.orekit.time.FieldAbsoluteDate)9 TimeScale (org.orekit.time.TimeScale)4 EOPHistory (org.orekit.frames.EOPHistory)3 Method (java.lang.reflect.Method)2 DSFactory (org.hipparchus.analysis.differentiation.DSFactory)2 DerivativeStructure (org.hipparchus.analysis.differentiation.DerivativeStructure)2 UT1Scale (org.orekit.time.UT1Scale)2 Field (java.lang.reflect.Field)1 RealFieldElement (org.hipparchus.RealFieldElement)1 BodiesElements (org.orekit.data.BodiesElements)1 FundamentalNutationArguments (org.orekit.data.FundamentalNutationArguments)1 PoissonSeries (org.orekit.data.PoissonSeries)1 PoissonSeriesParser (org.orekit.data.PoissonSeriesParser)1 OrekitInternalError (org.orekit.errors.OrekitInternalError)1 TimeVectorFunction (org.orekit.time.TimeVectorFunction)1 IERSConventions (org.orekit.utils.IERSConventions)1
About Me
UseOf

Java Tutorials :

Simple Java RabbitMQ Example with Spring
Simple Springboot JPA Eclipeslink Example
Simple SpringBoot JPA Hibernate Example
Jackson JSON @JsonIgnore and @JsonIgnoreProperties Examples
Java Infinite recursion (StackOverflowError) Jackson solutions
Simple Java Jackson Serialize Objects to JSON and convert them back To Object
ElasticSearch 5 - Upload files using Java API in binary field Example
Java JAXB marshall unmarshall Examples from String and Stream
Java 8 forEach List and Map examples
ElasticSearch 5 Java API SearchRequestBuilder examples
Java ElasticSearch 5 examples with node index and mapping - running local
Convert String to int or Integer Java 8
Java 9 in action - JShell - Ubuntu
Java - removing invalid characters from a file name
Hello world!? or Hello Tutorial