delete useless matrix transformer

Author: maki49

source/module_lr/utils/lr_util.cpp

@@ -102,105 +102,6 @@ out[i] = in[i];
         pztranc_(&n, &n, &alpha, tmp.data(), &i1, &i1, pmat.desc, &beta, inout, &i1, &i1, pmat.desc);
     }
 #endif
-    container::Tensor mat2ten_double(ModuleBase::matrix& m)
-    {
-        container::Tensor t(DAT::DT_DOUBLE, DEV::CpuDevice, { m.nr, m.nc });
-        for (int i = 0;i < t.NumElements();++i) {            
-t.data<double>()[i] = m.c[i];
-        }
-        return t;
-    }
-    std::vector<container::Tensor> mat2ten_double(std::vector<ModuleBase::matrix>& m)
-    {
-        std::vector<container::Tensor> t;
-        for (int i = 0;i < m.size();++i) {
-            t.push_back(mat2ten_double(m[i]));
-        }
-        return t;
-    }
-    ModuleBase::matrix ten2mat_double(container::Tensor& t)
-    {
-        ModuleBase::matrix m(t.shape().dims()[0], t.shape().dims()[1]);
-        for (int i = 0;i < t.NumElements();++i) {
-            m.c[i] = t.data<double>()[i];
-        }
-        return m;
-    }
-    std::vector<ModuleBase::matrix> ten2mat_double(std::vector<container::Tensor>& t)
-    {
-        std::vector<ModuleBase::matrix> m;
-        for (int i = 0;i < t.size();++i) {
-            m.push_back(ten2mat_double(t[i]));
-        }
-        return m;
-    }
-    container::Tensor mat2ten_complex(ModuleBase::ComplexMatrix& m)
-    {
-        container::Tensor t(DAT::DT_COMPLEX_DOUBLE, DEV::CpuDevice, { m.nr, m.nc });
-        for (int i = 0;i < t.NumElements();++i) {
-            t.data<std::complex<double>>()[i] = m.c[i];
-        }
-        return t;
-    }
-    std::vector<container::Tensor> mat2ten_complex(std::vector<ModuleBase::ComplexMatrix>& m)
-    {
-        std::vector<container::Tensor> t;
-        for (int i = 0;i < m.size();++i) {
-            t.push_back(mat2ten_complex(m[i]));
-        }
-        return t;
-    }
-    ModuleBase::ComplexMatrix ten2mat_complex(container::Tensor& t)
-    {
-        ModuleBase::ComplexMatrix m(t.shape().dims()[0], t.shape().dims()[1]);
-        for (int i = 0;i < t.NumElements();++i) {
-            m.c[i] = t.data<std::complex<double>>()[i];
-        }
-        return m;
-    }
-    std::vector<ModuleBase::ComplexMatrix> ten2mat_complex(std::vector<container::Tensor>& t)
-    {
-        std::vector<ModuleBase::ComplexMatrix> m;
-        for (int i = 0;i < t.size();++i) {
-            m.push_back(ten2mat_complex(t[i]));
-        }
-        return m;
-    }
-
-    ModuleBase::matrix vec2mat(const std::vector<double>& v, const int nr, const int nc)
-    {
-        assert(v.size() == nr * nc);
-        ModuleBase::matrix m(nr, nc, false);
-        for (int i = 0;i < v.size();++i) {
-            m.c[i] = v[i];
-        }
-        return m;
-    }
-    ModuleBase::ComplexMatrix vec2mat(const std::vector<std::complex<double>>& v, const int nr, const int nc)
-    {
-        assert(v.size() == nr * nc);
-        ModuleBase::ComplexMatrix m(nr, nc, false);
-        for (int i = 0;i < v.size();++i) {
-            m.c[i] = v[i];
-        }
-        return m;
-    }
-    std::vector<ModuleBase::matrix> vec2mat(const std::vector<std::vector<double>>& v, const int nr, const int nc)
-    {
-        std::vector<ModuleBase::matrix> m(v.size());
-        for (int i = 0;i < v.size();++i) {
-            m[i] = vec2mat(v[i], nr, nc);
-        }
-        return m;
-    }
-    std::vector<ModuleBase::ComplexMatrix> vec2mat(const std::vector<std::vector<std::complex<double>>>& v, const int nr, const int nc)
-    {
-        std::vector<ModuleBase::ComplexMatrix> m(v.size());
-        for (int i = 0;i < v.size();++i) {
-            m[i] = vec2mat(v[i], nr, nc);
-        }
-        return m;
-    }
 
     // for the first matrix in the commutator
     void setup_2d_division(Parallel_2D& pv, int nb, int gr, int gc)

source/module_lr/utils/lr_util.h

@@ -74,20 +74,6 @@ namespace LR_Util
     template<typename T>
     void matsym(T* inout, const int n, const Parallel_2D& pmat);
 #endif
-    ///======== Tensor - Matrix transformer==========
-    container::Tensor mat2ten_double(ModuleBase::matrix& m);
-    std::vector<container::Tensor> mat2ten_double(std::vector<ModuleBase::matrix>& m);
-    ModuleBase::matrix ten2mat_double(container::Tensor& t);
-    std::vector<ModuleBase::matrix> ten2mat_double(std::vector<container::Tensor>& t);
-    container::Tensor mat2ten_complex(ModuleBase::ComplexMatrix& m);
-    std::vector<container::Tensor> mat2ten_complex(std::vector<ModuleBase::ComplexMatrix>& m);
-    ModuleBase::ComplexMatrix ten2mat_complex(container::Tensor& t);
-    std::vector<ModuleBase::ComplexMatrix> ten2mat_complex(std::vector<container::Tensor>& t);
-
-    ModuleBase::matrix vec2mat(const std::vector<double>& v, const int nr, const int nc);
-    ModuleBase::ComplexMatrix vec2mat(const std::vector<std::complex<double>>& v, const int nr, const int nc);
-    std::vector<ModuleBase::matrix> vec2mat(const std::vector<std::vector<double>>& v, const int nr, const int nc);
-    std::vector<ModuleBase::ComplexMatrix> vec2mat(const std::vector<std::vector<std::complex<double>>>& v, const int nr, const int nc);
 
     ///===================Psi wrapper=================
     /// psi(nk=1, nbands=nb, nk * nbasis) -> psi(nb, nk, nbasis) without memory copy