ALPS MPS Codes
Reference documentation.
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
measurements::average< Matrix, SymmGroup > Class Template Reference

#include <average.h>

Inheritance diagram for measurements::average< Matrix, SymmGroup >:
measurement< Matrix, SymmGroup >

Public Types

typedef Matrix::value_type value_type
 

Public Member Functions

 average (std::string const &name_, const Lattice &lattice, op_vec const &identities, op_vec const &fillings, op_vec const &ops)
 
 average (std::string const &name_, const Lattice &lattice, op_vec const &identities, op_vec const &fillings, std::vector< bond_element > const &ops)
 
void evaluate (MPS< Matrix, SymmGroup > const &mps, boost::optional< reduced_mps< Matrix, SymmGroup > const & > rmps=boost::none)
 
template<class Archive >
void save (Archive &) const
 
void write_xml (alps::oxstream &) const
 
virtual void print (std::ostream &os) const
 
std::string const & name () const
 
int & eigenstate_index ()
 
int eigenstate_index () const
 
void set_super_meas (Index< SymmGroup > const &phys_psi_)
 
measurementclone () const
 

Protected Member Functions

measurement< Matrix, SymmGroup > * do_clone () const
 

Protected Attributes

bool cast_to_real
 
bool is_super_meas
 
std::vector< std::string > labels
 
value_type result
 
std::vector< typename MPS
< Matrix, SymmGroup >
::scalar_type > 
vector_results
 
Index< SymmGroup > phys_psi
 

Detailed Description

template<class Matrix, class SymmGroup>
class measurements::average< Matrix, SymmGroup >

Definition at line 39 of file average.h.

Member Typedef Documentation

template<class Matrix, class SymmGroup>
typedef Matrix::value_type measurement< Matrix, SymmGroup >::value_type
inherited

Definition at line 55 of file measurement.h.

Constructor & Destructor Documentation

template<class Matrix , class SymmGroup >
measurements::average< Matrix, SymmGroup >::average ( std::string const &  name_,
const Lattice lattice,
op_vec const &  identities,
op_vec const &  fillings,
op_vec const &  ops 
)
inline

Site term.

Definition at line 47 of file average.h.

51  : base(name_)
52  {
53  generator mpom(lattice, identities, fillings);
54 
55  for (std::size_t p = 0; p < lattice.size(); ++p)
56  {
58  term.operators.push_back( std::make_pair(p, ops[lattice.get_prop<int>("type", p)]) );
59  mpom.add_term(term);
60  }
61  mpo = mpom.create_mpo();
62 
63  this->cast_to_real = is_hermitian_meas(ops);
64  }
bool is_hermitian_meas(std::vector< block_matrix< Matrix, SymmGroup > > const &ops)
UTILITIES.
Definition: measurement.h:175
T get_prop(std::string property, pos_t site) const
Definition: lattice.h:103
std::vector< op_pair_t > operators
Definition: utils.hpp:118
bool cast_to_real
Definition: measurement.h:77
pos_t size() const
Definition: lattice.h:115
template<class Matrix , class SymmGroup >
measurements::average< Matrix, SymmGroup >::average ( std::string const &  name_,
const Lattice lattice,
op_vec const &  identities,
op_vec const &  fillings,
std::vector< bond_element > const &  ops 
)
inline

Bond term.

Definition at line 67 of file average.h.

71  : base(name_)
72  {
73  generator mpom(lattice, identities, fillings);
74 
75  for (std::size_t i = 0; i < ops.size(); ++i) {
76  for (std::size_t p = 0; p < lattice.size(); ++p) {
78  term.operators.push_back( std::make_pair(p, ops[i][0].first[lattice.get_prop<int>("type", p)]) );
79  term.with_sign = ops[i][0].second;
80  std::vector<Lattice::pos_t> neighs = lattice.forward(p);
81  for (typename std::vector<Lattice::pos_t>::const_iterator hopto = neighs.begin();
82  hopto != neighs.end(); ++hopto)
83  {
85  term2.operators.push_back( std::make_pair(*hopto, ops[i][1].first[lattice.get_prop<int>("type", p)]) );
86  mpom.add_term(term2);
87  }
88  }
89  }
90  mpo = mpom.create_mpo();
91 
92  /// TODO: this doesn't really work
93  this->cast_to_real = all_true(ops.begin(), ops.end(), static_cast<bool (*)(bond_element const&)>(&is_hermitian_meas));
94  }
bool all_true(InputIterator first, InputIterator last, Predicate pred)
Definition: utils.hpp:53
bool is_hermitian_meas(std::vector< block_matrix< Matrix, SymmGroup > > const &ops)
UTILITIES.
Definition: measurement.h:175
T get_prop(std::string property, pos_t site) const
Definition: lattice.h:103
std::vector< op_pair_t > operators
Definition: utils.hpp:118
std::vector< pos_t > forward(pos_t site) const
Definition: lattice.h:100
bool cast_to_real
Definition: measurement.h:77
pos_t size() const
Definition: lattice.h:115

