| [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 |  | 
|---|
|  | 21 | #include "Fragmentation/Converter/DataConverter.hpp" | 
|---|
|  | 22 | #include "Fragmentation/Summation/AllLevelOrthogonalSummator.hpp" | 
|---|
| [b8f0b25] | 23 | #include "Fragmentation/Summation/IndexSet.hpp" | 
|---|
| [0398bd] | 24 | #include "Fragmentation/Summation/IndexSetContainer.hpp" | 
|---|
|  | 25 |  | 
|---|
|  | 26 | #include "Jobs/MPQCData.hpp" | 
|---|
|  | 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_ */ | 
|---|