Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/rtie.ssp
There are 2 other files named rtie.ssp in the archive. Click here to see a list.
C                                                                       RTIE  10
C     ..................................................................RTIE  20
C                                                                       RTIE  30
C        SUBROUTINE RTIE                                                RTIE  40
C                                                                       RTIE  50
C        PURPOSE                                                        RTIE  60
C           ADJOIN TWO MATRICES WITH SAME COLUMN DIMENSION TO FORM ONE  RTIE  70
C           RESULTANT MATRIX (SEE METHOD)                               RTIE  80
C                                                                       RTIE  90
C        USAGE                                                          RTIE 100
C           CALL RTIE(A,B,R,N,M,MSA,MSB,L)                              RTIE 110
C                                                                       RTIE 120
C        DESCRIPTION OF PARAMETERS                                      RTIE 130
C           A - NAME OF FIRST INPUT MATRIX                              RTIE 140
C           B - NAME OF SECOND INPUT MATRIX                             RTIE 150
C           R - NAME OF OUTPUT MATRIX                                   RTIE 160
C           N - NUMBER OF ROWS IN A                                     RTIE 170
C           M - NUMBER OF COLUMNS IN A,B,R                              RTIE 180
C           MSA - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A         RTIE 190
C                  0 - GENERAL                                          RTIE 200
C                  1 - SYMMETRIC                                        RTIE 210
C                  2 - DIAGONAL                                         RTIE 220
C           MSB - SAME AS MSA EXCEPT FOR MATRIX B                       RTIE 230
C           L - NUMBER OF ROWS IN B                                     RTIE 240
C                                                                       RTIE 250
C        REMARKS                                                        RTIE 260
C           MATRIX R CANNOT BE IN THE SAME LOCATION AS MATRICES A OR B  RTIE 270
C           MATRIX R IS ALWAYS A GENERAL MATRIX                         RTIE 280
C           MATRIX A MUST HAVE THE SAME NUMBER OF COLUMNS AS MATRIX B   RTIE 290
C                                                                       RTIE 300
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  RTIE 310
C           LOC                                                         RTIE 320
C                                                                       RTIE 330
C        METHOD                                                         RTIE 340
C           MATRIX B IS ATTACHED TO THE BOTTOM OF MATRIX A .            RTIE 350
C           THE RESULTANT MATRIX R CONTAINS N+L ROWS AND M COLUMNS.     RTIE 360
C                                                                       RTIE 370
C     ..................................................................RTIE 380
C                                                                       RTIE 390
      SUBROUTINE RTIE(A,B,R,N,M,MSA,MSB,L)                              RTIE 400
      DIMENSION A(1),B(1),R(1)                                          RTIE 410
C                                                                       RTIE 420
      NN=N                                                              RTIE 430
      IR=0                                                              RTIE 440
      NX=NN                                                             RTIE 450
      MSX=MSA                                                           RTIE 460
      DO 9 J=1,M                                                        RTIE 470
      DO 8 II=1,2                                                       RTIE 480
      DO 7 I=1,NN                                                       RTIE 490
      IR=IR+1                                                           RTIE 500
      R(IR)=0.0                                                         RTIE 510
C                                                                       RTIE 520
C        LOCATE ELEMENT FOR ANY MATRIX STORAGE MODE                     RTIE 530
C                                                                       RTIE 540
      CALL LOC(I,J,IJ,NN,M,MSX)                                         RTIE 550
C                                                                       RTIE 560
C        TEST FOR ZERO ELEMENT IN DIAGONAL MATRIX                       RTIE 570
C                                                                       RTIE 580
      IF(IJ) 2,7,2                                                      RTIE 590
C                                                                       RTIE 600
C        MOVE ELEMENT TO MATRIX R                                       RTIE 610
C                                                                       RTIE 620
    2 GO TO(3,4),II                                                     RTIE 630
    3 R(IR)=A(IJ)                                                       RTIE 640
      GO TO 7                                                           RTIE 650
    4 R(IR)=B(IJ)                                                       RTIE 660
    7 CONTINUE                                                          RTIE 670
C                                                                       RTIE 680
C        REPEAT ABOVE FOR MATRIX B                                      RTIE 690
C                                                                       RTIE 700
      MSX=MSB                                                           RTIE 710
    8 NN=L                                                              RTIE 720
C                                                                       RTIE 730
C        RESET FOR NEXT COLUMN                                          RTIE 740
C                                                                       RTIE 750
      MSX=MSA                                                           RTIE 760
    9 NN=NX                                                             RTIE 770
      RETURN                                                            RTIE 780
      END                                                               RTIE 790