aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--src/main.F9043
1 files changed, 35 insertions, 8 deletions
diff --git a/src/main.F90 b/src/main.F90
index 97e8e2d..d6b8002 100644
--- a/src/main.F90
+++ b/src/main.F90
@@ -7,12 +7,11 @@ PROGRAM mul
integer, parameter :: seed = 123456
real :: time
- integer :: dim = 10, multype = 1, real_kind, stat
- character(5) :: kind_arg
+ integer :: dim = 10, multype, real_kind, stat
+ character(5) :: kind_arg, impl_arg
- IF (command_argument_count() < 1) THEN
- write (error_unit, *) "Usage: mul KIND", char(10), &
- "where KIND is one of: 4, 8, 16"
+ IF (command_argument_count() < 2) THEN
+ call print_usage()
STOP
END IF
@@ -21,10 +20,27 @@ PROGRAM mul
read (kind_arg, *, iostat = stat) real_kind
IF (stat .ne. 0) THEN
- write (error_unit, *) "Couldn't parse kind number argument"
+ write (error_unit, '(A)') "Couldn't parse kind number argument"
+ call print_usage()
STOP
END IF
-
+
+ call get_command_argument(2, impl_arg)
+
+ IF (trim(impl_arg) .eq. "naiv") THEN
+ multype = 1
+ ELSE IF (trim(impl_arg) .eq. "bett") THEN
+ multype = 2
+ ELSE IF (trim(impl_arg) .eq. "dot") THEN
+ multype = 3
+ ELSE IF (trim(impl_arg) .eq. "mat") THEN
+ multype = 4
+ ELSE
+ write (error_unit, '(A)') "Unrecognized implementation argument"
+ call print_usage()
+ STOP
+ END IF
+
DO WHILE (dim < 2000)
SELECT CASE(real_kind)
@@ -35,7 +51,8 @@ PROGRAM mul
CASE (16,3)
time = measure_16(dim)
CASE default
- write (error_unit, *) "wrong kind for real: ", real_kind
+ write (error_unit, '(A,I6)') &
+ "wrong kind for real:", real_kind
STOP
END SELECT
@@ -48,6 +65,16 @@ PROGRAM mul
CONTAINS
+ SUBROUTINE print_usage()
+
+ write (*, '(A)') &
+ "Usage: mul KIND IMPLEMENTATION" // char(10) // &
+ "where KIND is one of: 4, 8, 16" // char(10) // &
+ " IMPLEMENTATION is one of: naiv, bett, dot, mat"
+
+ END SUBROUTINE print_usage
+
+
SUBROUTINE init_random_seed()
integer :: i, n