Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/loc.ibm
There are 2 other files named loc.ibm in the archive. Click here to see a list.
C                                                                       LOC   10
C     ..................................................................LOC   20
C                                                                       LOC   30
C        SUBROUTINE LOC                                                 LOC   40
C                                                                       LOC   50
C        PURPOSE                                                        LOC   60
C           COMPUTE A VECTOR SUBSCRIPT FOR AN ELEMENT IN A MATRIX OF    LOC   70
C           SPECIFIED STORAGE MODE                                      LOC   80
C                                                                       LOC   90
C        USAGE                                                          LOC  100
C           CALL LOC (I,J,IR,N,M,MS)                                    LOC  110
C                                                                       LOC  120
C        DESCRIPTION OF PARAMETERS                                      LOC  130
C           I   - ROW NUMBER OF ELEMENT                                 LOC  140
C           J   - COLUMN NUMBER  OF ELEMENT                             LOC  150
C           IR  - RESULTANT VECTOR SUBSCRIPT                            LOC  160
C           N   - NUMBER OF ROWS IN MATRIX                              LOC  170
C           M   - NUMBER OF COLUMNS IN MATRIX                           LOC  180
C           MS  - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX           LOC  190
C                  0 - GENERAL                                          LOC  200
C                  1 - SYMMETRIC                                        LOC  210
C                  2 - DIAGONAL                                         LOC  220
C                                                                       LOC  230
C        REMARKS                                                        LOC  240
C           NONE                                                        LOC  250
C                                                                       LOC  260
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  LOC  270
C           NONE                                                        LOC  280
C                                                                       LOC  290
C        METHOD                                                         LOC  300
C           MS=0   SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N*M ELEMENTS LOC  310
C                  IN STORAGE (GENERAL MATRIX)                          LOC  320
C           MS=1   SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N*(N+1)/2 IN LOC  330
C                  STORAGE (UPPER TRIANGLE OF SYMMETRIC MATRIX). IF     LOC  340
C                  ELEMENT IS IN LOWER TRIANGULAR PORTION, SUBSCRIPT IS LOC  350
C                  CORRESPONDING ELEMENT IN UPPER TRIANGLE.             LOC  360
C           MS=2   SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N ELEMENTS   LOC  370
C                  IN STORAGE (DIAGONAL ELEMENTS OF DIAGONAL MATRIX).   LOC  380
C                  IF ELEMENT IS NOT ON DIAGONAL (AND THEREFORE NOT IN  LOC  390
C                  STORAGE), IR IS SET TO ZERO.                         LOC  400
C                                                                       LOC  410
C     ..................................................................LOC  420
C                                                                       LOC  430
      SUBROUTINE LOC(I,J,IR,N,M,MS)                                     LOC  440
C                                                                       LOC  450
      IX=I                                                              LOC  460
      JX=J                                                              LOC  470
      IF(MS-1) 10,20,30                                                 LOC  480
   10 IRX=N*(JX-1)+IX                                                   LOC  490
      GO TO 36                                                          LOC  500
   20 IF(IX-JX) 22,24,24                                                LOC  510
   22 IRX=IX+(JX*JX-JX)/2                                               LOC  520
      GO TO 36                                                          LOC  530
   24 IRX=JX+(IX*IX-IX)/2                                               LOC  540
      GO TO 36                                                          LOC  550
   30 IRX=0                                                             LOC  560
      IF(IX-JX) 36,32,36                                                LOC  570
   32 IRX=IX                                                            LOC  580
   36 IR=IRX                                                            LOC  590
      RETURN                                                            LOC  600
      END                                                               LOC  610