https://github.com/ITensor/ITensor
Tip revision: 7e0c480d71533289c95c67315b8c81cd1e2fed74 authored by Matthew Fishman on 18 October 2018, 18:22:10 UTC
Change *.ih -> *_impl.h
Change *.ih -> *_impl.h
Tip revision: 7e0c480
ExpIsing.h
//
// Distributed under the ITensor Library License, Version 1.1.
// (See accompanying LICENSE file.)
//
#ifndef __EXP_ISING_H
#define __EXP_ISING_H
#include "itensor/mps/mpo.h"
namespace itensor {
class ExpIsing
{
SpinHalf const& sites_;
int N_ = 0;
Cplx tau_ = 0;
Real h_ = 0;
bool initted_ = false;
MPO H_;
public:
ExpIsing(SpinHalf const& sites,
Cplx tau,
Args const& args = Global::args());
operator MPO() { init_(); return H_; }
private:
void
init_();
}; //class ExpIsing
inline ExpIsing::
ExpIsing(SpinHalf const& sites,
Cplx tau,
Args const& args)
: sites_(sites),
tau_(tau),
initted_(false)
{
N_ = sites_.N();
h_ = args.getReal("h",0.);
}
void inline ExpIsing::
init_()
{
if(initted_) return;
H_ = MPO(sites_);
auto links = std::vector<Index>(N_+1);
for(int l = 0; l <= N_; ++l)
{
links.at(l) = Index(nameint("hl",l),2);
}
for(int n = 1; n <= N_; ++n)
{
auto& W = H_.Aref(n);
auto row = dag(links.at(n-1));
auto col = links.at(n);
W = ITensor(sites_(n),prime(sites_(n)),row,col);
W += sites_.op("Id",n) * row(1) * col(1);
W += -tau_ * (-h_) * ITensor(sites_.op("Sx",n)) * row(1) * col(1);
W += -tau_ * sites_.op("Sz",n) * row(1) * col(2);
W += sites_.op("Sz",n) * row(2) * col(1);
}
H_.Aref(1) *= setElt(links.at(0)(1));
H_.Aref(N_) *= setElt(links.at(N_)(1));
initted_ = true;
}
} //namespace itensor
#endif