expec_cisajscktaltdc.cpp File Reference

File for calculating two-body green's functions. More...

#include "mltply.hpp"
#include "mltplyCommon.hpp"
#include "FileIO.hpp"
#include "bitcalc.hpp"
#include "expec_cisajscktaltdc.hpp"
#include "mltplySpinCore.hpp"
#include "mltplyHubbardCore.hpp"
#include "wrapperMPI.hpp"
#include "mltplyMPISpin.hpp"
#include "mltplyMPISpinCore.hpp"
#include "mltplyMPIHubbardCore.hpp"
#include "common/setmemory.hpp"

Go to the source code of this file.

Functions
int	Rearray_Interactions (int i, long int *org_isite1, long int *org_isite2, long int *org_isite3, long int *org_isite4, long int *org_sigma1, long int *org_sigma2, long int *org_sigma3, long int *org_sigma4, std::complex< double > *tmp_V, struct BindStruct *X)
	Rearray interactions. More...
int	expec_cisajscktalt_HubbardGC (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Child function to calculate two-body green's functions for Hubbard GC model. More...
int	expec_cisajscktalt_Hubbard (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Child function to calculate two-body green's functions for Hubbard model. More...
int	expec_cisajscktalt_SpinHalf (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Child function to calculate two-body green's functions for 1/2 Spin model. More...
int	expec_cisajscktalt_SpinGeneral (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Child function to calculate two-body green's functions for General Spin model. More...
int	expec_cisajscktalt_SpinGCHalf (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Child function to calculate two-body green's functions for 1/2 Spin GC model. More...
int	expec_cisajscktalt_SpinGCGeneral (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Child function to calculate two-body green's functions for General Spin GC model. More...
int	expec_cisajscktalt_Spin (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Parent function to calculate two-body green's functions for Spin model. More...
int	expec_cisajscktalt_SpinGC (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec, std::complex< double > **prod)
	Parent function to calculate two-body green's functions for Spin GC model. More...
int	expec_cisajscktaltdc (struct BindStruct *X, int nstate, std::complex< double > **Xvec, std::complex< double > **vec)
	Parent function to calculate two-body green's functions. More...

Detailed Description

File for calculating two-body green's functions.

Version

0.2

add function to treat the case of general spin

Version

0.1

Author

Takahiro Misawa (The University of Tokyo)

Kazuyoshi Yoshimi (The University of Tokyo)

Definition in file expec_cisajscktaltdc.cpp.

Function Documentation

expec_cisajscktalt_Hubbard()

int expec_cisajscktalt_Hubbard	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Child function to calculate two-body green's functions for Hubbard model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 258 of file expec_cisajscktaltdc.cpp.

References LargeList::A_spin, LargeList::B_spin, CheckPE(), child_CisAisCisAis_element(), child_CisAisCjtAku_element(), child_CisAjtCkuAku_element(), child_CisAjtCkuAlv_element(), child_general_int_GetInfo(), DefineList::CisAjtCkuAlvDC, BindStruct::Def, LargeList::i_max, DefineList::iFlgSzConserved, LargeList::is1_spin, LargeList::is2_spin, LargeList::is3_spin, LargeList::is4_spin, LargeList::isA_spin, LargeList::isB_spin, BindStruct::Large, MultiVecProdMPI(), DefineList::NCisAjtCkuAlvDC, DefineList::OrgTpow, X_child_CisAisCjtAjt_Hubbard_MPI(), X_child_CisAisCjtAku_Hubbard_MPI(), X_child_CisAjtCkuAku_Hubbard_MPI(), X_child_CisAjtCkuAlv_Hubbard_MPI(), and zclear().

Referenced by expec_cisajscktaltdc().

 {
  long int i, j;
  long int isite1, isite2, isite3, isite4;
  long int org_isite1, org_isite2, org_isite3, org_isite4;
  long int org_sigma1, org_sigma2, org_sigma3, org_sigma4;
  long int Asum, Bsum, Adiff, Bdiff;
  long int tmp_off = 0;
  long int tmp_off_2 = 0;
  std::complex<double> tmp_V;
  long int i_max;

  for (i = 0; i < X->Def.NCisAjtCkuAlvDC; i++) {
    zclear(X->Large.i_max*nstate, &Xvec[1][0]);
    org_isite1 = X->Def.CisAjtCkuAlvDC[i][0] + 1;
    org_sigma1 = X->Def.CisAjtCkuAlvDC[i][1];
    org_isite2 = X->Def.CisAjtCkuAlvDC[i][2] + 1;
    org_sigma2 = X->Def.CisAjtCkuAlvDC[i][3];
    org_isite3 = X->Def.CisAjtCkuAlvDC[i][4] + 1;
    org_sigma3 = X->Def.CisAjtCkuAlvDC[i][5];
    org_isite4 = X->Def.CisAjtCkuAlvDC[i][6] + 1;
    org_sigma4 = X->Def.CisAjtCkuAlvDC[i][7];
    tmp_V = 1.0;

    if (X->Def.iFlgSzConserved == TRUE) {
      if (org_sigma1 + org_sigma3 != org_sigma2 + org_sigma4) {
        zclear(nstate, prod[i]);
        continue;
      }
    }

    if (CheckPE(org_isite1 - 1, X) == TRUE || CheckPE(org_isite2 - 1, X) == TRUE ||
      CheckPE(org_isite3 - 1, X) == TRUE || CheckPE(org_isite4 - 1, X) == TRUE) {
      isite1 = X->Def.OrgTpow[2 * org_isite1 - 2 + org_sigma1];
      isite2 = X->Def.OrgTpow[2 * org_isite2 - 2 + org_sigma2];
      isite3 = X->Def.OrgTpow[2 * org_isite3 - 2 + org_sigma3];
      isite4 = X->Def.OrgTpow[2 * org_isite4 - 2 + org_sigma4];
      if (isite1 == isite2 && isite3 == isite4) {
        X_child_CisAisCjtAjt_Hubbard_MPI(org_isite1 - 1, org_sigma1,
          org_isite3 - 1, org_sigma3, 1.0, X, nstate, Xvec, vec);
      }
      else if (isite1 == isite2 && isite3 != isite4) {
        //printf("org_isite1=%d, org_isite2=%d, org_isite3=%d, org_isite4=%d\n", org_isite1, org_isite2, org_isite3, org_isite4);
        X_child_CisAisCjtAku_Hubbard_MPI(org_isite1 - 1, org_sigma1,
          org_isite3 - 1, org_sigma3, org_isite4 - 1, org_sigma4, 1.0, X, nstate, Xvec, vec);
      }
      else if (isite1 != isite2 && isite3 == isite4) {
        X_child_CisAjtCkuAku_Hubbard_MPI(org_isite1 - 1, org_sigma1, org_isite2 - 1, org_sigma2,
          org_isite3 - 1, org_sigma3, 1.0, X, nstate, Xvec, vec);

      }
      else if (isite1 != isite2 && isite3 != isite4) {
        X_child_CisAjtCkuAlv_Hubbard_MPI(org_isite1 - 1, org_sigma1, org_isite2 - 1, org_sigma2,
          org_isite3 - 1, org_sigma3, org_isite4 - 1, org_sigma4, 1.0, X, nstate, Xvec, vec);
      }
    }//InterPE
    else {
      child_general_int_GetInfo(
        X, org_isite1, org_isite2, org_isite3, org_isite4,
        org_sigma1, org_sigma2, org_sigma3, org_sigma4, tmp_V
      );

      i_max = X->Large.i_max;
      isite1 = X->Large.is1_spin;
      isite2 = X->Large.is2_spin;
      Asum = X->Large.isA_spin;
      Adiff = X->Large.A_spin;

      isite3 = X->Large.is3_spin;
      isite4 = X->Large.is4_spin;
      Bsum = X->Large.isB_spin;
      Bdiff = X->Large.B_spin;

      tmp_V = 1.0;
      if (isite1 == isite2 && isite3 == isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,tmp_V,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2)
        for (j = 1; j <= i_max; j++) {
          child_CisAisCisAis_element(j, isite1, isite3, tmp_V, nstate, Xvec, vec);
        }
      }
      else if (isite1 == isite2 && isite3 != isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,tmp_V,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2)
        for (j = 1; j <= i_max; j++) {
          child_CisAisCjtAku_element(j, isite1, isite3, isite4, Bsum, Bdiff, 
            tmp_V, nstate, Xvec, vec, X, &tmp_off);
        }
      }
      else if (isite1 != isite2 && isite3 == isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,tmp_V,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2)
        for (j = 1; j <= i_max; j++) {
          child_CisAjtCkuAku_element(j, isite1, isite2, isite3, Asum, Adiff, 
            tmp_V, nstate, Xvec, vec, X, &tmp_off);
        }
      }
      else if (isite1 != isite2 && isite3 != isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,tmp_V,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2)
        for (j = 1; j <= i_max; j++) {
          child_CisAjtCkuAlv_element(j, isite1, isite2, isite3, isite4, Asum, Adiff, Bsum, Bdiff, 
            tmp_V, nstate, Xvec, vec, X, &tmp_off_2);
        }
      }
    }
    MultiVecProdMPI(i_max, nstate, vec, Xvec, prod[i]);
  }
  return 0;
}

expec_cisajscktalt_HubbardGC()

int expec_cisajscktalt_HubbardGC	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Child function to calculate two-body green's functions for Hubbard GC model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 145 of file expec_cisajscktaltdc.cpp.

References LargeList::A_spin, LargeList::B_spin, CheckPE(), child_general_int_GetInfo(), DefineList::CisAjtCkuAlvDC, BindStruct::Def, GC_child_CisAisCisAis_element(), GC_child_CisAisCjtAku_element(), GC_child_CisAjtCkuAku_element(), GC_child_CisAjtCkuAlv_element(), I(), LargeList::i_max, LargeList::is1_spin, LargeList::is2_spin, LargeList::is3_spin, LargeList::is4_spin, LargeList::isA_spin, LargeList::isB_spin, BindStruct::Large, MultiVecProdMPI(), DefineList::NCisAjtCkuAlvDC, DefineList::OrgTpow, X_GC_child_CisAisCjtAjt_Hubbard_MPI(), X_GC_child_CisAisCjtAku_Hubbard_MPI(), X_GC_child_CisAjtCkuAku_Hubbard_MPI(), X_GC_child_CisAjtCkuAlv_Hubbard_MPI(), and zclear().

Referenced by expec_cisajscktaltdc().

  {
  long int i, j;
  long int isite1, isite2, isite3, isite4;
  long int org_isite1, org_isite2, org_isite3, org_isite4;
  long int org_sigma1, org_sigma2, org_sigma3, org_sigma4;
  long int Asum, Bsum, Adiff, Bdiff;
  long int tmp_off = 0;
  long int tmp_off_2 = 0;
  std::complex<double> tmp_V = 1.0 + 0.0*I;
  long int i_max;

  for (i = 0; i < X->Def.NCisAjtCkuAlvDC; i++) {
    zclear(X->Large.i_max*nstate, &Xvec[1][0]);
    org_isite1 = X->Def.CisAjtCkuAlvDC[i][0] + 1;
    org_sigma1 = X->Def.CisAjtCkuAlvDC[i][1];
    org_isite2 = X->Def.CisAjtCkuAlvDC[i][2] + 1;
    org_sigma2 = X->Def.CisAjtCkuAlvDC[i][3];
    org_isite3 = X->Def.CisAjtCkuAlvDC[i][4] + 1;
    org_sigma3 = X->Def.CisAjtCkuAlvDC[i][5];
    org_isite4 = X->Def.CisAjtCkuAlvDC[i][6] + 1;
    org_sigma4 = X->Def.CisAjtCkuAlvDC[i][7];

    if (CheckPE(org_isite1 - 1, X) == TRUE || CheckPE(org_isite2 - 1, X) == TRUE ||
      CheckPE(org_isite3 - 1, X) == TRUE || CheckPE(org_isite4 - 1, X) == TRUE) {
      isite1 = X->Def.OrgTpow[2 * org_isite1 - 2 + org_sigma1];
      isite2 = X->Def.OrgTpow[2 * org_isite2 - 2 + org_sigma2];
      isite3 = X->Def.OrgTpow[2 * org_isite3 - 2 + org_sigma3];
      isite4 = X->Def.OrgTpow[2 * org_isite4 - 2 + org_sigma4];
      if (isite1 == isite2 && isite3 == isite4) {
        X_GC_child_CisAisCjtAjt_Hubbard_MPI(org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3,
          1.0, X, nstate, Xvec, vec);
      }
      else if (isite1 == isite2 && isite3 != isite4) {
        X_GC_child_CisAisCjtAku_Hubbard_MPI(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, org_isite4 - 1, org_sigma4,
          1.0, X, nstate, Xvec, vec);
      }
      else if (isite1 != isite2 && isite3 == isite4) {
        X_GC_child_CisAjtCkuAku_Hubbard_MPI(org_isite1 - 1, org_sigma1, org_isite2 - 1, org_sigma2,
          org_isite3 - 1, org_sigma3, 1.0, X, nstate, Xvec, vec);
      }
      else if (isite1 != isite2 && isite3 != isite4) {
        X_GC_child_CisAjtCkuAlv_Hubbard_MPI(org_isite1 - 1, org_sigma1, org_isite2 - 1, org_sigma2,
          org_isite3 - 1, org_sigma3, org_isite4 - 1, org_sigma4, 1.0, X, nstate, Xvec, vec);
      }
    }//InterPE
    else {
      child_general_int_GetInfo(X, org_isite1, org_isite2, org_isite3, org_isite4,
        org_sigma1, org_sigma2, org_sigma3, org_sigma4, tmp_V);

      i_max = X->Large.i_max;
      isite1 = X->Large.is1_spin;
      isite2 = X->Large.is2_spin;
      Asum = X->Large.isA_spin;
      Adiff = X->Large.A_spin;

      isite3 = X->Large.is3_spin;
      isite4 = X->Large.is4_spin;
      Bsum = X->Large.isB_spin;
      Bdiff = X->Large.B_spin;

      if (isite1 == isite2 && isite3 == isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate)     \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2,tmp_V)
        for (j = 1; j <= i_max; j++) {
          GC_child_CisAisCisAis_element(j, isite1, isite3, tmp_V, nstate, Xvec, vec);
        }
      }
      else if (isite1 == isite2 && isite3 != isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2,tmp_V)
        for (j = 1; j <= i_max; j++) {
          GC_child_CisAisCjtAku_element(j, isite1, isite3, isite4, Bsum, Bdiff, 
            tmp_V, nstate, Xvec, vec, &tmp_off);
        }
      }
      else if (isite1 != isite2 && isite3 == isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2,tmp_V) 
        for (j = 1; j <= i_max; j++) {
          GC_child_CisAjtCkuAku_element(j, isite1, isite2, isite3, Asum, Adiff, 
            tmp_V, nstate, Xvec, vec, &tmp_off);
        }
      }
      else if (isite1 != isite2 && isite3 != isite4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isite1,isite2,isite4,isite3,Asum,Bsum,Adiff,Bdiff,tmp_off,tmp_off_2,tmp_V) 
        for (j = 1; j <= i_max; j++) {
          GC_child_CisAjtCkuAlv_element(j, isite1, isite2, isite3, isite4, Asum, Adiff, Bsum, Bdiff, 
            tmp_V, nstate, Xvec, vec, &tmp_off_2);
        }
      }
    }
    MultiVecProdMPI(i_max, nstate, vec, Xvec, prod[i]);
  }//Intra PE
  return 0;
}

