fix the waste of memory in matsym

Author: maki49

source/module_lr/utils/lr_util.cpp

@@ -62,37 +62,33 @@ namespace LR_Util
     template<>
     void matsym<double>(const double* in, const int n, const Parallel_2D& pmat, double* out)
     {
-        for (int i = 0;i < pmat.get_local_size();++i) {out[i] = in[i];
-}
+        for (int i = 0;i < pmat.get_local_size();++i) { out[i] = in[i]; }
         const double alpha = 0.5, beta = 0.5;
         const int i1 = 1;
         pdtran_(&n, &n, &alpha, in, &i1, &i1, pmat.desc, &beta, out, &i1, &i1, pmat.desc);
     }
     template<>
     void matsym<double>(double* inout, const int n, const Parallel_2D& pmat)
     {
-        std::vector<double> tmp(n * n);
-        for (int i = 0;i < pmat.get_local_size();++i) {tmp[i] = inout[i];
-}
+        std::vector<double> tmp(pmat.get_local_size());
+        std::copy(inout, inout + pmat.get_local_size(), tmp.begin());
         const double alpha = 0.5, beta = 0.5;
         const int i1 = 1;
         pdtran_(&n, &n, &alpha, tmp.data(), &i1, &i1, pmat.desc, &beta, inout, &i1, &i1, pmat.desc);
     }
     template<>
     void matsym<std::complex<double>>(const std::complex<double>* in, const int n, const Parallel_2D& pmat, std::complex<double>* out)
     {
-        for (int i = 0;i < pmat.get_local_size();++i) {out[i] = in[i];
-}
+        for (int i = 0;i < pmat.get_local_size();++i) { out[i] = in[i]; }
         const std::complex<double> alpha(0.5, 0.0), beta(0.5, 0.0);
         const int i1 = 1;
         pztranc_(&n, &n, &alpha, in, &i1, &i1, pmat.desc, &beta, out, &i1, &i1, pmat.desc);
     }
     template<>
     void matsym<std::complex<double>>(std::complex<double>* inout, const int n, const Parallel_2D& pmat)
     {
-        std::vector<std::complex<double>> tmp(n * n);
-        for (int i = 0;i < pmat.get_local_size();++i) {tmp[i] = inout[i];
-}
+        std::vector<std::complex<double>> tmp(pmat.get_local_size());
+        std::copy(inout, inout + pmat.get_local_size(), tmp.begin());
         const std::complex<double> alpha(0.5, 0.0), beta(0.5, 0.0);
         const int i1 = 1;
         pztranc_(&n, &n, &alpha, tmp.data(), &i1, &i1, pmat.desc, &beta, inout, &i1, &i1, pmat.desc);