source: src/Fragmentation/Summation/OrthogonalSummation_impl.hpp@ cb3363

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 Candidate_v1.7.0 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 cb3363 was 1e7dd4, checked in by Frederik Heber <heber@…>, 13 years ago

Added OrthogonalSummation and unit test.
  • Property mode set to 100644
File size: 2.4 KB
Line 
1/*
2 * OrthogonalSummation_impl.hpp
3 *
4 * Created on: Jun 25, 2012
5 * Author: heber
6 */
7
8#ifndef ORTHOGONALSUMMATION_IMPL_HPP_
9#define ORTHOGONALSUMMATION_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 OrthogonalSummation.
21 *
22 */
23template <class T>
24OrthogonalSummation<T>::OrthogonalSummation(
25 InputSets_t &indices,
26 InputValues_t& values,
27 SubsetMap::ptr _subsetmap) :
28 subsetmap(_subsetmap)
29{
30 ASSERT( indices.size() == values.size(),
31 "OrthogonalSummation<T>::OrthogonalSummation() - 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 OrthogonalSummation<T>::operator()() const
60{
61 return Sum();
62}
63
64template <class T>
65T OrthogonalSummation<T>::Sum() const
66{
67 typename SetValueMap<T>::const_iterator iter = setvalues.begin();
68 T sum = (iter->second)->getContribution();
69 LOG(1, "DEBUG: Contribution from subset "+toString(*(iter->second->getIndexSet()))
70 +" is "+toString(sum)+".");
71 for(++iter;iter != setvalues.end(); ++iter) {
72 const T tempvalue = (iter->second)->getContribution();
73 sum += tempvalue;
74 LOG(1, "DEBUG: Contribution from subset "+toString(*(iter->second->getIndexSet()))
75 +" is "+toString(tempvalue)+".");
76 }
77 return sum;
78}
79
80
81#endif /* ORTHOGONALSUMMATION_IMPL_HPP_ */
Note: See TracBrowser for help on using the repository browser.