/*
 * OrthogonalSumUpPerLevel.hpp
 *
 *  Created on: Aug 27, 2012
 *      Author: heber
 */

#ifndef ORTHOGONALSUMUPPERLEVEL_HPP_
#define ORTHOGONALSUMUPPERLEVEL_HPP_


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


#include <map>
#include <vector>

#include "Fragmentation/Summation/Converter/DataConverter.hpp"
#include "Fragmentation/Summation/AllLevelOrthogonalSummator.hpp"
#include "Fragmentation/Summation/IndexSet.hpp"
#include "Fragmentation/Summation/IndexSetContainer.hpp"

#include "Fragmentation/Summation/Containers/MPQCData.hpp"

#include <boost/mpl/for_each.hpp>

template <typename TypeMap, typename TypeData, typename TypeVector>
void OrthogonalSumUpPerLevel(
    const std::map<JobId_t, TypeData> &fragmentData,
    const std::map< JobId_t, size_t > &MatrixNrLookup,
    const IndexSetContainer::ptr &container,
    SubsetMap::ptr &subsetmap,
    std::vector<TypeMap> &levelresults,
    std::map<IndexSet::ptr, std::pair<TypeMap, TypeMap> > &keysetresults
    )
{
  // place data into boost::fusion::map instance
  std::map<JobId_t, TypeMap> Data_fused;
  convertDataTo<TypeData, TypeMap>(fragmentData, Data_fused);
  // instantiate summator
  levelresults.resize(subsetmap->getMaximumSetLevel());
  AllLevelOrthogonalSummator<TypeMap> Summer(
              subsetmap,
              Data_fused,
              container->getContainer(),
              MatrixNrLookup,
              levelresults,
              keysetresults);
  // sum up for each type key in TypeVector
  boost::mpl::for_each<TypeVector>(boost::ref(Summer));
}


#endif /* ORTHOGONALSUMUPPERLEVEL_HPP_ */