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/DArray.h>   // member
 
namespace Pscf {
 
   using namespace Util;

   template <class QT, class FT>
   class BlockTmpl : public Edge
   {
 
   public:
 
      // Public typename aliases

      using PropagatorT = QT;

      using FieldT = FT;
 
      // Public member functions
 
      // Protected constructor and destructor (see below).

      virtual void setKuhn(double kuhn);

      QT& propagator(int directionId);

      QT const & propagator(int directionId) const;

      FT& cField();

      FT const & cField() const;

      double kuhn() const;
 
   protected:

      BlockTmpl();

      ~BlockTmpl() = default;
 
   private:

      DArray<QT> propagators_;

      FT cField_;

      double kuhn_;
 
   };
 
   // Inline member functions
 
   /*
   * Get a Propagator indexed by direction.
   */
   template <class QT, class FT>
   inline
   QT& BlockTmpl<QT,FT>::propagator(int directionId)
   {  return propagators_[directionId]; }
 
   /*
   * Get a const Propagator indexed by direction.
   */
   template <class QT, class FT>
   inline
   QT const & BlockTmpl<QT,FT>::propagator(int directionId) const
   {  return propagators_[directionId]; }
 
   /*
   * Get the monomer concentration field.
   */
   template <class QT, class FT> inline
   FT& BlockTmpl<QT,FT>::cField()
   {  return cField_; }
 
   /*
   * Get the const monomer concentration field.
   */
   template <class QT, class FT> inline
   FT const & BlockTmpl<QT,FT>::cField() const
   {  return cField_; }
 
   /*
   * Get the monomer statistical segment length.
   */
   template <class QT, class FT> inline
   double BlockTmpl<QT,FT>::kuhn() const
   {  return kuhn_; }
 
}
#endif