Changeset 1728 in Sophya for trunk/Eval

Timestamp:

Oct 31, 2001, 6:55:43 PM (24 years ago)

Author:

ansari

Message:

extension prog test perfo f90 - Reza - 31/10/01

Location:

trunk/Eval/Speed

Files:

• makefile (modified) (1 diff)
• matrixf.f90 (modified) (4 diffs)

trunk/Eval/Speed/makefile

@@ -1 @@
-FC = f90 -O
-CC = cc -O -fast -tune host
-CXX = cxx -O -fast -tune host
+FC = f90 -O3
+CC = cc -O3
+CXX = cxx -O3
 LIBBLAS = -ldxml
 

trunk/Eval/Speed/matrixf.f90

@@ -4 +4 @@
 ! declaration d'interface de fonction 
       interface 
-        subroutine check_table(k, a, b, c)
+        subroutine check_table(k, a, b, c, ntot)
           integer, intent (in) :: k
+          integer, intent (inout) :: ntot
           double precision, intent (in), dimension(:,:) :: a,b,c
         end subroutine check_table
       end interface 
 
-      integer, parameter :: NL = 2, NC=2
+      integer, parameter :: NL = 4, NC=5
       double precision, dimension (NL, NC) :: m1,m2,m3,m4
 
-      integer :: i,j,k
-      
+      integer :: i,j,k,kk
+      integer :: ntot,nloop
+
       integer :: N = 5000000
 
       print *, '--- timing test : N=', N, ' NL=', NL, ' NC=',NC
+      ntot = 0
+      nloop = 0
+      do kk=1,1
       do k=1,N
         m1 = 0
@@ -29 +34 @@
         m3 = m1+m2;
         m4 = m1*m2;
-        call check_table(k, m2, m3, m4)
+        call check_table(k, m2, m3, m4, ntot)
 !        if (k .eq. 1) then
 !           print *, 'm1=', m1        
@@ -36 +41 @@
 !           print *, 'm4=', m4
 !        endif
+        nloop = nloop+1
       enddo
-
-      print *, '------- end of matrixf.f90 ------------'
+      enddo
+      
+      print *, ' nloop = ', nloop , ' ntot = ' , ntot 
+      print *, '------- end of matrixf.f90 -------- '
       end
 
 
-      subroutine check_table(k, a, b, c)
+      subroutine check_table(k, a, b, c, ntot)
         integer, intent (in) :: k
         double precision, intent (in), dimension(:,:) :: a,b,c
@@ -50 +58 @@
                ' size(a,2)= ', size(a,2)
       endif
+      ntot = ntot+1
       end subroutine check_table