From 033025c0eaa62ecb2ede97824ad31e38dca0f240 Mon Sep 17 00:00:00 2001
From: Anthony Scemama <scemama@irsamc.ups-tlse.fr>
Date: Fri, 3 Jun 2016 13:54:23 +0200
Subject: [PATCH] Prefetch in sparse-full-mv

---
 src/mo.irp.f | 8 ++++++++
 1 file changed, 8 insertions(+)

diff --git a/src/mo.irp.f b/src/mo.irp.f
index d1de6c2..f84378b 100644
--- a/src/mo.irp.f
+++ b/src/mo.irp.f
@@ -690,6 +690,14 @@ subroutine sparse_full_mv(A,LDA,                                     &
   
   ! LDC and LDA have to be factors of simd_sp
   
+  IRP_IF NO_PREFETCH
+  IRP_ELSE
+    call MM_PREFETCH (A(j,indices(1)),3)
+    call MM_PREFETCH (A(j,indices(2)),3)
+    call MM_PREFETCH (A(j,indices(3)),3)
+    call MM_PREFETCH (A(j,indices(4)),3)
+  IRP_ENDIF
+
   !DIR$ VECTOR ALIGNED
   do j=1,LDC
     C1(j) = 0.