d1/dfe/test__drotmg_8f_source.html

00001       PROGRAM TEST_DROTMG
00002
00003       DOUBLE PRECISION DD1,DD2,DX1,DY1
00004       DOUBLE PRECISION COPY_DD1,COPY_DD2,COPY_DX1,COPY_DY1
00005       DOUBLE PRECISION DPARAM(5),COPY_DPARAM(5)
00006       DOUBLE PRECISION RAND
00007       integer i                 ! Counts random numbers
00008       integer*4 timeArray(3)    ! Holds the hour, minute, and second
00009       call itime(timeArray)     ! Get the current time
00010       i = rand ( timeArray(1)+timeArray(2)+timeArray(3) )
00011       DD1=RAND(0)
00012       DD2=RAND(0)
00013       DX1=RAND(0)
00014       DY1=RAND(0)
00015       COPY_DD1=DD1
00016       COPY_DD2=DD2
00017       COPY_DX1=DX1
00018       COPY_DY1=DY1
00019       DPARAM(1)=0
00020       DPARAM(2)=0
00021       DPARAM(3)=0
00022       DPARAM(4)=0
00023       DPARAM(5)=0
00024       COPY_DPARAM(1)=0
00025       COPY_DPARAM(2)=0
00026       COPY_DPARAM(3)=0
00027       COPY_DPARAM(4)=0
00028       COPY_DPARAM(5)=0
00029       WRITE(*,*) "DD1=",DD1,"DD2=",DD2
00030       WRITE(*,*) "DX1=",DX1,"DY1=",DY1
00031
00032       CALL MYDROTMG (DD1,DD2,DX1,DY1,DPARAM)
00033       CALL DROTMG (COPY_DD1,COPY_DD2,COPY_DX1,COPY_DY1,COPY_DPARAM)
00034
00035       IF (DABS(DD1-COPY_DD1) .GT. 0.1 ) THEN
00036         WRITE(*,*) "DD1 DIFFERS"
00037       ENDIF
00038       IF (DABS(DD2-COPY_DD2) .GT. 0.1 ) THEN
00039         WRITE(*,*) "DD2 DIFFERS"
00040       ENDIF
00041       IF (DABS(DX1-COPY_DX1) .GT. 0.1 ) THEN
00042         WRITE(*,*) "DX1 DIFFERS"
00043       ENDIF
00044       IF (DABS(DY1-COPY_DY1) .GT. 0.1 ) THEN
00045         WRITE(*,*) "DX1 DIFFERS"
00046       ENDIF
00047       DO i=1,5
00048          IF (DABS(DPARAM(i)-COPY_DPARAM(i)) .GT. 0.1 ) THEN
00049             WRITE(*,*) "DPARAM DIFFERS"
00050          ENDIF
00051       END DO
00052       END PROGRAM