https://github.com/ITensor/ITensor
Tip revision: 12b397f26e1422cd5b2bfb325e691ce898caebc3 authored by Miles Stoudenmire on 13 August 2015, 19:52:28 UTC
Quieted most make output, except which files are currently compiling
Quieted most make output, except which files are currently compiling
Tip revision: 12b397f
spectrum.h
//
// Distributed under the ITensor Library License, Version 1.1.
// (See accompanying LICENSE file.)
//
#ifndef __ITENSOR_SPECTRUM_H
#define __ITENSOR_SPECTRUM_H
#include "iqtensor.h"
#include "iterpair.h"
namespace itensor {
//
// Spectrum
//
// Stores density matrix eigenvalue spectrum following a call to
// svd, denmatDecomp, etc.
//
class Spectrum
{
public:
using QNStorage = std::vector<QN>;
Spectrum(const Args& args = Global::args());
Spectrum(const ITensor& D, const Args& args = Global::args());
Spectrum(const IQTensor& D, const Args& args = Global::args());
Spectrum(const Vector& eigs, const Args& args = Global::args());
Spectrum(const Vector& eigs,
const QNStorage& qns,
const Args& args = Global::args());
QN
qn(int n) const;
const QNStorage&
qns() const { return qns_; }
Real
eig(int n) const { return eigs_(n); }
const Vector&
eigs() const { return eigs_; }
Real
truncerr() const { return truncerr_; }
int
size() const { return eigs_.Length(); }
bool
hasQNs() const { return !qns_.empty(); }
//
// Other Methods
//
const Vector&
eigsKept() const { return eigs_; }
int
numEigsKept() const { return eigs_.Length(); }
void
truncerr(Real val) { truncerr_ = val; }
void
eigsKept(const Vector& val) { eigs_ = val; }
void
read(std::istream& s);
void
write(std::ostream& s) const;
private:
/////////////////
Vector eigs_;
Real truncerr_;
std::vector<QN> qns_;
/////////////////
void
computeTruncerr(const Args& args);
}; //class Spectrum
std::ostream&
operator<<(std::ostream & s, const Spectrum& spec);
} //namespace itensor
#endif