Interaction.h

Go to the documentation of this file.

// Copyright András Vukics 2006–2014. Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE.txt)
// -*- C++ -*-
#ifndef CPPQEDCORE_STRUCTURE_INTERACTION_H_INCLUDED
#define CPPQEDCORE_STRUCTURE_INTERACTION_H_INCLUDED

#include "InteractionFwd.h"

#include "DynamicsBase.h"
#include "Free.h"

#include "BlitzTiny.h"
#include "SmartPtr.h"

#include <boost/bind.hpp>

#include <boost/range/adaptor/transformed.hpp>
#include <boost/range/algorithm/copy.hpp>


namespace structure {



template<int RANK>
class Interaction : public DynamicsBase, public DimensionsBookkeeper<RANK>
{
public:
  typedef boost::shared_ptr<const Interaction> Ptr;


  typedef blitz::TinyVector<Free::Ptr,RANK> Frees;

  template<typename F> using FreesTemplate=blitz::TinyVector<F,RANK>; 
  
  typedef typename DimensionsBookkeeper<RANK>::Dimensions Dimensions;
  
  explicit Interaction(const Frees& frees,
                       const    RealFreqs&    realFreqs=emptyRF, 
                       const ComplexFreqs& complexFreqs=emptyCF)
    : DynamicsBase(realFreqs,complexFreqs), DimensionsBookkeeper<RANK>(extractDimensions(frees)), frees_(frees) {}

  Interaction(const Frees& frees, const ComplexFreqs& complexFreqs) : Interaction(frees,emptyRF,complexFreqs) {}
  Interaction(const Frees& frees, RealFreqsInitializer rf, ComplexFreqsInitializer cf={}) : Interaction(frees,RealFreqs(rf),ComplexFreqs(cf)) {}
  Interaction(const Frees& frees, ComplexFreqsInitializer cf) : Interaction(frees,RealFreqsInitializer(),cf) {}
  Interaction(const Frees& frees, RF rf, CF cf=CF()) : Interaction(frees,RealFreqsInitializer{rf}, cf==CF() ? ComplexFreqsInitializer{} : ComplexFreqsInitializer{cf}) {}
  Interaction(const Frees& frees, CF cf) : Interaction(frees,ComplexFreqsInitializer{cf}) {}
  Interaction(const Frees& frees, RealFreqsInitializer rf, CF cf) : Interaction(frees,rf,{cf}) {}
  Interaction(const Frees& frees, RF rf, ComplexFreqsInitializer cf) : Interaction(frees,{rf},cf) {}

  template<typename F0, typename F1, typename... RCF_type>
  Interaction(const F0& f0, const F1& f1,  const RCF_type&... realAndComplexFreqs)
    :  Interaction(Frees(cpputils::sharedPointerize(f0),cpputils::sharedPointerize(f1)),realAndComplexFreqs...) {}

  const Frees& getFrees() const {return frees_;}

protected:
  class FreesProxy
  {
  public:
    template<typename... Pack>
    FreesProxy(Pack&&... pack) : frees_(cpputils::sharedPointerize(pack)...) {}
    
    operator const Frees&() const {return frees_;}
    
  private:
    const Frees frees_;
  };

private:
  static const Dimensions extractDimensions(const Frees& frees)
  {
    Dimensions res;
    boost::copy(frees | boost::adaptors::transformed(boost::bind(&Free::getTotalDimension,_1)), res.begin());
    return res;
  }
  
  const Frees frees_;

};


} // structure

#endif // CPPQEDCORE_STRUCTURE_INTERACTION_H_INCLUDED