From 6caf3521a423b80eed5bc9066f242e6e0682e17c Mon Sep 17 00:00:00 2001
From: Anthony Scemama <scemama@irsamc.ups-tlse.fr>
Date: Wed, 7 Feb 2024 17:06:27 +0100
Subject: [PATCH] commented loops

---
 org/qmckl_ao.org | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/org/qmckl_ao.org b/org/qmckl_ao.org
index c396bd8..2e2409c 100644
--- a/org/qmckl_ao.org
+++ b/org/qmckl_ao.org
@@ -5799,12 +5799,14 @@ qmckl_compute_ao_value_hpc_gaussian (const qmckl_context context,
 #endif
     for (int64_t ipoint=0 ; ipoint < point_num ; ++ipoint) {
 
+/*
       double* const ptr = &ao_value[ipoint*ao_num];
       for (size_t k=0 ; k<ao_num ; ++k) {
         ptr[k] = 0.;
       }
+*/
 
-//      memset(&ao_value[ipoint*ao_num], 0, ao_num*sizeof(double));
+      memset(&ao_value[ipoint*ao_num], 0, ao_num*sizeof(double));
 
       const double e_coord[3]  __attribute__((aligned(64))) =
         { coord[ipoint],
@@ -6595,12 +6597,14 @@ qmckl_compute_ao_vgl_hpc_gaussian (
 #endif
     for (int64_t ipoint=0 ; ipoint < point_num ; ++ipoint) {
 
+/*
       double* const ptr = &ao_vgl[ipoint*5*ao_num];
       for (size_t k=0 ; k<5*ao_num ; ++k) {
         ptr[k] = 0.;
       }
+*/
 
-//      memset(&ao_vgl[ipoint*ao_num*5], 0, 5*ao_num*sizeof(double));
+      memset(&ao_vgl[ipoint*ao_num*5], 0, 5*ao_num*sizeof(double));
 
       const double e_coord[3]  __attribute__((aligned(64))) =
         { coord[ipoint],