Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap2_198111 - decus/20-0026/pvsub.ssp
There are 2 other files named pvsub.ssp in the archive. Click here to see a list.
C                                                                       PVSU  10
C     ..................................................................PVSU  20
C                                                                       PVSU  30
C        SUBROUTINE PVSUB                                               PVSU  40
C                                                                       PVSU  50
C        PURPOSE                                                        PVSU  60
C           SUBSTITUTE VARIABLE OF A POLYNOMIAL BY ANOTHER POLYNOMIAL   PVSU  70
C                                                                       PVSU  80
C        USAGE                                                          PVSU  90
C           CALL PVSUB(Z,IDIMZ,X,IDIMX,Y,IDIMY,WORK1,WORK2)             PVSU 100
C                                                                       PVSU 110
C        DESCRIPTION OF PARAMETERS                                      PVSU 120
C           Z     - VECTOR OF COEFFICIENTS FOR RESULTANT POLYNOMIAL,    PVSU 130
C                   ORDERED FROM SMALLEST TO LARGEST POWER              PVSU 140
C           IDIMZ - DIMENSION OF Z                                      PVSU 150
C           X     - VECTOR OF COEFFICIENTS FOR ORIGINAL POLYNOMIAL,     PVSU 160
C                   ORDERED FROM SMALLEST TO LARGEST POWER              PVSU 170
C           IDIMX - DIMENSION OF X                                      PVSU 180
C           Y     - VECTOR OF COEFFICIENTS FOR POLYNOMIAL WHICH IS      PVSU 190
C                   SUBSTITUTED FOR VARIABLE, ORDERED FROM SMALLEST TO  PVSU 200
C                   LARGEST POWER                                       PVSU 210
C           IDIMY - DIMENSION OF Y                                      PVSU 220
C           WORK1 - WORKING STORAGE ARRAY (SAME DIMENSION AS Z)         PVSU 230
C           WORK2 - WORKING STORAGE ARRAY (SAME DIMENSION AS Z)         PVSU 240
C                                                                       PVSU 250
C        REMARKS                                                        PVSU 260
C           NONE                                                        PVSU 270
C                                                                       PVSU 280
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  PVSU 290
C           PMPY                                                        PVSU 300
C           PADDM                                                       PVSU 310
C           PCLA                                                        PVSU 320
C                                                                       PVSU 330
C        METHOD                                                         PVSU 340
C           VARIABLE OF POLYNOMIAL X IS SUBSTITUTED BY POLYNOMIAL Y     PVSU 350
C           TO FORM POLYNOMIAL Z. DIMENSION OF NEW POLYNOMIAL IS        PVSU 360
C           (IDIMX-1)*(IDIMY-1)+1. SUBROUTINE REQUIRES TWO WORK AREAS   PVSU 370
C                                                                       PVSU 380
C     ..................................................................PVSU 390
C                                                                       PVSU 400
      SUBROUTINE PVSUB(Z,IDIMZ,X,IDIMX,Y,IDIMY,WORK1,WORK2)             PVSU 410
      DIMENSION Z(1),X(1),Y(1),WORK1(1),WORK2(1)                        PVSU 420
C                                                                       PVSU 430
C     TEST OF DIMENSIONS                                                PVSU 440
C                                                                       PVSU 450
      IF (IDIMX-1) 1,3,3                                                PVSU 460
    1 IDIMZ=0                                                           PVSU 470
    2 RETURN                                                            PVSU 480
C                                                                       PVSU 490
    3 IDIMZ=1                                                           PVSU 500
      Z(1)=X(1)                                                         PVSU 510
      IF (IDIMY*IDIMX-IDIMY) 2,2,4                                      PVSU 520
    4 IW1=1                                                             PVSU 530
      WORK1(1)=1.                                                       PVSU 540
C                                                                       PVSU 550
      DO 5 I=2,IDIMX                                                    PVSU 560
      CALL PMPY(WORK2,IW2,Y,IDIMY,WORK1,IW1)                            PVSU 570
      CALL PCLA(WORK1,IW1,WORK2,IW2)                                    PVSU 580
      FACT=X(I)                                                         PVSU 590
      CALL PADDM(Z,IDIMR,Z,IDIMZ,FACT,WORK1,IW1)                        PVSU 600
      IDIMZ=IDIMR                                                       PVSU 605
    5 CONTINUE                                                          PVSU 610
      GO TO 2                                                           PVSU 620
      END                                                               PVSU 630