expec_cisajscktalt_Spin()

int expec_cisajscktalt_Spin	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Parent function to calculate two-body green's functions for Spin model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 884 of file expec_cisajscktaltdc.cpp.

References BindStruct::Def, expec_cisajscktalt_SpinGeneral(), expec_cisajscktalt_SpinHalf(), and DefineList::iFlgGeneralSpin.

Referenced by expec_cisajscktaltdc().

  {
  int info = 0;
  if (X->Def.iFlgGeneralSpin == FALSE) {
    info = expec_cisajscktalt_SpinHalf(X, nstate, Xvec, vec, prod);
  }
  else {
    info = expec_cisajscktalt_SpinGeneral(X, nstate, Xvec, vec, prod);
  }
  return info;
}

expec_cisajscktalt_SpinGC()

int expec_cisajscktalt_SpinGC	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Parent function to calculate two-body green's functions for Spin GC model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 910 of file expec_cisajscktaltdc.cpp.

References BindStruct::Def, expec_cisajscktalt_SpinGCGeneral(), expec_cisajscktalt_SpinGCHalf(), expec_cisajscktaltdc(), and DefineList::iFlgGeneralSpin.

Referenced by expec_cisajscktaltdc().

  {
  int info = 0;
  if (X->Def.iFlgGeneralSpin == FALSE) {
    info = expec_cisajscktalt_SpinGCHalf(X, nstate, Xvec, vec, prod);
  }
  else {
    info = expec_cisajscktalt_SpinGCGeneral(X, nstate, Xvec, vec, prod);
  }
  return info;
}

