Fix 2

Author: Critsium-xy

source/module_hsolver/diago_dav_subspace.cpp

@@ -25,6 +25,10 @@ Diago_DavSubspace<T, Device>::Diago_DavSubspace(const std::vector<Real>& precond
 {
     this->device = base_device::get_device_type<Device>(this->ctx);
 
+    this->one = &one_;
+    this->zero = &zero_;
+    this->neg_one = &neg_one_;
+
     assert(david_ndim_in > 1);
     assert(david_ndim_in * nband_in < nbasis_in * this->diag_comm.nproc);
 
@@ -193,7 +197,7 @@ int Diago_DavSubspace<T, Device>::diag_once(const HPsiFunc& hpsi_func,
                                  this->dim,
                                  this->vcc,
                                  this->nbase_x,
-                                 *this->zero,
+                                 this->zero,
                                  psi_in,
                                  psi_in_dmax);
 
@@ -279,7 +283,7 @@ void Diago_DavSubspace<T, Device>::cal_grad(const HPsiFunc& hpsi_func,
                          this->dim,
                          vcc,
                          this->nbase_x,
-                         *this->zero,
+                         this->zero,
                          psi_iter + (nbase) * this->dim,
                          this->dim);
 
@@ -413,7 +417,7 @@ void Diago_DavSubspace<T, Device>::cal_elem(const int& dim,
                          this->dim,
                          &hphi[nbase * this->dim],
                          this->dim,
-                         *this->zero,
+                         this->zero,
                          &hcc[nbase * this->nbase_x],
                          this->nbase_x);
 
@@ -433,7 +437,7 @@ void Diago_DavSubspace<T, Device>::cal_elem(const int& dim,
                          this->dim,
                          psi_iter + nbase * this->dim,
                          this->dim,
-                         *this->zero,
+                         this->zero,
                          &scc[nbase * this->nbase_x],
                          this->nbase_x);
 
@@ -640,7 +644,7 @@ void Diago_DavSubspace<T, Device>::refresh(const int& dim,
                          this->dim,
                          this->vcc,
                          this->nbase_x,
-                         *this->zero,
+                         this->zero,
                          psi_iter + nband * this->dim,
                          this->dim);