Proof ORBITAL_INTEGRATION #229 (#389)

* Proof FIXED_FIXED_CONSTANT_POOL (#207)

* wip(p1338/1354: DONE): Proof ORBITAL_INTEGRATION #229

* wip(p1354/1354: DONE): Proof ORBITAL_INTEGRATION #229
This commit is contained in:
Glenn Dwiyatcita 2019-05-11 20:46:35 +02:00 committed by James Harris
parent c1421dfa88
commit 80d781c20e

View File

@ -42,15 +42,15 @@
KEPPREP LXA,2 SETPD
PBODY
0
DLOAD* SQRT # SQRT(MU) (+18 OR +15) 0D PL 2D
DLOAD* SQRT # SQRT(MU) (+18 OR +15) 0D PL 2D
MUEARTH,2
PDVL UNIT # PL 8D
RCV
PDDL NORM # NORM R (+29 OR +27 - N1) 2D PL 4D
PDDL NORM # NORM R (+29 OR +27 - N1) 2D PL 4D
36D
X1
PDVL
DOT PDDL # F*SQRT(MU) (+7 OR +5) 4D PL 6D
DOT PDDL # F*SQRT(MU)(+7 OR +5) 4D PL 6D
VCV
TAU. # (+28)
DSU NORM
@ -59,25 +59,25 @@ KEPPREP LXA,2 SETPD
SR1
DDV PDDL
2D
DMP PUSH # FS (+6 +N1-N2) 6D PL 8D
DMP PUSH # FS(+6 +N1-N2) 6D PL 8D
4D
DSQ PDDL # (FS)SQ (+12 +2(N1-N2)) 8D PL 10D
DSQ PDDL # (FS)SQ(+12 +2(N1-N2)) 8D PL 10D
4D
DSQ PDDL* # SSQ/MU (-2 OR +2(N1-N2)) 10D PL 12D
DSQ PDDL* # SSQ/MU(-2OR +2(N1-N2)) 10D PL 12D
MUEARTH,2
SR3 SR4
PDVL VSQ # PREALIGN MU (+43 OR +37) 12D PL 14D
PDVL VSQ # PREALIGN MU (+43 OR +37) 12D PL 14D
VCV
DMP BDSU # PL 12D
DMP BDSU # PL 12D
36D
DDV DMP # PL 10D
2D # -(1/R-ALPHA) (+12 +3N1-2N2)
DDV DMP # PL 10D
2D # -(1/R-ALPHA)(+12 +3N1-2N2)
DMP SL*
DP2/3
0 -3,1 # 10L(1/R-ALPHA) (+13 +2(N1-N2))
XSU,1 DAD # 2(FS)SQ - ETCETERA PL 8D
0 -3,1 # 10L(1/R-ALPHA)(+13 +2(N1-N2))
XSU,1 DAD # 2(FS)SQ - ETCETERA PL 8D
S1 # X1 = N2-N1
SL* DSU # -FS+2(FS)SQ ETC (+6 +N1-N2) PL 6D
SL* DSU # -FS+2(FS)SQ ETC (+6 +N1-N2) PL 6D
8D,1
DMP DMP
0D
@ -85,8 +85,8 @@ KEPPREP LXA,2 SETPD
SL* SL*
# Page 1335
8D,1
0,1 # S(-FS(1-2FS)-1/6...) (+17 OR +16)
DAD PDDL # PL 6D
0,1 # S(-FS(1-2FS)-1/6...)(+17 OR +16)
DAD PDDL # PL 6D
XKEP
DMP SL* # S(+17 OR +16)
0D
@ -325,7 +325,6 @@ INT-ABRT EXIT
# THE OBLATE ROUTINE COMPUTES THE ACCELERATION DUE TO OBLATENESS. IT USES THE UNIT OF THE VEHICLE
# POSITION VECTOR FOUND IN ALPHAV AND THE DISTANCE TO THE CENTER IN ALPHAM. THIS IS ADDED TO THE SUM OF THE
# DISTURBING ACCELERATIONS IN FV AND THE PROPER DIFEQ STAGE IS CALLED VIA X1.
OBLATE LXA,2 DLOAD
PBODY
ALPHAM
@ -349,7 +348,7 @@ OBLATE LXA,2 DLOAD
VAD VXM
ZUNIT
MMATRIX
UNIT # PROBABLY UNNECESSARY.
UNIT # POSSIBLY UNNECESSARY
COMTERM STORE UZ
DLOAD DMPR
COSPHI/2
@ -378,7 +377,7 @@ COMTERM STORE UZ
DMP*
J4REQ/J3,2
DDV DAD # -3
ALPHAM # (((P5/256)B 2 /R+P4/32) /R+P3/8)ALPHAV
ALPHAM # (((P5/256)B 2 /R+P4/32) /R+P3/8)ALPHAV
4 # 4 3
DMPR* DDV
2J3RE/J2,2
@ -390,9 +389,9 @@ COMTERM STORE UZ
DMP* SR1
J4REQ/J3,2
DDV DAD
ALPHAM
ALPHAM # -3
DMPR* SR3
2J3RE/J2,2
2J3RE/J2,2 # 3 4
DDV DAD
ALPHAM
VXSC VSL1
@ -423,7 +422,7 @@ COMTERM STORE UZ
QUALITY1
QUALITY3 DSQ # J22 TERM X R**4 IN 2D. SCALED B61
# AS VECTOR.
PUSH DMP # STORE COSPHI**2 SCALED B2 IN 8D.
PUSH DMP # STORE COSPHI**2 SCALED B2 IN 8D
# Page 1343
5/8 # 5 SCALED B3
PDDL SR2 # PUT 5 COSPHI**2, D5, IN 8D. GET
@ -450,14 +449,14 @@ QUALITY3 DSQ # J22 TERM X R**4 IN 2D. SCALED B61
5/8 # 5 B3 ANSWER B5
SL1 DAD # FROM 12D FOR Z COMPONENT (SL1 GIVES 10
# INSTEAD OF 5 FOR COEFFICIENT)
PDDL NORM # BACK INTO 12D FOR Z COMPONENT
PDDL NORM # BACK INTO 12D FOR Z COMPONENT.
ALPHAM # SCALED B27 FOR MOON
X2
PUSH SLOAD # STORE IN 14D, DESTROYING URPV
# X COMPONENT
E32C31RM
DDV VXSC # IF X2 = 0, DIVISION GIVES B53, VXSC
# out of 8D B5 GIVES B58
# OUT OF 8D B5 GIVES B58
VSL* VAD # SHIFT MAKES B61, FOR ADDITION OF
# VECTOR IN 2D
0 -3,2
@ -471,7 +470,7 @@ QUALITY3 DSQ # J22 TERM X R**4 IN 2D. SCALED B61
# ALONE.
PBODY
RP-TO-R
VAD BOV # OVERFLOW INDICATOR RESET IN *RP-TO-R*
VAD BOV # OVERFLOW INDICATOR RESET IN "RP-TO-R"
FV
GOBAQUE
STORE FV
@ -617,6 +616,7 @@ ORIGCHNG STQ CALL
# THE RECTIFY SUBROUTINE IS CALLED BY THE INTEGRATION PROGRAM AND OCCASIONALLY BY THE MEASUREMENT INCORPORATION
# ROUTINES TO ESTABLISH A NEW CONIC.
RECTIFY LXA,2 VLOAD
PBODY
TDELTAV
@ -640,9 +640,8 @@ MINIRECT STORE VRECT
RVQ
# Page 1348
# THE THREE DIFEQ ROUTINES -- DIFEQ+0, DIFEQ+12, DIFEQ+24 -- ARE ENTERED TO PROCESS THE CONTRIBUTIONS AT THE
# THE THREE DIFEQ ROUTINES - DIFEQ+0, DIFEQ+12, DIFEQ+24 - ARE ENTEREDTO PROCESS THE CONTRIBUTIONS AT THE
# BEGINNING, MIDDLE, AND END OF THE TIMESTEP, RESPECTIVELY. THE UPDATING IS DONE BY THE NYSTROM METHOD.
DIFEQ+0 VLOAD VSR3
FV
STCALL PHIV
@ -782,7 +781,7 @@ DIFEQCOM DLOAD DAD # INCREMENT H AND DIFEQCNT.
FBR3
WMATEND CLEAR CLEAR
DIM0FLAG # DON'T INTEGRATE W THIS TIME
DIM0FLAG # DONT INTEGRATE W THIS TIME
ORBWFLAG # INVALIDATE W
CLEAR
RENDWFLG
@ -797,9 +796,8 @@ WMATEND CLEAR CLEAR
# Page 1352
# ORBITAL ROUTINE FOR EXTRAPOLATION OF THE W MATRIX. IT COMPUTES THE SECOND DERIVATIVE OF EACH COLUMN POSITION
# VECTOR OF THE MATRIX AND CALLS THE NYSTROM INTEGRATION ROUTINES TO SOLVE THE DIFFERENTIAL EQUATIONS. THE PROGRAM
# USES A TABLE OF VEHICLE POSITION VECTORS COMPUTED DURING THE INTEGRATION OF THE VEHICLE'S POSITION AND VELOCITY.
# VECTOR OF THE MATRIX AND CALLS THE NYSTROM INTEGRATION ROUTINES TO SOLVETHE DIFFERENTIAL EQUATIONS. THE PROGRAM
# USES A TABLE OF VEHICLE POSITION VECTORS COMPUTED DURING THE INTEGRATION OF THE VEHICLES POSITION AND VELOCITY.
DOW.. LXA,2 DLOAD*
PBODY
MUEARTH,2
@ -944,6 +942,3 @@ URPV EQUALS 14D
COSPHI/2 EQUALS URPV +4
UZ EQUALS 20D
TVEC EQUALS 26D