expec_cisajscktalt_SpinGCGeneral()

int expec_cisajscktalt_SpinGCGeneral	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Child function to calculate two-body green's functions for General Spin GC model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 742 of file expec_cisajscktaltdc.cpp.

References BitCheckGeneral(), BindStruct::Check, BindStruct::Def, GetOffCompGeneralSpin(), CheckList::idim_max, BindStruct::Large, LargeList::mode, MultiVecProdMPI(), DefineList::NCisAjtCkuAlvDC, DefineList::Nsite, Rearray_Interactions(), DefineList::SiteToBit, DefineList::Tpow, X_GC_child_CisAisCjuAju_GeneralSpin_MPIdouble(), X_GC_child_CisAisCjuAju_GeneralSpin_MPIsingle(), X_GC_child_CisAisCjuAjv_GeneralSpin_MPIdouble(), X_GC_child_CisAisCjuAjv_GeneralSpin_MPIsingle(), X_GC_child_CisAitCjuAju_GeneralSpin_MPIdouble(), X_GC_child_CisAitCjuAju_GeneralSpin_MPIsingle(), X_GC_child_CisAitCjuAjv_GeneralSpin_MPIdouble(), X_GC_child_CisAitCjuAjv_GeneralSpin_MPIsingle(), and zclear().

