obs_phip2_hist.hpp

   1 #ifndef OBS_PHIP2_HIST_HPP
   2 #define OBS_PHIP2_HIST_HPP
   3
   4 #include "latlib/o815/o815.h"
   5
   6 #include "latlib/writeout.h"
   7
   8 #include "latlib/obstat.hpp"
   9
  10 #include <iostream>
  11
  12 using namespace std;
  13
  14 class obs_phip2_hist : public o815::obs {
  15
  16 public:
  17   struct obsmem {
  18     double phip2[4];
  19   };
  20   obs_phip2_hist(o815 *_O815);
  21
  22 private:
  23   void _start() {}
  24   void _meas(bool loadedobs, const int& nthmeas);
  25   void _finish() {}
  26   obsmem *OM;
  27   sim *Sim;
  28 };
  29
  30 obs_phip2_hist::obs_phip2_hist(o815 *_O815) : o815::obs("phip2",
  31                                                         _O815->paraQ->getParaNames() + "nthstep:imode:kt:kx:phip2absdiffequi", _O815,
  32                                                         sizeof(obsmem),
  33                                                         "_hist") {
  34   OM = (obsmem*)obsMem;
  35   Sim = (sim*)O815->Sim;
  36 }
  37
  38 void obs_phip2_hist::_meas(bool loadedobs, const int& nthmeas) {
  39   for (int kpimode=0; kpimode<4; kpimode++) {
  40     const double ppt = 2.*M_PI/O815->comargs.lsize[0] * ( 0 + int( kpimode/4. * O815->comargs.lsize[0] ) ) ;
  41     const double ppx = 2.*M_PI/O815->comargs.lsize[1] * ( 0 + int( kpimode/4. * O815->comargs.lsize[1] ) );
  42
  43     if (!loadedobs) {
  44       OM->phip2[ kpimode ] = 0;
  45
  46       for (int ixt = 0; ixt < O815->comargs.lsize[0]; ixt++)
  47         for (int ixx = 0; ixx < O815->comargs.lsize[1]; ixx++)
  48           for (int ixpt = 0; ixpt < O815->comargs.lsize[0]; ixpt++)
  49             for (int ixpx = 0; ixpx < O815->comargs.lsize[1]; ixpx++)
  50               OM->phip2[ kpimode ] += real( conj( Sim->conf[ ixt*O815->comargs.lsize[1] + ixx ].phi )
  51                                             * Sim->conf[ ixpt*O815->comargs.lsize[1] + ixpx ].phi
  52                                             * exp ( _i_*(double)ppx*(double)(ixx-ixpx) + _i_*(double)ppt*(double)(ixt-ixpt) ) );
  53
  54       OM->phip2[ kpimode ] /= Sim->LSIZE2;
  55     }
  56
  57     *out << O815->paraQ->getParaVals();
  58     *out << "\t" << ( Sim->nequi + (nthmeas+1)*Sim->nskip );
  59     *out << "\t" << kpimode;
  60     *out << "\t" << 0 + int( kpimode/4. * O815->comargs.lsize[0] );
  61     *out << "\t" << 0 + int( kpimode/4. * O815->comargs.lsize[1] );
  62     *out << "\t" << abs( OM->phip2[ kpimode ] - 1./(ppt*ppt + ppx*ppx + Sim->m*Sim->m) ) << endl;
  63   }
  64 };
  65
  66 #endif