Commit 82bf85bc authored by Alvarez, Gonzalo's avatar Alvarez, Gonzalo
Browse files

CorrectionVectorAction: marking Matsubara and Complex Correction Vector cases

parent 0dc8c87e
......@@ -2,6 +2,7 @@
#define CORRECTIONVECTORACTION_H
#include "Vector.h"
#include "FreqEnum.h"
#include "AnsiColors.h"
namespace Dmrg {
......@@ -21,12 +22,15 @@ public:
RealType E0,
const VectorRealType& eigs)
: tstStruct_(tstStruct),E0_(E0),eigs_(eigs)
{}
ComplexOrRealType operator()(SizeType k) const
{
return (tstStruct_.omega().first == PsimagLite::FREQ_REAL) ? actionWhenFreqReal(k)
: actionWhenMatsubara(k);
if (tstStruct_.omega().first == PsimagLite::FREQ_REAL)
return; // <--- EARLY EXIT
std::cout<<PsimagLite::AnsiColor::red;
std::cerr<<PsimagLite::AnsiColor::red;
std::cout<<"CorrectionVectorActionBase:: Running matsubara\n";
std::cerr<<"CorrectionVectorActionBase:: Running matsubara\n";
std::cout<<PsimagLite::AnsiColor::reset;
std::cerr<<PsimagLite::AnsiColor::reset;
}
void setReal() const;
......@@ -37,8 +41,6 @@ public:
protected:
ComplexOrRealType actionWhenFreqReal(SizeType k) const;
RealType actionWhenMatsubara(SizeType k) const
{
RealType sign = (tstStruct_.type() == 0) ? -1.0 : 1.0;
......@@ -109,7 +111,15 @@ public:
CorrectionVectorAction(const TargetParamsType& tstStruct,
RealType E0,
const typename BaseType::VectorRealType& eigs)
: BaseType(tstStruct, E0, eigs) {}
: BaseType(tstStruct, E0, eigs)
{
std::cout<<PsimagLite::AnsiColor::red;
std::cerr<<PsimagLite::AnsiColor::red;
std::cout<<"CorrectionVectorAction:: Complex mode\n";
std::cerr<<"CorrectionVectorAction:: Complex mode\n";
std::cout<<PsimagLite::AnsiColor::reset;
std::cerr<<PsimagLite::AnsiColor::reset;
}
ComplexOrRealType operator()(SizeType k) const
{
......
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