TraceBuffer< T > Struct Template Reference

Inheritance diagram for TraceBuffer< T >:

Collaboration diagram for TraceBuffer< T >:

Public Member Functions
	TraceBuffer ()
void	set_verbose (bool v)
void	set_type (std::string stype)
void	reset ()
void	set_samples (TraceSamples< T > &s)
void	set_samples (TraceSamples< std::complex< T >> &s)
void	make_combined_trace (const std::string &name, std::vector< std::string > &names, std::vector< TraceReal > &weights)
void	order_and_resize ()
bool	same_as (TraceBuffer< T > &ref)
void	collect_sample ()
void	write ()
void	write_summary (std::string pad=" ")
void	user_report (const std::string &pad=" ")
void	register_hdf_data (hdf_archive &f)
void	write_hdf (hdf_archive &f)
void	write_hdf (hdf_archive &f, hsize_t &file_pointer)
void	test_buffer_write (int sample_size)
void	test_buffer_collect (int current_row)

Public Attributes
bool	has_complex
TraceSamples< T > *	samples
TraceSamples< std::complex< T > > *	complex_samples
std::string	type
Array< T, 2 >	buffer
bool	verbose
std::string	top
hsize_t	dims [2]
hsize_t	hdf_file_pointer

Detailed Description

template<typename T>
struct qmcplusplus::TraceBuffer< T >

Definition at line 1086 of file TraceManager.h.

Constructor & Destructor Documentation

TraceBuffer()

TraceBuffer ( )

inline

Definition at line 1101 of file TraceManager.h.

                : samples(0), complex_samples(0), verbose(false)
  {
    type        = "?";
    has_complex = false;
    reset();
  }

Member Function Documentation

collect_sample()

void collect_sample ( )

inline

Definition at line 1176 of file TraceManager.h.

Referenced by TraceManager::buffer_sample().

  {
    if (verbose)
      app_log() << " TraceBuffer<" << type << ">::collect_sample()" << std::endl;
    //make more room, if necessary
    int nrows    = buffer.size(0);
    int row_size = buffer.size(1);
    if (row_size > 0)
    {
      //make space for the row, if necessary
      int current_row = nrows;
      nrows++;
      // resizing buffer(type Array) doesn't preserve data. Thus keep old data and copy over
      auto buffer_old(buffer);
      buffer.resize(nrows, row_size);
      std::copy_n(buffer_old.data(), buffer_old.size(), buffer.data());
      if (verbose)
        app_log() << "  increasing # of rows to " << nrows << std::endl;
      //combine samples
      samples->combine_samples();
      if (has_complex)
        complex_samples->combine_samples();
      //collect data from all samples into the buffer row
      {
        std::vector<TraceSample<T>*>& ordered_samples = samples->ordered_samples;
        for (int s = 0; s < ordered_samples.size(); s++)
        {
          TraceSample<T>& tsample = *ordered_samples[s];
          if (tsample.write)
          {
            auto& sample = tsample.sample;
            for (int i = 0; i < sample.size(); ++i)
              buffer(current_row, tsample.buffer_start + i) = sample[i];
          }
        }
      }
      if (has_complex)
      {
        std::vector<TraceSample<std::complex<T>>*>& ordered_samples = complex_samples->ordered_samples;
        for (int s = 0; s < ordered_samples.size(); s++)
        {
          TraceSample<std::complex<T>>& tsample = *ordered_samples[s];
          if (tsample.write)
          {
            auto& sample = tsample.sample;
            for (int i = 0, ib = 0; i < sample.size(); ++i, ib += 2)
            {
              buffer(current_row, tsample.buffer_start + ib)     = sample[i].real();
              buffer(current_row, tsample.buffer_start + ib + 1) = sample[i].imag();
            }
          }
        }
      }
      //reset combined samples so they can be recombined on the next step
      samples->reset_combined_samples();
      if (has_complex)
        complex_samples->reset_combined_samples();
#if defined(TRACE_CHECK)
      test_buffer_collect(current_row);
#endif
    }
  }

make_combined_trace()

void make_combined_trace ( const std::string &  name, std::vector< std::string > &  names, std::vector< TraceReal > &  weights  )

inline

Definition at line 1132 of file TraceManager.h.

