ProductRule.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/MEAL/MEAL/ProductRule.h
 
#ifndef __MEAL_ProductRule_H
#define __MEAL_ProductRule_H
 
#include "MEAL/GroupRule.h"
 
namespace MEAL {

  template<class T>
  class ProductRule : public GroupRule<T>
  {
 
  public:
 
    typedef typename T::Result Result;

    void operator *= (T* model) { this->add_model (model); }
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Function implementation
    //
    // ///////////////////////////////////////////////////////////////////

    std::string get_name () const
    { return "ProductRule<" + std::string(T::Name)+ ">"; }
 
  protected:
    
    // ///////////////////////////////////////////////////////////////////
    //
    // GroupRule implementation
    //
    // ///////////////////////////////////////////////////////////////////

    const Result get_identity() const { return 1.0; }

    void operate (Result& total, const Result& element)
    { total = total * element; }

    const Result partial (const Result& element) const { return element; }
 
  };
  
  template<class T>
  ProductRule<T>* product (T* A, T* B)
  {
    Reference::To< ProductRule<T> > result = new ProductRule<T>;
    result->add_model (A);
    result->add_model (B);
    return result.release();
  }
 
  template<class T>
  ProductRule<T>* product (Reference::To<T>& A, T* B)
  { return product (A.get(), B); }
 
  template<class T>
  ProductRule<T>* product (T* A, Reference::To<T>& B)
  { return product (A, B.get()); }
 
  template<class T>
  ProductRule<T>* product (Reference::To<T>& A, Reference::To<T>& B)
  { return product (A.get(), B.get()); }
}
 
 
 
 
 
#endif