source: src/Fragmentation/Summation/Summation_impl.hpp@ 1b5a40

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
Last change on this file since 1b5a40 was 1b5a40, checked in by Frederik Heber <heber@…>, 12 years ago

VERBOSE: Changed verbosity of Summation process to make a lot less noise on default setting.
  • Property mode set to 100644
File size: 2.6 KB
Line 
1/*
2 * Summation_impl.hpp
3 *
4 * Created on: Jun 30, 2012
5 * Author: heber
6 */
7
8#ifndef SUMMATION_IMPL_HPP_
9#define SUMMATION_IMPL_HPP_
10
11
12// include config.h
13#ifdef HAVE_CONFIG_H
14#include <config.h>
15#endif
16
17#include "CodePatterns/IteratorAdaptors.hpp"
18#include "CodePatterns/Log.hpp"
19
20#include "Fragmentation/Summation/printKeyNames.hpp"
21
22/** Constructor of class Summation.
23 *
24 */
25template <class T>
26Summation<T>::Summation(
27 InputSets_t &indices,
28 InputValues_t& values,
29 SubsetMap::ptr _subsetmap) :
30 subsetmap(_subsetmap)
31{
32 ASSERT( indices.size() == values.size(),
33 "Summation<T>::Summation() - indices and values mismatch in size: "
34 +toString(indices.size())+" != "+toString(values.size())+".");
35 /// place each index
36 /// create own map if none is given
37 if (!subsetmap) {
38 typename InputSets_t::iterator iter = indices.begin();
39 IndexSetContainer container(*iter);
40 for (; iter != indices.end(); ++iter)
41 container.insert(*iter);
42 subsetmap.reset(new SubsetMap(container));
43 } else {
44 LOG(2, "DEBUG: Using given SubsetMap for " << printKeyNames::printName<T>() << ".");
45 }
46 /// instantiate all SubSetValue's by requesting the IndexSet from the Subsetmap
47 typename InputSets_t::iterator indexiter = indices.begin();
48 typename InputValues_t::iterator valueiter = values.begin();
49 for (;valueiter != values.end(); ++indexiter, ++valueiter) {
50 LOG(2, "DEBUG: Adding set " << **indexiter << " with value " << *valueiter << ".");
51 setvalues.addValue( *indexiter, *valueiter );
52 }
53 /// bind static lookup functions for SetValue<T>
54 SetValue<T>::lookupSubset =
55 boost::bind(&SubsetMap::getSubsets, boost::ref(*subsetmap), _1);
56 SetValue<T>::lookupValue =
57 boost::bind(&SetValueMap<T>::getValue, boost::ref(setvalues), _1);
58}
59
60template <class T>
61T Summation<T>::operator()(const size_t level) const
62{
63 return Sum(level);
64}
65
66template <class T>
67T Summation<T>::Sum(const size_t level) const
68{
69 typename SetValueMap<T>::const_iterator iter = setvalues.begin();
70 T sum; // we must assume here that cstor initializes to 0
71 if (iter->first->size() <= level)
72 sum = (iter->second)->getValue();
73 LOG(3, "DEBUG: Value from subset "+toString(*(iter->second->getIndexSet()))
74 +" is "+toString(sum)+".");
75 for(++iter;(iter != setvalues.end()) && (iter->first->size() <= level); ++iter) {
76 const T tempvalue = (iter->second)->getValue();
77 sum += tempvalue;
78 LOG(3, "DEBUG: Value from subset "+toString(*(iter->second->getIndexSet()))
79 +" is "+toString(tempvalue)+".");
80 }
81 return sum;
82}
83
84
85#endif /* SUMMATION_IMPL_HPP_ */
Note: See TracBrowser for help on using the repository browser.