Uploaded image for project: 'GMAT'
  1. GMAT
  2. GMT-2461

Major errors in state when changing specific coordinate systems.

    Details

      Description

      Created an attachment (id=1456)
      Load this file to convert to all coordinate systems

      Load the BS-AllCoords.script attached and open Default SC.

      Set the Coordinate system to EarthMJ2000Eq and set the Cartesian state to:

      X: 5071.29822692573
      Y: 8411.70980173605
      Z: 267.680557173554
      VX: -5.44396385662813
      VY: 3.33268919535591
      VZ: 0.727825646501984

      Then change the Coordinate System to SunMJ2000EC. The state should be:

      X: 125912803.1220470000
      Y: 163450023.4650460000
      Z: 66337826.4249301000
      VX: -33.5749844009
      VY: 0.0008503563
      VZ: 1.4931513809

      but GMAT outputs:

      X: -105610452.5292780000
      Y: 102950234.9720830000
      Z: 199.1551456824
      VX: -26.7090415138
      VY: -18.1104991712
      VZ: 0.7267681121

      There are large errors in every state element.

      Then change the state to CS_ESL2. The state should be:

      X: -560121.5615799140
      Y: 284125.8138175170
      Z: 13140.7216646326
      VX: 5.4669115787
      VY: 0.7694483179
      VZ: 0.6355953786

      but GMAT outputs:

      X: 1065146.9798736300
      Y: -940367.7992908120
      Z: -407406.8320355240
      VX: -5.2305237214
      VY: 2.9656336125
      VZ: 2.0790570669

      Again, lagre errors in every state element.

      Finally change the state to CS_SSBary. The state should be:

      X: 147428584.0468950000
      Y: -9578.3404675689
      Z: 266.8937442606
      VX: 6.6827948425
      VY: 1.3706720758
      VZ: 0.7269537851

      but GMAT outputs:

      X: -106261378.9198180000
      Y: 93752651.0989441000
      Z: 40671209.0779822000
      VX: -26.6947192601
      VY: -16.9105817466
      VZ: -6.5398538120

        Gliffy Diagrams

          Attachments

            Issue Links

              Activity

              Hide
              shoffman Shawn Hoffman added a comment -

              Attachment BS-AllCoords.script has been added with description: Load this file to convert to all coordinate systems

              Show
              shoffman Shawn Hoffman added a comment - Attachment BS-AllCoords.script has been added with description: Load this file to convert to all coordinate systems
              shoffman Shawn Hoffman created issue -
              Hide
              shughes Steven Hughes added a comment -

              I think there is a truth issue here as well as in other state conversion "bugs".
              Below is data from GMAT and STK for the SunMJ2000Ec conversion that agrees much better than the numbers in the previous comments. This data still does not agree well enough and is off by about 50 meters. I am reassigning to myself so I can regenerate truth data for this case.

              STK using TT for ephem
              X: -105610452.52926444 km
              Y: 102950234.97209756 km
              Z: 199.1531497495481600 km
              Vx: -26.7090415137802460 km/sec
              Vy: -18.1104991711882910 km/sec
              Vz: 0.7267681124253601 km/sec

              GMAT using TT for ephem
              GMAT DefaultSC.X = -105610452.5104495;
              GMAT DefaultSC.Y = 102950234.9910701;
              GMAT DefaultSC.Z = 199.1551466137171;
              GMAT DefaultSC.VX = -26.70904151765019;
              GMAT DefaultSC.VY = -18.11049916746633;
              GMAT DefaultSC.VZ = 0.7267681120753649;

              Show
              shughes Steven Hughes added a comment - I think there is a truth issue here as well as in other state conversion "bugs". Below is data from GMAT and STK for the SunMJ2000Ec conversion that agrees much better than the numbers in the previous comments. This data still does not agree well enough and is off by about 50 meters. I am reassigning to myself so I can regenerate truth data for this case. STK using TT for ephem X: -105610452.52926444 km Y: 102950234.97209756 km Z: 199.1531497495481600 km Vx: -26.7090415137802460 km/sec Vy: -18.1104991711882910 km/sec Vz: 0.7267681124253601 km/sec GMAT using TT for ephem GMAT DefaultSC.X = -105610452.5104495; GMAT DefaultSC.Y = 102950234.9910701; GMAT DefaultSC.Z = 199.1551466137171; GMAT DefaultSC.VX = -26.70904151765019; GMAT DefaultSC.VY = -18.11049916746633; GMAT DefaultSC.VZ = 0.7267681120753649;
              Hide
              shughes Steven Hughes added a comment -

              By hand in Matlab I get the data below for the SunMJ2000Ec output:

              105610452.510744
              -102950234.990774
              -199.155146621168

              This agrees w/ GMAT to the meter level. I suspect I do not have GMAT and STK configured the same.

              Show
              shughes Steven Hughes added a comment - By hand in Matlab I get the data below for the SunMJ2000Ec output: 105610452.510744 -102950234.990774 -199.155146621168 This agrees w/ GMAT to the meter level. I suspect I do not have GMAT and STK configured the same.
              jjkparker Joel Parker made changes -
              Field Original Value New Value
              Fix Version/s Dec 15 [ 10342 ]
              jjkparker Joel Parker made changes -
              Link This issue is duplicated by GMT-78 [ GMT-78 ]
              Hide
              jjkparker Joel Parker added a comment -

              Treating this as the master bug again. The last comment on GMT-78 was:

              Reassigning to Steve for investigation of possible truth data errors (as discussed at yesterday's meeting).

              Show
              jjkparker Joel Parker added a comment - Treating this as the master bug again. The last comment on GMT-78 was: Reassigning to Steve for investigation of possible truth data errors (as discussed at yesterday's meeting).
              jjkparker Joel Parker made changes -
              Priority P1 [ 1 ] P2 [ 3 ]
              Fix Version/s 2012a M3 [ 10103 ]
              jjkparker Joel Parker made changes -
              Assignee Steven Hughes [ shughes ] Joel Parker [ jjkparker ]
              Hide
              jjkparker Joel Parker added a comment -

              Moving all 2012a M3 items to 2012a.

              Show
              jjkparker Joel Parker added a comment - Moving all 2012a M3 items to 2012a.
              jjkparker Joel Parker made changes -
              Fix Version/s 2012a M3 [ 10103 ]
              Fix Version/s 2012a [ 10000 ]
              Hide
              jjkparker Joel Parker added a comment -

              Assigning to Steve as Orbit State feature lead.

              Show
              jjkparker Joel Parker added a comment - Assigning to Steve as Orbit State feature lead.
              jjkparker Joel Parker made changes -
              Assignee Joel Parker [ jjkparker ] Steven Hughes [ shughes ]
              jjkparker Joel Parker made changes -
              Fix Version/s 2012a [ 10000 ]
              Fix Version/s 2013a [ 10201 ]
              Hide
              shughes Steven Hughes added a comment -

              I verified that GMAT's output is correct. We need to simply use the data above as the truth data in the Test Complete tests.

              Show
              shughes Steven Hughes added a comment - I verified that GMAT's output is correct. We need to simply use the data above as the truth data in the Test Complete tests.
              shughes Steven Hughes made changes -
              Assignee Steven Hughes [ shughes ] Shawn Hoffman [ shoffman ]
              Hide
              shoffman Shawn Hoffman added a comment -

              Which data above should be used for truth?

              Show
              shoffman Shawn Hoffman added a comment - Which data above should be used for truth?
              Hide
              shughes Steven Hughes added a comment -

              For each case in the first comment, when it says "but GMAT outputs" I verified that what GMAT is outputting is correct and it should be used as truth

              Show
              shughes Steven Hughes added a comment - For each case in the first comment, when it says "but GMAT outputs" I verified that what GMAT is outputting is correct and it should be used as truth
              Hide
              shoffman Shawn Hoffman added a comment -

              I have added the correct truth data

              Show
              shoffman Shawn Hoffman added a comment - I have added the correct truth data
              shoffman Shawn Hoffman made changes -
              Resolution Fixed [ 1 ]
              Status Open [ 1 ] Resolved [ 5 ]
              Hide
              shoffman Shawn Hoffman added a comment -

              Tests now pass

              Show
              shoffman Shawn Hoffman added a comment - Tests now pass
              shoffman Shawn Hoffman made changes -
              Status Resolved [ 5 ] Verified [ 10002 ]
              shoffman Shawn Hoffman made changes -
              Status Verified [ 10002 ] Closed [ 6 ]

                People

                • Votes:
                  0 Vote for this issue
                  Watchers:
                  3 Start watching this issue

                  Dates

                  • Created:
                    Updated:
                    Resolved: