Only retrieve end thermodynamic states if we have not done so already
@@ -620,10 +620,12 @@ class PhaseAnalyzer(ABC):
The end thermodynamic states
"""
+ if self._end_thermodynamic_states is None:
self._end_thermodynamic_states = self._reporter.read_end_thermodynamic_states()
# Cache other useful informations since we have already read this.
# TODO should we read temperatures of all the states and let kT property depend on reference_states?
self._kT = self._end_thermodynamic_states[0].kT
+
return self._end_thermodynamic_states
@property
