add talk slides and source code for day2

Author: akohlmey

Diff for: day2/compile-link-make.pdf

@@ -0,0 +1,38 @@
+FC       = gfortran 
+#FC       = ifort
+#FC        = pgf77
+
+CFLAGS	= $(ABI) $(ARCH)
+
+TARGETS	= prog1 prog2a prog2b prog2c prog2d prog3
+
+default: all
+
+prog1: prog1.f
+	$(FC) -o $@ $(ABI) $(ARCH) \
+	      -O3 prog1.f
+
+prog2a: prog2.f
+	$(FC) -o $@ $(ABI) $(ARCH) \
+	      -O3 prog2.f
+
+prog2b: prog2.f
+	$(FC) -o $@ $(ABI) $(ARCH) \
+	      -O0 prog2.f
+
+prog2c: prog2.f
+	$(FC) -o $@ $(ABI) $(ARCH) \
+	      -O3 -ffast-math prog2.f
+
+prog2d: prog2.f
+	$(FC) -o $@ $(ABI) $(ARCH) \
+	      -O3 -ffast-math -funsafe-math-optimizations prog2.f
+
+prog3: prog3.f
+	$(FC) -o $@ $(ABI) $(ARCH) \
+	      -O3 prog3.f
+all: $(TARGETS)
+
+clean:
+	/bin/rm -f $(TARGETS)
+	/bin/rm -f *.s *.w2f.f

Diff for: day2/floating-point-code/mathopt/Makefile

@@ -0,0 +1,30 @@
+IEEE Floating Point Arithmetic Optimizations
+
+files: 
+prog1.f
+prog2.f
+prog3.f
+Makefile
+
+Objectives:
+
+    * understand the impact of IEEE-compliant compilation on performance 
+
+Instructions:
+
+   1. Look at the first program prog1.f: we calculate a two-point average, requiring a division by 2.
+   2. Make prog1 and measure the runtime. ( /bin/time ./prog1 ) 
+
+   3. Now we take the second program prog2.f, calculating a three-point average, requiring a division by 3.
+   4. Make prog2a and measure the runtime.
+   5. How do the runtimes compare to prog1?
+      Does the extra addition justify the increase in CPU time?
+
+   6. Look at the prog3.f: where is the difference between prog2.f and prog3.f? 
+   7, Make it and measure runtime
+
+   8. Now make prog2b, prog2c, and prog2d and take note of the extra compiler options used. 
+   9. Measure the runtime of prog2b, prog2c, and prog2d. Compare with prog3 and prog1
+      what are your conclusions?
+
+

Diff for: day2/floating-point-code/mathopt/README

@@ -0,0 +1,41 @@
+      program main1
+c
+c
+      implicit none
+      integer mx
+      parameter( mx = 100000 )
+      real*4 a(mx)
+
+      integer i
+
+      ! Initialize the array
+      
+      do i = 1, mx
+        a(i) = float(i)
+      end do
+      
+      ! Average the array
+      do i = 1, 10000
+        call xaver( a, mx )
+      end do
+      write(*,'(E22.14,2X,E22.14,2X,E22.14)')
+     &    a(10000), 15000.0e0,
+     &    (a(10000)-15000.0e0)
+      end
+
+
+      subroutine xaver( a, mx )
+c
+c  Average the data in i direction.
+c
+      implicit none
+      integer mx
+      real*4 a(mx)
+      integer i
+      
+      do i = 1, mx-1
+         a(i) = ( a(i) + a(i+1) ) / 2.0e0
+      end do
+
+      end
+

Diff for: day2/floating-point-code/mathopt/prog1

@@ -0,0 +1,42 @@
+      program main2
+c
+c  
+c
+      implicit none
+      integer mx
+      parameter( mx = 100000 )
+      real*4 a(mx)
+
+      integer i
+
+      ! Initialize the array
+      
+      do i = 1, mx
+         a(i) = float(i)
+      end do
+      
+      ! Average the array
+      do i = 1, 10000
+         call xaver( a, mx )
+      end do
+      write(*,'(E22.14,2X,E22.14,2X,E22.14)')
+     &    a(10000), 20000.0e0,
+     &    (a(10000)-20000.0e0)
+      end
+
+
+      subroutine xaver( a, mx )
+c
+c  Average the data in i direction.
+c
+      implicit none
+      integer mx
+      real a(mx)
+      integer i
+      
+      do i = 1, mx-2
+         a(i) = ( a(i) + a(i+1) + a(i+2) ) / 3.0e0
+      end do
+
+      end
+

Diff for: day2/floating-point-code/mathopt/prog1.f

@@ -0,0 +1,45 @@
+      program main3
+c
+c  
+c
+      implicit none
+      integer mx
+      parameter( mx = 100000 )
+      real a(mx)
+      real one_over3
+
+      integer i
+      ! define 1/3 
+      one_over3=1.0/3.0
+      ! Initialize the array
+      
+      do i = 1, mx
+         a(i) = float(i)
+      end do
+      
+      ! Average the array
+      do i = 1, 10000
+      	 call xaver( a, mx ,one_over3)
+      end do
+      write(*,'(E22.14,2X,E22.14,2X,E22.14)')
+     &    a(10000), 20000.0e0,
+     &    (a(10000)-20000.0e0)
+      end
+
+
+      subroutine xaver( a, mx,one_over3 )
+c
+c  Average the data in i direction.
+c
+      implicit none
+      integer mx
+      real a(mx)
+      real one_over3
+      integer i
+      
+      do i = 1, mx-2
+         a(i) = ( a(i) + a(i+1) + a(i+2) ) * one_over3
+      end do
+
+      end
+

Diff for: day2/floating-point-code/mathopt/prog2.f

@@ -0,0 +1,34 @@
+CC       = gcc
+
+CFLAGS	= $(ABI) $(ARCH)
+
+TARGETS	= prog1 prog2 prog3 prog4 prog5
+SRC = paranoia.c args.h
+
+default: all
+
+prog1:  $(SRC)
+	$(CC) -o $@ $(ABI) $(ARCH) \
+	      -O0 -mfpmath=387 paranoia.c -lm
+
+prog2:  $(SRC)
+	$(CC) -o $@ $(ABI) $(ARCH) \
+	      -O0 -mfpmath=sse paranoia.c -lm
+
+prog3:  $(SRC)
+	$(CC) -o $@ $(ABI) $(ARCH) \
+	      -O0 -mpc80 paranoia.c -lm
+
+prog4:  $(SRC)
+	$(CC) -o $@ $(ABI) $(ARCH) \
+	      -O2 -mpc64 -ansi -pedantic paranoia.c -lm
+
+prog5:  $(SRC)
+	$(CC) -o $@ $(ABI) $(ARCH) \
+	      -O3 -ffast-math -funsafe-math-optimizations paranoia.c -lm
+
+all: $(TARGETS)
+
+clean:
+	/bin/rm -f $(TARGETS)
+	/bin/rm -f *.s *.w2f.f

Diff for: day2/floating-point-code/mathopt/prog2a

@@ -0,0 +1,15 @@
+#if (defined(__cplusplus) || defined(__STDC__) || defined(c_plusplus))
+#define ARGS(parenthesized_list) parenthesized_list
+#define STDC		1
+#define VOID_ARG	void
+#else
+#define ARGS(parenthesized_list) ()
+#define STDC		0
+#define VOID_ARG
+#define const
+#endif
+
+#if !defined(EXIT_SUCCESS)
+#define EXIT_SUCCESS	0
+#define EXIT_FAILURE	1
+#endif