Referenced by TraceManager::make_combined_trace().

  {
    bool created_real = samples->make_combined_trace(name, names, weights);
    if (has_complex)
    {
      bool created_complex = complex_samples->make_combined_trace(name, names, weights);
      if (created_real && created_complex)
        APP_ABORT("TraceBuffer<" + type +
                  ">::make_combined_trace\n  cannot create real and complex combined traces for the same quantity\n  "
                  "attempted for quantity " +
                  name);
    }
  }

order_and_resize()

void order_and_resize ( )

inline

Definition at line 1149 of file TraceManager.h.

Referenced by TraceManager::initialize_traces().

  {
    //put the sample data in size order
    samples->set_unit_size(1);
    samples->order_by_size();
    if (has_complex)
    {
      complex_samples->set_unit_size(2);
      complex_samples->order_by_size();
    }
    //assign buffer ranges to each sample
    int sample_size = 0;
    samples->set_buffer_ranges(sample_size);
    if (has_complex)
      complex_samples->set_buffer_ranges(sample_size);
#if defined(TRACE_CHECK)
    test_buffer_write(sample_size);
#endif
    //resize the buffer
    int nsamples_init = 1;
    buffer.resize(nsamples_init, sample_size);
  }

register_hdf_data()

void register_hdf_data ( hdf_archive &  f )

inline

Definition at line 1264 of file TraceManager.h.

Referenced by TraceManager::open_hdf_file().

  {
    f.push(top);
    f.push("layout");
    samples->register_hdf_data(f);
    if (has_complex)
      complex_samples->register_hdf_data(f);
    f.pop();
    f.pop();
    if (!f.open_groups())
      APP_ABORT("TraceBuffer<" + type +
                ">::register_hdf_data() some hdf groups are still open at the end of registration");
    hdf_file_pointer = 0;
  }

reset()

void reset ( )

inline

Definition at line 1119 of file TraceManager.h.

Referenced by TraceManager::reset_buffers(), and TraceBuffer< TraceInt >::TraceBuffer().

{ buffer.resize(0, buffer.size(1)); }

same_as()

bool same_as ( TraceBuffer< T > &  ref )

inline

Definition at line 1173 of file TraceManager.h.

Referenced by TraceManager::check_clones().

{ return buffer.size(1) == ref.buffer.size(1); }

set_samples() [1/2]

void set_samples ( TraceSamples< T > &  s )

inline

Definition at line 1122 of file TraceManager.h.

Referenced by TraceManager::TraceManager().

{ samples = &s; }

set_samples() [2/2]

void set_samples ( TraceSamples< std::complex< T >> &  s )

inline

Definition at line 1125 of file TraceManager.h.

  {
    complex_samples = &s;
    has_complex     = true;
  }

set_type()

void set_type ( std::string  stype )

inline

Definition at line 1112 of file TraceManager.h.

Referenced by TraceManager::TraceManager().

  {
    type = stype;
    top  = type + "_data";
  }

set_verbose()

void set_verbose ( bool  v )

inline

Definition at line 1109 of file TraceManager.h.

Referenced by TraceManager::distribute().

{ verbose = v; }

test_buffer_collect()

void test_buffer_collect ( int  current_row )

inline

Definition at line 1376 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample().

  {
    if (verbose)
      app_log() << "TraceBuffer::test_buffer_collect" << std::endl;
    std::string scalars = "scalars";
    std::map<std::string, std::map<std::string, int>>::iterator dom;
    std::map<std::string, int>::iterator var;
    std::map<std::string, std::map<std::string, int>>& domains = samples->sample_indices;
    std::vector<TraceSample<T>*>& tsamples                     = samples->samples;
    std::map<std::string, int>& scalar_vars                    = domains[scalars];
    for (var = scalar_vars.begin(); var != scalar_vars.end(); var++)
    {
      const std::string& name = var->first;
      TraceSample<T>& sample  = *tsamples[var->second];
      T value                 = buffer(current_row, sample.buffer_start);
      T svalue                = 0;
      bool any_present        = false;
      for (dom = domains.begin(); dom != domains.end(); dom++)
      {
        const std::string& domain = dom->first;
        if (domain != scalars)
        {
          std::map<std::string, int>& vars = dom->second;
          if (vars.count(name) > 0)
          {
            any_present             = true;
            TraceSample<T>& ssample = *tsamples[vars[name]];
            int start               = ssample.buffer_start;
            int end                 = ssample.buffer_end;
            for (int i = start; i < end; i++)
              svalue += buffer(current_row, i);
          }
        }
      }
      if (any_present)
      {
        if (verbose)
          app_log() << "  " << name << " " << value << " " << svalue << std::endl;
      }
    }
  }

