https://github.com/ITensor/ITensor
Tip revision: 6280c1a7de2effbc23eb21feaa87812bde0ac91b authored by Miles Stoudenmire on 09 August 2017, 04:50:32 UTC
New exactApplyMPO algorithm with better scaling (m^3k^2+m^2k^3) and less ad hoc approach to compression of resulting MPS
New exactApplyMPO algorithm with better scaling (m^3k^2+m^2k^3) and less ad hoc approach to compression of resulting MPS
Tip revision: 6280c1a
model_test.cc
#include "test.h"
#include "model/spinhalf.h"
#include "model/spinone.h"
#include "model/hubbard.h"
#include "model/spinless.h"
#include "model/tj.h"
TEST_CASE("ModelTest")
{
const int N = 10;
SECTION("SpinHalfModel")
{
SpinHalf model(N);
model.op("Sz",2);
model.op("S+",2);
model.op("S-",2);
model.op("Sp",2);
model.op("Sm",2);
model.op("Sx",2);
model.op("Sy",2);
model.op("ISy",2);
}
SECTION("SpinOneModel")
{
SpinOne model(N);
model.op("Sz",2);
model.op("S+",2);
model.op("S-",2);
model.op("Sp",2);
model.op("Sm",2);
model.op("Sx",2);
model.op("Sy",2);
model.op("ISy",2);
}
SECTION("HubbardModel")
{
Hubbard model(N);
model.op("Nup",2);
model.op("Ndn",2);
model.op("Nupdn",2);
model.op("Ntot",2);
model.op("Sz",2);
model.op("Cup",2);
model.op("Cdn",2);
model.op("Aup",2);
model.op("Adn",2);
model.op("F",2);
}
SECTION("SpinlessModel")
{
Spinless model(N);
model.op("N",2);
model.op("A",2);
model.op("Adag",2);
model.op("F",2);
}
SECTION("tJModel")
{
tJ model(N);
model.op("Nup",2);
model.op("Ndn",2);
model.op("Aup",2);
model.op("Adn",2);
model.op("F",2);
}
}