Changeset 2366 in Sophya for trunk/Eval/JET/fmtx.f

Timestamp:

Apr 18, 2003, 5:48:16 PM (22 years ago)

Author:

ansari

Message:

Ajout fonctions Sin/Cos/... ds jet.h + prog en f90 - Reza 18/4/2003

File:

• trunk/Eval/JET/fmtx.f (modified) (3 diffs)

trunk/Eval/JET/fmtx.f

@@ -47 +47 @@
       CALL FPrtTim(4)
 
+      DO K=1,NTIME
+        C1 = K*0.135
+        C2 = COS(6.*C1)+0.04
+        C3 = SIN(C1*8.5)
+        CALL OPEFMTX(K, MX1, MX2, MX5, C1, C2, C3)
+      ENDDO
+      PRINT *, 'FIN  MX5 = C1*MX1 + MX3 + COS(C2*MX2*MX1) ' 
+      CALL FPrtTim(5)
+
       PRINT *, ' ------ FIN programme ---------- ' 
       END
@@ -82 +91 @@
       END
 
+      SUBROUTINE OPEFMTX(K, MX1, MX2, MX5, C1, C2, C3)
+      INTEGER I,J,K
+      DOUBLE PRECISION MX1(1000,500), MX2(1000,500)
+      DOUBLE PRECISION MX5(1000,500)
+      DOUBLE PRECISION C1,C2,C3
+      DO I=1,1000
+         DO J=1,500
+            MX5(I,J) = C1*MX1(I,J)+MX2(I,J)+
+     +                 C3*COS(C2*MX2(I,J)*MX1(I,J))
+         ENDDO
+      ENDDO
+      IF (K.EQ.5) PRINT *, 'OPEFMTX/COS - 5 -> ', MX5(20,5)
+      IF (K.EQ.15) PRINT *, 'OPEFMTX/COS - 15 -> ', MX5(20,15)
+      RETURN
+      END
+
       SUBROUTINE OPADDMTX(K, MX1, MX2, MX3, MX5, C1, C2, C3)
       INTEGER I,J,K
@@ -104 +129 @@
       DOUBLE PRECISION C1,C2,C3
       NTOT = 1000*500
-      CALL DVCAL(NTOT, C1, MX1, 1, MX5, 1)
-      CALL DZAXPY(NTOT, C2, MX2, 1, MX5, 1, MX4, 1)
-      CALL DZAXPY(NTOT, C3, MX3, 1, MX4, 1, MX5, 1)
+C  Sur OSF :
+C      CALL DVCAL(NTOT, C1, MX1, 1, MX5, 1)
+C      CALL DZAXPY(NTOT, C2, MX2, 1, MX5, 1, MX4, 1)
+C      CALL DZAXPY(NTOT, C3, MX3, 1, MX4, 1, MX5, 1)
+C  Sinon 
+      CALL DCOPY(NTOT, MX1, 1, MX5, 1)
+      CALL DSCAL(NTOT, C1, MX5, 1)
+      CALL DAXPY(NTOT, C2, MX2, 1, MX5, 1)
+      CALL DAXPY(NTOT, C3, MX3, 1, MX3, 1)
       IF (K.EQ.5) PRINT *, 'MTXBLAS - 5 -> ', MX5(20,5),
      + ' NTot=' , NTOT