source: src/Fragmentation/Summation/Summation_impl.hpp@ ff9963

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 ff9963 was ff9963, checked in by Frederik Heber <heber@…>, 12 years ago

(Orthogonal)Summation now goes level-wise.

  • print(Orthogonal)Sum go through all levels and call on (Orthogonal)Summator to get value which is stored in vector of result of MapType.
  • (Orthogonal)Summator hands on level to (Orthogonal)Summation
  • (Orthogonal)Summation::operator() and Sum() now require level as parameter.
  • NOTE: We cannot set the value to zero (this depends on the default cstor). If level==0, then we return an uninitialized value.
  • Renamed print(Orthogonal)Sum -> AllLevel(Orthogonal)Summator.
  • Property mode set to 100644
File size: 2.5 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/** Constructor of class Summation.
21 *
22 */
23template <class T>
24Summation<T>::Summation(
25 InputSets_t &indices,
26 InputValues_t& values,
27 SubsetMap::ptr _subsetmap) :
28 subsetmap(_subsetmap)
29{
30 ASSERT( indices.size() == values.size(),
31 "Summation<T>::Summation() - indices and values mismatch in size: "
32 +toString(indices.size())+" != "+toString(values.size())+".");
33 /// place each index
34 /// create own map if none is given
35 if (!subsetmap) {
36 typename InputSets_t::iterator iter = indices.begin();
37 IndexSetContainer container(*iter);
38 for (; iter != indices.end(); ++iter)
39 container.insert(*iter);
40 subsetmap.reset(new SubsetMap(container));
41 } else {
42 LOG(1, "INFO: Using given SubsetMap.");
43 }
44 /// instantiate all SubSetValue's by requesting the IndexSet from the Subsetmap
45 typename InputSets_t::iterator indexiter = indices.begin();
46 typename InputValues_t::iterator valueiter = values.begin();
47 for (;valueiter != values.end(); ++indexiter, ++valueiter) {
48 LOG(1, "INFO: Adding set " << **indexiter << " with value " << *valueiter << ".");
49 setvalues.addValue( *indexiter, *valueiter );
50 }
51 /// bind static lookup functions for SetValue<T>
52 SetValue<T>::lookupSubset =
53 boost::bind(&SubsetMap::getSubsets, boost::ref(*subsetmap), _1);
54 SetValue<T>::lookupValue =
55 boost::bind(&SetValueMap<T>::getValue, boost::ref(setvalues), _1);
56}
57
58template <class T>
59T Summation<T>::operator()(const size_t level) const
60{
61 return Sum(level);
62}
63
64template <class T>
65T Summation<T>::Sum(const size_t level) const
66{
67 typename SetValueMap<T>::const_iterator iter = setvalues.begin();
68 T sum; // we must assume here that cstor initializes to 0
69 if (iter->first->size() <= level)
70 sum = (iter->second)->getValue();
71 LOG(1, "DEBUG: Value from subset "+toString(*(iter->second->getIndexSet()))
72 +" is "+toString(sum)+".");
73 for(++iter;(iter != setvalues.end()) && (iter->first->size() <= level); ++iter) {
74 const T tempvalue = (iter->second)->getValue();
75 sum += tempvalue;
76 LOG(1, "DEBUG: Value from subset "+toString(*(iter->second->getIndexSet()))
77 +" is "+toString(tempvalue)+".");
78 }
79 return sum;
80}
81
82
83#endif /* SUMMATION_IMPL_HPP_ */
Note: See TracBrowser for help on using the repository browser.