Member Function Documentation

template<class Matrix, class SymmGroup>
measurement* measurement< Matrix, SymmGroup >::clone ( ) const
inlineinherited

Definition at line 72 of file measurement.h.

72 { return do_clone(); }
virtual measurement * do_clone() const =0
template<class Matrix , class SymmGroup >
measurement<Matrix, SymmGroup>* measurements::average< Matrix, SymmGroup >::do_clone ( ) const
inlineprotectedvirtual

Implements measurement< Matrix, SymmGroup >.

Definition at line 114 of file average.h.

115  {
116  return new average(*this);
117  }
average(std::string const &name_, const Lattice &lattice, op_vec const &identities, op_vec const &fillings, op_vec const &ops)
Site term.
Definition: average.h:47
template<class Matrix, class SymmGroup>
int& measurement< Matrix, SymmGroup >::eigenstate_index ( )
inlineinherited

Definition at line 67 of file measurement.h.

67 { return eigenstate; }
template<class Matrix, class SymmGroup>
int measurement< Matrix, SymmGroup >::eigenstate_index ( ) const
inlineinherited

Definition at line 68 of file measurement.h.

68 { return eigenstate; }
template<class Matrix , class SymmGroup >
void measurements::average< Matrix, SymmGroup >::evaluate ( MPS< Matrix, SymmGroup > const &  mps,
boost::optional< reduced_mps< Matrix, SymmGroup > const & >  rmps = boost::none 
)
inlinevirtual

Implements measurement< Matrix, SymmGroup >.

Definition at line 96 of file average.h.

97  {
98  if (!this->is_super_meas){
99  this->result = expval(mps, mpo);
100  } else {
101  typename MPS<Matrix, SymmGroup>::scalar_type nn = dm_trace(mps, this->phys_psi);
102  MPS<Matrix, SymmGroup> super_mpo = mpo_to_smps(mpo, this->phys_psi);
103  // static_cast needed for icpc 12.x
104 #ifdef __INTEL_COMPILER
105  typedef typename MPS<Matrix, SymmGroup>::scalar_type (*overlap_func)(MPS<Matrix, SymmGroup> const &, MPS<Matrix, SymmGroup> const &);
106  this->result = static_cast<overlap_func>(&overlap)(super_mpo, mps) / nn;
107 #else
108  this->result =::overlap<Matrix,SymmGroup>(super_mpo, mps) / nn;
109 #endif
110  }
111  }
double expval(MPS< Matrix, SymmGroup > const &mps, MPO< Matrix, SymmGroup > const &mpo, int d)
Definition: mps_mpo_ops.h:72
value_type result
Definition: measurement.h:79
bool is_super_meas
Definition: measurement.h:77
MPS< Matrix, SymmGroup >::scalar_type overlap(MPS< Matrix, SymmGroup > const &mps1, MPS< Matrix, SymmGroup > const &mps2)
Definition: mps_mpo_ops.h:136
Definition: mps.h:40
MPSTensor< Matrix, SymmGroup >::scalar_type scalar_type
Definition: mps.h:51
MPS< Matrix, SymmGroup > mpo_to_smps(MPO< Matrix, SymmGroup > const &mpo, Index< SymmGroup > const &phys_i)
Definition: super_mpo.h:60
Index< SymmGroup > phys_psi
Definition: measurement.h:82
MPS< Matrix, SymmGroup >::scalar_type dm_trace(MPS< Matrix, SymmGroup > const &mps, Index< SymmGroup > const &phys_psi)
Definition: mps_mpo_ops.h:258
template<class Matrix, class SymmGroup>
std::string const& measurement< Matrix, SymmGroup >::name ( ) const
inlineinherited

Definition at line 66 of file measurement.h.

66 { return name_; }
template<class Matrix , class SymmGroup >
void measurement< Matrix, SymmGroup >::print ( std::ostream &  os) const
virtualinherited

Definition at line 160 of file measurement.h.

161 {
162  os << "MEASURE[" << name_ << "]";
163 }
template<class Matrix , class SymmGroup >
template<class Archive >
void measurement< Matrix, SymmGroup >::save ( Archive &  ar) const
inherited

Definition at line 107 of file measurement.h.

