* OBSERVATION EQUATION (PARAMETRIC) LEAST SQUARES ADJUSTMENT 0020 * 0030 MAIN -ARADJ 0040 ON DATA,ASTOP 0050 L1 HEADING 0060 * 0070 * READ AND PRINT THE INPUT MATRICIES: A(COEFFICIENTS OF THE UNKNOWNS), 0080 * E (VARIANCE MATRIX OF THE OBSERVATIONS), AND L (CONSTANT TERMS) 0090 * 0100 DATA (SCALE,1.E-4) 0110 READ (NP,NN,MN,N) 0120 READ (A,E,L,LL) 0130 DOUBLE LL 0140 TEXT 'A MATRIX',C 0150 WRITE A,FORMAT=C5 0160 NEWPAGE 0170 TEXT 'VARIANCE-COVARIANCE MATRIX OF THE OBSERVATIONS',C 0180 WRITE E,FORMAT=C5 0190 FEED 5 0200 COPY N,NUA 0210 CANCEL N 0220 WRITE (NN,MN,NP,NUA),FORMAT=6 0230 NEWPAGE 0240 TEXT 'MISCLOSE MATRIX',C 0250 WRITE L,FORMAT=A5 0260 * 0270 * COMPUTATION AND PRINTING OF THE NORMAL EQUATIONS 0280 * 0290 INV E,P 0300 CANCEL E 0310 TRANS A,AT 0320 MULT AT,P,ATP 0330 MULT ATP,A,N 0340 MULT N,SCALE,N 0350 NEWPAGE 0360 TEXT 'NORMAL EQUATIONS',C 0370 WRITE N,FORMAT=C5 0380 * 0390 * COMPUTATION AND PRINTING OF THE SOLUTION VECTOR (X) 0400 * 0410 MULT ATP,L,U 0420 MULT U,SCALE,U 0430 TEXT 'TIME BEFORE INVERSION IS ',T 0440 FOR (MN,LT,100),LAB4 0450 COPY MN,RD 0460 SUB RD,100,C1 0470 DIV 4,100,INC 0480 CALL INVERT,(C1,RD,INC,N,NI) 0490 GOTO LAB5 0500 LAB4 INV N,NI 0510 LAB5 TEXT 'TIME AFTER INVERSION ',T 0520 COPYC NI,NIN 0530 MULT NIN,U,X 0540 NEWPAGE 0550 TEXT 'SOLUTION VECTOR (X MATRIX)',C 0560 WRITE X,FORMAT=A10 0570 * 0580 * COMPUTATION AND PRINTING OF THE ADJUSTED CO-ORDINATES (E,N) 0590 * 0600 MULT 2,NUA,NNE 0610 EXSUBM X,(1,1),(NNE,1),XL 0620 DOUBLE XL 0630 NEWPAGE 0640 TEXT 'APPROXIAMATE CO-ORDINATES (BY SEQUENCE NUMBER)',C 0650 WRITE LL,FORMAT=A10 0660 ADD XL,LL,LL 0670 TEXT 'ADJUSTED CO-ORDINATES (BY SEQUENCE NUMBER)',C 0680 WRITE LL,FORMAT=A10 0690 * 0700 * COMPUTATION AND PRINTING OF THE RESIDUAL VECTOR (V) 0710 * 0720 MULT A,X,AX 0730 CANCEL A 0740 CANCEL AT 0750 CANCEL ATP 0760 ADD AX,L,V 0770 NEWPAGE 0780 TEXT 'RESIDUAL VECTOR (V MATRIX)',C 0790 WRITE V,FORMAT=A10 0800 * 0810 * TWO SOLUTIONS FOR V'PV AND THE DIFFERENCE BETWEEN THEM 0820 * 0830 TRANS V,VT 0840 MULT VT,P,VTP 0850 MULT VTP,V,VPV1 0860 TRANS X,XT 0870 TRANS L,LT 0880 DIV SCALE,U,U 0890 MULT XT,U,XTU 0900 MULT LT,P,LTP 0910 MULT LTP,L,LPL 0920 ADD XTU,LPL,VPV2 0930 SUB VPV2,VPV1,DIFF 0940 SUB NN,MN,DF 0950 NEWPAGE 0960 * 0970 * PRINTING OF THE TWO SOLUTIONS OF V'PV,DIFFERENCE BETWEEN THE TWO SOLUTION0980 * FOR THE MINIMUM(V'PV1-V'PV2),AND THE DEGREES OF FREEDOM 0990 * 1000 FEED 5 1010 TEXT 'V''PV..DIRECT SOLUTION FROM THE RESIDUAL VECTOR V',5 1020 WRITE VPV1,FORMAT=A10 1030 FEED 7 1040 TEXT 'V''PV FROM THE FORMULA: V''PV=X''U+L''PL',5 1050 WRITE VPV2,FORMAT=A10 1060 FEED 7 1070 TEXT 'DIFFERENCE IN SOLUTIONS FOR V''PV',5 1080 WRITE DIFF,FORMAT=A10 1090 FEED 5 1100 TEXT 'DEGREES OF FREEDOM=#DF ',5 1110 * 1120 * COMPUTATION AND PRINTING OF THE ESTIMATED VARIANCE FACTOR(EVF) 1130 * 1140 DIV DF,VPV2,EVF 1150 FEED 5 1160 TEXT 'ESTIMATED VARIANCE FACTOR (EVF)',5 1170 WRITE EVF,FORMAT=A10 1180 * 1190 * COMPUTATION AND PRINTING OF THE VAR.-COVAR. OF THE ADJ. PARAMS. 1200 * 1210 MULT SCALE,NI,NI 1220 MULT EVF,NI,Q 1230 CANCEL NI 1240 NEWPAGE 1250 TEXT 'VARIANCE-COVARIANCE OF THE ADJUSTED COORDINATES',C 1260 EXSUBM Q,(1,1),(NNE,NNE),QX 1270 CANCEL Q 1280 WRITE QX,FORMAT=C6 1290 * 1300 * COMPUTATION OF STANDARD ELLIPSES FOR ALL ADJUSTED STATIONS 1310 * 1320 DIV 2,NNE,NCA 1330 SUB NP,NCA,CA 1340 ADD -1,NNE,NNE 1350 COPY -1,I 1360 LAB1 ADD 2,I,I 1370 EXSUBM QX,(I,I),(2,2),QXQ 1380 DET QXQ,DQ 1390 EXSUBM QXQ,(1,1),(1,1),AQ 1400 EXSUBM QXQ,(2,2),(1,1),BQ 1410 EXSUBM QXQ,(1,2),(1,1),CQ 1420 ADD AQ,BQ,TRC 1430 MULT -4,DQ,DQ 1440 MULT TRC,TRC,TRCS 1450 ADD TRCS,DQ,SRT 1460 ESQRT SRT,SRT 1470 ADD TRC,SRT,TSP 1480 DIV 2,TSP,LAI 1490 SUB TRC,SRT,TSM 1500 DIV 2,TSM,LAII 1510 SUB AQ,LAI,XQX 1520 FOR (CQ,EQ,0),LAB2 1530 DIV CQ,XQX,TAN 1540 COPYC TAN,TAN 1550 EARCTAN TAN,FI 1560 MULT 57.296,FI,FI 1570 GOTO LAB3 1580 LAB2 ADD CQ,CQ,FI 1590 LAB3 ESQRT LAI,LAI 1600 ESQRT LAII,LAII 1610 ADD 1,CA,CA 1620 NEWPAGE 1630 TEXT 'STANDARD ELLIPSE FOR POINT NO. =#CA ',40 1640 FEED 5 1650 TEXT 'LENGTH OF SEMI MAJOR AXIS',4 1660 WRITE LAI 1670 FEED 3 1680 TEXT 'LENGTH OF SEMI MINOR AXIS',4 1690 WRITE LAII 1700 FEED 3 1710 TEXT 'ROTATION OF MAJOR AXIS WITH RESPECT TO',4 1720 TEXT 'THE NORTH DIRECTION ( IN DEGREES )',4 1730 WRITE FI 1740 FOR (I,LT,NNE),LAB1 1750 GOTO L1 1760 ASTOP STOP 1770 END 1780 * 1790 SUBPRO INVERT,(C1,RD,INC,NS,VI) 1800 ATTRIB (NS,VI),SYM=1 1810 EXSUBM NS,(1,1),(C1,C1),V11 1820 INV V11,VI 1830 ADD C1,INC,D1 1840 ADD 1,C1,C2 1850 COPY D1,II 1860 START LOOP END1,D1,II,RD,INC 1870 EXSUBM NS,(1,1),(D1,D1),V 1880 EXSUBM V,(1,C2),(C1,INC),V12 1890 EXSUBM V,(C2,1),(INC,C1),V21 1900 EXSUBM V,(C2,C2),(INC,INC),V22 1910 MULT VI,V12,VIV 1920 MULT V21,VI,VVI 1930 MULT V21,VIV,V2VIV1 1940 SUB V22,V2VIV1,UV 1950 INV UV,U22 1960 COPYC U22,NU22 1970 MULT NU22,VVI,U21 1980 MULT VIV,NU22,U12 1990 MULT VIV,U21,VUV 2000 SUB VI,VUV,U11 2010 NULLMAT UI,(D1,D1) 2020 INSUBM U11,UI,(1,1) 2030 INSUBM U12,UI,(1,C2) 2040 INSUBM U21,UI,(C2,1) 2050 INSUBM U22,UI,(C2,C2) 2060 RENAME (UI,VI) 2070 ADD C1,INC,C1 2080 ADD 1,C1,C2 2090 END1 LOOPEND 2100 RETURN 2110 END 2120