Referenced by expec_cisajscktalt_SpinGC().

 {
  long int i, j;
  long int org_isite1, org_isite2, org_isite3, org_isite4;
  long int org_sigma1, org_sigma2, org_sigma3, org_sigma4;
  long int tmp_off = 0;
  long int tmp_off_2 = 0;
  int  num1, one = 1;
  std::complex<double> tmp_V;
  long int i_max;
  i_max = X->Check.idim_max;
  X->Large.mode = M_CORR;

  for(i=0;i<X->Def.NCisAjtCkuAlvDC;i++){
    zclear(i_max*nstate, &Xvec[1][0]);

    if (Rearray_Interactions(i, &org_isite1, &org_isite2, &org_isite3, &org_isite4, &org_sigma1, &org_sigma2, &org_sigma3, &org_sigma4, &tmp_V, X) != 0) {
      //error message will be added
      zclear(nstate, prod[i]);
      continue;
    }

    if (org_isite1 > X->Def.Nsite && org_isite3 > X->Def.Nsite) {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        X_GC_child_CisAisCjuAju_GeneralSpin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 == org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAisCjuAjv_GeneralSpin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 == org_sigma4) {
        X_GC_child_CisAitCjuAju_GeneralSpin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAitCjuAjv_GeneralSpin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4,
          tmp_V, X, nstate, Xvec, vec);
      }
    }
    else if (org_isite3 > X->Def.Nsite || org_isite1 > X->Def.Nsite) {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        X_GC_child_CisAisCjuAju_GeneralSpin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 == org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAisCjuAjv_GeneralSpin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 == org_sigma4) {
        X_GC_child_CisAitCjuAju_GeneralSpin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAitCjuAjv_GeneralSpin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
    }
    else {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
#pragma omp parallel for default(none) private(j, num1) shared(vec,Xvec,nstate,one) \
firstprivate(i_max,X,org_isite1, org_sigma1,org_isite3, org_sigma3, tmp_V)
        for (j = 1; j <= i_max; j++) {
          num1 = BitCheckGeneral(j - 1, org_isite1, org_sigma1, X->Def.SiteToBit, X->Def.Tpow);
          if (num1 != FALSE) {
            num1 = BitCheckGeneral(j - 1, org_isite3, org_sigma3, X->Def.SiteToBit, X->Def.Tpow);
            if (num1 != FALSE) {
              zaxpy_(&nstate, &tmp_V, &vec[j][0], &one, &Xvec[j][0], &one);
            }
          }
        }
      }
      else if (org_sigma1 == org_sigma2 && org_sigma3 != org_sigma4) {
#pragma omp parallel for default(none) private(j, num1) shared(vec,Xvec,nstate,one) \
firstprivate(i_max,X, org_isite1, org_isite3, org_sigma1,org_sigma3,org_sigma4, tmp_off, tmp_V)
        for (j = 1; j <= i_max; j++) {
          num1 = GetOffCompGeneralSpin(j - 1, org_isite3, org_sigma4, org_sigma3,
            &tmp_off, X->Def.SiteToBit, X->Def.Tpow);
          if (num1 != FALSE) {
            num1 = BitCheckGeneral(tmp_off, org_isite1, org_sigma1, X->Def.SiteToBit, X->Def.Tpow);
            if (num1 != FALSE) {
              zaxpy_(&nstate, &tmp_V, &vec[j][0], &one, &Xvec[tmp_off + 1][0], &one);
            }
          }
        }
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 == org_sigma4) {
#pragma omp parallel for default(none) private(j, num1) shared(vec,Xvec,nstate,one) \
firstprivate(i_max,X, org_isite1, org_isite3, org_sigma1,org_sigma2, org_sigma3, tmp_off, tmp_V)
        for (j = 1; j <= i_max; j++) {
          num1 = BitCheckGeneral(j - 1, org_isite3, org_sigma3, X->Def.SiteToBit, X->Def.Tpow);
          if (num1 != FALSE) {
            num1 = GetOffCompGeneralSpin(j - 1, org_isite1, org_sigma2, org_sigma1,
              &tmp_off, X->Def.SiteToBit, X->Def.Tpow);
            if (num1 != FALSE) {
              zaxpy_(&nstate, &tmp_V, &vec[j][0], &one, &Xvec[tmp_off + 1][0], &one);
            }
          }
        }
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
#pragma omp parallel for default(none) private(j, num1) \
firstprivate(i_max,X,org_isite1,org_isite3,org_sigma1,org_sigma2,org_sigma3,org_sigma4,tmp_off,tmp_off_2,tmp_V) \
  shared(vec,Xvec,nstate,one)
        for (j = 1; j <= i_max; j++) {
          num1 = GetOffCompGeneralSpin(j - 1, org_isite3, org_sigma4, org_sigma3,
            &tmp_off, X->Def.SiteToBit, X->Def.Tpow);
          if (num1 != FALSE) {
            num1 = GetOffCompGeneralSpin(tmp_off, org_isite1, org_sigma2, org_sigma1, 
              &tmp_off_2, X->Def.SiteToBit, X->Def.Tpow);
            if (num1 != FALSE) {
              zaxpy_(&nstate, &tmp_V, &vec[j][0], &one, &Xvec[tmp_off_2 + 1][0], &one);
            }
          }
        }
      }
    }
    MultiVecProdMPI(i_max, nstate, vec, Xvec, prod[i]);
  }
  return 0;
}

expec_cisajscktalt_SpinGCHalf()

int expec_cisajscktalt_SpinGCHalf	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Child function to calculate two-body green's functions for 1/2 Spin GC model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 616 of file expec_cisajscktaltdc.cpp.

References BindStruct::Check, BindStruct::Def, GC_child_CisAisCisAis_spin_element(), GC_child_CisAisCitAiu_spin_element(), GC_child_CisAitCiuAiu_spin_element(), GC_child_CisAitCiuAiv_spin_element(), CheckList::idim_max, MultiVecProdMPI(), DefineList::NCisAjtCkuAlvDC, DefineList::Nsite, Rearray_Interactions(), DefineList::Tpow, X_GC_child_CisAis_spin_MPIdouble(), X_GC_child_CisAisCjuAju_spin_MPIdouble(), X_GC_child_CisAisCjuAju_spin_MPIsingle(), X_GC_child_CisAisCjuAjv_spin_MPIdouble(), X_GC_child_CisAisCjuAjv_spin_MPIsingle(), X_GC_child_CisAitCiuAiv_spin_MPIdouble(), X_GC_child_CisAitCiuAiv_spin_MPIsingle(), X_GC_child_CisAitCjuAju_spin_MPIdouble(), X_GC_child_CisAitCjuAju_spin_MPIsingle(), and zclear().

