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],