Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/pmpy.ssp
There are 2 other files named pmpy.ssp in the archive. Click here to see a list.
C                                                                       PMPY  10
C     ..................................................................PMPY  20
C                                                                       PMPY  30
C        SUBROUTINE PMPY                                                PMPY  40
C                                                                       PMPY  50
C        PURPOSE                                                        PMPY  60
C           MULTIPLY TWO POLYNOMIALS                                    PMPY  70
C                                                                       PMPY  80
C        USAGE                                                          PMPY  90
C           CALL PMPY(Z,IDIMZ,X,IDIMX,Y,IDIMY)                          PMPY 100
C                                                                       PMPY 110
C        DESCRIPTION OF PARAMETERS                                      PMPY 120
C           Z     - VECTOR OF RESULTANT COEFFICIENTS, ORDERED FROM      PMPY 130
C                   SMALLEST TO LARGEST POWER                           PMPY 140
C           IDIMZ - DIMENSION OF Z (CALCULATED)                         PMPY 150
C           X     - VECTOR OF COEFFICIENTS FOR FIRST POLYNOMIAL, ORDEREDPMPY 160
C                   FROM SMALLEST TO LARGEST POWER                      PMPY 170
C           IDIMX - DIMENSION OF X (DEGREE IS IDIMX-1)                  PMPY 180
C           Y     - VECTOR OF COEFFICIENTS FOR SECOND POLYNOMIAL,       PMPY 190
C                   ORDERED FROM SMALLEST TO LARGEST POWER              PMPY 200
C           IDIMY - DIMENSION OF Y (DEGREE IS IDIMY-1)                  PMPY 210
C                                                                       PMPY 220
C        REMARKS                                                        PMPY 230
C           Z CANNOT BE IN THE SAME LOCATION AS X                       PMPY 240
C           Z CANNOT BE IN THE SAME LOCATION AS Y                       PMPY 250
C                                                                       PMPY 260
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  PMPY 270
C           NONE                                                        PMPY 280
C                                                                       PMPY 290
C        METHOD                                                         PMPY 300
C           DIMENSION OF Z IS CALCULATED AS IDIMX+IDIMY-1               PMPY 310
C           THE COEFFICIENTS OF Z ARE CALCULATED AS SUM OF PRODUCTS     PMPY 320
C           OF COEFFICIENTS OF X AND Y , WHOSE EXPONENTS ADD UP TO THE  PMPY 330
C           CORRESPONDING EXPONENT OF Z.                                PMPY 340
C                                                                       PMPY 350
C     ..................................................................PMPY 360
C                                                                       PMPY 370
      SUBROUTINE PMPY(Z,IDIMZ,X,IDIMX,Y,IDIMY)                          PMPY 380
      DIMENSION Z(1),X(1),Y(1)                                          PMPY 390
C                                                                       PMPY 400
      IF(IDIMX*IDIMY)10,10,20                                           PMPY 410
   10 IDIMZ=0                                                           PMPY 420
      GO TO 50                                                          PMPY 430
   20 IDIMZ=IDIMX+IDIMY-1                                               PMPY 440
      DO 30 I=1,IDIMZ                                                   PMPY 450
   30 Z(I)=0.                                                           PMPY 460
      DO 40 I=1,IDIMX                                                   PMPY 470
      DO 40 J=1,IDIMY                                                   PMPY 480
      K=I+J-1                                                           PMPY 490
   40 Z(K)=X(I)*Y(J)+Z(K)                                               PMPY 500
   50 RETURN                                                            PMPY 510
      END                                                               PMPY 520