Commit 6d606120 by Francois Gygi

Compute polarization using ElectricEnthalpy whenever polarization is

not OFF. This can be used to compute the polarization in zero field.

git-svn-id: cba15fb0-1239-40c8-b417-11db7ca47a34
parent 1191e431
......@@ -162,7 +162,7 @@ EnergyFunctional::EnergyFunctional( Sample& s, const ChargeDensity& cd)
// Electric enthalpy
el_enth_ = 0;
if ( norm2(s_.ctrl.e_field) != 0.0 )
if ( s_.ctrl.polarization != "OFF" )
el_enth_ = new ElectricEnthalpy(s_);
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment