BlockTmpl.h

#ifndef PSCF_BLOCK_TMPL_H
#define PSCF_BLOCK_TMPL_H
 
/*
* PSCF - Polymer Self-Consistent Field
*
* Copyright 2015 - 2025, 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 QT>
   class BlockTmpl : public Edge
   {
 
   public:
 
      // Public typename alias

      using PropagatorT = QT;
 
      // Public member functions

      BlockTmpl();
 
      virtual ~BlockTmpl();

      virtual void setKuhn(double kuhn);

      QT& propagator(int directionId);
   
      QT const & propagator(int directionId) const;
   
      typename QT::FieldT& cField();

      typename QT::FieldT const & cField() const;
   
      double kuhn() const;
  
   private:

      Pair<PropagatorT> propagators_;

      typename QT::FieldT cField_;

      double kuhn_;
 
   };
 
   // Inline member functions
 
   /*
   * Get a Propagator indexed by direction.
   */
   template <class QT>
   inline 
   QT& BlockTmpl<QT>::propagator(int directionId)
   {  return propagators_[directionId]; }
 
   /*
   * Get a const Propagator indexed by direction.
   */
   template <class QT>
   inline 
   QT const & BlockTmpl<QT>::propagator(int directionId) const
   {  return propagators_[directionId]; }
 
   /*
   * Get the monomer concentration field.
   */
   template <class QT>
   inline
   typename QT::FieldT& BlockTmpl<QT>::cField()
   {  return cField_; }
 
   /*
   * Get the const monomer concentration field.
   */
   template <class QT>
   inline
   typename QT::FieldT const & BlockTmpl<QT>::cField() const
   {  return cField_; }
 
   /*
   * Get the monomer statistical segment length. 
   */
   template <class QT>
   inline double BlockTmpl<QT>::kuhn() const
   {  return kuhn_; }
 
   // Non-inline functions
 
   /*
   * Constructor.
   */
   template <class QT>
   BlockTmpl<QT>::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 QT>
   BlockTmpl<QT>::~BlockTmpl()
   {}
 
   /*
   * Set the monomer statistical segment length.
   */
   template <class QT>
   void BlockTmpl<QT>::setKuhn(double kuhn)
   {  kuhn_ = kuhn; }
 
}
#endif