BlockTmpl.h

#ifndef PSCF_BLOCK_TMPL_H
#define PSCF_BLOCK_TMPL_H
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <pscf/chem/Edge.h>         // base class
#include <util/containers/Pair.h>   // member template
 
#include <cmath>
 
namespace Pscf
{ 
 
   using namespace Util;
 
   template <class TP>
   class BlockTmpl : public Edge
   {
 
   public:
 
      // Modified diffusion equation propagator for one block direction.
      typedef TP Propagator;
 
      // Monomer concentration field.
      typedef typename TP::CField CField;
 
      // Chemical potential field.
      typedef typename TP::WField WField;
 
      BlockTmpl();
 
      virtual ~BlockTmpl();
 
      virtual void setKuhn(double kuhn);
 
      TP& propagator(int directionId);
   
      TP const & propagator(int directionId) const;
   
      typename TP::CField& cField();
 
      typename TP::CField const & cField() const;
   
      double kuhn() const;
  
   private:
 
      Pair<Propagator> propagators_;
 
      CField cField_;
 
      double kuhn_;
 
   };
 
   // Inline member functions
 
   /*
   * Get a Propagator indexed by direction.
   */
   template <class TP>
   inline 
   TP& BlockTmpl<TP>::propagator(int directionId)
   {  return propagators_[directionId]; }
 
   /*
   * Get a const Propagator indexed by direction.
   */
   template <class TP>
   inline 
   TP const & BlockTmpl<TP>::propagator(int directionId) const
   {  return propagators_[directionId]; }
 
   /*
   * Get the monomer concentration field.
   */
   template <class TP>
   inline
   typename TP::CField& BlockTmpl<TP>::cField()
   {  return cField_; }
 
   /*
   * Get the const monomer concentration field.
   */
   template <class TP>
   inline
   typename TP::CField const & BlockTmpl<TP>::cField() const
   {  return cField_; }
 
   /*
   * Get the monomer statistical segment length. 
   */
   template <class TP>
   inline double BlockTmpl<TP>::kuhn() const
   {  return kuhn_; }
 
   // Non-inline functions
 
   /*
   * Constructor.
   */
   template <class TP>
   BlockTmpl<TP>::BlockTmpl()
    : propagators_(),
      cField_(),
      kuhn_(0.0)
   {
      propagator(0).setDirectionId(0);
      propagator(1).setDirectionId(1);
      propagator(0).setPartner(propagator(1));
      propagator(1).setPartner(propagator(0));
   }
 
   /*
   * Destructor.
   */
   template <class TP>
   BlockTmpl<TP>::~BlockTmpl()
   {}
 
   /*
   * Set the monomer statistical segment length.
   */
   template <class TP>
   void BlockTmpl<TP>::setKuhn(double kuhn)
   {  kuhn_ = kuhn; }
 
}
#endif