programs/libtest/src/matrix-test.f

   1 C
   2 C
   3 C
   4       COMMON /DDD/Z
   5
   6       REAL A,B,C,Z
   7       DIMENSION A(4),B(4),C(4)
   8       INTEGER J
   9       DIMENSION J(2)
  10
  11 C      A(1)=10.0
  12 C      A(2)=1.22
  13 C      A(3)=23.0
  14 C      A(4)=1.5
  15
  16       Z=1.7070
  17       Z=1.557
  18       CALL M$INIT(A)
  19       CALL M$INIT(B)
  20       CALL M$INIT(C)
  21  10   FORMAT (9HMATRIX A:)
  22  11   FORMAT (9HMATRIX B:)
  23  12   FORMAT (9HMATRIX C:)
  24
  25
  26 C      A(1)=1.0
  27 C      A(2)=2.0
  28 C      A(3)=3.0
  29 C      A(4)=4.0
  30
  31       WRITE (1,10)
  32       CALL PMAT(A)
  33       CALL M$ROT(A,Z)
  34       CALL PMAT (A)
  35       PAUSE 4
  36 C      A(1)=1.0
  37 C      A(2)=2.0
  38 C      A(3)=3.0
  39 C      A(4)=4.0
  40
  41       B(1)=1.7
  42       B(4)=1.7
  43 C      B(2)=6.0
  44 C      B(3)=7.0
  45 C      B(4)=-8.0
  46
  47       CALL M$MUL(C,A,B)
  48
  49       WRITE (1,10)
  50       CALL PMAT(A)
  51       WRITE (1,11)
  52       CALL PMAT(B)
  53       WRITE (1,12)
  54       CALL PMAT(C)
  55
  56       J(1)=4
  57       J(2)=10
  58
  59  13   FORMAT (8HPOINT J:)
  60       WRITE (1,13)
  61       CALL PPOINT(J)
  62
  63       CALL M$APLI(C,J)
  64
  65  14   FORMAT (30HPOINT J AFTER TRANSFORMATION:   )
  66       WRITE (1,14)
  67       CALL PPOINT(J)
  68
  69
  70  22   FORMAT (20HPROGRAMMENDE            )
  71
  72
  73       WRITE (1,22)
  74
  75       PAUSE 2
  76  99   GO TO 99
  77       END
  78
  79 CCCCCCCCCCCCCCCCCCCCCCC
  80
  81       SUBROUTINE PMAT (MATRIX)
  82
  83       DIMENSION MATRIX(4)
  84       REAL MATRIX
  85       REAL M1,M2,M3,M4
  86       M1=MATRIX(0)
  87       M2=MATRIX(1)
  88       M3=MATRIX(2)
  89       M4=MATRIX(3)
  90
  91       WRITE (1,343) MATRIX
  92       WRITE (1,344)
  93
  94  343  FORMAT(2H[ ,F18.10,2X,F18.10,2H ])
  95  344  FORMAT(X)
  96
  97       RETURN
  98       END
  99
 100 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
 101
 102       SUBROUTINE PPOINT(POINT)
 103       DIMENSION POINT(2)
 104       INTEGER POINT
 105
 106       WRITE (1,200) POINT
 107  200  FORMAT(2H( ,I4,2H, I4,2H ))
 108       RETURN
 109       END
 110
 111
 112 CCCCCCCCCCCCCCCCCCC
 113 $0