Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/mtra.ssp
There are 2 other files named mtra.ssp in the archive. Click here to see a list.
C                                                                       MTRA  10
C     ..................................................................MTRA  20
C                                                                       MTRA  30
C        SUBROUTINE MTRA                                                MTRA  40
C                                                                       MTRA  50
C        PURPOSE                                                        MTRA  60
C           TRANSPOSE A MATRIX                                          MTRA  70
C                                                                       MTRA  80
C        USAGE                                                          MTRA  90
C           CALL MTRA(A,R,N,M,MS)                                       MTRA 100
C                                                                       MTRA 110
C        DESCRIPTION OF PARAMETERS                                      MTRA 120
C           A - NAME OF MATRIX TO BE TRANSPOSED                         MTRA 130
C           R - NAME OF OUTPUT MATRIX                                   MTRA 140
C           N - NUMBER OF ROWS OF A AND COLUMNS OF R                    MTRA 150
C           M - NUMBER OF COLUMNS OF A AND ROWS OF R                    MTRA 160
C           MS  - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A (AND R) MTRA 170
C                  0 - GENERAL                                          MTRA 180
C                  1 - SYMMETRIC                                        MTRA 190
C                  2 - DIAGONAL                                         MTRA 200
C                                                                       MTRA 210
C        REMARKS                                                        MTRA 220
C           MATRIX R CANNOT BE IN THE SAME LOCATION AS MATRIX A         MTRA 230
C                                                                       MTRA 240
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  MTRA 250
C           MCPY                                                        MTRA 260
C                                                                       MTRA 270
C        METHOD                                                         MTRA 280
C           TRANSPOSE N BY M MATRIX A TO FORM M BY N MATRIX R BY MOVING MTRA 290
C           EACH ROW OF A INTO THE CORRESPONDING COLUMN OF R. IF MATRIX MTRA 300
C           A IS SYMMETRIC OR DIAGONAL, MATRIX R IS THE SAME AS A.      MTRA 310
C                                                                       MTRA 320
C     ..................................................................MTRA 330
C                                                                       MTRA 340
      SUBROUTINE MTRA(A,R,N,M,MS)                                       MTRA 350
      DIMENSION A(1),R(1)                                               MTRA 360
C                                                                       MTRA 370
C        IF MS IS 1 OR 2, COPY A                                        MTRA 380
C                                                                       MTRA 390
      IF(MS) 10,20,10                                                   MTRA 400
   10 CALL MCPY(A,R,N,N,MS)                                             MTRA 410
      RETURN                                                            MTRA 420
C                                                                       MTRA 430
C        TRANSPOSE GENERAL MATRIX                                       MTRA 440
C                                                                       MTRA 450
   20 IR=0                                                              MTRA 460
      DO 30 I=1,N                                                       MTRA 470
      IJ=I-N                                                            MTRA 480
      DO 30 J=1,M                                                       MTRA 490
      IJ=IJ+N                                                           MTRA 500
      IR=IR+1                                                           MTRA 510
   30 R(IR)=A(IJ)                                                       MTRA 520
      RETURN                                                            MTRA 530
      END                                                               MTRA 540