Modify pointer usage

Author: Critsium-xy

source/module_hsolver/diag_const_nums.cpp

@@ -13,13 +13,13 @@ const_nums<double>::const_nums()
 {
     base_device::memory::resize_memory_op<double, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->zero, 1);
-    *this->zero = 0.0;
+    this->zero[0] = 0.0;
     base_device::memory::resize_memory_op<double, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->one, 1);
-    *this->one = 1.0;
+    this->one[0] = 1.0;
     base_device::memory::resize_memory_op<double, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->neg_one, 1);
-    *this->neg_one = -1.0;
+    this->neg_one[0] = -1.0;
 }
 
 // Specialize templates to support double types
@@ -28,13 +28,13 @@ const_nums<float>::const_nums()
 {
     base_device::memory::resize_memory_op<float, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->zero, 1);
-    *this->zero = 0.0;
+    this->zero[0] = 0.0;
     base_device::memory::resize_memory_op<float, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->one, 1);
-    *this->one = 1.0;
+    this->one[0] = 1.0;
     base_device::memory::resize_memory_op<float, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->neg_one, 1);
-    *this->neg_one = -1.0;
+    this->neg_one[0] = -1.0;
 }
 
 // Specialized templates to support std:: complex<double>types
@@ -43,13 +43,13 @@ const_nums<std::complex<double>>::const_nums()
 {
     base_device::memory::resize_memory_op<std::complex<double>, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->zero, 1);
-    *this->zero = std::complex<double>(0.0, 0.0);
+    this->zero[0] = std::complex<double>(0.0, 0.0);
     base_device::memory::resize_memory_op<std::complex<double>, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->one, 1);
-    *this->one = std::complex<double>(1.0, 0.0);
+    this->one[0] = std::complex<double>(1.0, 0.0);
     base_device::memory::resize_memory_op<std::complex<double>, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->neg_one, 1);
-    *this->neg_one = std::complex<double>(-1.0, 0.0);
+    this->neg_one[0] = std::complex<double>(-1.0, 0.0);
 }
 
 // Specialized templates to support std:: complex<float>types
@@ -58,11 +58,11 @@ const_nums<std::complex<float>>::const_nums()
 {
     base_device::memory::resize_memory_op<std::complex<float>, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->zero, 1);
-    *this->zero = std::complex<float>(0.0, 0.0);
+    this->zero[0] = std::complex<float>(0.0, 0.0);
     base_device::memory::resize_memory_op<std::complex<float>, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->one, 1);
-    *this->one = std::complex<float>(1.0, 0.0);
+    this->one[0] = std::complex<float>(1.0, 0.0);
     base_device::memory::resize_memory_op<std::complex<float>, base_device::DEVICE_CPU>()(
                         this->cpu_ctx, this->neg_one, 1);
-    *this->neg_one = std::complex<float>(-1.0, 0.0);
+    this->neg_one[0] = std::complex<float>(-1.0, 0.0);
 }

source/module_hsolver/diago_dav_subspace.cpp

@@ -25,9 +25,9 @@ Diago_DavSubspace<T, Device>::Diago_DavSubspace(const std::vector<Real>& precond
 {
     this->device = base_device::get_device_type<Device>(this->ctx);
 
-    this->one = this->cs.one;
-    this->zero = this->cs.zero;
-    this->neg_one = this->cs.neg_one;
+    this->one = this->cs.one[0];
+    this->zero = this->cs.zero[0];
+    this->neg_one = this->cs.neg_one[0];
 
     assert(david_ndim_in > 1);
     assert(david_ndim_in * nband_in < nbasis_in * this->diag_comm.nproc);
@@ -534,8 +534,8 @@ void Diago_DavSubspace<T, Device>::diag_zhegvx(const int& nbase,
         }
         else
         {
-            std::vector<std::vector<T>> h_diag(nbase, std::vector<T>(nbase, *cs.zero));
-            std::vector<std::vector<T>> s_diag(nbase, std::vector<T>(nbase, *cs.zero));
+            std::vector<std::vector<T>> h_diag(nbase, std::vector<T>(nbase, cs.zero[0]));
+            std::vector<std::vector<T>> s_diag(nbase, std::vector<T>(nbase, cs.zero[0]));
 
             for (size_t i = 0; i < nbase; i++)
             {
@@ -564,10 +564,10 @@ void Diago_DavSubspace<T, Device>::diag_zhegvx(const int& nbase,
 
                 for (size_t j = nbase; j < this->nbase_x; j++)
                 {
-                    hcc[i * this->nbase_x + j] = *cs.zero;
-                    hcc[j * this->nbase_x + i] = *cs.zero;
-                    scc[i * this->nbase_x + j] = *cs.zero;
-                    scc[j * this->nbase_x + i] = *cs.zero;
+                    hcc[i * this->nbase_x + j] = cs.zero[0];
+                    hcc[j * this->nbase_x + i] = cs.zero[0];
+                    scc[i * this->nbase_x + j] = cs.zero[0];
+                    scc[j * this->nbase_x + i] = cs.zero[0];
                 }
             }
         }