Referenced by expec_cisajscktalt_SpinGC().

 {
  long int i, j;
  long int org_isite1, org_isite2, org_isite3, org_isite4;
  long int org_sigma1, org_sigma2, org_sigma3, org_sigma4;
  long int isA_up, isB_up;
  long int tmp_off = 0;
  std::complex<double> tmp_V;
  long int i_max;
  i_max = X->Check.idim_max;

  for (i = 0; i < X->Def.NCisAjtCkuAlvDC; i++) {
    zclear(i_max*nstate, &Xvec[1][0]);

    if (Rearray_Interactions(i, &org_isite1, &org_isite2, &org_isite3, &org_isite4, &org_sigma1, &org_sigma2, &org_sigma3, &org_sigma4, &tmp_V, X) != 0) {
      //error message will be added
      zclear(nstate, prod[i]);
      continue;
    }

    if (org_isite1 > X->Def.Nsite && org_isite3 > X->Def.Nsite) { //org_isite3 >= org_isite1 > Nsite

      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        X_GC_child_CisAisCjuAju_spin_MPIdouble(
          org_isite1 - 1, org_sigma1, (org_isite3 - 1), org_sigma3, tmp_V, X, nstate, Xvec, vec);

      }
      else if (org_isite1 == org_isite3 && org_sigma1 == org_sigma4 && org_sigma2 == org_sigma3) { //diagonal (for spin: cuadcdau=cuau)
        X_GC_child_CisAis_spin_MPIdouble(
          org_isite1 - 1, org_sigma1, tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 == org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAisCjuAjv_spin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 == org_sigma4) {
        X_GC_child_CisAitCjuAju_spin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3,
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAitCiuAiv_spin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
    }
    else if (org_isite3 > X->Def.Nsite || org_isite1 > X->Def.Nsite) { //org_isite3 > Nsite >= org_isite1
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        X_GC_child_CisAisCjuAju_spin_MPIsingle(
          org_isite1 - 1, org_sigma1, (org_isite3 - 1), org_sigma3, tmp_V, X, nstate, Xvec, vec);

      }
      else if (org_sigma1 == org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAisCjuAjv_spin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, org_sigma4,
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 == org_sigma4) {
        X_GC_child_CisAitCjuAju_spin_MPIsingle(
          org_isite1 - 1, org_sigma2, org_isite3 - 1, org_sigma3, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
        X_GC_child_CisAitCiuAiv_spin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4,
          tmp_V, X, nstate, Xvec, vec);
      }
    }
    else {
      if (org_isite1 == org_isite2 && org_isite3 == org_isite4) {
        isA_up = X->Def.Tpow[org_isite2 - 1];
        isB_up = X->Def.Tpow[org_isite4 - 1];
        if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isA_up,isB_up,org_sigma2,org_sigma4,tmp_off,tmp_V)
          for (j = 1; j <= i_max; j++) {
            GC_child_CisAisCisAis_spin_element(j, isA_up, isB_up, org_sigma2, org_sigma4,
              tmp_V, nstate, Xvec, vec);
          }
        }
        else if (org_sigma1 == org_sigma2 && org_sigma3 != org_sigma4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isA_up,isB_up,org_sigma2,org_sigma4,tmp_off,tmp_V)
          for (j = 1; j <= i_max; j++) {
            GC_child_CisAisCitAiu_spin_element(j, org_sigma2, org_sigma4, isA_up, isB_up,
              tmp_V, nstate, Xvec, vec, &tmp_off);
          }
        }
        else if (org_sigma1 != org_sigma2 && org_sigma3 == org_sigma4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isA_up,isB_up,org_sigma2,org_sigma4,tmp_off,tmp_V)
          for (j = 1; j <= i_max; j++) {
            GC_child_CisAitCiuAiu_spin_element(j, org_sigma2, org_sigma4, isA_up, isB_up, 
              tmp_V, nstate, Xvec, vec, &tmp_off);
          }
        }
        else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isA_up,isB_up,org_sigma2,org_sigma4,tmp_off,tmp_V)
          for (j = 1; j <= i_max; j++) {
            GC_child_CisAitCiuAiv_spin_element(j, org_sigma2, org_sigma4, isA_up, isB_up,
              tmp_V, nstate, Xvec, vec, &tmp_off);
          }
        }
      }
    }
    MultiVecProdMPI(i_max, nstate, vec, Xvec, prod[i]);
  }
  return 0;
}

expec_cisajscktalt_SpinGeneral()

int expec_cisajscktalt_SpinGeneral	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Child function to calculate two-body green's functions for General Spin model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 503 of file expec_cisajscktaltdc.cpp.

References BitCheckGeneral(), BindStruct::Check, DefineList::CisAjtCkuAlvDC, ConvertToList1GeneralSpin(), BindStruct::Def, GetLocal2Sz(), GetOffCompGeneralSpin(), CheckList::idim_max, BindStruct::Large, list_1, LargeList::mode, MultiVecProdMPI(), DefineList::NCisAjtCkuAlvDC, DefineList::Nsite, Rearray_Interactions(), CheckList::sdim, DefineList::SiteToBit, DefineList::Tpow, X_child_CisAisCjuAju_GeneralSpin_MPIdouble(), X_child_CisAisCjuAju_GeneralSpin_MPIsingle(), X_child_CisAitCjuAjv_GeneralSpin_MPIdouble(), X_child_CisAitCjuAjv_GeneralSpin_MPIsingle(), and zclear().

