source: src/Fragmentation/Summation/printOrthogonalSum.hpp@ 8b58ac

/*
 * printOrthogonalSum.hpp
 *
 *  Created on: 29.07.2012
 *      Author: heber
 */

#ifndef PRINTORTHOGONALSUM_HPP_
#define PRINTORTHOGONALSUM_HPP_

// include config.h
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include <vector>

#include "Fragmentation/Summation/IndexSetContainer.hpp"
#include "Fragmentation/Summation/SubsetMap.hpp"
#include "Fragmentation/Summation/OrthogonalSummator.hpp"

#include "Fragmentation/Summation/printKeyNames.hpp"

/** Tiny template functor to use OrthogonalSummation, sum up and print the result.
 *
 */
template <typename MapType>
struct printOrthogonalSum {
  /** Constructor takes the arguments that \a OrthogonalSummator also needs and stores
   * them internally.
   *
   * \param _subsetmap map with hierarchy of IndexSet's
   * \param _data MPQCData converted to MPQCDataMap_t type
   * \param _jobids job ids to sum data in correct order
   * \param _container container of IndexSet's such that each set has correct order
   *        to job id and hence to _data.
   * \param _MatrixNrLookup lookup from job id to ordering in above vectors
   */
  printOrthogonalSum(
      SubsetMap::ptr &_subsetmap,
      const std::vector<MapType> &_data,
      const std::vector<JobId_t> &_jobids,
      const IndexSetContainer::Container_t &_container,
      std::map< JobId_t, size_t > &_MatrixNrLookup) : /* cannot make this const due to operator[] */
    subsetmap(_subsetmap),
    data(_data),
    jobids(_jobids),
    container(_container),
    MatrixNrLookup(_MatrixNrLookup)
  {}

  /** Operator that calls on OrthogonalSummator and prints the value.
   *
   * \note the parameter is needed for boost::mpl::for_each but is not
   * used here.
   */
  template <typename MapKey>
  void operator()(MapKey &) {
    // We retrieve the type of the MPQCData member variable from the boost::fusion::map.
    typedef typename boost::fusion::result_of::value_at_key<MapType, MapKey>::type MapValue;
    // create OrthogonalSummator instance
    OrthogonalSummator<MapType, MapKey> sum_value(
        subsetmap, data, jobids, container, MatrixNrLookup
        );
    // sum up
    MapValue value = sum_value();
    // print value
    LOG(0, "STATUS: Resulting " << printKeyNames::printName<MapKey>() << " is " << value << ".");
  }

private:
  //!> Hierarchy of IndexSet's
  SubsetMap::ptr &subsetmap;
  //!> vector of data converted from MPQCData
  const std::vector<MapType> &data;
  //!> vector of jobids
  const std::vector<JobId_t> &jobids;
  //!> container with all IndexSet's
  const IndexSetContainer::Container_t &container;
  //!> lookup map from job ids to ordering in above vectors
  std::map< JobId_t, size_t > MatrixNrLookup;
};

#endif /* PRINTORTHOGONALSUM_HPP_ */