diff options
| -rw-r--r-- | src/main.F90 | 80 |
1 files changed, 75 insertions, 5 deletions
diff --git a/src/main.F90 b/src/main.F90 index 43eeabb..be9dd58 100644 --- a/src/main.F90 +++ b/src/main.F90 @@ -6,17 +6,17 @@ PROGRAM mul integer, parameter :: seed = 123456 real :: time - integer :: dim = 10, multype = 1, realtype = 4 + integer :: dim = 10, multype = 1, realtype = 8 DO WHILE (dim < 2000) SELECT CASE(realtype) CASE (4) time = measure_4(dim) - ! CASE (8) - ! time = measure_8(dim) - ! CASE (16) - ! time = measure_16(dim) + CASE (8) + time = measure_8(dim) + CASE (16) + time = measure_16(dim) CASE default write (*, *) "wrong kind for real: ", realtype @@ -73,5 +73,75 @@ CONTAINS time = end - start END FUNCTION measure_4 + + real FUNCTION measure_8(dim) result(time) + integer, intent(in) :: dim + real(kind=8), dimension(:,:), allocatable :: mat1, mat2, res + real :: start, end + + call init_random_seed() + + allocate(mat1(dim,dim)) + allocate(mat2(dim,dim)) + allocate(res(dim,dim)) + + call random_number(mat1) + call random_number(mat2) + + call cpu_time(start) + + SELECT CASE(multype) + CASE (1) + res = naivmull(mat1, mat2) + CASE (2) + res = bettmull(mat1, mat2) + CASE (3) + res = dotmull(mat1, mat2) + CASE default + res = matmul(mat1, mat2) + + END SELECT + + call cpu_time(end) + + time = end - start + + END FUNCTION measure_8 + + + real FUNCTION measure_16(dim) result(time) + integer, intent(in) :: dim + real(kind=16), dimension(:,:), allocatable :: mat1, mat2, res + real :: start, end + + call init_random_seed() + + allocate(mat1(dim,dim)) + allocate(mat2(dim,dim)) + allocate(res(dim,dim)) + + call random_number(mat1) + call random_number(mat2) + + call cpu_time(start) + + SELECT CASE(multype) + CASE (1) + res = naivmull(mat1, mat2) + CASE (2) + res = bettmull(mat1, mat2) + CASE (3) + res = dotmull(mat1, mat2) + CASE default + res = matmul(mat1, mat2) + + END SELECT + + call cpu_time(end) + + time = end - start + + END FUNCTION measure_16 + END PROGRAM mul |