Referenced by expec_cisajscktalt_Spin().

 {
  long int i, j;
  long int org_isite1, org_isite2, org_isite3, org_isite4;
  long int org_sigma1, org_sigma2, org_sigma3, org_sigma4;
  long int tmp_off = 0;
  long int tmp_off_2 = 0;
  long int list1_off = 0;
  int num1, one = 1;
  std::complex<double> tmp_V;
  long int i_max;
  int tmp_Sz;
  long int tmp_org = 0;
  i_max = X->Check.idim_max;
  X->Large.mode = M_CORR;

  for (i = 0; i < X->Def.NCisAjtCkuAlvDC; i++) {
    zclear(i_max*nstate, &Xvec[1][0]);

    if (Rearray_Interactions(i, &org_isite1, &org_isite2, &org_isite3, &org_isite4, &org_sigma1, &org_sigma2, &org_sigma3, &org_sigma4, &tmp_V, X) != 0) {
      zclear(nstate, prod[i]);
      continue;
    }
    tmp_Sz = 0;

    for (j = 0; j < 2; j++) {
      tmp_org = X->Def.CisAjtCkuAlvDC[i][4 * j + 1] * X->Def.Tpow[X->Def.CisAjtCkuAlvDC[i][4 * j]];
      tmp_Sz += GetLocal2Sz(X->Def.CisAjtCkuAlvDC[i][4 * j] + 1, tmp_org, X->Def.SiteToBit, X->Def.Tpow);
      tmp_org = X->Def.CisAjtCkuAlvDC[i][4 * j + 3] * X->Def.Tpow[X->Def.CisAjtCkuAlvDC[i][4 * j + 2]];
      tmp_Sz -= GetLocal2Sz(X->Def.CisAjtCkuAlvDC[i][4 * j + 2] + 1, tmp_org, X->Def.SiteToBit, X->Def.Tpow);
    }
    if (tmp_Sz != 0) { // not Sz conserved
      zclear(nstate, prod[i]);
      continue;
    }

    if (org_isite1 > X->Def.Nsite && org_isite3 > X->Def.Nsite) {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        X_child_CisAisCjuAju_GeneralSpin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
        X_child_CisAitCjuAjv_GeneralSpin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4,
          tmp_V, X, nstate, Xvec, vec);
      }
      else {
      }
    }
    else if (org_isite3 > X->Def.Nsite || org_isite1 > X->Def.Nsite) {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        X_child_CisAisCjuAju_GeneralSpin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_isite3 - 1, org_sigma3, tmp_V, X, nstate, Xvec, vec);
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
        X_child_CisAitCjuAjv_GeneralSpin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4,
          tmp_V, X, nstate, Xvec, vec);
      }
      else {
      }
    }
    else {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
#pragma omp parallel for default(none) private(j, num1) shared(vec,list_1,Xvec,nstate,one) \
firstprivate(i_max,X,org_isite1, org_sigma1,org_isite3, org_sigma3, tmp_V)
        for (j = 1; j <= i_max; j++) {
          num1 = BitCheckGeneral(list_1[j], org_isite1, org_sigma1, X->Def.SiteToBit, X->Def.Tpow);
          if (num1 != FALSE) {
            num1 = BitCheckGeneral(list_1[j], org_isite3, org_sigma3, X->Def.SiteToBit, X->Def.Tpow);
            if (num1 != FALSE) {
              zaxpy_(&nstate, &tmp_V, &vec[j][0], &one, &Xvec[j][0], &one);
            }
          }
        }
      }
      else if (org_sigma1 != org_sigma2 && org_sigma3 != org_sigma4) {
#pragma omp parallel for default(none) private(j,num1) \
firstprivate(i_max,X,org_isite1,org_isite3,org_sigma1,org_sigma2,org_sigma3,org_sigma4,tmp_off,tmp_off_2,list1_off,myrank,tmp_V) \
  shared(vec,list_1,Xvec,nstate,one)
        for (j = 1; j <= i_max; j++) {
          num1 = GetOffCompGeneralSpin(list_1[j], org_isite3, org_sigma4, org_sigma3, &tmp_off, X->Def.SiteToBit, X->Def.Tpow);
          if (num1 != FALSE) {
            num1 = GetOffCompGeneralSpin(tmp_off, org_isite1, org_sigma2, org_sigma1, &tmp_off_2,
              X->Def.SiteToBit, X->Def.Tpow);
            if (num1 != FALSE) {
              ConvertToList1GeneralSpin(tmp_off_2, X->Check.sdim, &list1_off);
              zaxpy_(&nstate, &tmp_V, &vec[j][0], &one, &Xvec[list1_off][0], &one);
            }
          }
        }
      }
    }
    MultiVecProdMPI(i_max, nstate, vec, Xvec, prod[i]);
  }
  return 0;
}

expec_cisajscktalt_SpinHalf()

int expec_cisajscktalt_SpinHalf	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec,
		std::complex< double > **	prod
	)

Child function to calculate two-body green's functions for 1/2 Spin model.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors
_fp	[in] output file name

Return values

0	normally finished
-1	abnormally finished

Definition at line 383 of file expec_cisajscktaltdc.cpp.

References BindStruct::Check, child_CisAisCisAis_spin_element(), BindStruct::Def, CheckList::idim_max, BindStruct::Large, LargeList::mode, MultiVecProdMPI(), myrank, DefineList::NCisAjtCkuAlvDC, DefineList::Nsite, Rearray_Interactions(), DefineList::Tpow, X_child_exchange_spin_element(), X_child_general_int_spin_MPIdouble(), X_child_general_int_spin_MPIsingle(), X_Spin_CisAis(), X_SpinGC_CisAis(), zaxpy_long(), and zclear().

Referenced by expec_cisajscktalt_Spin().

 {
  long int i, j;
  long int org_isite1, org_isite2, org_isite3, org_isite4;
  long int org_sigma1, org_sigma2, org_sigma3, org_sigma4;
  long int isA_up, isB_up;
  long int is1_up, is2_up;
  long int tmp_off = 0;
  int tmp_sgn, num1, num2, one = 1;
  std::complex<double> tmp_V;
  long int i_max;
  std::complex<double> dmv;

  i_max = X->Check.idim_max;
  X->Large.mode = M_CORR;

  for (i = 0; i < X->Def.NCisAjtCkuAlvDC; i++) {
    zclear(i_max*nstate, &Xvec[1][0]);
    if (Rearray_Interactions(i, &org_isite1, &org_isite2, &org_isite3, &org_isite4, &org_sigma1, &org_sigma2, &org_sigma3, &org_sigma4, &tmp_V, X) != 0) {
      //error message will be added
      zclear(nstate, prod[i]);
      continue;
    }

    if (org_isite1 > X->Def.Nsite && org_isite3 > X->Def.Nsite) {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        is1_up = X->Def.Tpow[org_isite1 - 1];
        is2_up = X->Def.Tpow[org_isite3 - 1];
        num1 = X_SpinGC_CisAis((long int)myrank + 1, is1_up, org_sigma1);
        num2 = X_SpinGC_CisAis((long int)myrank + 1, is2_up, org_sigma3);
        zaxpy_long(i_max*nstate, tmp_V * (std::complex<double>)(num1*num2),
          &vec[1][0], &Xvec[1][0]);
      }
      else if (org_isite1 == org_isite3 && org_sigma1 == org_sigma4 && org_sigma2 == org_sigma3) {
        is1_up = X->Def.Tpow[org_isite1 - 1];
        num1 = X_SpinGC_CisAis((long int)myrank + 1, is1_up, org_sigma1);
        zaxpy_long(i_max*nstate, tmp_V * (std::complex<double>)num1, &vec[1][0], &Xvec[1][0]);
      }
      else if (org_sigma1 == org_sigma4 && org_sigma2 == org_sigma3) {//exchange
        X_child_general_int_spin_MPIdouble(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else {  // other process is not allowed
                // error message will be added
      }
    }
    else if (org_isite1 > X->Def.Nsite || org_isite3 > X->Def.Nsite) {
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
        is1_up = X->Def.Tpow[org_isite1 - 1];
        is2_up = X->Def.Tpow[org_isite3 - 1];
        num2 = X_SpinGC_CisAis((long int)myrank + 1, is2_up, org_sigma3);
#pragma omp parallel for default(none)shared(vec,Xvec,nstate,one)      \
  firstprivate(i_max, tmp_V, is1_up, org_sigma1, X, num2) private(j, num1,dmv)
        for (j = 1; j <= i_max; j++) {
          num1 = X_Spin_CisAis(j, is1_up, org_sigma1);
          dmv = tmp_V * (std::complex<double>)(num1*num2);
          zaxpy_(&nstate, &dmv, &vec[j][0], &one, &Xvec[j][0], &one);
        }
      }
      else if (org_sigma1 == org_sigma4 && org_sigma2 == org_sigma3) {//exchange
        X_child_general_int_spin_MPIsingle(
          org_isite1 - 1, org_sigma1, org_sigma2, org_isite3 - 1, org_sigma3, org_sigma4, 
          tmp_V, X, nstate, Xvec, vec);
      }
      else {  // other process is not allowed
                // error message will be added
      }
    }
    else {
      isA_up = X->Def.Tpow[org_isite1 - 1];
      isB_up = X->Def.Tpow[org_isite3 - 1];
      if (org_sigma1 == org_sigma2 && org_sigma3 == org_sigma4) { //diagonal
#pragma omp parallel for default(none) private(j) shared(vec,Xvec,nstate) \
firstprivate(i_max,X,isA_up,isB_up,org_sigma2,org_sigma4,tmp_off, tmp_V)
        for (j = 1; j <= i_max; j++) {
          child_CisAisCisAis_spin_element(j, isA_up, isB_up, org_sigma2, org_sigma4,
            tmp_V, nstate, Xvec, vec);
        }
      }
      else if (org_isite1 == org_isite3 && org_sigma1 == org_sigma4 && org_sigma3 == org_sigma2) {
#pragma omp parallel for default(none) private(j, dmv) \
  firstprivate(i_max,X,isA_up,org_sigma1, tmp_V) shared(vec, list_1,Xvec,nstate,one)
        for (j = 1; j <= i_max; j++) {
          dmv = tmp_V * (std::complex<double>)X_Spin_CisAis(j, isA_up, org_sigma1);
          zaxpy_(&nstate, &dmv, &vec[j][0], &one, &Xvec[j][0], &one);
        }
      }
      else if (org_sigma1 == org_sigma4 && org_sigma2 == org_sigma3) { // exchange
#pragma omp parallel for default(none) private(j, tmp_sgn, dmv) shared(vec,Xvec,nstate,one) \
firstprivate(i_max,X,isA_up,isB_up,org_sigma2,org_sigma4,tmp_off,tmp_V)
        for (j = 1; j <= i_max; j++) {
          tmp_sgn = X_child_exchange_spin_element(j, X, isA_up, isB_up, org_sigma2, org_sigma4, &tmp_off);
          dmv = tmp_sgn;
          zaxpy_(&nstate, &dmv, &vec[j][0], &one, &Xvec[tmp_off][0], &one);
        }
      }
      else {  // other process is not allowed
                // error message will be added
      }
    }
    MultiVecProdMPI(i_max, nstate, vec, Xvec, prod[i]);
  }
  return 0;
}

