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

/*
 * OrthogonalSummation_impl.hpp
 *
 *  Created on: Jun 25, 2012
 *      Author: heber
 */

#ifndef ORTHOGONALSUMMATION_IMPL_HPP_
#define ORTHOGONALSUMMATION_IMPL_HPP_


// include config.h
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include "CodePatterns/IteratorAdaptors.hpp"
#include "CodePatterns/Log.hpp"

#include "Fragmentation/Summation/OrthogonalSummation.hpp"

#include "Fragmentation/Summation/ZeroInstance.hpp"
#include "Fragmentation/Summation/printKeyNames.hpp"

/** Constructor of class OrthogonalSummation.
 *
 */
template <class T>
OrthogonalSummation<T>::OrthogonalSummation(
    InputSets_t &indices,
    InputValues_t& values,
    SubsetMap::ptr _subsetmap) :
  subsetmap(_subsetmap),
  zeroinstance(ZeroInstance<T>())
{
  ASSERT( indices.size() == values.size(),
      "OrthogonalSummation<T>::OrthogonalSummation() - indices and values mismatch in size: "
      +toString(indices.size())+" != "+toString(values.size())+".");
  /// place each index
  /// create own map if none is given
  if (!subsetmap) {
    typename InputSets_t::iterator iter = indices.begin();
    IndexSetContainer container(*iter);
    for (; iter != indices.end(); ++iter)
      container.insert(*iter);
    subsetmap.reset(new SubsetMap(container));
  } else {
    LOG(2, "DEBUG: Using given SubsetMap.");
  }
  /// instantiate all SubSetValue's by requesting the IndexSet from the Subsetmap
  typename InputSets_t::iterator indexiter = indices.begin();
  typename InputValues_t::iterator valueiter = values.begin();
  for (;valueiter != values.end(); ++indexiter, ++valueiter) {
    LOG(2, "DEBUG: Adding set " << **indexiter << " with value " << *valueiter << ".");
    setvalues.addValue( *indexiter, *valueiter );
  }
  /// bind static lookup functions for SetValue<T>
  SetValue<T>::lookupSubset =
      boost::bind(&SubsetMap::getSubsets, boost::ref(*subsetmap), _1);
  SetValue<T>::lookupValue =
      boost::bind(&SetValueMap<T>::getValue, boost::ref(setvalues), _1);
}

template <class T>
T OrthogonalSummation<T>::operator()(const size_t level) const
{
  return Sum(level);
}

template <class T>
T OrthogonalSummation<T>::Sum(const size_t level) const
{
  T sum(zeroinstance);
  for(typename SetValueMap<T>::const_iterator iter = setvalues.begin();
      (iter != setvalues.end()) && (iter->first->size() <= level); ++iter) {
    const T tempvalue = (iter->second)->getContribution();
    sum += tempvalue;
    LOG(3, "DEBUG: Contribution from subset "+toString(*(iter->second->getIndexSet()))
        +" is "+toString(tempvalue)+".");
  }
  return sum;
}


#endif /* ORTHOGONALSUMMATION_IMPL_HPP_ */