From 80d781c20e6bc18b427818579c6e21dfb6ef00a9 Mon Sep 17 00:00:00 2001 From: Glenn Dwiyatcita Date: Sat, 11 May 2019 20:46:35 +0200 Subject: [PATCH] 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 --- Comanche055/ORBITAL_INTEGRATION.agc | 61 +++++++++++++---------------- 1 file changed, 28 insertions(+), 33 deletions(-) diff --git a/Comanche055/ORBITAL_INTEGRATION.agc b/Comanche055/ORBITAL_INTEGRATION.agc index 4086ccb..06151e5 100644 --- a/Comanche055/ORBITAL_INTEGRATION.agc +++ b/Comanche055/ORBITAL_INTEGRATION.agc @@ -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 - - -