Fast MKL on AMD

src/ezfio_files/NEED

@@ -1,2 +1,3 @@
 mpi
 zmq
+utils

src/ezfio_files/ezfio.irp.f

@@ -7,6 +7,8 @@ BEGIN_PROVIDER [ character*(1024), ezfio_filename ]
 
   PROVIDE mpi_initialized
 
+  integer :: i
+
   ! Get the QPACKAGE_INPUT environment variable
   call getenv('QPACKAGE_INPUT',ezfio_filename)
   if (ezfio_filename == '') then
@@ -44,11 +46,14 @@ BEGIN_PROVIDER [ character*(1024), ezfio_filename ]
 END_PROVIDER
 
 BEGIN_PROVIDER [ character*(1024), ezfio_work_dir ]
+ use c_functions
  implicit none
  BEGIN_DOC
  ! EZFIO/work/
  END_DOC
- call ezfio_set_work_empty(.False.)
+ logical :: b
+ b = mkl_serv_intel_cpu_true() /= 1
+ call ezfio_set_work_empty(b)
  ezfio_work_dir = trim(ezfio_filename)//'/work/'
 END_PROVIDER
 

src/utils/c_functions.f90

@@ -57,6 +57,12 @@ module c_functions
      end subroutine sscanf_sd_c
   end interface
 
+  interface
+    integer(kind=c_int) function mkl_serv_intel_cpu_true() bind(C)
+      use iso_c_binding
+    end function
+  end interface
+
 contains
 
   integer function atoi(a)
@@ -131,4 +137,3 @@ subroutine usleep(us)
   call usleep_c(u)
 end subroutine usleep
 
-