expec_cisajscktaltdc()

int expec_cisajscktaltdc	(	struct BindStruct *	X,
		int	nstate,
		std::complex< double > **	Xvec,
		std::complex< double > **	vec
	)

Parent function to calculate two-body green's functions.

Parameters

X	[in] data list for calculation
vec	[in] eigenvectors

Return values

0	normally finished
-1	abnormally finished

Note

The origin of function's name cisajscktalt comes from c=creation, i=ith site, s=spin, a=annihiration, j=jth site and so on.

Version

0.2

add function to treat the case of general spin

Version

0.1

Author

Takahiro Misawa (The University of Tokyo)

Kazuyoshi Yoshimi (The University of Tokyo)

Definition at line 944 of file expec_cisajscktaltdc.cpp.

References DefineList::CDataFileHead, childfopenMPI(), DefineList::CisAjtCkuAlvDC, BindStruct::Def, expec_cisajscktalt_Hubbard(), expec_cisajscktalt_HubbardGC(), expec_cisajscktalt_Spin(), expec_cisajscktalt_SpinGC(), GetSplitBitByModel(), DefineList::iCalcModel, DefineList::iCalcType, DefineList::istep, BindStruct::Large, LargeList::mode, DefineList::NCisAjtCkuAlvDC, DefineList::Nsite, TimeKeeperWithRandAndStep(), and TimeKeeperWithStep().

Referenced by CalcByTEM(), CalcByTPQ(), expec_cisajscktalt_SpinGC(), and phys().

  {
  FILE *fp;
  char sdt[D_FileNameMax];
  long int irght, ilft, ihfbit, icaca;
  std::complex<double> **prod;
  //For TPQ
  int step = 0, rand_i = 0, istate;

  if (X->Def.NCisAjtCkuAlvDC < 1) return 0;
  X->Large.mode = M_CORR;

  if (GetSplitBitByModel(X->Def.Nsite, X->Def.iCalcModel, &irght, &ilft, &ihfbit) != 0) {
    return -1;
  }

  //Make File Name for output
  prod = cd_2d_allocate(X->Def.NCisAjtCkuAlvDC, nstate);
  switch (X->Def.iCalcType) {
  case TPQCalc:
    step = X->Def.istep;
    TimeKeeperWithRandAndStep(X, "%s_TimeKeeper.dat", "set %d step %d:expec_cisajscktaltdc finishes: %s", "a", 0, step);
    break;
  case TimeEvolution:
    step = X->Def.istep;
    TimeKeeperWithStep(X, "%s_TimeKeeper.dat", "step %d:expec_cisajscktaltdc finishes: %s", "a", step);
    break;
  case FullDiag:
  case CG:
    break;
  }

  switch (X->Def.iCalcModel) {
  case HubbardGC:
    if (expec_cisajscktalt_HubbardGC(X, nstate, Xvec, vec, prod) != 0) {
      return -1;
    }
    break;

  case KondoGC:
  case Hubbard:
  case Kondo:
    if (expec_cisajscktalt_Hubbard(X, nstate, Xvec, vec, prod) != 0) {
      return -1;
    }
    break;

  case Spin:
    if (expec_cisajscktalt_Spin(X, nstate, Xvec, vec, prod) != 0) {
      return -1;
    }
    break;

  case SpinGC:
    if (expec_cisajscktalt_SpinGC(X, nstate, Xvec, vec, prod) != 0) {
      return -1;
    }
    break;

  default:
    return -1;
  }

  for (istate = 0; istate < nstate; istate++) {
    switch (X->Def.iCalcType) {
    case TPQCalc:
      step = X->Def.istep;
      sprintf(sdt, "%s_cisajscktalt_set%dstep%d.dat", X->Def.CDataFileHead, istate, step);
      break;
    case TimeEvolution:
      step = X->Def.istep;
      sprintf(sdt, "%s_cisajscktalt_step%d.dat", X->Def.CDataFileHead, step);
      break;
    case FullDiag:
    case CG:
      sprintf(sdt, "%s_cisajscktalt_eigen%d.dat", X->Def.CDataFileHead, istate);
      break;
    }
    if (childfopenMPI(sdt, "w", &fp) == 0) {
      for (icaca = 0; icaca < X->Def.NCisAjtCkuAlvDC; icaca++) {
        fprintf(fp, " %4d %4d %4d %4d %4d %4d %4d %4d %.10lf %.10lf\n",
          X->Def.CisAjtCkuAlvDC[icaca][0], X->Def.CisAjtCkuAlvDC[icaca][1],
          X->Def.CisAjtCkuAlvDC[icaca][2], X->Def.CisAjtCkuAlvDC[icaca][3],
          X->Def.CisAjtCkuAlvDC[icaca][4], X->Def.CisAjtCkuAlvDC[icaca][5],
          X->Def.CisAjtCkuAlvDC[icaca][6], X->Def.CisAjtCkuAlvDC[icaca][7],
          real(prod[icaca][istate]), imag(prod[icaca][istate]));
      }
      fclose(fp);
    }
    else return -1;
  }/*for (istate = 0; istate < nstate; istate++)*/

  if (X->Def.iCalcType == TPQCalc) {
    TimeKeeperWithRandAndStep(X, "%s_TimeKeeper.dat", "set %d step %d:expec_cisajscktaltdc finishes: %s", "a", rand_i, step);
  }
  else if (X->Def.iCalcType == TimeEvolution) {
    TimeKeeperWithStep(X, "%s_TimeKeeper.dat", "step %d:expec_cisajscktaltdc finishes: %s", "a", step);
  }
  //[s] this part will be added
  /* For FullDiag, it is convinient to calculate the total spin for each vector.
     Such functions will be added
     if(X->Def.iCalcType==FullDiag){
     if(X->Def.iCalcModel==Spin){
     expec_cisajscktaltdc_alldiag_spin(X,vec);
     }else if(X->Def.iCalcModel==Hubbard || X->Def.iCalcModel==Kondo){
     expec_cisajscktaltdc_alldiag(X,vec);
     }else{//
     X->Phys.s2=0.0;
     }
     }
  */
  //[e]
  free_cd_2d_allocate(prod);
  return 0;
}