108 {
109  if (vector_results.size() > 0) {
110  if (cast_to_real) {
111  save_val_at_index(ar, storage::encode(name()) + std::string("/mean/value"), maquis::real(vector_results), eigenstate_index());
112  } else {
113  save_val_at_index(ar, storage::encode(name()) + std::string("/mean/value"), vector_results, eigenstate_index());
114  }
115  if (labels.size() > 0)
116  ar[storage::encode(name()) + std::string("/labels")] << labels;
117  } else {
118  if (cast_to_real) {
119  save_val_at_index(ar, storage::encode(name()) + std::string("/mean/value"), maquis::real(result), eigenstate_index());
120  } else {
121  save_val_at_index(ar, storage::encode(name()) + std::string("/mean/value"), result, eigenstate_index());
122  }
123  }
124 }
std::vector< std::string > labels
Definition: measurement.h:78
value_type result
Definition: measurement.h:79
std::string encode(std::string const &s)
Definition: archive.h:82
void save_val_at_index(Archive &ar, std::string const &archive_path, T const &val, std::size_t eig)
Definition: measurement.h:96
std::vector< typename MPS< Matrix, SymmGroup >::scalar_type > vector_results
Definition: measurement.h:80
int & eigenstate_index()
Definition: measurement.h:67
bool cast_to_real
Definition: measurement.h:77
alps::numeric::real_type< T >::type real(T f)
Definition: bindings.hpp:38
std::string const & name() const
Definition: measurement.h:66
template<class Matrix , class SymmGroup >
void measurement< Matrix, SymmGroup >::set_super_meas ( Index< SymmGroup > const &  phys_psi_)
inherited

Definition at line 153 of file measurement.h.

154 {
155  phys_psi = phys_psi_;
156  is_super_meas = true;
157 }
bool is_super_meas
Definition: measurement.h:77
Index< SymmGroup > phys_psi
Definition: measurement.h:82
template<class Matrix , class SymmGroup >
void measurement< Matrix, SymmGroup >::write_xml ( alps::oxstream &  out) const
inherited

Definition at line 127 of file measurement.h.

128 {
129  if (labels.size() > 0) {
130  out << alps::start_tag("VECTOR_AVERAGE") << alps::attribute("name", name());
131  for (int i=0; i<labels.size(); ++i) {
132  out << alps::start_tag("SCALAR_AVERAGE");
133  out << alps::attribute("indexvalue",labels[i]) << alps::no_linebreak;
134 
135  if (cast_to_real) {
136  out << alps::start_tag("MEAN") << alps::no_linebreak << maquis::real(vector_results[i]) << alps::end_tag("MEAN");
137  } else {
138  out << alps::start_tag("MEAN") << alps::no_linebreak << vector_results[i] << alps::end_tag("MEAN");
139  }
140  out << alps::end_tag("SCALAR_AVERAGE");
141 
142  }
143  out << alps::end_tag("VECTOR_AVERAGE");
144  } else {
145  out << alps::start_tag("SCALAR_AVERAGE") << alps::attribute("name", name()) << alps::no_linebreak;
146  out << alps::start_tag("MEAN") << alps::no_linebreak << ((cast_to_real) ? maquis::real(result) : result) << alps::end_tag("MEAN");
147  out << alps::end_tag("SCALAR_AVERAGE");
148  }
149 }
std::vector< std::string > labels
Definition: measurement.h:78
value_type result
Definition: measurement.h:79
std::vector< typename MPS< Matrix, SymmGroup >::scalar_type > vector_results
Definition: measurement.h:80
bool cast_to_real
Definition: measurement.h:77
alps::numeric::real_type< T >::type real(T f)
Definition: bindings.hpp:38
std::string const & name() const
Definition: measurement.h:66

Member Data Documentation

template<class Matrix, class SymmGroup>
bool measurement< Matrix, SymmGroup >::cast_to_real
protectedinherited

Definition at line 77 of file measurement.h.

template<class Matrix, class SymmGroup>
bool measurement< Matrix, SymmGroup >::is_super_meas
protectedinherited

Definition at line 77 of file measurement.h.

template<class Matrix, class SymmGroup>
std::vector<std::string> measurement< Matrix, SymmGroup >::labels
protectedinherited

Definition at line 78 of file measurement.h.

template<class Matrix, class SymmGroup>
Index<SymmGroup> measurement< Matrix, SymmGroup >::phys_psi
protectedinherited

Definition at line 82 of file measurement.h.

template<class Matrix, class SymmGroup>
value_type measurement< Matrix, SymmGroup >::result
protectedinherited

Definition at line 79 of file measurement.h.

template<class Matrix, class SymmGroup>
std::vector<typename MPS<Matrix, SymmGroup>::scalar_type> measurement< Matrix, SymmGroup >::vector_results
protectedinherited

Definition at line 80 of file measurement.h.


The documentation for this class was generated from the following file: