| [0398bd] | 1 | /*
 | 
|---|
 | 2 |  * OrthogonalSumUpPerLevel.hpp
 | 
|---|
 | 3 |  *
 | 
|---|
 | 4 |  *  Created on: Aug 27, 2012
 | 
|---|
 | 5 |  *      Author: heber
 | 
|---|
 | 6 |  */
 | 
|---|
 | 7 | 
 | 
|---|
 | 8 | #ifndef ORTHOGONALSUMUPPERLEVEL_HPP_
 | 
|---|
 | 9 | #define ORTHOGONALSUMUPPERLEVEL_HPP_
 | 
|---|
 | 10 | 
 | 
|---|
 | 11 | 
 | 
|---|
 | 12 | // include config.h
 | 
|---|
 | 13 | #ifdef HAVE_CONFIG_H
 | 
|---|
 | 14 | #include <config.h>
 | 
|---|
 | 15 | #endif
 | 
|---|
 | 16 | 
 | 
|---|
 | 17 | 
 | 
|---|
 | 18 | #include <map>
 | 
|---|
 | 19 | #include <vector>
 | 
|---|
 | 20 | 
 | 
|---|
| [fbf143] | 21 | #include "Fragmentation/Summation/Converter/DataConverter.hpp"
 | 
|---|
| [0398bd] | 22 | #include "Fragmentation/Summation/AllLevelOrthogonalSummator.hpp"
 | 
|---|
| [b8f0b25] | 23 | #include "Fragmentation/Summation/IndexSet.hpp"
 | 
|---|
| [0398bd] | 24 | #include "Fragmentation/Summation/IndexSetContainer.hpp"
 | 
|---|
 | 25 | 
 | 
|---|
| [fbf143] | 26 | #include "Fragmentation/Summation/Containers/MPQCData.hpp"
 | 
|---|
| [0398bd] | 27 | 
 | 
|---|
 | 28 | #include <boost/mpl/for_each.hpp>
 | 
|---|
 | 29 | 
 | 
|---|
| [a3112d] | 30 | template <typename TypeMap, typename TypeData, typename TypeVector>
 | 
|---|
| [b8f0b25] | 31 | void OrthogonalSumUpPerLevel(
 | 
|---|
| [a3112d] | 32 |     const std::map<JobId_t, TypeData> &fragmentData,
 | 
|---|
| [3dd32f] | 33 |     const std::map< JobId_t, size_t > &MatrixNrLookup,
 | 
|---|
| [0398bd] | 34 |     const IndexSetContainer::ptr &container,
 | 
|---|
| [b8f0b25] | 35 |     SubsetMap::ptr &subsetmap,
 | 
|---|
 | 36 |     std::vector<TypeMap> &levelresults,
 | 
|---|
| [79398d] | 37 |     std::map<IndexSet::ptr, std::pair<TypeMap, TypeMap> > &keysetresults
 | 
|---|
| [0398bd] | 38 |     )
 | 
|---|
 | 39 | {
 | 
|---|
 | 40 |   // place data into boost::fusion::map instance
 | 
|---|
| [a3112d] | 41 |   std::map<JobId_t, TypeMap> Data_fused;
 | 
|---|
 | 42 |   convertDataTo<TypeData, TypeMap>(fragmentData, Data_fused);
 | 
|---|
| [0398bd] | 43 |   // instantiate summator
 | 
|---|
| [b8f0b25] | 44 |   levelresults.resize(subsetmap->getMaximumSetLevel());
 | 
|---|
| [0398bd] | 45 |   AllLevelOrthogonalSummator<TypeMap> Summer(
 | 
|---|
 | 46 |               subsetmap,
 | 
|---|
| [a3112d] | 47 |               Data_fused,
 | 
|---|
| [0398bd] | 48 |               container->getContainer(),
 | 
|---|
 | 49 |               MatrixNrLookup,
 | 
|---|
| [b8f0b25] | 50 |               levelresults,
 | 
|---|
 | 51 |               keysetresults);
 | 
|---|
| [0398bd] | 52 |   // sum up for each type key in TypeVector
 | 
|---|
 | 53 |   boost::mpl::for_each<TypeVector>(boost::ref(Summer));
 | 
|---|
 | 54 | }
 | 
|---|
 | 55 | 
 | 
|---|
 | 56 | 
 | 
|---|
 | 57 | #endif /* ORTHOGONALSUMUPPERLEVEL_HPP_ */
 | 
|---|