Added missing imaginary number in ucorr

ttadano · Nov 29, 2018 · 1d1e034 · 1d1e034
1 parent fd47005
commit 1d1e034
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/anphon/thermodynamics.cpp b/anphon/thermodynamics.cpp
@@ -324,7 +324,7 @@ double Thermodynamics::disp_corrfunc(const double T_in,
     unsigned int ik, is;
     double omega;
     double phase;
-
+    const std::complex<double> im(0.0, 1.0);
     double ret = 0.0;
 
     if (classical) {
@@ -362,7 +362,7 @@ double Thermodynamics::disp_corrfunc(const double T_in,
 
             ret += real(std::conj(evec_in[ik][is][ncrd1])
                     * evec_in[ik][is][ncrd2]
-                    * std::exp(phase))
+                    * std::exp(im * phase))
                 * (fB(omega, T_in) + 0.5) / omega;
         }
     }

diff --git a/anphon/write_phonons.h b/anphon/write_phonons.h
@@ -40,7 +40,7 @@ namespace PHON_NS
 
         unsigned int anime_cellsize[3];
         double anime_kpoint[3];
-        unsigned int shift_ucorr[3];
+        int shift_ucorr[3];
 
         double in_kayser(const double) const;