PairEx.cpp

Go to the documentation of this file.

/* HPhi  -  Quantum Lattice Model Simulator */
/* Copyright (C) 2015 The University of Tokyo */

/* This program is free software: you can redistribute it and/or modify */
/* it under the terms of the GNU General Public License as published by */
/* the Free Software Foundation, either version 3 of the License, or */
/* (at your option) any later version. */

/* This program is distributed in the hope that it will be useful, */
/* but WITHOUT ANY WARRANTY; without even the implied warranty of */
/* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the */
/* GNU General Public License for more details. */

/* You should have received a copy of the GNU General Public License */
/* along with this program.  If not, see <http://www.gnu.org/licenses/>. */
/*-------------------------------------------------------------*/
#include "bitcalc.hpp"
#include "mltplyCommon.hpp"
#include "PairEx.hpp"
#include "PairExHubbard.hpp"
#include "PairExSpin.hpp"

int GetPairExcitedState
(
  struct BindStruct *X,
  int nstate, std::complex<double> **tmp_v0,
  std::complex<double> **tmp_v1,
  int iEx
)
{
  int iret;
  long int irght, ilft, ihfbit;

  //  i_max = X->Check.idim_max;
  if (X->Def.iFlgGeneralSpin == FALSE) {
    if (GetSplitBitByModel(X->Def.Nsite, X->Def.iCalcModel, &irght, &ilft, &ihfbit) != 0) {
      return -1;
    }
  }
  else {
    if (GetSplitBitForGeneralSpin(X->Def.Nsite, &ihfbit, X->Def.SiteToBit) != 0) {
      return -1;
    }
  }

  X->Large.i_max = X->Check.idim_maxOrg;
  X->Large.irght = irght;
  X->Large.ilft = ilft;
  X->Large.ihfbit = ihfbit;
  X->Large.mode = M_CALCSPEC;

  switch (X->Def.iCalcModel) {
  case HubbardGC:
    iret = GetPairExcitedStateHubbardGC(X, nstate, tmp_v0, tmp_v1, iEx);
    break;

  case KondoGC:
  case Hubbard:
  case Kondo:
    iret = GetPairExcitedStateHubbard(X, nstate, tmp_v0, tmp_v1, iEx);
    break;

  case Spin: // for the Sz-conserved spin system
    iret = GetPairExcitedStateSpin(X, nstate, tmp_v0, tmp_v1, iEx);
    break;

  case SpinGC:
    iret = GetPairExcitedStateSpinGC(X, nstate, tmp_v0, tmp_v1, iEx);
    break;

  default:
    iret = FALSE;
    break;
  }

  return iret;
}