StdFace_main.c

Go to the documentation of this file.

/*
HPhi-mVMC-StdFace - Common input generator
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 <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <ctype.h>
#include <math.h>
#include "StdFace_vals.h"
#include "StdFace_ModelUtil.h"
#include <complex.h>

#if defined(_HPhi)

static void StdFace_LargeValue(struct StdIntList *StdI) {
  int ktrans, kintr;
  double LargeValue0;

  LargeValue0 = 0.0;
  for (ktrans = 0; ktrans < StdI->ntrans; ktrans++) {
    LargeValue0 += cabs(StdI->trans[ktrans]);
  }
  for (kintr = 0; kintr < StdI->nintr; kintr++) {
    LargeValue0 += cabs(StdI->intr[kintr]);
  }
  for (kintr = 0; kintr < StdI->NCintra; kintr++) {
    LargeValue0 += fabs(StdI->Cintra[kintr]);
  }
  for (kintr = 0; kintr < StdI->NCinter; kintr++) {
    LargeValue0 += fabs(StdI->Cinter[kintr]);
  }
  for (kintr = 0; kintr < StdI->NEx; kintr++) {
    LargeValue0 += 2.0 * fabs(StdI->Ex[kintr]);
  }
  for (kintr = 0; kintr < StdI->NPairLift; kintr++) {
    LargeValue0 += 2.0 * fabs(StdI->PairLift[kintr]);
  }
  for (kintr = 0; kintr < StdI->NHund; kintr++) {
    LargeValue0 += 2.0 * fabs(StdI->Hund[kintr]);
  }
  LargeValue0 /= (double)StdI->nsite;
  StdFace_PrintVal_d("LargeValue", &StdI->LargeValue, LargeValue0);
}/*static void StdFace_LargeValue*/
static void PrintCalcMod(struct StdIntList *StdI)
{
  FILE *fp;
  int iCalcType, iCalcModel, iRestart, iCalcSpec, 
    iCalcEigenvec, iInitialVecTpye, InputEigenVec, OutputEigenVec,
    iInputHam, iOutputHam, iOutputExVec;
  /*
  First, check all parameters and exit if invalid parameters
  */
  fprintf(stdout, "\n  @ CalcMod\n\n");
  /*
   Method
  */
  iCalcEigenvec = 0;
  if (strcmp(StdI->method, "****") == 0){
    fprintf(stdout, "ERROR ! Method is NOT specified !\n");
    StdFace_exit(-1);
  }
  else if (strcmp(StdI->method, "lanczos") == 0) iCalcType = 0;
  else if (strcmp(StdI->method, "lanczosenergy") == 0) { 
    iCalcType = 0; 
    iCalcEigenvec = 1;
  }
  else if (strcmp(StdI->method, "tpq") == 0) iCalcType = 1;
  else if (strcmp(StdI->method, "fulldiag") == 0 ) iCalcType = 2;
  else if (strcmp(StdI->method, "cg") == 0) iCalcType = 3;
  else if (strcmp(StdI->method, "timeevolution") == 0) iCalcType = 4;
  else{
    fprintf(stdout, "\n ERROR ! Unsupported Solver : %s\n", StdI->method);
    StdFace_exit(-1);
  }/*if (strcmp(StdI->method, METHODS) != 0*/
  if (iCalcType != 4) StdI->PumpBody = 0;
  /*
   Model
  */
  if (strcmp(StdI->model, "hubbard") == 0) {
    if (StdI->lGC == 0)iCalcModel = 0;
    else iCalcModel = 3;
  }/*if (strcmp(StdI->model, "hubbard") == 0)*/
  else if (strcmp(StdI->model, "spin") == 0) {
    if (StdI->lGC == 0)iCalcModel = 1;
    else iCalcModel = 4;
  }/*if (strcmp(StdI->model, "spin") == 0)*/
  else if (strcmp(StdI->model, "kondo") == 0) {
    if (StdI->lGC == 0)iCalcModel = 2;
    else iCalcModel = 5;
  }/*if (strcmp(StdI->model, "kondo") == 0)*/
  /*
  Restart
  */
  if (strcmp(StdI->Restart, "****") == 0) {
    strcpy(StdI->Restart, "none\0");
    fprintf(stdout, "          Restart = none        ######  DEFAULT VALUE IS USED  ######\n");
    iRestart = 0;
  }/*if (strcmp(StdI->Restart, "****") == 0)*/
  else {
    fprintf(stdout, "          Restart = %s\n", StdI->Restart);
    if (strcmp(StdI->Restart, "none") == 0) iRestart = 0;
    else if (strcmp(StdI->Restart, "restart_out") == 0 ||
             strcmp(StdI->Restart, "save") == 0) iRestart = 1;
    else if (strcmp(StdI->Restart, "restartsave") == 0 ||
             strcmp(StdI->Restart, "restart")     == 0) iRestart = 2;
    else if (strcmp(StdI->Restart, "restart_in") == 0) iRestart = 3;
    else {
      fprintf(stdout, "\n ERROR ! Restart Mode : %s\n", StdI->Restart);
      StdFace_exit(-1);
    }
  }/*if (strcmp(StdI->Restart, "****") != 0)*/
  /*
  InitialVecType
  */
  if (strcmp(StdI->InitialVecType, "****") == 0) {
    strcpy(StdI->InitialVecType, "c\0");
    fprintf(stdout, "   InitialVecType = c           ######  DEFAULT VALUE IS USED  ######\n");
    iInitialVecTpye = 0;
  }/*if (strcmp(StdI->InitialVecType, "****") == 0)*/
  else {
    fprintf(stdout, "   InitialVecType = %s\n", StdI->InitialVecType);
    if (strcmp(StdI->InitialVecType, "c") == 0) iInitialVecTpye = 0;
    else if (strcmp(StdI->InitialVecType, "r") == 0) iInitialVecTpye = 1;
    else {
      fprintf(stdout, "\n ERROR ! Restart Mode : %s\n", StdI->Restart);
      StdFace_exit(-1);
    }
  }/*if (strcmp(StdI->InitialVecType, "****") != 0)*/
  /*
  EigenVecIO
  */
  InputEigenVec = 0;
  OutputEigenVec = 0;
  if (strcmp(StdI->EigenVecIO, "****") == 0) {
    strcpy(StdI->EigenVecIO, "none\0");
    fprintf(stdout, "       EigenVecIO = none        ######  DEFAULT VALUE IS USED  ######\n");
  }/*if (strcmp(StdI->EigenVecIO, "****") == 0)*/
  else {
    fprintf(stdout, "       EigenVecIO = %s\n", StdI->EigenVecIO);
    if (strcmp(StdI->EigenVecIO, "none") == 0) InputEigenVec = 0;
    else if (strcmp(StdI->EigenVecIO, "in") == 0) InputEigenVec = 1;
    else if (strcmp(StdI->EigenVecIO, "out") == 0) OutputEigenVec = 1;
    else if (strcmp(StdI->EigenVecIO, "inout") == 0) {
      InputEigenVec = 1;
      OutputEigenVec = 1;
    }/*if (strcmp(StdI->EigenVecIO, "inout") == 0)*/
    else {
      fprintf(stdout, "\n ERROR ! EigenVecIO Mode : %s\n", StdI->Restart);
      StdFace_exit(-1);
    }
  }/*if (strcmp(StdI->EigenVecIO, "****") != 0)*/
  if (strcmp(StdI->method, "timeevolution") == 0) InputEigenVec = 1;
  /*
   * HamIO
   */
  iOutputHam = 0;
  iInputHam = 0;
  if (strcmp(StdI->HamIO, "****") == 0) {
    strcpy(StdI->HamIO, "none\0");
    fprintf(stdout, "         HamIO = none        ######  DEFAULT VALUE IS USED  ######\n");
  }/*if (strcmp(StdI->HamIO, "****") == 0)*/
  else {
    fprintf(stdout, "         HamIO = %s\n", StdI->HamIO);
    if (strcmp(StdI->HamIO, "none") == 0){ iOutputHam = 0; iInputHam=0;}
    else if (strcmp(StdI->HamIO, "out") == 0) iOutputHam = 1;
    else if (strcmp(StdI->HamIO, "in") == 0) iInputHam = 1;
    else {
      fprintf(stdout, "\n ERROR ! HamIO mode : %s\n", StdI->HamIO);
      StdFace_exit(-1);
    }
  }
  /*
  CalcSpec
  */
  if (strcmp(StdI->CalcSpec, "****") == 0) {
    strcpy(StdI->CalcSpec, "none\0");
    fprintf(stdout, "         CalcSpec = none        ######  DEFAULT VALUE IS USED  ######\n");
    iCalcSpec = 0;
  }/*if (strcmp(StdI->CalcSpec, "****") == 0)*/
  else {
    fprintf(stdout, "         CalcSpec = %s\n", StdI->CalcSpec);
    if (strcmp(StdI->CalcSpec, "none") == 0) iCalcSpec = 0;
    else if (strcmp(StdI->CalcSpec, "normal") == 0) iCalcSpec = 1;
    else if (strcmp(StdI->CalcSpec, "noiteration") == 0) iCalcSpec = 2;
    else if (strcmp(StdI->CalcSpec, "restart_out") == 0) iCalcSpec = 3;
    else if (strcmp(StdI->CalcSpec, "restart_in") == 0) iCalcSpec = 4;
    else if (strcmp(StdI->CalcSpec, "restartsave") == 0 ||
             strcmp(StdI->CalcSpec, "restart")     == 0) iCalcSpec = 5;
    else {
      fprintf(stdout, "\n ERROR ! CalcSpec : %s\n", StdI->CalcSpec);
      StdFace_exit(-1);
    }
  }/*if (strcmp(StdI->CalcSpec, "****") != 0)*/
  /*
  OutputExcitedVec
  */
  iOutputExVec = 0;
  if (strcmp(StdI->OutputExVec, "****") == 0) {
    strcpy(StdI->OutputExVec, "none\0");
    fprintf(stdout, "         OutputExcitedVec = none        ######  DEFAULT VALUE IS USED  ######\n");
  }/*if (strcmp(StdI->OutputExVec, "****") == 0)*/
  else {
    fprintf(stdout, "         OutputExcitedVec = %s\n", StdI->OutputExVec);
    if (strcmp(StdI->OutputExVec, "none") == 0) iOutputExVec = 0;
    else if (strcmp(StdI->OutputExVec, "out") == 0) iOutputExVec = 1;
    else {
      fprintf(stdout, "\n ERROR ! OutputExcitedVec : %s\n", StdI->OutputExVec);
      StdFace_exit(-1);
    }
  }/*if (strcmp(StdI->OutputExVec, "****") != 0)*/

  fp = fopen("calcmod.def", "w");
  fprintf(fp, "#CalcType = 0:Lanczos, 1:TPQCalc, 2:FullDiag, 3:CG, 4:Time-evolution\n");
  fprintf(fp, "#CalcModel = 0:Hubbard, 1:Spin, 2:Kondo, 3:HubbardGC, 4:SpinGC, 5:KondoGC\n");
  fprintf(fp, "#Restart = 0:None, 1:Save, 2:Restart&Save, 3:Restart\n");
  fprintf(fp, "#CalcSpec = 0:None, 1:Normal, 2:No H*Phi, 3:Save, 4:Restart, 5:Restart&Save\n");
  fprintf(fp, "CalcType %3d\n", iCalcType);
  fprintf(fp, "CalcModel %3d\n", iCalcModel);
  fprintf(fp, "ReStart %3d\n", iRestart);
  fprintf(fp, "CalcSpec %3d\n", iCalcSpec);
  fprintf(fp, "CalcEigenVec %3d\n", iCalcEigenvec);
  fprintf(fp, "InitialVecType %3d\n", iInitialVecTpye);
  fprintf(fp, "InputEigenVec %3d\n", InputEigenVec);
  fprintf(fp, "OutputEigenVec %3d\n", OutputEigenVec);
  fprintf(fp, "InputHam %3d\n", iInputHam);
  fprintf(fp, "OutputHam %3d\n", iOutputHam);
  fprintf(fp, "OutputExVec %3d\n", iOutputExVec);
  fflush(fp);
  fclose(fp);
  fprintf(stdout, "     calcmod.def is written.\n\n");
}/*static void PrintCalcMod*/
static void PrintExcitation(struct StdIntList *StdI) {
  FILE *fp;
  int NumOp, **spin, isite, ispin, icell, itau;
  double *coef, pi, Cphase, S, Sz;
  double *fourier_r, *fourier_i;

  if (strcmp(StdI->model, "spin") == 0 && StdI->S2 > 1) {
    coef = (double *)malloc(sizeof(double) * (StdI->S2 + 1));
    spin = (int **)malloc(sizeof(int*) * (StdI->S2 + 1));
    for (ispin = 0; ispin < StdI->S2 + 1; ispin++) spin[ispin] = (int *)malloc(sizeof(int) * 2);
  }
  else {
    coef = (double *)malloc(sizeof(double) * 2);
    spin = (int **)malloc(sizeof(int*) * 2);
    for (ispin = 0; ispin < 2; ispin++) spin[ispin] = (int *)malloc(sizeof(int) * 2);
  }

  fourier_r = (double *)malloc(sizeof(double) * StdI->nsite);
  fourier_i = (double *)malloc(sizeof(double) * StdI->nsite);
  
  fprintf(stdout, "\n  @ Spectrum\n\n");

  StdFace_PrintVal_d("SpectrumQW", &StdI->SpectrumQ[0], 0.0);
  StdFace_PrintVal_d("SpectrumQL", &StdI->SpectrumQ[1], 0.0);
  StdFace_PrintVal_d("SpectrumQH", &StdI->SpectrumQ[2], 0.0);

  if (strcmp(StdI->SpectrumType, "****") == 0) {
    strcpy(StdI->SpectrumType, "szsz\0");
    fprintf(stdout, "     SpectrumType = szsz        ######  DEFAULT VALUE IS USED  ######\n");
    if (strcmp(StdI->model, "spin") == 0) {
      NumOp = StdI->S2 + 1;
      for (ispin = 0; ispin <= StdI->S2; ispin++) {
        Sz = (double)ispin - (double)StdI->S2 * 0.5;
        coef[ispin] = Sz;
        spin[ispin][0] = ispin;
        spin[ispin][1] = ispin;
      }
    }
    else {
      NumOp = 2;
      coef[0] = 0.5;
      coef[1] = -0.5;
      spin[0][0] = 0;
      spin[0][1] = 0;
      spin[1][0] = 1;
      spin[1][1] = 1;
    }
    StdI->SpectrumBody = 2;
  }
  else {
    fprintf(stdout, "     SpectrumType = %s\n", StdI->SpectrumType);
    if (strcmp(StdI->SpectrumType, "szsz") == 0) {
      if (strcmp(StdI->model, "spin") == 0) {
        NumOp = StdI->S2 + 1;
        for (ispin = 0; ispin <= StdI->S2; ispin++) {
          Sz = (double)ispin - (double)StdI->S2 * 0.5;
          coef[ispin] = Sz;
          spin[ispin][0] = ispin;
          spin[ispin][1] = ispin;
        }
      }
      else {
        NumOp = 2;
        coef[0] = 0.5;
        coef[1] = -0.5;
        spin[0][0] = 0;
        spin[0][1] = 0;
        spin[1][0] = 1;
        spin[1][1] = 1;
      }
      StdI->SpectrumBody = 2;
    }
    else if (strcmp(StdI->SpectrumType, "s+s-") == 0) {
      if (strcmp(StdI->model, "spin") == 0 && StdI->S2 > 1) {
        NumOp = StdI->S2;
        S = (double)StdI->S2 * 0.5;
        for (ispin = 0; ispin < StdI->S2; ispin++) {
          Sz = (double)ispin - (double)StdI->S2 * 0.5;
          coef[ispin] = sqrt(S*(S + 1.0) - Sz*(Sz + 1.0));
          spin[ispin][0] = ispin;
          spin[ispin][1] = ispin + 1;
        }
      }
      else {
        NumOp = 1;
        coef[0] = 1.0;
        spin[0][0] = 0;
        spin[0][1] = 1;
      }
      StdI->SpectrumBody = 2;
    }
    else if (strcmp(StdI->SpectrumType, "density") == 0) {
      NumOp = 2;
      coef[0] = 1,0;
      coef[1] = 1.0;
      spin[0][0] = 0;
      spin[0][1] = 0;
      spin[1][0] = 1;
      spin[1][1] = 1;
      StdI->SpectrumBody = 2;
    }
    else if (strcmp(StdI->SpectrumType, "up") == 0) {
      NumOp = 1;
      coef[0] = 1.0;
      spin[0][0] = 0;
      StdI->SpectrumBody = 1;
    }
    else if (strcmp(StdI->SpectrumType, "down") == 0) {
      NumOp = 1;
      coef[0] = 1.0;
      spin[0][0] = 1;
      StdI->SpectrumBody = 1;
    }
    else {
      fprintf(stdout, "\n ERROR ! SpectrumType : %s\n", StdI->SpectrumType);
      StdFace_exit(-1);
    }
  }

  isite = 0;
  for (icell = 0; icell < StdI->NCell; icell++) {
    for (itau = 0; itau < StdI->NsiteUC; itau++) {
      Cphase = (StdI->Cell[icell][0] + StdI->tau[itau][0])*StdI->SpectrumQ[0]
             + (StdI->Cell[icell][1] + StdI->tau[itau][1])*StdI->SpectrumQ[1]
             + (StdI->Cell[icell][2] + StdI->tau[itau][2])*StdI->SpectrumQ[2];
      fourier_r[isite] = cos(2.0*StdI->pi*Cphase);
      fourier_i[isite] = sin(2.0*StdI->pi*Cphase);
      isite += 1;
    }
  }
  if (strcmp(StdI->model, "kondo") == 0) {
    for (isite = 0; isite < StdI->nsite / 2; isite++) {
      fourier_r[isite + StdI->nsite / 2] = fourier_r[isite];
      fourier_i[isite + StdI->nsite / 2] = fourier_i[isite];
    }/*for (isite = 0; isite < StdI->nsite; isite++)*/
  }/*if (strcmp(StdI->model, "kondo") == 0)*/

  if (StdI->SpectrumBody == 1) {
    fp = fopen("single.def", "w");
    fprintf(fp, "=============================================\n");
    if (strcmp(StdI->model, "kondo") == 0) {
      fprintf(fp, "NSingle %d\n", StdI->nsite / 2 * NumOp);
    }
    else {
      fprintf(fp, "NSingle %d\n", StdI->nsite * NumOp);
    }
    fprintf(fp, "=============================================\n");
    fprintf(fp, "============== Single Excitation ============\n");
    fprintf(fp, "=============================================\n");
    if (strcmp(StdI->model, "kondo") == 0) {
      for (isite = StdI->nsite / 2; isite < StdI->nsite; isite++) {
        fprintf(fp, "%d %d 0 %25.15f %25.15f\n", isite, spin[0][0],
          fourier_r[isite] * coef[0], fourier_i[isite] * coef[0]);
      }/*for (isite = 0; isite < StdI->nsite; isite++)*/
    }/*if (strcmp(StdI->model, "kondo") == 0)*/
    else {
      for (isite = 0; isite < StdI->nsite; isite++) {
        fprintf(fp, "%d %d 0 %25.15f %25.15f\n", isite, spin[0][0],
          fourier_r[isite] * coef[0], fourier_i[isite] * coef[0]);
      }/*for (isite = 0; isite < StdI->nsite; isite++)*/
    }
    fprintf(stdout, "      single.def is written.\n\n");
  }
  else {
    fp = fopen("pair.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NPair %d\n", StdI->nsite * NumOp);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "=============== Pair Excitation =============\n");
    fprintf(fp, "=============================================\n");
    for (isite = 0; isite < StdI->nsite; isite++) {
      for (ispin = 0; ispin < NumOp; ispin++) {
        fprintf(fp, "%d %d %d %d 1 %25.15f %25.15f\n", 
          isite, spin[ispin][0], isite, spin[ispin][1],
          fourier_r[isite] * coef[ispin], fourier_i[isite] * coef[ispin]);
      }
    }
    fprintf(stdout, "        pair.def is written.\n\n");
  }
  fflush(fp);
  fclose(fp);

  free(fourier_r);
  free(fourier_i);
  if (strcmp(StdI->model, "spin") == 0) 
    for (ispin = 0; ispin < StdI->S2 + 1; ispin++) free(spin[ispin]);
  else 
    for (ispin = 0; ispin < 2; ispin++) free(spin[ispin]);
  free(spin);
  free(coef);

}/*static void PrintExcitation()*/
/*
@brief Compute vectorpotential
*/
static void VectorPotential(struct StdIntList *StdI) {
  FILE *fp;
  int it, ii, isite, icell, itau, itrans, jsite, jcell, jtau, ntrans0;
  double Cphase, time, dR[3];
  double **Et;
  double complex coef;

  fprintf(stdout, "\n  @ Time-evolution\n\n");

  StdFace_PrintVal_d("VecPotW", &StdI->VecPot[0], 0.0);
  StdFace_PrintVal_d("VecPotL", &StdI->VecPot[1], 0.0);
  StdFace_PrintVal_d("VecPotH", &StdI->VecPot[2], 0.0);
  StdFace_PrintVal_i("Lanczos_max", &StdI->Lanczos_max, 1000);
  StdFace_PrintVal_d("dt", &StdI->dt, 0.1);
  StdFace_PrintVal_d("freq", &StdI->freq, 0.1);
  StdFace_PrintVal_d("tshift", &StdI->tshift, 0.0);
  StdFace_PrintVal_d("tdump", &StdI->tdump, 0.1);
  StdFace_PrintVal_d("Uquench", &StdI->Uquench, 0.0);
  StdFace_PrintVal_i("ExpandCoef", &StdI->ExpandCoef, 10);
  StdI->At = (double **)malloc(sizeof(double*) * StdI->Lanczos_max);
  Et = (double **)malloc(sizeof(double*) * StdI->Lanczos_max);
  for (it = 0; it < StdI->Lanczos_max; it++) {
    StdI->At[it] = (double *)malloc(sizeof(double) * 3);
    Et[it] = (double *)malloc(sizeof(double) * 3);
  }

  if (strcmp(StdI->PumpType, "****") == 0) {
    strcpy(StdI->PumpType, "quench\0");
    fprintf(stdout, "     PumpType = quench        ######  DEFAULT VALUE IS USED  ######\n");
    StdI->PumpBody = 2;
  }/*if (strcmp(StdI->PumpType, "****")*/
  else {
    fprintf(stdout, "     PumpType = %s\n", StdI->PumpType);
    if (strcmp(StdI->PumpType, "quench") == 0) {
      StdI->PumpBody = 2;
    }/*if (strcmp(StdI->PumpType, "quench")*/
    else if (strcmp(StdI->PumpType, "pulselaser") == 0) {
      for (it = 0; it < StdI->Lanczos_max; it++) {
        time = StdI->dt*(double)it;
        for (ii = 0; ii < 3; ii++) {
          StdI->At[it][ii] = StdI->VecPot[ii] * cos(StdI->freq*(time - StdI->tshift))
            * exp(-0.5* (time - StdI->tshift)*(time - StdI->tshift) / (StdI->tdump*StdI->tdump));
          Et[it][ii] = -StdI->VecPot[ii]
            * (
            (StdI->tshift - time) / (StdI->tdump*StdI->tdump) * cos(StdI->freq*(time - StdI->tshift))
              - StdI->freq* sin(StdI->freq*(time - StdI->tshift))
              )
            * exp(-0.5* (time - StdI->tshift)*(time - StdI->tshift) / (StdI->tdump*StdI->tdump));
        }
      }/*for (it = 0; it < StdI->Lanczos_max; it++)*/
      StdI->PumpBody = 1;
    }/*if (strcmp(StdI->PumpType, "pulselaser") == 0)*/
    else if (strcmp(StdI->PumpType, "aclaser") == 0) {
      for (it = 0; it < StdI->Lanczos_max; it++) {
        time = StdI->dt*(double)it;
        for (ii = 0; ii < 3; ii++) {
          StdI->At[it][ii] = StdI->VecPot[ii] * sin(StdI->freq*(time - StdI->tshift));
          Et[it][ii] = StdI->VecPot[ii] * cos(StdI->freq*(time - StdI->tshift)) * StdI->freq;
        }
      }/*for (it = 0; it < StdI->Lanczos_max; it++)*/
      StdI->PumpBody = 1;
    }/*if (strcmp(StdI->PumpType, "aclaser") == 0)*/
    else if (strcmp(StdI->PumpType, "dclaser") == 0) {
      for (it = 0; it < StdI->Lanczos_max; it++) {
        time = StdI->dt*(double)it;
        for (ii = 0; ii < 3; ii++) {
          StdI->At[it][ii] = StdI->VecPot[ii] * time;
          Et[it][ii] = -StdI->VecPot[ii];
        }
      }/*for (it = 0; it < StdI->Lanczos_max; it++)*/
      StdI->PumpBody = 1;
    }/* if (strcmp(StdI->PumpType, "dclaser") == 0)*/
    else {
      fprintf(stdout, "\n ERROR ! PumpType : %s\n", StdI->PumpType);
      StdFace_exit(-1);
    }
  }/*if (! strcmp(StdI->PumpType, "****"))*/

  if (StdI->PumpBody == 1) {
    fp = fopen("potential.dat", "w");
    fprintf(fp, "# Time A_W A_L A_H E_W E_L E_H\n");
    for (it = 0; it < StdI->Lanczos_max; it++) {
      time = StdI->dt*(double)it;
      fprintf(fp, "%f %f %f %f %f %f %f\n",
        time, StdI->At[it][0], StdI->At[it][1], StdI->At[it][2], Et[it][0], Et[it][1], Et[it][2]);
    }
    fflush(fp);
    fclose(fp);
  }/*if (StdI->PumpBody == 1)*/

  for (it = 0; it < StdI->Lanczos_max; it++) free(Et[it]);
  free(Et);
}/*static void VectorPotential(struct StdIntList *StdI)*/
static void PrintPump(struct StdIntList *StdI) {
  FILE *fp;
  int it, ii, isite, ipump, jpump, npump0;

  if (StdI->PumpBody == 1) {

    fp = fopen("teone.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "AllTimeStep %d\n", StdI->Lanczos_max);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "=========  OneBody Time Evolution  ==========\n");
    fprintf(fp, "=============================================\n");
    for (it = 0; it < StdI->Lanczos_max; it++) {
      /*
      Sum equivalent pumping
      */
      for (ipump = 0; ipump < StdI->npump[it]; ipump++) {
        for (jpump = ipump + 1; jpump < StdI->npump[it]; jpump++) {
          if (StdI->pumpindx[it][ipump][0] == StdI->pumpindx[it][jpump][0]
            && StdI->pumpindx[it][ipump][1] == StdI->pumpindx[it][jpump][1]
            && StdI->pumpindx[it][ipump][2] == StdI->pumpindx[it][jpump][2]
            && StdI->pumpindx[it][ipump][3] == StdI->pumpindx[it][jpump][3]) {
            StdI->pump[it][ipump] = StdI->pump[it][ipump] + StdI->pump[it][jpump];
            StdI->pump[it][jpump] = 0.0;
          }
        }/*for (ktrans = jtrans + 1; ktrans < StdI->ntrans; ktrans++)*/
      }/*for (jtrans = 0; jtrans < StdI->ntrans; jtrans++)*/
      /*
      Count the number of finite pumping
      */
      npump0 = 0;
      for (ipump = 0; ipump < StdI->npump[it]; ipump++) 
        if (cabs(StdI->pump[it][ipump]) > 0.000001) npump0 += 1;

      fprintf(fp, "%f  %d\n", StdI->dt*(double)it, npump0);
      for (ipump = 0; ipump < StdI->npump[it]; ipump++) {

        if (cabs(StdI->pump[it][ipump]) <= 0.000001) continue;

        fprintf(fp, "%5d %5d %5d %5d %25.15f %25.15f\n",
          StdI->pumpindx[it][ipump][0], StdI->pumpindx[it][ipump][1],
          StdI->pumpindx[it][ipump][2], StdI->pumpindx[it][ipump][3],
          creal(StdI->pump[it][ipump]), cimag(StdI->pump[it][ipump]));
      }/*for (itrans = 0; itrans < StdI->ntrans; itrans++)*/
    }/*for (it = 0; it < StdI->Lanczos_max; it++)*/
    fprintf(stdout, "      teone.def is written.\n\n");
  }
  else {
    fp = fopen("tetwo.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "AllTimeStep %d\n", StdI->Lanczos_max);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "========== TwoBody Time Evolution ===========\n");
    fprintf(fp, "=============================================\n");
    for (it = 0; it < StdI->Lanczos_max; it++) {
      fprintf(fp, "%f  %d\n", StdI->dt*(double)it, StdI->nsite);
      for (isite = 0; isite < StdI->nsite; isite++) {
        fprintf(fp, "%5d %5d %5d %5d %5d %5d %5d %5d %25.15f  %25.15f\n",
          isite, 0, isite, 0, isite, 1, isite, 1, StdI->Uquench, 0.0);
      }/*for (isite = 0; isite < StdI->nsite; isite++)*/
    }/*for (it = 0; it < StdI->Lanczos_max; it++)*/
    fprintf(stdout, "        tetwo.def is written.\n\n");
  }
  fflush(fp);
  fclose(fp);
}/*tatic void PrintPump*/
#elif defined(_mVMC)

static void PrintOrb(struct StdIntList *StdI) {
  FILE *fp;
  int isite, jsite, iOrb;

  fp = fopen("orbitalidx.def", "w");
  fprintf(fp, "=============================================\n");
  fprintf(fp, "NOrbitalIdx %10d\n", StdI->NOrb);
  fprintf(fp, "ComplexType %10d\n", StdI->ComplexType);
  fprintf(fp, "=============================================\n");
  fprintf(fp, "=============================================\n");

  for (isite = 0; isite < StdI->nsite; isite++) {
    for (jsite = 0; jsite < StdI->nsite; jsite++) {
      if (StdI->AntiPeriod[0] == 1 || StdI->AntiPeriod[1] == 1 || StdI->AntiPeriod[2] == 1) {
        fprintf(fp, "%5d  %5d  %5d  %5d\n", isite, jsite, StdI->Orb[isite][jsite], StdI->AntiOrb[isite][jsite]);
      }
      else {
        fprintf(fp, "%5d  %5d  %5d\n", isite, jsite, StdI->Orb[isite][jsite]);
      }
    }/*for (jsite = 0; jsite < isite; jsite++)*/
  }/*for (isite = 0; isite < StdI->nsite; isite++)*/

  for (iOrb = 0; iOrb < StdI->NOrb; iOrb++)
    fprintf(fp, "%5d  %5d\n", iOrb, 1);

  fflush(fp);
  fclose(fp);
  fprintf(stdout, "    orbitalidx.def is written.\n");

  for (isite = 0; isite < StdI->nsite; isite++) free(StdI->Orb[isite]);
  free(StdI->Orb);
}/*void PrintOrb*/
static void PrintOrbPara(struct StdIntList *StdI) {
  FILE *fp;
  int isite, jsite, NOrbGC, iOrbGC, isite1, jsite1, iorb;
  int **OrbGC, **AntiOrbGC;
  OrbGC = (int **)malloc(sizeof(int*) * StdI->nsite);
  AntiOrbGC = (int **)malloc(sizeof(int*) * StdI->nsite);
  for (isite = 0; isite < StdI->nsite; isite++) {
    OrbGC[isite] = (int *)malloc(sizeof(int) * StdI->nsite);
    AntiOrbGC[isite] = (int *)malloc(sizeof(int) * StdI->nsite);
    for (jsite = 0; jsite < StdI->nsite; jsite++) {
      OrbGC[isite][jsite] = StdI->Orb[isite][jsite];
      AntiOrbGC[isite][jsite] = StdI->AntiOrb[isite][jsite];
    }/*for (jsite = 0; jsite < isite; jsite++)*/
  }/*for (isite = 0; isite < StdI->nsite; isite++)*/
  for (iorb = 0; iorb < StdI->NOrb; iorb++) {
    for (isite = 0; isite < StdI->nsite; isite++) {
      for (jsite = 0; jsite < StdI->nsite; jsite++) {
        if (OrbGC[isite][jsite] == iorb) {
          OrbGC[jsite][isite] = OrbGC[isite][jsite];
        }
      }/*for (jsite = 0; jsite < isite; jsite++)*/
    }/*for (isite = 0; isite < StdI->nsite; isite++)*/
  }/*for (iorb = 0; iorb < StdI->NOrb; iorb++)*/
   
  NOrbGC = 0;
  for (isite = 0; isite < StdI->nsite; isite++) {
    for (jsite = 0; jsite < isite; jsite++) {
      if (OrbGC[isite][jsite] >= 0) {
        iOrbGC = OrbGC[isite][jsite];
        NOrbGC -= 1;
        for (isite1 = 0; isite1 < StdI->nsite; isite1++) {
          for (jsite1 = 0; jsite1 < StdI->nsite; jsite1++) {
            if (OrbGC[isite1][jsite1] == iOrbGC)
              OrbGC[isite1][jsite1] = NOrbGC;
          }/*for (jsite1 = 0; jsite1 < StdI->nsite; jsite1++)*/
        }/*for (isite1 = 0; isite1 < StdI->nsite; isite1++)*/
      }/*if (OrbGC[isite][jsite] >= 0)*/
    }/*for (jsite = 0; jsite < isite; jsite++)*/
  }/*for (isite = 0; isite < StdI->nsite; isite++)*/
   
  NOrbGC = -NOrbGC;
  for (isite = 0; isite < StdI->nsite; isite++) {
    for (jsite = 0; jsite < StdI->nsite; jsite++) {
      OrbGC[isite][jsite] = -1 - OrbGC[isite][jsite];
    }/*for (jsite = 0; jsite < isite; jsite++)*/
  }/*for (isite = 0; isite < StdI->nsite; isite++)*/

  fp = fopen("orbitalidxpara.def", "w");
  fprintf(fp, "=============================================\n");
  fprintf(fp, "NOrbitalIdx %10d\n", NOrbGC);
  fprintf(fp, "ComplexType %10d\n", StdI->ComplexType);
  fprintf(fp, "=============================================\n");
  fprintf(fp, "=============================================\n");

  for (isite = 0; isite < StdI->nsite; isite++) {
    for (jsite = 0; jsite < StdI->nsite; jsite++) {
      if (isite >= jsite) continue;
      if (StdI->AntiPeriod[0] == 1 || StdI->AntiPeriod[1] == 1 || StdI->AntiPeriod[2] == 1)
        fprintf(fp, "%5d  %5d  %5d  %5d\n", isite, jsite, OrbGC[isite][jsite], AntiOrbGC[isite][jsite]);
      else
        fprintf(fp, "%5d  %5d  %5d\n", isite, jsite, OrbGC[isite][jsite]);
    }/*for (jsite = 0; jsite < isite; jsite++)*/
  }/*for (isite = 0; isite < StdI->nsite; isite++)*/

  for (iOrbGC = 0; iOrbGC < NOrbGC; iOrbGC++)
    fprintf(fp, "%5d  %5d\n", iOrbGC, 1);

  fflush(fp);
  fclose(fp);
  fprintf(stdout, "    orbitalidxpara.def is written.\n");

  for (isite = 0; isite < StdI->nsite; isite++) {
    free(OrbGC[isite]);
    free(AntiOrbGC[isite]);
  }
  free(OrbGC);
  free(AntiOrbGC);
}/*static void PrintOrbPara*/
static void PrintGutzwiller(struct StdIntList *StdI)
{
  FILE *fp;
  int iCell, isite, jsite, NGutzwiller, iGutz;
  int *Gutz;

  Gutz = (int *)malloc(sizeof(int) * StdI->nsite);

  if (abs(StdI->NMPTrans) == 1 || StdI->NMPTrans == StdI->NaN_i) {
    if (strcmp(StdI->model, "hubbard") == 0) NGutzwiller = 0;
    else NGutzwiller = -1;

    for (isite = 0; isite < StdI->nsite; isite++) Gutz[isite] = StdI->Orb[isite][isite];

    for (isite = 0; isite < StdI->nsite; isite++) {
      /*
      For Local spin
      */
      if (StdI->locspinflag[isite] != 0) {
        Gutz[isite] = -1;
        continue;
      }
      
      if (Gutz[isite] >= 0) {
        iGutz = Gutz[isite];
        NGutzwiller -= 1;
        for (jsite = 0; jsite < StdI->nsite; jsite++) {
          if (Gutz[jsite] == iGutz)
            Gutz[jsite] = NGutzwiller;
        }/*for (jsite = 0; jsite < StdI->nsite; jsite++)*/
      }/*if (Gutz[isite] >= 0)*/
    }/*for (isite = 0; isite < StdI->nsite; isite++)*/
     
    NGutzwiller = -NGutzwiller;
    for (isite = 0; isite < StdI->nsite; isite++) {
      Gutz[isite] = -1 - Gutz[isite];
    }/*for (isite = 0; isite < StdI->nsite; isite++)*/
  }/*if (abs(StdI->NMPTrans) == 1)*/
  else {
    if (strcmp(StdI->model, "hubbard") == 0) NGutzwiller = StdI->NsiteUC;
    else if (strcmp(StdI->model, "spin") == 0) NGutzwiller = 1;
    else NGutzwiller = StdI->NsiteUC + 1;

    for (iCell = 0; iCell < StdI->NCell; iCell++) {
      for (isite = 0; isite < StdI->NsiteUC; isite++) {
        if (strcmp(StdI->model, "hubbard") == 0)
          Gutz[isite + StdI->NsiteUC*iCell] = isite;
        else if (strcmp(StdI->model, "spin") == 0)
          Gutz[isite + StdI->NsiteUC*iCell] = 0;
        else {
          Gutz[isite + StdI->NsiteUC*iCell] = 0;
          Gutz[isite + StdI->NsiteUC*(iCell + StdI->NCell)] = isite + 1;
        }
      }/*for (isite = 0; isite < StdI->NsiteUC; isite++)*/
    }/*for (iCell = 0; iCell < StdI->NCell; iCell++)*/
  }/*if (abs(StdI->NMPTrans) != 1)*/

  fp = fopen("gutzwilleridx.def", "w");
  fprintf(fp, "=============================================\n");
  fprintf(fp, "NGutzwillerIdx %10d\n", NGutzwiller);
  fprintf(fp, "ComplexType %10d\n", 0);
  fprintf(fp, "=============================================\n");
  fprintf(fp, "=============================================\n");

  for (isite = 0; isite < StdI->nsite; isite++)
    fprintf(fp, "%5d  %5d\n", isite, Gutz[isite]);

  for (iGutz = 0; iGutz < NGutzwiller; iGutz++) {
    if (strcmp(StdI->model, "hubbard") == 0 || iGutz > 0)
      fprintf(fp, "%5d  %5d\n", iGutz, 1);
    else
      fprintf(fp, "%5d  %5d\n", iGutz, 0);
  }/*for (iGutz = 0; iGutz < NGutzwiller; iGutz++)*/
  fflush(fp);
  fclose(fp);
  fprintf(stdout, "    gutzwilleridx.def is written.\n");

  free(Gutz);
}/*static void PrintGutzwiller*/
#endif

static void StdFace_ResetVals(struct StdIntList *StdI) {
  int i, j;
  double NaN_d;
  /*
  NaN is used for not inputed variable
  */
  NaN_d = 0.0 / 0.0;
  StdI->NaN_i = 2147483647;
  StdI->pi = acos(-1.0);
  
  StdI->a = NaN_d;
  for (i = 0; i < 3; i++) StdI->length[i] = NaN_d;
  for (i = 0; i < 3; i++)
    for (j = 0; j < 3; j++)
      StdI->box[i][j] = StdI->NaN_i;
  StdI->Gamma = NaN_d;
  StdI->h = NaN_d;
  StdI->Height = StdI->NaN_i;
  StdI->JAll = NaN_d;
  StdI->JpAll = NaN_d;
  StdI->JpAll = NaN_d;
  StdI->JppAll = NaN_d;
  StdI->J0All = NaN_d;
  StdI->J0pAll = NaN_d;
  StdI->J0ppAll = NaN_d;
  StdI->J1All = NaN_d;
  StdI->J1pAll = NaN_d;
  StdI->J1ppAll = NaN_d;
  StdI->J2All = NaN_d;
  StdI->J2pAll = NaN_d;
  StdI->J2ppAll = NaN_d;
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      StdI->J[i][j] = NaN_d;
      StdI->Jp[i][j] = NaN_d;
      StdI->Jpp[i][j] = NaN_d;
      StdI->J0[i][j] = NaN_d;
      StdI->J0p[i][j] = NaN_d;
      StdI->J0pp[i][j] = NaN_d;
      StdI->J1[i][j] = NaN_d;
      StdI->J1p[i][j] = NaN_d;
      StdI->J1pp[i][j] = NaN_d;
      StdI->J2[i][j] = NaN_d;
      StdI->J2p[i][j] = NaN_d;
      StdI->J2pp[i][j] = NaN_d;
      StdI->D[i][j] = 0.0;
    }
  }
  StdI->D[2][2] = NaN_d;
  StdI->K = NaN_d;
  StdI->L = StdI->NaN_i;
  for (i = 0; i < 3; i++) 
    for (j = 0; j < 3; j++)
      StdI->direct[i][j] = NaN_d;
  StdI->mu = NaN_d;
  StdI->S2 = StdI->NaN_i;
  StdI->t = NaN_d;
  StdI->tp = NaN_d;
  StdI->tpp = NaN_d;
  StdI->t0 = NaN_d;
  StdI->t0p = NaN_d;
  StdI->t0pp = NaN_d;
  StdI->t1 = NaN_d;
  StdI->t1p = NaN_d;
  StdI->t1pp = NaN_d;
  StdI->t2 = NaN_d;
  StdI->t2p = NaN_d;
  StdI->t2pp = NaN_d;
  StdI->U = NaN_d;
  StdI->V = NaN_d;
  StdI->Vp = NaN_d;
  StdI->Vpp = NaN_d;
  StdI->V0 = NaN_d;
  StdI->V0p = NaN_d;
  StdI->V0pp = NaN_d;
  StdI->V1 = NaN_d;
  StdI->V1p = NaN_d;
  StdI->V1pp = NaN_d;
  StdI->V2 = NaN_d;
  StdI->V2p = NaN_d;
  StdI->V2pp = NaN_d;
  StdI->W = StdI->NaN_i;
  for (i = 0; i < 3; i++)StdI->phase[i] = NaN_d;
  StdI->pi180 = StdI->pi / 180.0;

  StdI->nelec = StdI->NaN_i;
  StdI->Sz2 = StdI->NaN_i;
  strcpy(StdI->model, "****\0");
  strcpy(StdI->lattice, "****\0");
  strcpy(StdI->outputmode, "****\0");
  strcpy(StdI->CDataFileHead, "****\0");
  StdI->cutoff_t = NaN_d;
  StdI->cutoff_u = NaN_d;
  StdI->cutoff_j = NaN_d;
  StdI->cutoff_length_t = NaN_d;
  StdI->cutoff_length_U = NaN_d;
  StdI->cutoff_length_J = NaN_d;
  for (i = 0; i < 3; i++)StdI->cutoff_tR[i] = StdI->NaN_i;
  for (i = 0; i < 3; i++)StdI->cutoff_UR[i] = StdI->NaN_i;
  for (i = 0; i < 3; i++)StdI->cutoff_JR[i] = StdI->NaN_i;
#if defined(_HPhi)
  StdI->LargeValue = NaN_d;
  StdI->OmegaMax = NaN_d;
  StdI->OmegaMin = NaN_d;
  StdI->OmegaIm = NaN_d;
  StdI->Nomega = StdI->NaN_i;
  for (i = 0; i < 3; i++)StdI->SpectrumQ[i] = NaN_d;
  strcpy(StdI->method, "****\0");
  strcpy(StdI->Restart, "****\0");
  strcpy(StdI->EigenVecIO, "****\0");
  strcpy(StdI->InitialVecType, "****\0");
  strcpy(StdI->HamIO, "****\0");
  strcpy(StdI->CalcSpec, "****\0");
  strcpy(StdI->SpectrumType, "****\0");
  strcpy(StdI->OutputExVec, "****\0");
  StdI->FlgTemp = 1;
  StdI->Lanczos_max = StdI->NaN_i;
  StdI->initial_iv = StdI->NaN_i;
  StdI->nvec = StdI->NaN_i;
  StdI->exct = StdI->NaN_i;
  StdI->LanczosEps = StdI->NaN_i;
  StdI->LanczosTarget = StdI->NaN_i;
  StdI->NumAve = StdI->NaN_i;
  StdI->ExpecInterval = StdI->NaN_i;
  StdI->dt = NaN_d;
  StdI->tdump = NaN_d;
  StdI->tshift = NaN_d;
  StdI->freq = NaN_d;
  StdI->Uquench = NaN_d;
  for (i = 0; i < 3; i++)StdI->VecPot[i] = NaN_d;;
  strcpy(StdI->PumpType, "****\0");
  StdI->ExpandCoef = StdI->NaN_i;
#elif defined(_mVMC)
  strcpy(StdI->CParaFileHead, "****\0");
  StdI->NVMCCalMode = StdI->NaN_i;
  StdI->NLanczosMode = StdI->NaN_i;
  StdI->NDataIdxStart = StdI->NaN_i;
  StdI->NDataQtySmp = StdI->NaN_i;
  StdI->NSPGaussLeg = StdI->NaN_i;
  StdI->NSPStot = StdI->NaN_i;
  StdI->NMPTrans = StdI->NaN_i;
  StdI->NSROptItrStep = StdI->NaN_i;
  StdI->NSROptItrSmp = StdI->NaN_i;
  StdI->DSROptRedCut = NaN_d;
  StdI->DSROptStaDel = NaN_d;
  StdI->DSROptStepDt = NaN_d;
  StdI->NVMCWarmUp = StdI->NaN_i;
  StdI->NVMCInterval = StdI->NaN_i;
  StdI->NVMCSample = StdI->NaN_i;
  StdI->NExUpdatePath = StdI->NaN_i;
  StdI->RndSeed = StdI->NaN_i;
  StdI->NSplitSize = StdI->NaN_i;
  StdI->NStore = StdI->NaN_i;
  StdI->NSRCG = StdI->NaN_i;
  StdI->ComplexType = StdI->NaN_i;
  for (i = 0; i < 3; i++)
    for (j = 0; j < 3; j++)
      StdI->boxsub[i][j] = StdI->NaN_i;
  StdI->Hsub = StdI->NaN_i;
  StdI->Lsub = StdI->NaN_i;
  StdI->Wsub = StdI->NaN_i;
#endif
}/*static void StdFace_ResetVals*/
/*
@brief Make all characters lower
@author Mitsuaki Kawamura (The University of Tokyo)
*/
static void Text2Lower(char *value 
){
  char value2;
  int valuelen, ii;

  valuelen = strlen(value);
  for (ii = 0; ii < valuelen; ii++) {
    value2 = tolower(value[ii]);
    value[ii] = value2;
  }
}/*static void Text2Lower*/
static void TrimSpaceQuote(char *value 
){
  char value2[256];
  int valuelen, valuelen2, ii;

  valuelen = strlen(value);
  valuelen2 = 0;
  for (ii = 0; ii < valuelen; ii++){
    if (value[ii] != ' ' &&
      value[ii] != ':' &&
      value[ii] != ';' &&
      value[ii] != '\"' &&
      value[ii] != '\b' &&
      value[ii] != '\\' &&
      value[ii] != '\v' &&
      value[ii] != '\n' &&
      value[ii] != '\0'){
      value2[valuelen2] = value[ii];
      valuelen2++;
    }
  }

  strncpy(value, value2, valuelen2);
  value[valuelen2] = '\0';

}/*static void TrimSpaceQuote*/
static void StoreWithCheckDup_s(
  char *keyword,
  char *valuestring,
  char *value
)
{
  if (strcmp(value, "****") != 0){
    fprintf(stdout, "ERROR !  Keyword %s is duplicated ! \n", keyword);
    StdFace_exit(-1);
  }
  else{
    strcpy(value, valuestring);
  }
}/*static void StoreWithCheckDup_s*/
static void StoreWithCheckDup_sl(
  char *keyword,
  char *valuestring,
  char *value
)
{
  if (strcmp(value, "****") != 0) {
    fprintf(stdout, "ERROR !  Keyword %s is duplicated ! \n", keyword);
    StdFace_exit(-1);
  }
  else {
    strcpy(value, valuestring);
    Text2Lower(value);
  }
}/*static void StoreWithCheckDup_sl*/
static void StoreWithCheckDup_i(
  char *keyword,
  char *valuestring,
  int *value
)
{
  int NaN_i = 2147483647;

  if (*value != NaN_i){
    fprintf(stdout, "ERROR !  Keyword %s is duplicated ! \n", keyword);
    StdFace_exit(-1);
  }
  else{
    sscanf(valuestring, "%d", value);
  }
}/*static void StoreWithCheckDup_i*/
static void StoreWithCheckDup_d(
  char *keyword,
  char *valuestring,
  double *value
)
{
  if (isnan(*value) == 0){
    fprintf(stdout, "ERROR !  Keyword %s is duplicated ! \n", keyword);
    StdFace_exit(-1);
  }
  else{
    sscanf(valuestring, "%lf", value);
  }
}/*static void StoreWithCheckDup_d*/
static void StoreWithCheckDup_c(
  char *keyword,
  char *valuestring,
  double complex *value
)
{
  int num;
  char *valuestring_r, *valuestring_i;
  double value_r, value_i;

  if (isnan(creal(*value)) == 0) {
    fprintf(stdout, "ERROR !  Keyword %s is duplicated ! \n", keyword);
    StdFace_exit(-1);
  }
  else {

    if (valuestring[0] == ',') {
      valuestring_r = NULL;
      valuestring_i = strtok(valuestring, ",");
    }
    else {
      valuestring_r = strtok(valuestring, ",");
      valuestring_i = strtok(NULL, ",");
    }
    
    if (valuestring_r == NULL) {
      *value = 0.0;
    }
    else {
      num = sscanf(valuestring_r, "%lf", &value_r);
      if (num == 1) *value = value_r;
      else *value = 0.0;
    }

    if (valuestring_i == NULL) {
      *value += I * 0.0;
    }
    else {
        num = sscanf(valuestring_i, "%lf", &value_i);
      if (num == 1) *value += I * value_i;
      else *value += I * 0.0;
    }
  }
}/*static void StoreWithCheckDup_c*/
static void PrintLocSpin(struct StdIntList *StdI) {
  FILE *fp;
  int isite, nlocspin;

  nlocspin = 0;
  for (isite = 0; isite < StdI->nsite; isite++)
    if (StdI->locspinflag[isite] != 0) nlocspin = nlocspin + 1;

  fp = fopen("locspn.def", "w");
  fprintf(fp, "================================ \n");
  fprintf(fp, "NlocalSpin %5d  \n", nlocspin);
  fprintf(fp, "================================ \n");
  fprintf(fp, "========i_0LocSpn_1IteElc ====== \n");
  fprintf(fp, "================================ \n");

  for (isite = 0; isite < StdI->nsite; isite++)
    fprintf(fp, "%5d  %5d\n", isite, StdI->locspinflag[isite]);

  fflush(fp);
  fclose(fp);
  fprintf(stdout, "    locspn.def is written.\n");
}/*static void PrintLocSpin*/
static void PrintTrans(struct StdIntList *StdI){
  FILE *fp;
  int jtrans, ktrans, ntrans0;

  for (jtrans = 0; jtrans < StdI->ntrans; jtrans++){
    for (ktrans = jtrans + 1; ktrans < StdI->ntrans; ktrans++){
      if (StdI->transindx[jtrans][0] == StdI->transindx[ktrans][0]
        && StdI->transindx[jtrans][1] == StdI->transindx[ktrans][1]
        && StdI->transindx[jtrans][2] == StdI->transindx[ktrans][2]
        && StdI->transindx[jtrans][3] == StdI->transindx[ktrans][3]){
        StdI->trans[jtrans] = StdI->trans[jtrans] + StdI->trans[ktrans];
        StdI->trans[ktrans] = 0.0;
      }
    }/*for (ktrans = jtrans + 1; ktrans < StdI->ntrans; ktrans++)*/
  }/*for (jtrans = 0; jtrans < StdI->ntrans; jtrans++)*/

  ntrans0 = 0;
  for (ktrans = 0; ktrans < StdI->ntrans; ktrans++){
    if (cabs(StdI->trans[ktrans]) > 0.000001) ntrans0 = ntrans0 + 1;
  }

  fp = fopen("trans.def", "w");
  fprintf(fp, "======================== \n");
  fprintf(fp, "NTransfer %7d  \n", ntrans0);
  fprintf(fp, "======================== \n");
  fprintf(fp, "========i_j_s_tijs====== \n");
  fprintf(fp, "======================== \n");

  ntrans0 = 0;
  for (ktrans = 0; ktrans < StdI->ntrans; ktrans++) {
    if (cabs(StdI->trans[ktrans]) > 0.000001)
      fprintf(fp, "%5d %5d %5d %5d %25.15f %25.15f\n",
        StdI->transindx[ktrans][0], StdI->transindx[ktrans][1],
        StdI->transindx[ktrans][2], StdI->transindx[ktrans][3],
        creal(StdI->trans[ktrans]), cimag(StdI->trans[ktrans]));
  }

  fflush(fp);
  fclose(fp);
  fprintf(stdout, "      trans.def is written.\n");
}/*static void PrintTrans*/
static void PrintNamelist(struct StdIntList *StdI){
  FILE *fp;

  fp = fopen("namelist.def", "w");
  fprintf(                         fp, "         ModPara  modpara.def\n");
  fprintf(                         fp, "         LocSpin  locspn.def\n");
  fprintf(                         fp, "           Trans  trans.def\n");
  if (StdI->LCintra == 1) fprintf( fp, "    CoulombIntra  coulombintra.def\n");
  if (StdI->LCinter == 1) fprintf( fp, "    CoulombInter  coulombinter.def\n");
  if (StdI->LHund == 1)fprintf(    fp, "            Hund  hund.def\n");
  if (StdI->LEx == 1)fprintf(      fp, "        Exchange  exchange.def\n");
  if (StdI->LPairLift == 1)fprintf(fp, "        PairLift  pairlift.def\n");
  if (StdI->LPairHopp == 1)fprintf(fp, "         PairHop  pairhopp.def\n");
  if (StdI->Lintr == 1)fprintf(    fp, "        InterAll  interall.def\n");
  if (StdI->ioutputmode != 0) {
    fprintf(                       fp, "        OneBodyG  greenone.def\n");
    fprintf(                       fp, "        TwoBodyG  greentwo.def\n");
  }
#if defined(_HPhi)
  fprintf(                         fp, "         CalcMod  calcmod.def\n");
  if(StdI->SpectrumBody == 1) 
    fprintf(                       fp, "SingleExcitation  single.def\n");
  else fprintf(                    fp, "  PairExcitation  pair.def\n");
  if (strcmp(StdI->method, "timeevolution") == 0) {
    if (StdI->PumpBody == 1)
      fprintf(fp, "       TEOneBody  teone.def\n");
    else if (StdI->PumpBody == 2)
      fprintf(fp, "       TETwoBody  tetwo.def\n");
  }/*if (strcmp(StdI->method, "timeevolution") == 0)*/
  fprintf(                         fp, "     SpectrumVec  %s_eigenvec_0\n",
                                   StdI->CDataFileHead);
  if (StdI->lBoost == 1) fprintf(  fp, "           Boost  boost.def\n");
#elif defined(_mVMC)
  fprintf(                         fp, "      Gutzwiller  gutzwilleridx.def\n");
  fprintf(                         fp, "         Jastrow  jastrowidx.def\n");
  fprintf(                         fp, "         Orbital  orbitalidx.def\n");
  if (StdI->lGC == 1 || (StdI->Sz2 != 0 && StdI->Sz2 != StdI->NaN_i))
    fprintf(fp, " OrbitalParallel  orbitalidxpara.def\n");
  fprintf(                         fp, "        TransSym  qptransidx.def\n");
#endif
  
  fflush(fp);
  fclose(fp);
  fprintf(stdout, "    namelist.def is written.\n");
}/*static void PrintNamelist*/
static void PrintModPara(struct StdIntList *StdI)
{
  FILE *fp;

  fp = fopen("modpara.def", "w");
  fprintf(fp, "--------------------\n");
  fprintf(fp, "Model_Parameters   0\n");
  fprintf(fp, "--------------------\n");
#if defined(_HPhi)
  fprintf(fp, "HPhi_Cal_Parameters\n");
  fprintf(fp, "--------------------\n");
  fprintf(fp, "CDataFileHead  %s\n", StdI->CDataFileHead);
  fprintf(fp, "CParaFileHead  zqp\n");
  fprintf(fp, "--------------------\n");
  fprintf(fp, "Nsite          %-5d\n", StdI->nsite);
  if (StdI->Sz2 != StdI->NaN_i) fprintf(fp, "2Sz            %-5d\n", StdI->Sz2);
  if (StdI->nelec != StdI->NaN_i) fprintf(fp, "Ncond          %-5d\n", StdI->nelec);
  fprintf(fp, "Lanczos_max    %-5d\n", StdI->Lanczos_max);
  fprintf(fp, "initial_iv     %-5d\n", StdI->initial_iv);
  if(StdI->nvec != StdI->NaN_i) fprintf(fp, "nvec           %-5d\n", StdI->nvec);
  fprintf(fp, "exct           %-5d\n", StdI->exct);
  fprintf(fp, "LanczosEps     %-5d\n", StdI->LanczosEps);
  fprintf(fp, "LanczosTarget  %-5d\n", StdI->LanczosTarget);
  fprintf(fp, "LargeValue     %-25.15e\n", StdI->LargeValue);
  fprintf(fp, "NumAve         %-5d\n", StdI->NumAve);
  fprintf(fp, "ExpecInterval  %-5d\n", StdI->ExpecInterval);
  fprintf(fp, "NOmega         %-5d\n", StdI->Nomega);
  fprintf(fp, "OmegaMax       %-25.15e %-25.15e\n", StdI->OmegaMax, StdI->OmegaIm);
  fprintf(fp, "OmegaMin       %-25.15e %-25.15e\n", StdI->OmegaMin, StdI->OmegaIm);
  fprintf(fp, "OmegaOrg       0.0 0.0\n");
  if (strcmp(StdI->method, "timeevolution") == 0)
    fprintf(fp, "ExpandCoef     %-5d\n", StdI->ExpandCoef);
#elif defined(_mVMC)
  fprintf(fp, "VMC_Cal_Parameters\n");
  fprintf(fp, "--------------------\n");
  fprintf(fp, "CDataFileHead  %s\n", StdI->CDataFileHead);
  fprintf(fp, "CParaFileHead  %s\n", StdI->CParaFileHead);
  fprintf(fp, "--------------------\n");
  fprintf(fp, "NVMCCalMode    %d\n", StdI->NVMCCalMode);
  /*fprintf(fp, "NLanczosMode   %d\n", StdI->NLanczosMode);*/
  fprintf(fp, "--------------------\n");
  fprintf(fp, "NDataIdxStart  %d\n", StdI->NDataIdxStart);
  fprintf(fp, "NDataQtySmp    %d\n", StdI->NDataQtySmp);
  fprintf(fp, "--------------------\n");
  fprintf(fp, "Nsite          %d\n", StdI->nsite);
  fprintf(fp, "Ncond          %-5d\n", StdI->nelec);
  if (StdI->Sz2 != StdI->NaN_i)
    fprintf(fp, "2Sz            %d\n", StdI->Sz2);
  if (StdI->NSPGaussLeg != StdI->NaN_i)
    fprintf(fp, "NSPGaussLeg    %d\n", StdI->NSPGaussLeg);
  if (StdI->NSPStot != StdI->NaN_i)
    fprintf(fp, "NSPStot        %d\n", StdI->NSPStot);
  fprintf(fp, "NMPTrans       %d\n", StdI->NMPTrans);
  fprintf(fp, "NSROptItrStep  %d\n", StdI->NSROptItrStep);
  fprintf(fp, "NSROptItrSmp   %d\n", StdI->NSROptItrSmp);
  fprintf(fp, "DSROptRedCut   %.10f\n", StdI->DSROptRedCut);
  fprintf(fp, "DSROptStaDel   %.10f\n", StdI->DSROptStaDel);
  fprintf(fp, "DSROptStepDt   %.10f\n", StdI->DSROptStepDt);
  fprintf(fp, "NVMCWarmUp     %d\n", StdI->NVMCWarmUp);
  fprintf(fp, "NVMCInterval   %d\n", StdI->NVMCInterval);
  fprintf(fp, "NVMCSample     %d\n", StdI->NVMCSample);
  fprintf(fp, "NExUpdatePath  %d\n", StdI->NExUpdatePath);
  fprintf(fp, "RndSeed        %d\n", StdI->RndSeed);
  fprintf(fp, "NSplitSize     %d\n", StdI->NSplitSize);
  fprintf(fp, "NStore         %d\n", StdI->NStore);
  fprintf(fp, "NSRCG          %d\n", StdI->NSRCG);
#endif

  fflush(fp);
  fclose(fp);
  fprintf(stdout, "     modpara.def is written.\n");
}/*static void PrintModPara*/
static void Print1Green(struct StdIntList *StdI)
{
  FILE *fp;
  int ngreen, igreen, store, xkondo;
  int isite, jsite, ispin, jspin, SiMax, SjMax;
  int **greenindx;
  /*
   Set Indices of correlation functions
  */
  ngreen = 0;
  if (StdI->ioutputmode != 0) {
    for (store = 0; store < 2; store++) {

      if (store == 1) {
        greenindx = (int **)malloc(sizeof(int*) * (ngreen + 1));
        for (igreen = 0; igreen < ngreen; igreen++) {
          greenindx[igreen] = (int *)malloc(sizeof(int) * 4);
        }
        ngreen = 0;
      }/*if (store == 1)*/

      if (strcmp(StdI->model, "kondo") == 0) xkondo = 2;
      else xkondo = 1;

      if (StdI->ioutputmode == 1) {
        for (isite = 0; isite < StdI->NsiteUC*xkondo; isite++) {

          if (isite >= StdI->NsiteUC) isite += StdI->nsite / 2;

          if (StdI->locspinflag[isite] == 0) SiMax = 1;
          else SiMax = StdI->locspinflag[isite];

          for (ispin = 0; ispin <= SiMax; ispin++) {
            for (jsite = 0; jsite < StdI->nsite; jsite++) {

              if (StdI->locspinflag[jsite] == 0) SjMax = 1;
              else SjMax = StdI->locspinflag[jsite];

              for (jspin = 0; jspin <= SjMax; jspin++) {

                if (isite != jsite &&
                  (StdI->locspinflag[isite] != 0 && StdI->locspinflag[jsite] != 0)) continue;

                if (ispin == jspin){
                  if (store == 1) {
                    greenindx[ngreen][0] = isite;
                    greenindx[ngreen][1] = ispin;
                    greenindx[ngreen][2] = jsite;
                    greenindx[ngreen][3] = jspin;
                  }
                  ngreen++;
                }

              }/*for (jspin = 0; jspin <= SjMax; jspin++)*/
            }/*for (jsite = 0; jsite < StdI->nsite; jsite++)*/
          }/*for (ispin = 0; ispin <= SiMax; ispin++)*/
        }/*for (isite = 0; isite < StdI->nsite; isite++)*/
      }/*if (StdI->ioutputmode == 1)*/
      else {
        for (isite = 0; isite < StdI->nsite; isite++) {

          if (StdI->locspinflag[isite] == 0) SiMax = 1;
          else SiMax = StdI->locspinflag[isite];

          for (ispin = 0; ispin <= SiMax; ispin++) {
            for (jsite = 0; jsite < StdI->nsite; jsite++) {

              if (StdI->locspinflag[jsite] == 0) SjMax = 1;
              else SjMax = StdI->locspinflag[jsite];

              for (jspin = 0; jspin <= SjMax; jspin++) {

                if (isite != jsite &&
                  (StdI->locspinflag[isite] != 0 && StdI->locspinflag[jsite] != 0)) continue;

                if (store == 1) {
                  greenindx[ngreen][0] = isite;
                  greenindx[ngreen][1] = ispin;
                  greenindx[ngreen][2] = jsite;
                  greenindx[ngreen][3] = jspin;
                }
                ngreen++;

              }/*for (jspin = 0; jspin <= SjMax; jspin++)*/
            }/*for (jsite = 0; jsite < StdI->nsite; jsite++)*/
          }/*for (ispin = 0; ispin <= SiMax; ispin++)*/
        }/*for (isite = 0; isite < StdI->nsite; isite++)*/
      }/*if (StdI->ioutputmode == 2)*/
    }/*if (StdI->ioutputmode != 0)*/

    fp = fopen("greenone.def", "w");
    fprintf(fp, "===============================\n");
    fprintf(fp, "NCisAjs %10d\n", ngreen);
    fprintf(fp, "===============================\n");
    fprintf(fp, "======== Green functions ======\n");
    fprintf(fp, "===============================\n");
    for (igreen = 0; igreen < ngreen; igreen++) {
    fprintf(fp,   "%5d %5d %5d %5d\n",
      greenindx[igreen][0], greenindx[igreen][1], greenindx[igreen][2], greenindx[igreen][3]);
    }
    fflush(fp);
    fclose(fp);

    fprintf(stdout, "    greenone.def is written.\n");

    for (igreen = 0; igreen < ngreen; igreen++) {
      free(greenindx[igreen]);
    }
    free(greenindx);

  }/*if (StdI->ioutputmode != 0) */
}/*static void Print1Green*/
static void Print2Green(struct StdIntList *StdI) {
  FILE *fp;
  int ngreen, store, igreen, xkondo;
  int site1, site2, site3, site4;
  int spin1, spin2, spin3, spin4;
  int S1Max, S2Max, S3Max, S4Max;
  int **greenindx;
  /*
   Set Indices of correlation functions
  */
  ngreen = 0;
  if (StdI->ioutputmode == 1) {
    for (store = 0; store < 2; store++) {

      if (store == 1) {
        greenindx = (int **)malloc(sizeof(int*) * (ngreen + 1));
        for (igreen = 0; igreen < ngreen; igreen++)
          greenindx[igreen] = (int *)malloc(sizeof(int) * 8);
        ngreen = 0;
      }/*if (store == 1)*/

      if (strcmp(StdI->model, "kondo") == 0) xkondo = 2;
      else xkondo = 1;

      for (site1 = 0; site1 < StdI->NsiteUC*xkondo; site1++) {

        if (site1 >= StdI->NsiteUC) site1 += StdI->nsite / 2;

        if (StdI->locspinflag[site1] == 0) S1Max = 1;
        else S1Max = StdI->locspinflag[site1];
        for (spin1 = 0; spin1 <= S1Max; spin1++) {
          for (spin2 = 0; spin2 <= S1Max; spin2++) {

            for (site3 = 0; site3 < StdI->nsite; site3++) {

              if (StdI->locspinflag[site3] == 0) S3Max = 1;
              else S3Max = StdI->locspinflag[site3];
              for (spin3 = 0; spin3 <= S3Max; spin3++) {
                for (spin4 = 0; spin4 <= S3Max; spin4++) {

                  if (spin1 - spin2 + spin3 - spin4 == 0) {
                    if (store == 1) {
#if defined(_mVMC)
                      if (spin1 != spin2 || spin3 != spin4)
                      {
                        greenindx[ngreen][0] = site1;
                        greenindx[ngreen][1] = spin1;
                        greenindx[ngreen][2] = site3;
                        greenindx[ngreen][3] = spin4;
                        greenindx[ngreen][4] = site3;
                        greenindx[ngreen][5] = spin3;
                        greenindx[ngreen][6] = site1;
                        greenindx[ngreen][7] = spin2;
                      }
                      else
#endif
                      {
                        greenindx[ngreen][0] = site1;
                        greenindx[ngreen][1] = spin1;
                        greenindx[ngreen][2] = site1;
                        greenindx[ngreen][3] = spin2;
                        greenindx[ngreen][4] = site3;
                        greenindx[ngreen][5] = spin3;
                        greenindx[ngreen][6] = site3;
                        greenindx[ngreen][7] = spin4;
                      }
                    }/*if (store == 1)*/
                    ngreen++;
                  }/*if (spin1 - spin2 + spin3 - spin4 == 0)*/

                }/*for (spin4 = 0; spin4 <= S3Max; spin4++)*/
              }/*for (spin3 = 0; spin3 <= S3Max; spin3++*/
            }/*for (site3 = 0; site3 < StdI->nsite; site3++)*/
          }/*for (spin2 = 0; spin2 <= S1Max; spin2++)*/
        }/*for (spin1 = 0; spin1 <= S1Max; spin1++)*/
      }/*for (site1 = 0; site1 < StdI->nsite; site1++)*/

    }/*for (store = 0; store < 2; store++)*/
  }/*if (StdI->ioutputmode == 1)*/
  else if (StdI->ioutputmode == 2) {
    for (store = 0; store < 2; store++) {

      if (store == 1) {
        greenindx = (int **)malloc(sizeof(int*) * (ngreen + 1));
        for (igreen = 0; igreen < ngreen; igreen++)
          greenindx[igreen] = (int *)malloc(sizeof(int) * 8);
        ngreen = 0;
      }/*if (store == 1)*/

      for (site1 = 0; site1 < StdI->nsite; site1++) {

        if (StdI->locspinflag[site1] == 0) S1Max = 1;
        else S1Max = StdI->locspinflag[site1];
        for (spin1 = 0; spin1 <= S1Max; spin1++) {

          for (site2 = 0; site2 < StdI->nsite; site2++) {

            if (StdI->locspinflag[site1] != 0 && StdI->locspinflag[site2] != 0
              && site1 != site2) continue;

            if (StdI->locspinflag[site2] == 0) S2Max = 1;
            else S2Max = StdI->locspinflag[site2];
            for (spin2 = 0; spin2 <= S2Max; spin2++) {

              for (site3 = 0; site3 < StdI->nsite; site3++) {

                if (StdI->locspinflag[site3] == 0) S3Max = 1;
                else S3Max = StdI->locspinflag[site3];
                for (spin3 = 0; spin3 <= S3Max; spin3++) {

                  for (site4 = 0; site4 < StdI->nsite; site4++) {

                    if (StdI->locspinflag[site3] != 0 && StdI->locspinflag[site4] != 0
                      && site3 != site4) continue;

                    if (StdI->locspinflag[site4] == 0) S4Max = 1;
                    else S4Max = StdI->locspinflag[site4];
                    for (spin4 = 0; spin4 <= S4Max; spin4++) {

                      if (store == 1) {
                        greenindx[ngreen][0] = site1;
                        greenindx[ngreen][1] = spin1;
                        greenindx[ngreen][2] = site2;
                        greenindx[ngreen][3] = spin2;
                        greenindx[ngreen][4] = site3;
                        greenindx[ngreen][5] = spin3;
                        greenindx[ngreen][6] = site4;
                        greenindx[ngreen][7] = spin4;
                      }/*if (store == 1)*/
                      ngreen++;

                    }/*for (spin4 = 0; spin4 <= S4Max; spin4++)*/
                  }/*for (site4 = 0; site4 < StdI->nsite; site4++)*/
                }/*for (spin3 = 0; spin3 <= S3Max; spin3++*/
              }/*for (site3 = 0; site3 < StdI->nsite; site3++)*/
            }/*for (spin2 = 0; spin2 <= S2Max; spin2++)*/
          }/*for (site2 = 0; site2 < StdI->nsite; site2++)*/
        }/*for (spin1 = 0; spin1 <= S1Max; spin1++)*/
      }/*for (site1 = 0; site1 < StdI->nsite; site1++)*/

    }/*for (store = 0; store < 2; store++)*/
  }/*if (StdI->ioutputmode == 2)*/
  if (StdI->ioutputmode != 0) {
    fp = fopen("greentwo.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NCisAjsCktAltDC %10d\n", ngreen);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "======== Green functions for Sq AND Nq ======\n");
    fprintf(fp, "=============================================\n");
    for (igreen = 0; igreen < ngreen; igreen++) {
      fprintf(fp, "%5d %5d %5d %5d %5d %5d %5d %5d\n",
        greenindx[igreen][0], greenindx[igreen][1], greenindx[igreen][2], greenindx[igreen][3],
        greenindx[igreen][4], greenindx[igreen][5], greenindx[igreen][6], greenindx[igreen][7]);
    }
    fflush(fp);
    fclose(fp);

    fprintf(stdout, "    greentwo.def is written.\n");

    for (igreen = 0; igreen < ngreen; igreen++) {
      free(greenindx[igreen]);
    }
    free(greenindx);
  }/*if (StdI->ioutputmode != 0)*/
}/*static void Print2Green(struct StdIntList *StdI)*/
static void UnsupportedSystem(
  char *model,
  char *lattice
)
{
  fprintf(stdout, "\nSorry, specified combination, \n");
  fprintf(stdout, "    MODEL : %s  \n", model);
  fprintf(stdout, "  LATTICE : %s, \n", lattice);
  fprintf(stdout, "is unsupported in the STANDARD MODE...\n");
  fprintf(stdout, "Please use the EXPART MODE, or write a NEW FUNCTION and post us.\n");
  StdFace_exit(-1);
}/*static void UnsupportedSystem*/
static void CheckOutputMode(struct StdIntList *StdI)
{
  /*
  Form for Correlation function
  */
  if (strcmp(StdI->outputmode, "non") == 0
    || strcmp(StdI->outputmode, "none") == 0
    || strcmp(StdI->outputmode, "off") == 0) {
    StdI->ioutputmode = 0;
    fprintf(stdout, "      ioutputmode = %-10d\n", StdI->ioutputmode);
  }
  else if (strcmp(StdI->outputmode, "cor") == 0
    || strcmp(StdI->outputmode, "corr") == 0
    || strcmp(StdI->outputmode, "correlation") == 0) {
    StdI->ioutputmode = 1;
    fprintf(stdout, "      ioutputmode = %-10d\n", StdI->ioutputmode);
  }
  else if (strcmp(StdI->outputmode, "****") == 0) {
    StdI->ioutputmode = 1;
    fprintf(stdout, "      ioutputmode = %-10d  ######  DEFAULT VALUE IS USED  ######\n", StdI->ioutputmode);
  }
  else if (strcmp(StdI->outputmode, "raw") == 0
    || strcmp(StdI->outputmode, "all") == 0
    || strcmp(StdI->outputmode, "full") == 0) {
    StdI->ioutputmode = 2;
    fprintf(stdout, "      ioutputmode = %-10d\n", StdI->ioutputmode);
  }
  else{
    fprintf(stdout, "\n ERROR ! Unsupported OutPutMode : %s\n", StdI->outputmode);
    StdFace_exit(-1);
  }
}/*static void CheckOutputMode*/
static void CheckModPara(struct StdIntList *StdI)
{

  
#if defined(_HPhi)
  StdFace_PrintVal_i("Lanczos_max", &StdI->Lanczos_max, 2000);
  StdFace_PrintVal_i("initial_iv", &StdI->initial_iv, -1);
  /*StdFace_PrintVal_i("nvec", &StdI->nvec, 1);*/
  StdFace_PrintVal_i("exct", &StdI->exct, 1);
  StdFace_PrintVal_i("LanczosEps", &StdI->LanczosEps, 14);
  StdFace_PrintVal_i("LanczosTarget", &StdI->LanczosTarget, 2);
  if(StdI->LanczosTarget < StdI->exct) StdI->LanczosTarget = StdI->exct;
  StdFace_PrintVal_i("NumAve", &StdI->NumAve, 5);
  StdFace_PrintVal_i("ExpecInterval", &StdI->ExpecInterval, 20);
  StdFace_PrintVal_i("NOmega", &StdI->Nomega, 200);
  StdFace_PrintVal_d("OmegaMax", &StdI->OmegaMax, StdI->LargeValue*StdI->nsite);
  StdFace_PrintVal_d("OmegaMin", &StdI->OmegaMin, -StdI->LargeValue*StdI->nsite);
  StdFace_PrintVal_d("OmegaIm", &StdI->OmegaIm, 0.01* (int)StdI->LargeValue);
#elif defined(_mVMC)
  if (strcmp(StdI->CParaFileHead, "****") == 0) {
    strcpy(StdI->CParaFileHead, "zqp\0");
    fprintf(stdout, "    CParaFileHead = %-12s######  DEFAULT VALUE IS USED  ######\n", StdI->CParaFileHead);
  }
  else fprintf(stdout, "    CParaFileHead = %-s\n", StdI->CParaFileHead);
  
  StdFace_PrintVal_i("NVMCCalMode", &StdI->NVMCCalMode, 0);
  StdFace_PrintVal_i("NLanczosMode", &StdI->NLanczosMode, 0);
  StdFace_PrintVal_i("NDataIdxStart", &StdI->NDataIdxStart, 1);

  if (StdI->NVMCCalMode == 0) StdFace_NotUsed_i("NDataQtySmp", StdI->NDataQtySmp);
  /*else*/StdFace_PrintVal_i("NDataQtySmp", &StdI->NDataQtySmp, 1);

  if (StdI->lGC == 0 && (StdI->Sz2 == 0 || StdI->Sz2 == StdI->NaN_i)) {
    StdFace_PrintVal_i("NSPGaussLeg", &StdI->NSPGaussLeg, 8);
    StdFace_PrintVal_i("NSPStot", &StdI->NSPStot, 0);
  }
  else {
    StdFace_NotUsed_i("NSPGaussLeg", StdI->NSPGaussLeg);
    StdFace_NotUsed_i("NSPStot", StdI->NSPStot);
  }
 
  if (StdI->AntiPeriod[0] == 1 || StdI->AntiPeriod[1] == 1 || StdI->AntiPeriod[2] == 2)
    StdFace_PrintVal_i("NMPTrans", &StdI->NMPTrans, -1);
  else StdFace_PrintVal_i("NMPTrans", &StdI->NMPTrans, 1);

  StdFace_PrintVal_i("NSROptItrStep", &StdI->NSROptItrStep, 1000);
  
  if (StdI->NVMCCalMode == 1) StdFace_NotUsed_i("NSROptItrSmp", StdI->NSROptItrSmp);
  /*else*/ StdFace_PrintVal_i("NSROptItrSmp", &StdI->NSROptItrSmp, StdI->NSROptItrStep/10);

  StdFace_PrintVal_i("NVMCWarmUp", &StdI->NVMCWarmUp, 10);
  StdFace_PrintVal_i("NVMCInterval", &StdI->NVMCInterval, 1);
  StdFace_PrintVal_i("NVMCSample", &StdI->NVMCSample, 1000);

  if (strcmp(StdI->model, "hubbard") == 0) StdI->NExUpdatePath = 0;
  else if (strcmp(StdI->model, "spin") == 0) StdI->NExUpdatePath = 2;
  else if (strcmp(StdI->model, "kondo") == 0) { 
    if(StdI->lGC==0) StdI->NExUpdatePath = 1; 
    else StdI->NExUpdatePath = 3;
  }
  fprintf(stdout, "  %15s = %-10d\n", "NExUpdatePath", StdI->NExUpdatePath);

  StdFace_PrintVal_i("RndSeed", &StdI->RndSeed, 123456789);
  StdFace_PrintVal_i("NSplitSize", &StdI->NSplitSize, 1);
  StdFace_PrintVal_i("NStore", &StdI->NStore, 1);
  StdFace_PrintVal_i("NSRCG", &StdI->NSRCG, 0);

  StdFace_PrintVal_d("DSROptRedCut", &StdI->DSROptRedCut, 0.001);
  StdFace_PrintVal_d("DSROptStaDel", &StdI->DSROptStaDel, 0.02);
  StdFace_PrintVal_d("DSROptStepDt", &StdI->DSROptStepDt, 0.02);
#endif
  /*
   (Un)Conserved variables (Number of electrons, total Sz)
  */
  if (strcmp(StdI->model, "hubbard") == 0){
#if defined(_HPhi)
    if (StdI->lGC == 0) StdFace_RequiredVal_i("nelec", StdI->nelec);
    else {
      StdFace_NotUsed_i("nelec", StdI->nelec);
      StdFace_NotUsed_i("2Sz", StdI->Sz2);
    }
#else
    StdFace_RequiredVal_i("nelec", StdI->nelec);
    if (StdI->lGC == 0) StdFace_PrintVal_i("2Sz", &StdI->Sz2, 0);
    else StdFace_NotUsed_i("2Sz", StdI->Sz2);
#endif
  }
  else if (strcmp(StdI->model, "spin") == 0) {
    StdFace_NotUsed_i("nelec", StdI->nelec);
#if defined(_mVMC)
    StdI->nelec = 0;
#endif
    if (StdI->lGC == 0) StdFace_RequiredVal_i("2Sz", StdI->Sz2);
    else StdFace_NotUsed_i("2Sz", StdI->Sz2);
  }/*else if (strcmp(StdI->model, "spin") == 0)*/
  else if (strcmp(StdI->model, "kondo") == 0) {
#if defined(_HPhi)
    if (StdI->lGC == 0) StdFace_RequiredVal_i("nelec", StdI->nelec);
    else {
      StdFace_NotUsed_i("nelec", StdI->nelec);
      StdFace_NotUsed_i("2Sz", StdI->Sz2);
    }
#else
    StdFace_RequiredVal_i("nelec", StdI->nelec);
    if (StdI->lGC == 0) StdFace_PrintVal_i("2Sz", &StdI->Sz2, 0);
    else StdFace_NotUsed_i("2Sz", StdI->Sz2);
#endif
  }/*else if (strcmp(StdI->model, "kondo") == 0)*/
}/*static void CheckModPara*/
static void PrintInteractions(struct StdIntList *StdI)
{
  FILE *fp;
  int nintr0, kintr, jintr;
  /*
   Coulomb INTRA
  */
  for (kintr = 0; kintr < StdI->NCintra; kintr++) {
    for (jintr = kintr + 1; jintr < StdI->NCintra; jintr++) 
      if(StdI->CintraIndx[jintr][0] == StdI->CintraIndx[kintr][0])
      {
        StdI->Cintra[kintr] += StdI->Cintra[jintr];
        StdI->Cintra[jintr] = 0.0;
      }
  }
  nintr0 = 0;
  for (kintr = 0; kintr < StdI->NCintra; kintr++) {
    if (fabs(StdI->Cintra[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->LCintra = 0;
  else StdI->LCintra = 1;

  if (StdI->LCintra == 1) {
    fp = fopen("coulombintra.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NCoulombIntra %10d\n", nintr0);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "================== CoulombIntra ================\n");
    fprintf(fp, "=============================================\n");
    for (kintr = 0; kintr < StdI->NCintra; kintr++) {
      if (fabs(StdI->Cintra[kintr]) > 0.000001)
        fprintf(fp, "%5d %25.15f\n",
          StdI->CintraIndx[kintr][0], StdI->Cintra[kintr]);
    }
    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    coulombintra.def is written.\n");
  }/*if (StdI->LCintra == 1)*/
  /*
  Coulomb INTER
  */
  for (kintr = 0; kintr < StdI->NCinter; kintr++) {
    for (jintr = kintr + 1; jintr < StdI->NCinter; jintr++)
      if (
        (    StdI->CinterIndx[jintr][0] == StdI->CinterIndx[kintr][0]
          && StdI->CinterIndx[jintr][1] == StdI->CinterIndx[kintr][1])
        ||
        (    StdI->CinterIndx[jintr][0] == StdI->CinterIndx[kintr][1]
          && StdI->CinterIndx[jintr][1] == StdI->CinterIndx[kintr][0])
        )
      {
        StdI->Cinter[kintr] += StdI->Cinter[jintr];
        StdI->Cinter[jintr] = 0.0;
      }
  }/*for (kintr = 0; kintr < StdI->NCinter; kintr++)*/
  nintr0 = 0;
  for (kintr = 0; kintr < StdI->NCinter; kintr++) {
    if (fabs(StdI->Cinter[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->LCinter = 0;
  else StdI->LCinter = 1;

  if (StdI->LCinter == 1) {
    fp = fopen("coulombinter.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NCoulombInter %10d\n", nintr0);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "================== CoulombInter ================\n");
    fprintf(fp, "=============================================\n");
    for (kintr = 0; kintr < StdI->NCinter; kintr++) {
      if (fabs(StdI->Cinter[kintr]) > 0.000001)
        fprintf(fp, "%5d %5d %25.15f\n",
          StdI->CinterIndx[kintr][0], StdI->CinterIndx[kintr][1], StdI->Cinter[kintr]);
    }
    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    coulombinter.def is written.\n");
  }/*if (StdI->LCinter == 1)*/
  /*
  Hund
  */
  for (kintr = 0; kintr < StdI->NHund; kintr++) {
    for (jintr = kintr + 1; jintr < StdI->NHund; jintr++)
      if (
        (StdI->HundIndx[jintr][0] == StdI->HundIndx[kintr][0]
          && StdI->HundIndx[jintr][1] == StdI->HundIndx[kintr][1])
        ||
        (StdI->HundIndx[jintr][0] == StdI->HundIndx[kintr][1]
          && StdI->HundIndx[jintr][1] == StdI->HundIndx[kintr][0])
        )
      {
        StdI->Hund[kintr] += StdI->Hund[jintr];
        StdI->Hund[jintr] = 0.0;
      }
  }/*for (kintr = 0; kintr < StdI->NHund; kintr++)*/
  nintr0 = 0;
  for (kintr = 0; kintr < StdI->NHund; kintr++) {
    if (fabs(StdI->Hund[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->LHund = 0;
  else StdI->LHund = 1;

  if (StdI->LHund == 1) {
    fp = fopen("hund.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NHund %10d\n", nintr0);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "=============== Hund coupling ===============\n");
    fprintf(fp, "=============================================\n");
    for (kintr = 0; kintr < StdI->NHund; kintr++) {
      if (fabs(StdI->Hund[kintr]) > 0.000001)
        fprintf(fp, "%5d %5d %25.15f\n",
          StdI->HundIndx[kintr][0], StdI->HundIndx[kintr][1], StdI->Hund[kintr]);
    }
    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    hund.def is written.\n");
  }/*if (StdI->LHund == 1)*/
  /*
  Exchange
  */
  for (kintr = 0; kintr < StdI->NEx; kintr++) {
    for (jintr = kintr + 1; jintr < StdI->NEx; jintr++)
      if (
        (StdI->ExIndx[jintr][0] == StdI->ExIndx[kintr][0]
          && StdI->ExIndx[jintr][1] == StdI->ExIndx[kintr][1])
        ||
        (StdI->ExIndx[jintr][0] == StdI->ExIndx[kintr][1]
          && StdI->ExIndx[jintr][1] == StdI->ExIndx[kintr][0])
        )
      {
        StdI->Ex[kintr] += StdI->Ex[jintr];
        StdI->Ex[jintr] = 0.0;
      }
  }/*for (kintr = 0; kintr < StdI->NEx; kintr++)*/
  nintr0 = 0;
  for (kintr = 0; kintr < StdI->NEx; kintr++) {
    if (fabs(StdI->Ex[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->LEx = 0;
  else StdI->LEx = 1;

  if (StdI->LEx == 1) {
    fp = fopen("exchange.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NExchange %10d\n", nintr0);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "====== ExchangeCoupling coupling ============\n");
    fprintf(fp, "=============================================\n");
    for (kintr = 0; kintr < StdI->NEx; kintr++) {
      if (fabs(StdI->Ex[kintr]) > 0.000001)
        fprintf(fp, "%5d %5d %25.15f\n",
          StdI->ExIndx[kintr][0], StdI->ExIndx[kintr][1], StdI->Ex[kintr]);
    }
    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    exchange.def is written.\n");
  }
  /*
    PairLift
  */
  for (kintr = 0; kintr < StdI->NPairLift; kintr++) {
    for (jintr = kintr + 1; jintr < StdI->NPairLift; jintr++)
      if (
        (StdI->PLIndx[jintr][0] == StdI->PLIndx[kintr][0]
          && StdI->PLIndx[jintr][1] == StdI->PLIndx[kintr][1])
        ||
        (StdI->PLIndx[jintr][0] == StdI->PLIndx[kintr][1]
          && StdI->PLIndx[jintr][1] == StdI->PLIndx[kintr][0])
        )
      {
        StdI->PairLift[kintr] += StdI->PairLift[jintr];
        StdI->PairLift[jintr] = 0.0;
      }
  }/*for (kintr = 0; kintr < StdI->NPairLift; kintr++)*/
  nintr0 = 0;
  for (kintr = 0; kintr < StdI->NPairLift; kintr++) {
    if (fabs(StdI->PairLift[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->LPairLift = 0;
  else StdI->LPairLift = 1;

  if (StdI->LPairLift == 1) {
    fp = fopen("pairlift.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NPairLift %10d\n", nintr0);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "====== Pair-Lift term ============\n");
    fprintf(fp, "=============================================\n");
    for (kintr = 0; kintr < StdI->NPairLift; kintr++) {
      if (fabs(StdI->PairLift[kintr]) > 0.000001)
        fprintf(fp, "%5d %5d %25.15f\n",
          StdI->PLIndx[kintr][0], StdI->PLIndx[kintr][1], StdI->PairLift[kintr]);
    }
    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    pairlift.def is written.\n");
  }
  /*
  PairHopp
  */
  for (kintr = 0; kintr < StdI->NPairHopp; kintr++) {
    for (jintr = kintr + 1; jintr < StdI->NPairHopp; jintr++)
      if (
        (StdI->PHIndx[jintr][0] == StdI->PHIndx[kintr][0]
          && StdI->PHIndx[jintr][1] == StdI->PHIndx[kintr][1])
        ||
        (StdI->PHIndx[jintr][0] == StdI->PHIndx[kintr][1]
          && StdI->PHIndx[jintr][1] == StdI->PHIndx[kintr][0])
        )
      {
        StdI->PairHopp[kintr] += StdI->PairHopp[jintr];
        StdI->PairHopp[jintr] = 0.0;
      }
  }/*for (kintr = 0; kintr < StdI->NPairHopp; kintr++)*/
  nintr0 = 0;
  for (kintr = 0; kintr < StdI->NPairHopp; kintr++) {
    if (fabs(StdI->PairHopp[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->LPairHopp = 0;
  else StdI->LPairHopp = 1;

  if (StdI->LPairHopp == 1) {
    fp = fopen("pairhopp.def", "w");
    fprintf(fp, "=============================================\n");
    fprintf(fp, "NPairHopp %10d\n", nintr0);
    fprintf(fp, "=============================================\n");
    fprintf(fp, "====== Pair-Hopping term ============\n");
    fprintf(fp, "=============================================\n");
    for (kintr = 0; kintr < StdI->NPairHopp; kintr++) {
      if (fabs(StdI->PairHopp[kintr]) > 0.000001)
        fprintf(fp, "%5d %5d %25.15f\n",
          StdI->PHIndx[kintr][0], StdI->PHIndx[kintr][1], StdI->PairHopp[kintr]);
    }
    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    pairhopp.def is written.\n");
  }
  /*
   InterAll
  */
  for (jintr = 0; jintr < StdI->nintr; jintr++) {
    for (kintr = jintr + 1; kintr < StdI->nintr; kintr++) {
      if (
        (StdI->intrindx[jintr][0] == StdI->intrindx[kintr][0]
          && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][1]
          && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][2]
          && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][3]
          && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][4]
          && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][5]
          && StdI->intrindx[jintr][6] == StdI->intrindx[kintr][6]
          && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][7])
        ||
        (StdI->intrindx[jintr][0] == StdI->intrindx[kintr][4]
          && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][5]
          && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][6]
          && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][7]
          && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][0]
          && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][1]
          && StdI->intrindx[jintr][6] == StdI->intrindx[kintr][2]
          && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][3])
        ) {
        StdI->intr[jintr] = StdI->intr[jintr] + StdI->intr[kintr];
        StdI->intr[kintr] = 0.0;
      }
      else if (
        (StdI->intrindx[jintr][0] == StdI->intrindx[kintr][4]
          && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][5]
          && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][2]
          && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][3]
          && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][0]
          && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][1]
          && StdI->intrindx[jintr][6] == StdI->intrindx[kintr][6]
          && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][7])
        ||
        (StdI->intrindx[jintr][0] == StdI->intrindx[kintr][0]
          && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][1]
          && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][6]
          && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][7]
          && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][4]
          && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][5]
          && StdI->intrindx[jintr][6] == StdI->intrindx[kintr][2]
          && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][3])
        ) {
        StdI->intr[jintr] = StdI->intr[jintr] - StdI->intr[kintr];
        StdI->intr[kintr] = 0.0;
      }
    }/*for (kintr = jintr + 1; kintr < StdI->nintr; kintr++)*/
  }/*for (jintr = 0; jintr < StdI->nintr; jintr++)*/

  for (jintr = 0; jintr < StdI->nintr; jintr++) {
    for (kintr = jintr + 1; kintr < StdI->nintr; kintr++) {
      if (StdI->intrindx[jintr][6] == StdI->intrindx[kintr][4]
        && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][5]
        && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][6]
        && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][7]
        && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][0]
        && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][1]
        && StdI->intrindx[jintr][0] == StdI->intrindx[kintr][2]
        && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][3]
        ) {
        StdI->intrindx[kintr][0] = StdI->intrindx[jintr][6];
        StdI->intrindx[kintr][1] = StdI->intrindx[jintr][7];
        StdI->intrindx[kintr][2] = StdI->intrindx[jintr][4];
        StdI->intrindx[kintr][3] = StdI->intrindx[jintr][5];
        StdI->intrindx[kintr][4] = StdI->intrindx[jintr][2];
        StdI->intrindx[kintr][5] = StdI->intrindx[jintr][3];
        StdI->intrindx[kintr][6] = StdI->intrindx[jintr][0];
        StdI->intrindx[kintr][7] = StdI->intrindx[jintr][1];
      }
      else if (
        (StdI->intrindx[jintr][6] == StdI->intrindx[kintr][4]
          && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][5]
          && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][2]
          && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][3]
          && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][0]
          && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][1]
          && StdI->intrindx[jintr][0] == StdI->intrindx[kintr][6]
          && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][7])
        ||
        (StdI->intrindx[jintr][6] == StdI->intrindx[kintr][0]
          && StdI->intrindx[jintr][7] == StdI->intrindx[kintr][1]
          && StdI->intrindx[jintr][4] == StdI->intrindx[kintr][6]
          && StdI->intrindx[jintr][5] == StdI->intrindx[kintr][7]
          && StdI->intrindx[jintr][2] == StdI->intrindx[kintr][4]
          && StdI->intrindx[jintr][3] == StdI->intrindx[kintr][5]
          && StdI->intrindx[jintr][0] == StdI->intrindx[kintr][2]
          && StdI->intrindx[jintr][1] == StdI->intrindx[kintr][3])
        ) {
        StdI->intrindx[kintr][0] = StdI->intrindx[jintr][6];
        StdI->intrindx[kintr][1] = StdI->intrindx[jintr][7];
        StdI->intrindx[kintr][2] = StdI->intrindx[jintr][4];
        StdI->intrindx[kintr][3] = StdI->intrindx[jintr][5];
        StdI->intrindx[kintr][4] = StdI->intrindx[jintr][2];
        StdI->intrindx[kintr][5] = StdI->intrindx[jintr][3];
        StdI->intrindx[kintr][6] = StdI->intrindx[jintr][0];
        StdI->intrindx[kintr][7] = StdI->intrindx[jintr][1];

        StdI->intr[kintr] = -StdI->intr[kintr];
      }
    }/*for (kintr = jintr + 1; kintr < StdI->nintr; kintr++)*/
  }/*for (jintr = 0; jintr < StdI->nintr; jintr++)*/

  for (jintr = 0; jintr < StdI->nintr; jintr++) {

    if (
      (StdI->intrindx[jintr][0] == StdI->intrindx[jintr][4]
        && StdI->intrindx[jintr][1] == StdI->intrindx[jintr][5]) ||
        (StdI->intrindx[jintr][2] == StdI->intrindx[jintr][6]
          && StdI->intrindx[jintr][3] == StdI->intrindx[jintr][7])
      ) {

      if (!(
        (StdI->intrindx[jintr][0] == StdI->intrindx[jintr][2]
          && StdI->intrindx[jintr][1] == StdI->intrindx[jintr][3])
        ||
        (StdI->intrindx[jintr][0] == StdI->intrindx[jintr][6]
          && StdI->intrindx[jintr][1] == StdI->intrindx[jintr][7])
        ||
        (StdI->intrindx[jintr][4] == StdI->intrindx[jintr][2]
          && StdI->intrindx[jintr][5] == StdI->intrindx[jintr][3])
        ||
        (StdI->intrindx[jintr][4] == StdI->intrindx[jintr][6]
          && StdI->intrindx[jintr][5] == StdI->intrindx[jintr][7])
        ))
        StdI->intr[jintr] = 0.0;
    }
  }/*for (jintr = 0; jintr < StdI->nintr; jintr++)*/

  nintr0 = 0;
  for (kintr = 0; kintr < StdI->nintr; kintr++) {
    if (cabs(StdI->intr[kintr]) > 0.000001) nintr0 = nintr0 + 1;
  }
  if (nintr0 == 0 || StdI->lBoost == 1) StdI->Lintr = 0;
  else StdI->Lintr = 1;

  if (StdI->Lintr == 1) {
    fp = fopen("interall.def", "w");
    fprintf(fp, "====================== \n");
    fprintf(fp, "NInterAll %7d  \n", nintr0);
    fprintf(fp, "====================== \n");
    fprintf(fp, "========zInterAll===== \n");
    fprintf(fp, "====================== \n");

    if (StdI->lBoost == 0) {
      nintr0 = 0;
      for (kintr = 0; kintr < StdI->nintr; kintr++) {
        if (cabs(StdI->intr[kintr]) > 0.000001)
          fprintf(fp, "%5d %5d %5d %5d %5d %5d %5d %5d %25.15f  %25.15f\n",
            StdI->intrindx[kintr][0], StdI->intrindx[kintr][1],
            StdI->intrindx[kintr][2], StdI->intrindx[kintr][3],
            StdI->intrindx[kintr][4], StdI->intrindx[kintr][5],
            StdI->intrindx[kintr][6], StdI->intrindx[kintr][7],
            creal(StdI->intr[kintr]), cimag(StdI->intr[kintr]));
      }/*for (kintr = 0; kintr < StdI->nintr; kintr++)*/
    }/* if (StdI->lBoost == 0)*/

    fflush(fp);
    fclose(fp);
    fprintf(stdout, "    interall.def is written.\n");
  }
}/*static void PrintInteractions*/
void StdFace_main(
  char *fname
)
{
  struct StdIntList *StdI;
  FILE *fp;
  int ktrans, kintr;
  char ctmpline[256];
  char *keyword, *value;

  StdI = (struct StdIntList *)malloc(sizeof(struct StdIntList));

  fprintf(stdout, "\n######  Input Parameter of Standard Intarface  ######\n");
  if ((fp = fopen(fname, "r")) == NULL) {
    fprintf(stdout, "\n  ERROR !  Cannot open input file %s !\n\n", fname);
    StdFace_exit(-1);
  }
  else {
    fprintf(stdout, "\n  Open Standard-Mode Inputfile %s \n\n", fname);
  }

  StdFace_ResetVals(StdI);

  while (fgets(ctmpline, 256, fp) != NULL) {

    TrimSpaceQuote(ctmpline);
    if (strncmp(ctmpline, "//", 2) == 0) {
      fprintf(stdout, "  Skipping a line.\n");
      continue;
    }
    else if (ctmpline[0] == '\0') {
      fprintf(stdout, "  Skipping a line.\n");
      continue;
    }
    keyword = strtok(ctmpline, "=");
    value = strtok(NULL, "=");
    if (value == NULL) {
      fprintf(stdout, "\n  ERROR !  \"=\" is NOT found !\n\n");
      StdFace_exit(-1);
    }
    Text2Lower(keyword);
    fprintf(stdout, "  KEYWORD : %-20s | VALUE : %s \n", keyword, value);

    if (strcmp(keyword, "a") == 0) StoreWithCheckDup_d(keyword, value, &StdI->a);
    else if (strcmp(keyword, "a0h") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[0][2]);
    else if (strcmp(keyword, "a0l") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[0][1]);
    else if (strcmp(keyword, "a0w") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[0][0]);
    else if (strcmp(keyword, "a1h") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[1][2]);
    else if (strcmp(keyword, "a1l") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[1][1]);
    else if (strcmp(keyword, "a1w") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[1][0]);
    else if (strcmp(keyword, "a2h") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[2][2]);
    else if (strcmp(keyword, "a2l") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[2][1]);
    else if (strcmp(keyword, "a2w") == 0) StoreWithCheckDup_i(keyword, value, &StdI->box[2][0]);
    else if (strcmp(keyword, "cutoff_j") == 0) StoreWithCheckDup_d(keyword, value, &StdI->cutoff_j);
    else if (strcmp(keyword, "cutoff_jh") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_JR[2]);
    else if (strcmp(keyword, "cutoff_jl") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_JR[1]);
    else if (strcmp(keyword, "cutoff_jw") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_JR[0]);
    else if (strcmp(keyword, "cutoff_length_j") == 0) StoreWithCheckDup_d(keyword, value, &StdI->cutoff_length_J);
    else if (strcmp(keyword, "cutoff_length_u") == 0) StoreWithCheckDup_d(keyword, value, &StdI->cutoff_length_U);
    else if (strcmp(keyword, "cutoff_length_t") == 0) StoreWithCheckDup_d(keyword, value, &StdI->cutoff_length_t);
    else if (strcmp(keyword, "cutoff_t") == 0) StoreWithCheckDup_d(keyword, value, &StdI->cutoff_t);
    else if (strcmp(keyword, "cutoff_th") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_tR[2]);
    else if (strcmp(keyword, "cutoff_tl") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_tR[1]);
    else if (strcmp(keyword, "cutoff_tw") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_tR[0]);
    else if (strcmp(keyword, "cutoff_u") == 0) StoreWithCheckDup_d(keyword, value, &StdI->cutoff_u);
    else if (strcmp(keyword, "cutoff_uh") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_UR[2]);
    else if (strcmp(keyword, "cutoff_ul") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_UR[1]);
    else if (strcmp(keyword, "cutoff_uw") == 0) StoreWithCheckDup_i(keyword, value, &StdI->cutoff_UR[0]);
    else if (strcmp(keyword, "d") == 0) StoreWithCheckDup_d(keyword, value, &StdI->D[2][2]);
    else if (strcmp(keyword, "gamma") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Gamma);
    else if (strcmp(keyword, "h") == 0) StoreWithCheckDup_d(keyword, value, &StdI->h);
    else if (strcmp(keyword, "height") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Height);
    else if (strcmp(keyword, "hlength") == 0) StoreWithCheckDup_d(keyword, value, &StdI->length[2]);
    else if (strcmp(keyword, "hx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[2][0]);
    else if (strcmp(keyword, "hy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[2][1]);
    else if (strcmp(keyword, "hz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[2][2]);
    else if (strcmp(keyword, "j") == 0) StoreWithCheckDup_d(keyword, value, &StdI->JAll);
    else if (strcmp(keyword, "jx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[0][0]);
    else if (strcmp(keyword, "jxy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[0][1]);
    else if (strcmp(keyword, "jxz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[0][2]);
    else if (strcmp(keyword, "jy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[1][1]);
    else if (strcmp(keyword, "jyx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[1][0]);
    else if (strcmp(keyword, "jyz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[1][2]);
    else if (strcmp(keyword, "jz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[2][2]);
    else if (strcmp(keyword, "jzx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[2][0]);
    else if (strcmp(keyword, "jzy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J[2][1]);
    else if (strcmp(keyword, "j0") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0All);
    else if (strcmp(keyword, "j0x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[0][0]);
    else if (strcmp(keyword, "j0xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[0][1]);
    else if (strcmp(keyword, "j0xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[0][2]);
    else if (strcmp(keyword, "j0y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[1][1]);
    else if (strcmp(keyword, "j0yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[1][0]);
    else if (strcmp(keyword, "j0yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[1][2]);
    else if (strcmp(keyword, "j0z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[2][2]);
    else if (strcmp(keyword, "j0zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[2][0]);
    else if (strcmp(keyword, "j0zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0[2][1]);
    else if (strcmp(keyword, "j0'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pAll);
    else if (strcmp(keyword, "j0'x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[0][0]);
    else if (strcmp(keyword, "j0'xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[0][1]);
    else if (strcmp(keyword, "j0'xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[0][2]);
    else if (strcmp(keyword, "j0'y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[1][1]);
    else if (strcmp(keyword, "j0'yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[1][0]);
    else if (strcmp(keyword, "j0'yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[1][2]);
    else if (strcmp(keyword, "j0'z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[2][2]);
    else if (strcmp(keyword, "j0'zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[2][0]);
    else if (strcmp(keyword, "j0'zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0p[2][1]);
    else if (strcmp(keyword, "j0''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0ppAll);
    else if (strcmp(keyword, "j0''x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[0][0]);
    else if (strcmp(keyword, "j0''xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[0][1]);
    else if (strcmp(keyword, "j0''xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[0][2]);
    else if (strcmp(keyword, "j0''y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[1][1]);
    else if (strcmp(keyword, "j0''yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[1][0]);
    else if (strcmp(keyword, "j0''yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[1][2]);
    else if (strcmp(keyword, "j0''z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[2][2]);
    else if (strcmp(keyword, "j0''zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[2][0]);
    else if (strcmp(keyword, "j0''zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J0pp[2][1]);
    else if (strcmp(keyword, "j1") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1All);
    else if (strcmp(keyword, "j1x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[0][0]);
    else if (strcmp(keyword, "j1xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[0][1]);
    else if (strcmp(keyword, "j1xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[0][2]);
    else if (strcmp(keyword, "j1y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[1][1]);
    else if (strcmp(keyword, "j1yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[1][0]);
    else if (strcmp(keyword, "j1yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[1][2]);
    else if (strcmp(keyword, "j1z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[2][2]);
    else if (strcmp(keyword, "j1zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[2][0]);
    else if (strcmp(keyword, "j1zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1[2][1]);
    else if (strcmp(keyword, "j1'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pAll);
    else if (strcmp(keyword, "j1'x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[0][0]);
    else if (strcmp(keyword, "j1'xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[0][1]);
    else if (strcmp(keyword, "j1'xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[0][2]);
    else if (strcmp(keyword, "j1'y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[1][1]);
    else if (strcmp(keyword, "j1'yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[1][0]);
    else if (strcmp(keyword, "j1'yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[1][2]);
    else if (strcmp(keyword, "j1'z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[2][2]);
    else if (strcmp(keyword, "j1'zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[2][0]);
    else if (strcmp(keyword, "j1'zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1p[2][1]);
    else if (strcmp(keyword, "j1''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1ppAll);
    else if (strcmp(keyword, "j1''x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[0][0]);
    else if (strcmp(keyword, "j1''xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[0][1]);
    else if (strcmp(keyword, "j1''xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[0][2]);
    else if (strcmp(keyword, "j1''y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[1][1]);
    else if (strcmp(keyword, "j1''yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[1][0]);
    else if (strcmp(keyword, "j1''yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[1][2]);
    else if (strcmp(keyword, "j1''z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[2][2]);
    else if (strcmp(keyword, "j1''zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[2][0]);
    else if (strcmp(keyword, "j1''zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J1pp[2][1]);
    else if (strcmp(keyword, "j2") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2All);
    else if (strcmp(keyword, "j2x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[0][0]);
    else if (strcmp(keyword, "j2xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[0][1]);
    else if (strcmp(keyword, "j2xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[0][2]);
    else if (strcmp(keyword, "j2y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[1][1]);
    else if (strcmp(keyword, "j2yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[1][0]);
    else if (strcmp(keyword, "j2yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[1][2]);
    else if (strcmp(keyword, "j2z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[2][2]);
    else if (strcmp(keyword, "j2zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[2][0]);
    else if (strcmp(keyword, "j2zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2[2][1]);
    else if (strcmp(keyword, "j2'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pAll);
    else if (strcmp(keyword, "j2'x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[0][0]);
    else if (strcmp(keyword, "j2'xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[0][1]);
    else if (strcmp(keyword, "j2'xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[0][2]);
    else if (strcmp(keyword, "j2'y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[1][1]);
    else if (strcmp(keyword, "j2'yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[1][0]);
    else if (strcmp(keyword, "j2'yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[1][2]);
    else if (strcmp(keyword, "j2'z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[2][2]);
    else if (strcmp(keyword, "j2'zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[2][0]);
    else if (strcmp(keyword, "j2'zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2p[2][1]);
    else if (strcmp(keyword, "j2''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2ppAll);
    else if (strcmp(keyword, "j2''x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[0][0]);
    else if (strcmp(keyword, "j2''xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[0][1]);
    else if (strcmp(keyword, "j2''xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[0][2]);
    else if (strcmp(keyword, "j2''y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[1][1]);
    else if (strcmp(keyword, "j2''yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[1][0]);
    else if (strcmp(keyword, "j2''yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[1][2]);
    else if (strcmp(keyword, "j2''z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[2][2]);
    else if (strcmp(keyword, "j2''zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[2][0]);
    else if (strcmp(keyword, "j2''zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->J2pp[2][1]);
    else if (strcmp(keyword, "j'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->JpAll);
    else if (strcmp(keyword, "j'x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[0][0]);
    else if (strcmp(keyword, "j'xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[0][1]);
    else if (strcmp(keyword, "j'xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[0][2]);
    else if (strcmp(keyword, "j'y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[1][1]);
    else if (strcmp(keyword, "j'yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[1][0]);
    else if (strcmp(keyword, "j'yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[1][2]);
    else if (strcmp(keyword, "j'z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[2][2]);
    else if (strcmp(keyword, "j'zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[2][0]);
    else if (strcmp(keyword, "j'zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jp[2][1]);
    else if (strcmp(keyword, "j''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->JppAll);
    else if (strcmp(keyword, "j''x") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[0][0]);
    else if (strcmp(keyword, "j''xy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[0][1]);
    else if (strcmp(keyword, "j''xz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[0][2]);
    else if (strcmp(keyword, "j''y") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[1][1]);
    else if (strcmp(keyword, "j''yx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[1][0]);
    else if (strcmp(keyword, "j''yz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[1][2]);
    else if (strcmp(keyword, "j''z") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[2][2]);
    else if (strcmp(keyword, "j''zx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[2][0]);
    else if (strcmp(keyword, "j''zy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Jpp[2][1]);
    else if (strcmp(keyword, "k") == 0) StoreWithCheckDup_d(keyword, value, &StdI->K);
    else if (strcmp(keyword, "l") == 0) StoreWithCheckDup_i(keyword, value, &StdI->L);
    else if (strcmp(keyword, "lattice") == 0) StoreWithCheckDup_sl(keyword, value, StdI->lattice);
    else if (strcmp(keyword, "llength") == 0) StoreWithCheckDup_d(keyword, value, &StdI->length[1]);
    else if (strcmp(keyword, "lx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[1][0]);
    else if (strcmp(keyword, "ly") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[1][1]);
    else if (strcmp(keyword, "lz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[1][2]);
    else if (strcmp(keyword, "model") == 0) StoreWithCheckDup_sl(keyword, value, StdI->model);
    else if (strcmp(keyword, "mu") == 0) StoreWithCheckDup_d(keyword, value, &StdI->mu);
    else if (strcmp(keyword, "nelec") == 0) StoreWithCheckDup_i(keyword, value, &StdI->nelec);
    else if (strcmp(keyword, "outputmode") == 0) StoreWithCheckDup_sl(keyword, value, StdI->outputmode);
    else if (strcmp(keyword, "phase0") == 0) StoreWithCheckDup_d(keyword, value, &StdI->phase[0]);
    else if (strcmp(keyword, "phase1") == 0) StoreWithCheckDup_d(keyword, value, &StdI->phase[1]);
    else if (strcmp(keyword, "phase2") == 0) StoreWithCheckDup_d(keyword, value, &StdI->phase[2]);
    else if (strcmp(keyword, "t") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t);
    else if (strcmp(keyword, "t0") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t0);
    else if (strcmp(keyword, "t0'") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t0p);
    else if (strcmp(keyword, "t0''") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t0pp);
    else if (strcmp(keyword, "t1") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t1);
    else if (strcmp(keyword, "t1'") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t1p);
    else if (strcmp(keyword, "t1''") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t1pp);
    else if (strcmp(keyword, "t2") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t2);
    else if (strcmp(keyword, "t2'") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t2p);
    else if (strcmp(keyword, "t2''") == 0) StoreWithCheckDup_c(keyword, value, &StdI->t2pp);
    else if (strcmp(keyword, "t'") == 0) StoreWithCheckDup_c(keyword, value, &StdI->tp);
    else if (strcmp(keyword, "t''") == 0) StoreWithCheckDup_c(keyword, value, &StdI->tpp);
    else if (strcmp(keyword, "u") == 0) StoreWithCheckDup_d(keyword, value, &StdI->U);
    else if (strcmp(keyword, "v") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V);
    else if (strcmp(keyword, "v0") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V0);
    else if (strcmp(keyword, "v0'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V0p);
    else if (strcmp(keyword, "v0''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V0pp);
    else if (strcmp(keyword, "v1") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V1);
    else if (strcmp(keyword, "v1'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V1p);
    else if (strcmp(keyword, "v1''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V1pp);
    else if (strcmp(keyword, "v2") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V2);
    else if (strcmp(keyword, "v2'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V2p);
    else if (strcmp(keyword, "v2''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->V2pp);
    else if (strcmp(keyword, "v'") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Vp);
    else if (strcmp(keyword, "v''") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Vpp);
    else if (strcmp(keyword, "w") == 0) StoreWithCheckDup_i(keyword, value, &StdI->W);
    else if (strcmp(keyword, "wlength") == 0) StoreWithCheckDup_d(keyword, value, &StdI->length[0]);
    else if (strcmp(keyword, "wx") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[0][0]);
    else if (strcmp(keyword, "wy") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[0][1]);
    else if (strcmp(keyword, "wz") == 0) StoreWithCheckDup_d(keyword, value, &StdI->direct[0][2]);
    else if (strcmp(keyword, "2sz") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Sz2);

#if defined(_HPhi)
    else if (strcmp(keyword, "calcspec") == 0) StoreWithCheckDup_sl(keyword, value, StdI->CalcSpec);
    else if (strcmp(keyword, "exct") == 0) StoreWithCheckDup_i(keyword, value, &StdI->exct);
    else if (strcmp(keyword, "eigenvecio") == 0) StoreWithCheckDup_sl(keyword, value, StdI->EigenVecIO);
    else if (strcmp(keyword, "expandcoef") == 0) StoreWithCheckDup_i(keyword, value, &StdI->ExpandCoef);
    else if (strcmp(keyword, "expecinterval") == 0) StoreWithCheckDup_i(keyword, value, &StdI->ExpecInterval);
    else if (strcmp(keyword, "cdatafilehead") == 0) StoreWithCheckDup_s(keyword, value, StdI->CDataFileHead);
    else if (strcmp(keyword, "dt") == 0) StoreWithCheckDup_d(keyword, value, &StdI->dt);
    else if (strcmp(keyword, "flgtemp") == 0) StoreWithCheckDup_i(keyword, value, &StdI->FlgTemp);
    else if (strcmp(keyword, "freq") == 0) StoreWithCheckDup_d(keyword, value, &StdI->freq);
    else if (strcmp(keyword, "hamio") == 0) StoreWithCheckDup_sl(keyword, value, StdI->HamIO);
    else if (strcmp(keyword, "initialvectype") == 0) StoreWithCheckDup_sl(keyword, value, StdI->InitialVecType);
    else if (strcmp(keyword, "initial_iv") == 0) StoreWithCheckDup_i(keyword, value, &StdI->initial_iv);
    else if (strcmp(keyword, "lanczoseps") == 0) StoreWithCheckDup_i(keyword, value, &StdI->LanczosEps);
    else if (strcmp(keyword, "lanczostarget") == 0) StoreWithCheckDup_i(keyword, value, &StdI->LanczosTarget);
    else if (strcmp(keyword, "lanczos_max") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Lanczos_max);
    else if (strcmp(keyword, "largevalue") == 0) StoreWithCheckDup_d(keyword, value, &StdI->LargeValue);
    else if (strcmp(keyword, "method") == 0) StoreWithCheckDup_sl(keyword, value, StdI->method);
    else if (strcmp(keyword, "nomega") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Nomega);
    else if (strcmp(keyword, "numave") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NumAve);
    else if (strcmp(keyword, "nvec") == 0) StoreWithCheckDup_i(keyword, value, &StdI->nvec);
    else if (strcmp(keyword, "omegamax") == 0) StoreWithCheckDup_d(keyword, value, &StdI->OmegaMax);
    else if (strcmp(keyword, "omegamin") == 0) StoreWithCheckDup_d(keyword, value, &StdI->OmegaMin);
    else if (strcmp(keyword, "omegaim") == 0) StoreWithCheckDup_d(keyword, value, &StdI->OmegaIm);
    else if (strcmp(keyword, "outputexcitedvec") == 0) StoreWithCheckDup_sl(keyword, value, StdI->OutputExVec);
    else if (strcmp(keyword, "pumptype") == 0) StoreWithCheckDup_sl(keyword, value, StdI->PumpType);
    else if (strcmp(keyword, "restart") == 0) StoreWithCheckDup_sl(keyword, value, StdI->Restart);
    else if (strcmp(keyword, "spectrumqh") == 0) StoreWithCheckDup_d(keyword, value, &StdI->SpectrumQ[2]);
    else if (strcmp(keyword, "spectrumql") == 0) StoreWithCheckDup_d(keyword, value, &StdI->SpectrumQ[1]);
    else if (strcmp(keyword, "spectrumqw") == 0) StoreWithCheckDup_d(keyword, value, &StdI->SpectrumQ[0]);
    else if (strcmp(keyword, "spectrumtype") == 0) StoreWithCheckDup_sl(keyword, value, StdI->SpectrumType);
    else if (strcmp(keyword, "tdump") == 0) StoreWithCheckDup_d(keyword, value, &StdI->tdump);
    else if (strcmp(keyword, "tshift") == 0) StoreWithCheckDup_d(keyword, value, &StdI->tshift);
    else if (strcmp(keyword, "uquench") == 0) StoreWithCheckDup_d(keyword, value, &StdI->Uquench);
    else if (strcmp(keyword, "vecpoth") == 0) StoreWithCheckDup_d(keyword, value, &StdI->VecPot[2]);
    else if (strcmp(keyword, "vecpotl") == 0) StoreWithCheckDup_d(keyword, value, &StdI->VecPot[1]);
    else if (strcmp(keyword, "vecpotw") == 0) StoreWithCheckDup_d(keyword, value, &StdI->VecPot[0]);
    else if (strcmp(keyword, "2s") == 0) StoreWithCheckDup_i(keyword, value, &StdI->S2);
#elif defined(_mVMC)
    else if (strcmp(keyword, "a0hsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[0][2]);
    else if (strcmp(keyword, "a0lsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[0][1]);
    else if (strcmp(keyword, "a0wsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[0][0]);
    else if (strcmp(keyword, "a1hsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[1][2]);
    else if (strcmp(keyword, "a1lsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[1][1]);
    else if (strcmp(keyword, "a1wsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[1][0]);
    else if (strcmp(keyword, "a2hsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[2][2]);
    else if (strcmp(keyword, "a2lsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[2][1]);
    else if (strcmp(keyword, "a2wsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->boxsub[2][0]);
    else if (strcmp(keyword, "complextype") == 0) StoreWithCheckDup_i(keyword, value, &StdI->ComplexType);
    else if (strcmp(keyword, "cparafilehead") == 0) StoreWithCheckDup_s(keyword, value, StdI->CParaFileHead);
    else if (strcmp(keyword, "dsroptredcut") == 0) StoreWithCheckDup_d(keyword, value, &StdI->DSROptRedCut);
    else if (strcmp(keyword, "dsroptstadel") == 0) StoreWithCheckDup_d(keyword, value, &StdI->DSROptStaDel);
    else if (strcmp(keyword, "dsroptstepdt") == 0) StoreWithCheckDup_d(keyword, value, &StdI->DSROptStepDt);
    else if (strcmp(keyword, "hsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Hsub);
    else if (strcmp(keyword, "lsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Lsub);
    else if (strcmp(keyword, "nvmccalmode") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NVMCCalMode);
    else if (strcmp(keyword, "ndataidxstart") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NDataIdxStart);
    else if (strcmp(keyword, "ndataqtysmp") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NDataQtySmp);
    else if (strcmp(keyword, "nlanczosmode") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NLanczosMode);
    else if (strcmp(keyword, "nmptrans") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NMPTrans);
    else if (strcmp(keyword, "nspgaussleg") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NSPGaussLeg);
    else if (strcmp(keyword, "nsplitsize") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NSplitSize);
    else if (strcmp(keyword, "nspstot") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NSPStot);
    else if (strcmp(keyword, "nsroptitrsmp") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NSROptItrSmp);
    else if (strcmp(keyword, "nsroptitrstep") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NSROptItrStep);
    else if (strcmp(keyword, "nstore") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NStore);
    else if (strcmp(keyword, "nsrcg") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NSRCG);
    else if (strcmp(keyword, "nvmcinterval") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NVMCInterval);
    else if (strcmp(keyword, "nvmcsample") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NVMCSample);
    else if (strcmp(keyword, "nvmcwarmup") == 0) StoreWithCheckDup_i(keyword, value, &StdI->NVMCWarmUp);
    else if (strcmp(keyword, "rndseed") == 0) StoreWithCheckDup_i(keyword, value, &StdI->RndSeed);
    else if (strcmp(keyword, "wsub") == 0) StoreWithCheckDup_i(keyword, value, &StdI->Wsub);
#endif
    else {
      fprintf(stdout, "ERROR ! Unsupported Keyword in Standard mode!\n");
      StdFace_exit(-1);
    }
  }
  fflush(fp);
  fclose(fp);
  fprintf(stdout, "\n");
  fprintf(stdout, "#######  Construct Model  #######\n");
  fprintf(stdout, "\n");
  /*
  Check the model
  */
  if (strcmp(StdI->CDataFileHead, "****") == 0) {
    strcpy(StdI->CDataFileHead, "zvo\0");
    fprintf(stdout, "    CDataFileHead = %-12s######  DEFAULT VALUE IS USED  ######\n", StdI->CDataFileHead);
  }
  else fprintf(stdout, "    CDataFileHead = %-s\n", StdI->CDataFileHead);
  
  StdI->lGC = 0;
  StdI->lBoost = 0;
  if (strcmp(StdI->model, "fermionhubbard") == 0
    || strcmp(StdI->model, "hubbard") == 0)
    strcpy(StdI->model, "hubbard\0");
  else if(strcmp(StdI->model, "fermionhubbardgc") == 0
    || strcmp(StdI->model, "hubbardgc") == 0) {
    strcpy(StdI->model, "hubbard\0");
    StdI->lGC = 1;
  }
  else if (strcmp(StdI->model, "spin") == 0)
    strcpy(StdI->model, "spin\0");
  else if (strcmp(StdI->model, "spingc") == 0) {
    strcpy(StdI->model, "spin\0");
    StdI->lGC = 1;
  }
#if defined(_HPhi)
  else if(strcmp(StdI->model, "spingcboost") == 0 ||
    strcmp(StdI->model, "spingccma") == 0) {
    strcpy(StdI->model, "spin\0");
    StdI->lGC = 1;
    StdI->lBoost = 1;
  }
#endif
  else if (strcmp(StdI->model, "kondolattice") == 0
    || strcmp(StdI->model, "kondo") == 0) {
    strcpy(StdI->model, "kondo\0");
  }
  else if(strcmp(StdI->model, "kondolatticegc") == 0
    || strcmp(StdI->model, "kondogc") == 0) {
    strcpy(StdI->model, "kondo\0");
    StdI->lGC = 1;
  }
  else UnsupportedSystem(StdI->model, StdI->lattice);
#if defined(_HPhi)
  /*
  Check the method
  */
  if (strcmp(StdI->method, "direct") == 0
    || strcmp(StdI->method, "alldiag") == 0)
    strcpy(StdI->method, "fulldiag\0");
  else if (strcmp(StdI->method, "te") == 0
    || strcmp(StdI->method, "time-evolution") == 0) {
    strcpy(StdI->method, "timeevolution\0");
  }
  /*
  Compute vector potential and electrical field
  */
  if (strcmp(StdI->method, "timeevolution") == 0) VectorPotential(StdI);
#endif
  /*>>
  Generate Hamiltonian definition files
  */
  if (strcmp(StdI->lattice, "chain") == 0
    || strcmp(StdI->lattice, "chainlattice") == 0) StdFace_Chain(StdI);
  else if (strcmp(StdI->lattice, "face-centeredorthorhombic") == 0
    || strcmp(StdI->lattice, "fcorthorhombic") == 0
    || strcmp(StdI->lattice, "fco") == 0) StdFace_FCOrtho(StdI);
  else if (strcmp(StdI->lattice, "honeycomb") == 0
    || strcmp(StdI->lattice, "honeycomblattice") == 0) StdFace_Honeycomb(StdI);
  else if (strcmp(StdI->lattice, "kagome") == 0
    || strcmp(StdI->lattice, "kagomelattice") == 0) StdFace_Kagome(StdI);
  else if (strcmp(StdI->lattice, "ladder") == 0
    || strcmp(StdI->lattice, "ladderlattice") == 0) StdFace_Ladder(StdI);
  else if (strcmp(StdI->lattice, "orthorhombic") == 0
    || strcmp(StdI->lattice, "simpleorthorhombic") == 0) StdFace_Orthorhombic(StdI);
  else if (strcmp(StdI->lattice, "pyrochlore") == 0) StdFace_Pyrochlore(StdI);
  else if (strcmp(StdI->lattice, "tetragonal") == 0
    || strcmp(StdI->lattice, "tetragonallattice") == 0
    || strcmp(StdI->lattice, "square") == 0
    || strcmp(StdI->lattice, "squarelattice") == 0) StdFace_Tetragonal(StdI);
  else if (strcmp(StdI->lattice, "triangular") == 0
    || strcmp(StdI->lattice, "triangularlattice") == 0) StdFace_Triangular(StdI);
  else if (strcmp(StdI->lattice, "wannier90") == 0) StdFace_Wannier90(StdI);
  else UnsupportedSystem(StdI->model, StdI->lattice);//<<
  
#if defined(_HPhi)
  StdFace_LargeValue(StdI);
  /*
  Generate Hamiltonian for Boost
  */
  if (StdI->lBoost == 1) {
    if (strcmp(StdI->lattice, "chain") == 0
      || strcmp(StdI->lattice, "chainlattice") == 0) StdFace_Chain_Boost(StdI);
    else if (strcmp(StdI->lattice, "honeycomb") == 0
      || strcmp(StdI->lattice, "honeycomblattice") == 0) StdFace_Honeycomb_Boost(StdI);
    else if (strcmp(StdI->lattice, "kagome") == 0
      || strcmp(StdI->lattice, "kagomelattice") == 0) StdFace_Kagome_Boost(StdI);
    else if (strcmp(StdI->lattice, "ladder") == 0
      || strcmp(StdI->lattice, "ladderlattice") == 0) StdFace_Ladder_Boost(StdI);
    else UnsupportedSystem(StdI->model, StdI->lattice);
  }
#endif
  
  fprintf(stdout, "\n");
  fprintf(stdout, "######  Print Expert input files  ######\n");
  fprintf(stdout, "\n");
  PrintLocSpin(StdI);
  PrintTrans(StdI);
  PrintInteractions(StdI);
  CheckModPara(StdI);
  PrintModPara(StdI); 
#if defined(_HPhi)
  PrintExcitation(StdI);
  if (strcmp(StdI->method, "timeevolution") == 0) PrintPump(StdI);
  PrintCalcMod(StdI);
#elif defined(_mVMC)

  if(StdI->lGC == 0 && (StdI->Sz2 == 0 || StdI->Sz2 == StdI->NaN_i)) 
    StdFace_PrintVal_i("ComplexType", &StdI->ComplexType, 0);
  else StdFace_PrintVal_i("ComplexType", &StdI->ComplexType, 1);

  StdFace_generate_orb(StdI);
  StdFace_Proj(StdI);
  PrintJastrow(StdI);
  if(StdI->lGC == 1 || (StdI->Sz2 != 0 && StdI->Sz2 != StdI->NaN_i) )
    PrintOrbPara(StdI);
  PrintGutzwiller(StdI);
  PrintOrb(StdI);
#endif
  CheckOutputMode(StdI);
  Print1Green(StdI);
  Print2Green(StdI);
  PrintNamelist(StdI);
  /*
  Finalize All
  */
  free(StdI->locspinflag);
  for (ktrans = 0; ktrans < StdI->ntrans; ktrans++) {
    free(StdI->transindx[ktrans]);
  }
  free(StdI->transindx);
  free(StdI->trans);
  for (kintr = 0; kintr < StdI->nintr; kintr++) {
    free(StdI->intrindx[kintr]);
  }
  free(StdI->intrindx);
  free(StdI->intr);

  fprintf(stdout, "\n######  Input files are generated.  ######\n\n");
  free(StdI);
}/*void StdFace_main*/