Added a new speedup while computing the harmonic free energy

Author: mesonepigreco

Modules/SchaMinimizer.py

@@ -1136,7 +1136,7 @@ def run(self, verbose = 1, custom_function_pre = None, custom_function_post = No
             self.__fe__.append(np.real(fe))
             self.__fe_err__.append(np.real(err))
 
-            harm_fe = self.dyn.GetHarmonicFreeEnergy(self.ensemble.current_T) / np.prod(self.ensemble.supercell)
+            harm_fe = self.dyn.GetHarmonicFreeEnergy(self.ensemble.current_T, w_pols = (self.ensemble.current_w, self.ensemble.current_pols) / np.prod(self.ensemble.supercell)
             anharm_fe = np.real(fe - harm_fe)
 
             # Compute the KL ratio