Optimize ompGemm_m multiply kernel

.github/workflows/makefile.yml

@@ -70,8 +70,6 @@ jobs:
             multiply_kernel: ompDojk
           - test_all_multiply_kernels: false
             multiply_kernel: ompGemm
-          - test_all_multiply_kernels: false
-            multiply_kernel: ompGemm_m
 
     steps:
     - uses: actions/checkout@v3

benchmarks/.gitignore

@@ -0,0 +1,8 @@
+__pycache__
+Conquest_o*
+Conquest_w*
+*.i00*
+*.bib
+*.dat
+*.log
+fort.*

src/Makefile

@@ -102,7 +102,7 @@ initial_read_module.o:initial_read_module.f90 datestamp.o
 # for GCC13, possibly only on Mac.  It doesn't need any other
 # compiler flags (libraries or communications) so should be OK like this
 pseudo_tm_info.o:pseudo_tm_info.f90
-	$(FC) -c $<
+	$(FC) $(XC_COMPFLAGS) -c $<
 
 #datestamp.f90: $(COMMENT)
 #	$(ECHOSTR) "module datestamp\n" > datestamp.f90

src/matrix_module.f90

@@ -121,16 +121,16 @@ module matrix_module
   type matrix_halo
      integer :: np_in_halo ! Partns in halo
      integer :: ni_in_halo ! Atoms in halo
-     integer,pointer :: nh_part(:)    ! No of halo atoms in halo part
-     integer,pointer :: j_beg(:)      ! accumulator of nh_part
-     integer,pointer :: lab_hcell(:)  ! sim cell part (CC) = halo part
-     integer,pointer :: lab_hcover(:) ! CS part (CC) = halo part
-     integer,pointer :: j_seq(:)      ! Part seq of halo atom
-     integer,pointer :: i_h2d(:)      ! halo->neigh trans for halo atom
-     integer,pointer :: i_halo(:)     ! halo seq of atom in CS
-     integer,pointer :: i_hbeg(:)     ! Where CS part starts in i_halo
-     integer,pointer :: ndimi(:)
-     integer,pointer :: ndimj(:)
+     integer,allocatable :: nh_part(:)    ! No of halo atoms in halo part
+     integer,allocatable :: j_beg(:)      ! accumulator of nh_part
+     integer,allocatable :: lab_hcell(:)  ! sim cell part (CC) = halo part
+     integer,allocatable :: lab_hcover(:) ! CS part (CC) = halo part
+     integer,allocatable :: j_seq(:)      ! Part seq of halo atom
+     integer,allocatable :: i_h2d(:)      ! halo->neigh trans for halo atom
+     integer,allocatable :: i_halo(:)     ! halo seq of atom in CS
+     integer,allocatable :: i_hbeg(:)     ! Where CS part starts in i_halo
+     integer,allocatable :: ndimi(:)
+     integer,allocatable :: ndimj(:)
      ! This type's values for maxima
      integer :: mx_part,mx_halo
   end type matrix_halo