Enable PT for all frustmag sims, incl. Heisenberg

frustmag/config_frustmag_sim.hpp

@@ -84,15 +84,6 @@ namespace {
 
 template <typename Hamiltonian>
 struct update_t {
-    template <typename LatticeH>
-    using type = update::muxer<
-        update::single_flip
-        , update::global_trafo
-        >::type<LatticeH>;
-};
-
-template <template <typename> typename Lat>
-struct update_t<hamiltonian::ising<Lat>> {
     template <typename LatticeH>
     using type = update::muxer<
         update::single_flip
@@ -108,6 +99,7 @@ struct update_t<hamiltonian::heisenberg<Lat>> {
         update::single_flip
         , update::overrelaxation
         , update::global_trafo
+        , update::parallel_tempering
         >::type<LatticeH>;
 };
 

frustmag/hamiltonian/heisenberg.hpp

@@ -22,6 +22,7 @@
 #include "update/single_flip.hpp"
 #include "update/overrelaxation.hpp"
 #include "update/global_trafo.hpp"
+#include "update/parallel_tempering.hpp"
 
 #include <alps/params.hpp>
 #include <alps/hdf5.hpp>
@@ -163,6 +164,10 @@ struct heisenberg {
         return true;
     }
 
+    double log_weight(phase_point const& other) const {
+        return (ppoint.temp - other.temp) * energy();
+    }
+
     virtual void save(alps::hdf5::archive & ar) const {
         ar["phase_point"] << std::vector<double>{ppoint.begin(), ppoint.end()};
         ar["lattice"] << lattice_;