source: src/Fragmentation/Summation/AllLevelSummator.hpp@ c8f12e

Action_Thermostats Add_AtomRandomPerturbation Add_FitFragmentPartialChargesAction Add_RotateAroundBondAction Add_SelectAtomByNameAction Added_ParseSaveFragmentResults AddingActions_SaveParseParticleParameters Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_ParticleName_to_Atom Adding_StructOpt_integration_tests AtomFragments Automaking_mpqc_open AutomationFragmentation_failures Candidate_v1.5.4 Candidate_v1.6.0 Candidate_v1.6.1 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator CombiningParticlePotentialParsing Combining_Subpackages Debian_Package_split Debian_package_split_molecuildergui_only Disabling_MemDebug Docu_Python_wait EmpiricalPotential_contain_HomologyGraph EmpiricalPotential_contain_HomologyGraph_documentation Enable_parallel_make_install Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph FitPartialCharges_GlobalError Fix_BoundInBox_CenterInBox_MoleculeActions Fix_ChargeSampling_PBC Fix_ChronosMutex Fix_FitPartialCharges Fix_FitPotential_needs_atomicnumbers Fix_ForceAnnealing Fix_IndependentFragmentGrids Fix_ParseParticles Fix_ParseParticles_split_forward_backward_Actions Fix_PopActions Fix_QtFragmentList_sorted_selection Fix_Restrictedkeyset_FragmentMolecule Fix_StatusMsg Fix_StepWorldTime_single_argument Fix_Verbose_Codepatterns Fix_fitting_potentials Fixes ForceAnnealing_goodresults ForceAnnealing_oldresults ForceAnnealing_tocheck ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion FragmentAction_writes_AtomFragments FragmentMolecule_checks_bonddegrees GeometryObjects Gui_Fixes Gui_displays_atomic_force_velocity ImplicitCharges IndependentFragmentGrids IndependentFragmentGrids_IndividualZeroInstances IndependentFragmentGrids_IntegrationTest IndependentFragmentGrids_Sole_NN_Calculation JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool JobMarket_unresolvable_hostname_fix MoreRobust_FragmentAutomation ODR_violation_mpqc_open PartialCharges_OrthogonalSummation PdbParser_setsAtomName PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks Rewrite_FitPartialCharges RotateToPrincipalAxisSystem_UndoRedo SaturateAtoms_findBestMatching SaturateAtoms_singleDegree StoppableMakroAction Subpackage_CodePatterns Subpackage_JobMarket Subpackage_LinearAlgebra Subpackage_levmar Subpackage_mpqc_open Subpackage_vmg Switchable_LogView ThirdParty_MPQC_rebuilt_buildsystem TrajectoryDependenant_MaxOrder TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps TremoloParser_setsAtomName Ubuntu_1604_changes stable v1.3.7
Last change on this file since c8f12e was 19c50e, checked in by Frederik Heber <heber@…>, 12 years ago

IndexSetContainer now treats super set specially.

  • The super set must not gather its subsets via the gatherSubsets() as by construction all other sets are its subsets! As the super set is very large the power set way is no good idea.
  • added default cstor for SortedVector
  • removed SubsetMap::getMaximumSubsetLevel() as is replaced by ::getMaximumSetLevel() which is the level to sum up to.
  • changed all uses of getMaximumSubsetLevel() to getMaximumSetLevel().
  • TESTFIX: Changed unit test function on getMaximumSubsetLevel() to check on getMaximumSetLevel()
  • removed OrthogonalFullSummator as is fully replacable by OrthogonalSummator.
  • changed IndexSetContainer::createSuperSet a bit.
  • IndexSetContainer::AllIndices is now no more static convenience entity but truely contains the super set (non-statically). ::createSuperSet() is for convenience to be called in cstor for AllIndices.
  • Property mode set to 100644
File size: 3.0 KB
Line 
1/*
2 * AllLevelSummator.hpp
3 *
4 * Created on: 29.07.2012
5 * Author: heber
6 */
7
8#ifndef ALLLEVELSUMMATOR_HPP_
9#define ALLLEVELSUMMATOR_HPP_
10
11// include config.h
12#ifdef HAVE_CONFIG_H
13#include <config.h>
14#endif
15
16#include <vector>
17
18#include "CodePatterns/Assert.hpp"
19
20#include "Fragmentation/Summation/IndexSetContainer.hpp"
21#include "Fragmentation/Summation/SubsetMap.hpp"
22#include "Fragmentation/Summation/Summator.hpp"
23
24#include "Fragmentation/Summation/printKeyNames.hpp"
25
26/** Tiny template functor to use Summation, sum up each level, and store the
27 * result in a given vector.
28 *
29 */
30template <typename MapType>
31struct AllLevelSummator {
32 /** Constructor takes the arguments that \a Summator also needs and stores
33 * them internally.
34 *
35 * \param _subsetmap map with hierarchy of IndexSet's
36 * \param _data MPQCData converted to MPQCDataMap_t type, associated to job id
37 * \param _container container of IndexSet's such that each set has correct order
38 * to job id and hence to _data.
39 * \param _MatrixNrLookup lookup from job id to ordering in above vectors
40 * \param _results vector place results into
41 */
42 AllLevelSummator(
43 SubsetMap::ptr &_subsetmap,
44 const std::map<JobId_t, MapType> &_data,
45 const IndexSetContainer::Container_t &_container,
46 const std::map< JobId_t, size_t > &_MatrixNrLookup,
47 std::vector<MapType> &_results) :
48 subsetmap(_subsetmap),
49 data(_data),
50 container(_container),
51 MatrixNrLookup(_MatrixNrLookup),
52 results(_results)
53 {
54 ASSERT( results.size() >= subsetmap->getMaximumSetLevel(),
55 "AllLevelSummator() - result vector is not large enough.");
56 }
57
58 /** Operator that calls on Summator and prints the value.
59 *
60 * \note the parameter is needed for boost::mpl::for_each but is not
61 * used here.
62 */
63 template <typename MapKey>
64 void operator()(MapKey &) {
65 // We retrieve the type of the MPQCData member variable from the boost::fusion::map.
66 typedef typename boost::fusion::result_of::value_at_key<MapType, MapKey>::type MapValue;
67 // create Summator instance
68 Summator<MapType, MapKey> sum_value(
69 subsetmap, data, container, MatrixNrLookup
70 );
71 const size_t MaxLevel = subsetmap->getMaximumSetLevel();
72 for (size_t level=1; level <= MaxLevel; ++level) {
73 MapType &LevelResults = results[level-1];
74 // sum up and store in results
75 const MapValue value = sum_value(level);
76 boost::fusion::at_key<MapKey>(LevelResults) = value;
77 // print value
78 //LOG(0, "STATUS: Level " << level << " resulting " << printKeyNames::printName<MapKey>() << " is " << value << ".");
79 }
80 }
81
82private:
83 //!> Hierarchy of IndexSet's
84 SubsetMap::ptr &subsetmap;
85 //!> vector of data converted from MPQCData
86 const std::map<JobId_t, MapType> &data;
87 //!> container with all IndexSet's
88 const IndexSetContainer::Container_t &container;
89 //!> lookup map from job ids to ordering in above vectors
90 const std::map< JobId_t, size_t > &MatrixNrLookup;
91 //!> vector of results
92 std::vector<MapType> &results;
93};
94
95#endif /* ALLLEVELSUMMATOR_HPP_ */
Note: See TracBrowser for help on using the repository browser.