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[2];
  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:iruv: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   if (!loadedobs) {
  40     for (int kpiruv=0; kpiruv<2; kpiruv++) {
  41       double ppt = 2.*M_PI/O815->comargs.lsize[0] * ( 0 + kpiruv * O815->comargs.lsize[0]/2 ) ;
  42       double ppx = 2.*M_PI/O815->comargs.lsize[1] * ( 0 + kpiruv * O815->comargs.lsize[1]/2 );
  43
  44       OM->phip2[ kpiruv ] = 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[ kpiruv ] += 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[ kpiruv ] /= Sim->LSIZE2;
  55     }
  56   }
  57
  58   for (int kpiruv=0; kpiruv<2; kpiruv++) {
  59     double ppt = 2.*M_PI/O815->comargs.lsize[0] * ( 0 + kpiruv * O815->comargs.lsize[0]/2 ) ;
  60     double ppx = 2.*M_PI/O815->comargs.lsize[1] * ( 0 + kpiruv * O815->comargs.lsize[1]/2 );
  61
  62     *out << O815->paraQ->getParaVals();
  63     *out << "\t" << ( Sim->nequi + (nthmeas+1)*Sim->nskip );
  64     *out << "\t" << kpiruv;
  65     *out << "\t" << abs( OM->phip2[ kpiruv ] - 1./(ppt*ppt + ppx*ppx + Sim->m*Sim->m) ) << endl;
  66   }
  67 };
  68
  69 #endif