test_buffer_write()

void test_buffer_write ( int  sample_size )

inline

Definition at line 1300 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::order_and_resize().

  {
    //check that the size is correct
    int ssize = samples->total_size();
    if (has_complex)
      ssize += complex_samples->total_size();
    if (sample_size != ssize)
    {
      app_log() << "sample_size = " << sample_size << "\ntotal_size = " << ssize << std::endl;
      APP_ABORT("TraceBuffer::test_buffer_write sample_size and total_size do not match");
    }
    //check that buffer indices fall in the expected range
    int nsamples  = samples->size();
    int min_index = samples->min_buffer_index();
    int max_index = samples->max_buffer_index();
    if (has_complex)
    {
      nsamples += complex_samples->size();
      min_index = std::min(min_index, complex_samples->min_buffer_index());
      max_index = std::max(max_index, complex_samples->max_buffer_index());
    }
    if (nsamples > 0)
    {
      if (min_index != 0)
        APP_ABORT("TraceBuffer::test_buffer_write min_index!=0\n  min_index=" << min_index);
      if (max_index != sample_size)
        APP_ABORT("TraceBuffer::test_buffer_write max_index!=sample_size");
      //check that no overlap exists in writes to buffer
      Array<int, 2> test_buffer;
      test_buffer.resize(1, sample_size);
      std::fill(test_buffer.begin(), test_buffer.end(), 0);
      int row      = 0;
      int row_size = test_buffer.size(1);
      int offset   = row * row_size;
      int* loc1    = &test_buffer(offset);
      int* loc2    = &test_buffer(row, 0);
      if (loc1 != loc2)
        APP_ABORT("TraceBuffer::test_buffer_write serialized buffer offset improperly computed");
      {
        int boffset;
        std::vector<TraceSample<T>*>& ordered_samples = samples->ordered_samples;
        for (int s = 0; s < ordered_samples.size(); s++)
        {
          TraceSample<T>& tsample = *ordered_samples[s];
          std::vector<T>& sample  = tsample.sample;
          boffset                 = offset + tsample.buffer_start;
          for (int i = 0; i < sample.size(); ++i)
            test_buffer(boffset + i) = 1;
        }
      }
      if (has_complex)
      {
        int boffset;
        std::vector<TraceSample<std::complex<T>>*>& ordered_samples = complex_samples->ordered_samples;
        for (int s = 0; s < ordered_samples.size(); s++)
        {
          TraceSample<std::complex<T>>& tsample = *ordered_samples[s];
          std::vector<std::complex<T>>& sample  = tsample.sample;
          boffset                               = offset + tsample.buffer_start;
          for (int i = 0, ib = 0; i < sample.size(); ++i, ib += tsample.unit_size)
          {
            test_buffer(boffset + ib)     = 1;
            test_buffer(boffset + ib + 1) = 1;
          }
        }
      }
      //app_log()<<"test_buffer:"<< std::endl;
      //for(int i=0;i<row_size;++i)
      //  app_log()<<"  "<<i<<"  "<<test_buffer(row,i)<< std::endl;
      for (int i = 0; i < row_size; ++i)
        if (!test_buffer(row, i))
          APP_ABORT("TraceBuffer::test_buffer_write write to row is not contiguous");
    }
  }

user_report()

void user_report ( const std::string &  pad = "  " )

inline

Definition at line 1257 of file TraceManager.h.

Referenced by TraceManager::user_report().

  {
    samples->user_report(type, pad);
    if (has_complex)
      complex_samples->user_report("complex " + type, pad);
  }

write()

void write ( )

inline

Definition at line 1240 of file TraceManager.h.

{ APP_ABORT("TraceBuffer::write has not yet been implemented"); }

write_hdf() [1/2]

void write_hdf ( hdf_archive &  f )

inline

Definition at line 1280 of file TraceManager.h.

Referenced by TraceManager::write_buffers_hdf(), and TraceBuffer< TraceInt >::write_hdf().

{ write_hdf(f, hdf_file_pointer); }

write_hdf() [2/2]

void write_hdf ( hdf_archive &  f, hsize_t &  file_pointer  )

inline

