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.