summaryrefslogtreecommitdiff
path: root/doc/measurement.rst
blob: 1ca41ccb488a2af4177deb32f18ea70516317953 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67

.. default-domain:: cpp

************
Measurements
************

One of the most complicated parts of using this library is setting up measurements to take during a run, but once understood they provide a powerful way to quickly measure arbitrary properties. The class is defined in the header :file:`lib/wolff.hpp`.

.. class:: template<R_t, X_t, G_t = graph<>> measurement

   This class defines several virtual member functions. To use the library, you must supply your own measurement class that inherits this one and defines those functions, which may be trivial.

   Each member function defines a hook that is run at a specific time during execution. These hooks can be used to modify member objects of your inheritor measurement class, and thereby extract information from the simulation.

   .. function:: void measurement::pre_cluster(unsigned n, unsigned N, const system<R_t, X_t, G_t>& S, const G_t::vertex& v, const R_t& r)

      This hook is run prior to cluster formation.

      :param unsigned n: The number of clusters already flipped.
      :param unsigned N: The total number of clusters to flip.
      :param const system<R_t, X_t, G_t>& S: The current system state.
      :param const G_t\:\:vertex& v: The seed vertex for the cluster.
      :param const R_t& r: The transformation by which the cluster will be flipped.

   .. function:: void measurement::plain_bond_visited(const system<R_t, X_t, G_t>& S, const G_t::halfedge& e, const X_t& s1_new, double dE)

         This hook is run when an ordinary edge is visited during cluster formation, whether it is activated or not.

         :param const system<R_t, X_t, G_t>& S: The current system state.
         :param const G_t\:\:halfedge& e:   The edge being considered, with e.self referencing the vertex soon to be flipped.
         :param const X_t& s1_new:         The state the vertex will be flipped to.
         :param double dE:                 The change in energy that will result when the spin is flipped.

   .. cpp:function:: void measurement::plain_site_transformed(const system<R_t, X_t, G_t>& S, const G_t::vertex& v, const X_t& si_new)

         This hook is run when an ordinary site is about to be flipped.

         :param const system<R_t, X_t, G_t>& S: The current system state.
         :param const G_t\:\:vertex& v:                     The vertex soon to be flipped.
         :param const X_t& si_new:          The state that vertex will be flipped to.

   .. cpp:function:: void measurement::ghost_bond_visited(const system<R_t, X_t, G_t>& S, const G_t::vertex& v, const X_t& s0si_old, const X_t& s0si_new, double dE)

         This hook is run when an edge containing the ghost site is visited during cluster formation, whether activated or not.

         :param const system<R_t, X_t, G_t>& S: The current system state.
         :param const G_t\:\:vertex& v:                     The non-ghost site in this edge.
         :param const X_t& s0si_old:       The state that the spin on the non-ghost site has before the transformation is applied, rotated by the inverse action of the ghost site.
         :param const X_t& s0si_new:       The state that the spin on the non-ghost site will have after the transformation is applied, rotated by the inverse action of the ghost site.
         :param double dE:                 The change in energy that will result when one site is transformed.

   .. cpp:function:: void measurement::ghost_site_transformed(const system<R_t, X_t, G_t>& S, const R_t& R_new)

         This hook is run when the ghost site is about to be flipped.

         :param const system<R_t, X_t, G_t>& S: The current system state.
         :param const R_t& R_new:          The state the ghost site will be flipped to.

   .. cpp:function:: void measurement::post_cluster(unsigned n, unsigned N, const system<R_t, X_t, G_t>& S) 

         This hook is run after a cluster has been flipped.

         :param unsigned n:                     The number of clusters already flipped.
         :param unsigned N:                     The total number of clusters to flip.
         :param const system<R_t, X_t, G_t>& S: The current system state.