Rearray_Interactions()

int Rearray_Interactions	(	int	i,
		long int *	org_isite1,
		long int *	org_isite2,
		long int *	org_isite3,
		long int *	org_isite4,
		long int *	org_sigma1,
		long int *	org_sigma2,
		long int *	org_sigma3,
		long int *	org_sigma4,
		std::complex< double > *	tmp_V,
		struct BindStruct *	X
	)

Rearray interactions.

Parameters

i
org_isite1	a site number on the site 1.
org_isite2	a site number on the site 2.
org_isite3	a site number on the site 3.
org_isite4	a site number on the site 4.
org_sigma1	a spin index on the site 1.
org_sigma2	a spin index on the site 2.
org_sigma3	a spin index on the site 3.
org_sigma4	a spin index on the site 4.
tmp_V	a value of interaction
X	data list for calculation

Returns

0 normally finished

-1 unnormally finished

Definition at line 57 of file expec_cisajscktaltdc.cpp.

Referenced by expec_cisajscktalt_SpinGCGeneral(), expec_cisajscktalt_SpinGCHalf(), expec_cisajscktalt_SpinGeneral(), and expec_cisajscktalt_SpinHalf().

{
  long int tmp_org_isite1,tmp_org_isite2,tmp_org_isite3,tmp_org_isite4;
  long int tmp_org_sigma1,tmp_org_sigma2,tmp_org_sigma3,tmp_org_sigma4;

  tmp_org_isite1   = X->Def.CisAjtCkuAlvDC[i][0]+1;
  tmp_org_sigma1   = X->Def.CisAjtCkuAlvDC[i][1];
  tmp_org_isite2   = X->Def.CisAjtCkuAlvDC[i][2]+1;
  tmp_org_sigma2   = X->Def.CisAjtCkuAlvDC[i][3];
  tmp_org_isite3   = X->Def.CisAjtCkuAlvDC[i][4]+1;
  tmp_org_sigma3   = X->Def.CisAjtCkuAlvDC[i][5];
  tmp_org_isite4   = X->Def.CisAjtCkuAlvDC[i][6]+1;
  tmp_org_sigma4   = X->Def.CisAjtCkuAlvDC[i][7];

  if(tmp_org_isite1==tmp_org_isite2 && tmp_org_isite3==tmp_org_isite4){
    if(tmp_org_isite1 > tmp_org_isite3){
      *org_isite1   = tmp_org_isite3;
      *org_sigma1   = tmp_org_sigma3;
      *org_isite2   = tmp_org_isite4;
      *org_sigma2   = tmp_org_sigma4;
      *org_isite3   = tmp_org_isite1;
      *org_sigma3   = tmp_org_sigma1;
      *org_isite4   = tmp_org_isite2;
      *org_sigma4   = tmp_org_sigma2;
    }
    else{
      *org_isite1   = tmp_org_isite1;
      *org_sigma1   = tmp_org_sigma1;
      *org_isite2   = tmp_org_isite2;
      *org_sigma2   = tmp_org_sigma2;
      *org_isite3   = tmp_org_isite3;
      *org_sigma3   = tmp_org_sigma3;
      *org_isite4   = tmp_org_isite4;
      *org_sigma4   = tmp_org_sigma4;
    }
    *tmp_V = 1.0;

  }
  else if(tmp_org_isite1==tmp_org_isite4 && tmp_org_isite3==tmp_org_isite2){
    if(tmp_org_isite1 > tmp_org_isite3){
      *org_isite1   = tmp_org_isite3;
      *org_sigma1   = tmp_org_sigma3;
      *org_isite2   = tmp_org_isite2;
      *org_sigma2   = tmp_org_sigma2;
      *org_isite3   = tmp_org_isite1;
      *org_sigma3   = tmp_org_sigma1;
      *org_isite4   = tmp_org_isite4;
      *org_sigma4   = tmp_org_sigma4;
    }
    else{
      *org_isite1   = tmp_org_isite1;
      *org_sigma1   = tmp_org_sigma1;
      *org_isite2   = tmp_org_isite4;
      *org_sigma2   = tmp_org_sigma4;
      *org_isite3   = tmp_org_isite3;
      *org_sigma3   = tmp_org_sigma3;
      *org_isite4   = tmp_org_isite2;
      *org_sigma4   = tmp_org_sigma2;
    }
    *tmp_V =-1.0;
  }
  else{
    return -1;
  }
  return 0;
}