https://github.com/ITensor/ITensor
Tip revision: e4f80002b0aeb8034681b1ed07666709dc803adf authored by Miles Stoudenmire on 18 May 2015, 23:21:44 UTC
Removed static keywords from workspace vectors in lapack_wrap.h. Removed a variable sizes array from utility.cc. Fixes #79 - thanks kyungminlee
Removed static keywords from workspace vectors in lapack_wrap.h. Removed a variable sizes array from utility.cc. Fixes #79 - thanks kyungminlee
Tip revision: e4f8000
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