From 0e767b06e541debcbb87fbddbf8175b315c3e72a Mon Sep 17 00:00:00 2001 From: Wojtek Kosior Date: Thu, 25 Apr 2019 19:29:12 +0200 Subject: adjust blocks size for different real types --- src/blockmath.F90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/blockmath.F90 b/src/blockmath.F90 index 247d0b0..2a05b94 100644 --- a/src/blockmath.F90 +++ b/src/blockmath.F90 @@ -28,8 +28,6 @@ MODULE blockmat IMPLICIT none PRIVATE - - integer, public :: blocksize = 25 PUBLIC :: blockmull PRIVATE :: blockmull_4, blockmull_8, blockmull_16 @@ -45,6 +43,7 @@ CONTAINS real(kind=4), intent(in), dimension(1:,1:) :: A, B real(kind=4), dimension(size(A, 1), size(B, 2)) :: C integer :: i, ii, ii_max, iiend, j, k, kk, kk_max, kkend + integer, parameter :: blocksize = 400 C = 0 @@ -81,6 +80,7 @@ CONTAINS real(kind=8), intent(in), dimension(1:,1:) :: A, B real(kind=8), dimension(size(A, 1), size(B, 2)) :: C integer :: i, ii, ii_max, iiend, j, k, kk, kk_max, kkend + integer, parameter :: blocksize = 44 C = 0 @@ -117,6 +117,7 @@ CONTAINS real(kind=16), intent(in), dimension(1:,1:) :: A, B real(kind=16), dimension(size(A, 1), size(B, 2)) :: C integer :: i, ii, ii_max, iiend, j, k, kk, kk_max, kkend + integer, parameter :: blocksize = 28 C = 0 -- cgit v1.2.3