Definition at line 1283 of file TraceManager.h.

  {
    if (verbose)
      app_log() << "TraceBuffer<" << type << ">::write_hdf() " << file_pointer << " " << buffer.size(0) << " "
                << buffer.size(1) << std::endl;
    dims[0] = buffer.size(0);
    dims[1] = buffer.size(1);
    if (dims[0] > 0)
    {
      f.push(top);
      h5d_append(f.top(), "traces", file_pointer, buffer.dim(), dims, buffer.data());
      f.pop();
    }
    f.flush();
  }

write_summary()

void write_summary ( std::string  pad = "  " )

inline

Definition at line 1243 of file TraceManager.h.

Referenced by TraceManager::write_summary().

  {
    std::string pad2 = pad + "  ";
    app_log() << pad << "TraceBuffer<" << type << ">" << std::endl;
    app_log() << pad2 << "nrows       = " << buffer.size(0) << std::endl;
    app_log() << pad2 << "row_size    = " << buffer.size(1) << std::endl;
    app_log() << pad2 << "has_complex = " << has_complex << std::endl;
    samples->write_summary(type, pad2);
    if (has_complex)
      complex_samples->write_summary("complex " + type, pad2);
    app_log() << pad << "end TraceBuffer<" << type << ">" << std::endl;
  }

Member Data Documentation

buffer

Array<T, 2> buffer

Definition at line 1092 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample(), TraceBuffer< TraceInt >::order_and_resize(), TraceBuffer< TraceInt >::reset(), TraceBuffer< TraceInt >::same_as(), TraceBuffer< TraceInt >::test_buffer_collect(), TraceBuffer< TraceInt >::write_hdf(), and TraceBuffer< TraceInt >::write_summary().

complex_samples

TraceSamples<std::complex<T> >* complex_samples

Definition at line 1090 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample(), TraceBuffer< TraceInt >::make_combined_trace(), TraceBuffer< TraceInt >::order_and_resize(), TraceBuffer< TraceInt >::register_hdf_data(), TraceBuffer< TraceInt >::set_samples(), TraceBuffer< TraceInt >::test_buffer_write(), TraceBuffer< TraceInt >::user_report(), and TraceBuffer< TraceInt >::write_summary().

dims

hsize_t dims[2]

Definition at line 1097 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::write_hdf().

has_complex

bool has_complex

Definition at line 1088 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample(), TraceBuffer< TraceInt >::make_combined_trace(), TraceBuffer< TraceInt >::order_and_resize(), TraceBuffer< TraceInt >::register_hdf_data(), TraceBuffer< TraceInt >::set_samples(), TraceBuffer< TraceInt >::test_buffer_write(), TraceBuffer< TraceInt >::TraceBuffer(), TraceBuffer< TraceInt >::user_report(), and TraceBuffer< TraceInt >::write_summary().

hdf_file_pointer

hsize_t hdf_file_pointer

Definition at line 1098 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::register_hdf_data(), TraceManager::write_buffers_hdf(), and TraceBuffer< TraceInt >::write_hdf().

samples

TraceSamples<T>* samples

Definition at line 1089 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample(), TraceBuffer< TraceInt >::make_combined_trace(), TraceBuffer< TraceInt >::order_and_resize(), TraceBuffer< TraceInt >::register_hdf_data(), TraceBuffer< TraceInt >::set_samples(), TraceBuffer< TraceInt >::test_buffer_collect(), TraceBuffer< TraceInt >::test_buffer_write(), TraceBuffer< TraceInt >::user_report(), and TraceBuffer< TraceInt >::write_summary().

top

std::string top

Definition at line 1096 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::register_hdf_data(), TraceBuffer< TraceInt >::set_type(), and TraceBuffer< TraceInt >::write_hdf().

type

std::string type

Definition at line 1091 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample(), TraceBuffer< TraceInt >::make_combined_trace(), TraceBuffer< TraceInt >::register_hdf_data(), TraceBuffer< TraceInt >::set_type(), TraceBuffer< TraceInt >::TraceBuffer(), TraceBuffer< TraceInt >::user_report(), TraceBuffer< TraceInt >::write_hdf(), and TraceBuffer< TraceInt >::write_summary().

verbose

bool verbose

Definition at line 1093 of file TraceManager.h.

Referenced by TraceBuffer< TraceInt >::collect_sample(), TraceBuffer< TraceInt >::set_verbose(), TraceBuffer< TraceInt >::test_buffer_collect(), and TraceBuffer< TraceInt >::write_hdf().

---

The documentation for this struct was generated from the following file:

• /home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/Estimators/TraceManager.h