Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-02 - 43,50145/load.ssp
There are 2 other files named load.ssp in the archive. Click here to see a list.
C                                                                       LOAD  10
C     ..................................................................LOAD  20
C                                                                       LOAD  30
C        SUBROUTINE LOAD                                                LOAD  40
C                                                                       LOAD  50
C        PURPOSE                                                        LOAD  60
C           COMPUTE A FACTOR MATRIX (LOADING) FROM EIGENVALUES AND      LOAD  70
C           ASSOCIATED EIGENVECTORS.  THIS SUBROUTINE NORMALLY OCCURS   LOAD  80
C           IN A SEQUENCE OF CALLS TO SUBROUTINES CORRE, EIGEN, TRACE,  LOAD  90
C           LOAD, AND VARMX IN THE PERFORMANCE OF A FACTOR ANALYSIS.    LOAD 100
C                                                                       LOAD 110
C        USAGE                                                          LOAD 120
C           CALL LOAD (M,K,R,V)                                         LOAD 130
C                                                                       LOAD 140
C        DESCRIPTION OF PARAMETERS                                      LOAD 150
C           M     - NUMBER OF VARIABLES.                                LOAD 160
C           K     - NUMBER OF FACTORS. K MUST BE GREATER THAN OR EQUAL  LOAD 170
C                   TO 1 AND LESS THAN OR EQUAL TO M.                   LOAD 171
C           R     - A MATRIX (SYMMETRIC AND STORED IN COMPRESSED FORM   LOAD 180
C                   WITH ONLY UPPER TRIANGLE BY COLUMN IN CORE) CON-    LOAD 190
C                   TAINING EIGENVALUES IN DIAGONAL.  EIGENVALUES ARE   LOAD 200
C                   ARRANGED IN DESCENDING ORDER, AND FIRST K           LOAD 210
C                   EIGENVALUES ARE USED BY THIS SUBROUTINE.  THE ORDER LOAD 220
C                   OF MATRIX R IS M BY M.  ONLY M*(M+1)/2 ELEMENTS ARE LOAD 230
C                   IN STORAGE.  (STORAGE MODE OF 1)                    LOAD 240
C           V     - WHEN THIS SUBROUTINE IS CALLED, MATRIX V (M X M)    LOAD 250
C                   CONTAINS EIGENVECTORS COLUMNWISE.  UPON RETURNING TOLOAD 260
C                   THE CALLING PROGRAM, MATRIX V CONTAINS A FACTOR     LOAD 270
C                   MATRIX (M X K).                                     LOAD 280
C                                                                       LOAD 290
C        REMARKS                                                        LOAD 300
C           NONE                                                        LOAD 310
C                                                                       LOAD 320
C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED                  LOAD 330
C           NONE                                                        LOAD 340
C                                                                       LOAD 350
C        METHOD                                                         LOAD 360
C           NORMALIZED EIGENVECTORS ARE CONVERTED TO THE FACTOR PATTERN LOAD 370
C           BY MULTIPLYING THE ELEMENTS OF EACH VECTOR BY THE SQUARE    LOAD 380
C           ROOT OF THE CORRESPONDING EIGENVALUE.                       LOAD 390
C                                                                       LOAD 400
C     ..................................................................LOAD 410
C                                                                       LOAD 420
      SUBROUTINE LOAD (M,K,R,V)                                         LOAD 430
      DIMENSION R(1),V(1)                                               LOAD 440
C                                                                       LOAD 450
C        ...............................................................LOAD 460
C                                                                       LOAD 470
C        IF A DOUBLE PRECISION VERSION OF THIS ROUTINE IS DESIRED, THE  LOAD 480
C        C IN COLUMN 1 SHOULD BE REMOVED FROM THE DOUBLE PRECISION      LOAD 490
C        STATEMENT WHICH FOLLOWS.                                       LOAD 500
C                                                                       LOAD 510
C     DOUBLE PRECISION R,V,SQ                                           LOAD 520
C                                                                       LOAD 530
C        THE C MUST ALSO BE REMOVED FROM DOUBLE PRECISION STATEMENTS    LOAD 540
C        APPEARING IN OTHER ROUTINES USED IN CONJUNCTION WITH THIS      LOAD 550
C        ROUTINE.                                                       LOAD 560
C                                                                       LOAD 570
C        THE DOUBLE PRECISION VERSION OF THIS SUBROUTINE MUST ALSO      LOAD 580
C        CONTAIN DOUBLE PRECISION FORTRAN FUNCTIONS.  SQRT IN STATEMENT LOAD 590
C        150 MUST BE CHANGED TO DSQRT.                                  LOAD 600
C                                                                       LOAD 610
C        ...............................................................LOAD 620
C                                                                       LOAD 630
      L=0                                                               LOAD 640
      JJ=0                                                              LOAD 650
      DO 160 J=1,K                                                      LOAD 660
      JJ=JJ+J                                                           LOAD 670
  150 SQ= SQRT(R(JJ))                                                   LOAD 680
      DO 160 I=1,M                                                      LOAD 690
      L=L+1                                                             LOAD 700
  160 V(L)=SQ*V(L)                                                      LOAD 710
      RETURN                                                            LOAD 720
      END                                                               LOAD 730