Commit 76e263d6 authored by Alvarez, Gonzalo's avatar Alvarez, Gonzalo
Browse files

dyndmrg: type n+1 or n-1 added

parent 30f9ccd7
%% File CHANGES
What's new in Version 2
* Engine: SU(2) symmetry now supported and integrated
* Engine: Checkpointing
* Engine: DiskStack to support checkpointing
* Engine: Faster WaveFunctionTransformation for local symmetries
* User Interface: Customizable finite sweeps
* New Model: FeBasedSc for Fe-based Superconductors
* New Geometry: GeometryLadderFeAs to go with FeBasedSc
* Tests: TestSuite added (under /TestSuite)
* Documentation: Manual added (under /doc)
%% Extracted from the DMRG++ Manual
%% by ../scripts/extractFromManual.pl for CHANGES
See
https://github.com/g1257/dmrgpp/commits/master
......@@ -195,11 +195,6 @@ namespace Dmrg {
{
if (direction==INFINITE) return;
stage_ = ENABLED;
// block.size()==1 or throw
// size_t i = block[0]; // center site
// FermionSign fs(lrs_.left(),correctionStruct_.electrons);
// const BasisWithOperatorsType *basis = (direction==EXPAND_SYSTEM) ?
// &(lrs_.left()) : &(lrs_.right());
// operators in the one-site basis:
std::vector<OperatorType> creationMatrix;
......
......@@ -109,7 +109,9 @@ namespace Dmrg {
{
io.rewind();
this->concatenation = SUM;
io.readline(type,"DynamicDmrgType=");
}
size_t type;
}; // class DynamicDmrgParams
......@@ -121,6 +123,7 @@ namespace Dmrg {
const typename TimeStepParams<ModelType>::TargetParamsCommonType&
tp = t;
os<<tp;
os<<"DynamicDmrgType="<<t.type<<"\n";
return os;
}
} // namespace Dmrg
......
......@@ -266,8 +266,11 @@ namespace Dmrg {
{
if (block.size()!=1) throw std::runtime_error(
"DynamicTargetting only supports blocks of size 1\n");
ContinuedFractionType cf(ab_,Eg_,weightForContinuedFraction_);
size_t type = tstStruct_.type;
int s = (type&1) ? -1 : 1;
int s2 = (type>1) ? -1 : 1;
ContinuedFractionType cf(ab_,Eg_,s2*weightForContinuedFraction_,
s);
DynamicSerializerType dynS(cf,block[0],targetVectors_);
dynS.save(io);
psi_.save(io,"PSI");
......
......@@ -366,8 +366,8 @@ namespace Dmrg {
//weStack_.pop();
//std::cerr<<"CHANGED-COUNTER0 We taken from stack\n";
} else {
std::cerr<<"PUSHING-COUNTER0 STACK ERROR E\n";
throw std::runtime_error("Environ Stack is empty\n");
// std::cerr<<"PUSHING-COUNTER0 STACK ERROR E\n";
// throw std::runtime_error("Environ Stack is empty\n");
}
}
......@@ -386,8 +386,8 @@ namespace Dmrg {
//weStack_.pop();
//std::cerr<<"CHANGED-COUNTER0 We taken from stack\n";
} else {
std::cerr<<"PUSHING-COUNTER0 STACK ERROR E\n";
throw std::runtime_error("System Stack is empty\n");
// std::cerr<<"PUSHING-COUNTER0 STACK ERROR E\n";
// throw std::runtime_error("System Stack is empty\n");
}
}
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment