MixtureTmpl.h

#ifndef PSCF_MIXTURE_TMPL_H
#define PSCF_MIXTURE_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/MixtureBase.h>       // base class
#include <util/param/ParamComposite.h>   // base class
#include <util/containers/DArray.h>      // member template
 
namespace Pscf
{
 
   using namespace Util;

   template <class PT, class ST>
   class MixtureTmpl : public MixtureBase, public ParamComposite
   {
   public:
 
      // Public typedefs

      using SolventT = ST;

      using PolymerT = PT;
 
      // Public member functions

      MixtureTmpl();

      ~MixtureTmpl();

      virtual void readParameters(std::istream& in);

 
      /*
      * Public member functions inherited from MixtureBase:
      *
      *   int nMonomer() const;
      *   int nPolymer() const;
      *   int nSolvent() const;
      *   int nBlock() const;
      *   double vMonomer() const;
      *   Monomer const & monomer(int id) const;
      *
      *   void setVmonomer(double);
      */

      PolymerT& polymer(int id);

      PolymerT const & polymer(int id) const;

      PolymerSpecies const & polymerSpecies(int id) const final;

      SolventT& solvent(int id);

      SolventT const & solvent(int id) const;

      SolventSpecies const & solventSpecies(int id) const final;

 
   private:

      DArray<PolymerT> polymers_;

      DArray<SolventT> solvents_;
 
      // Restrict access to inherited protected data
      using MixtureBase::monomers_;
      using MixtureBase::nMonomer_;
      using MixtureBase::nPolymer_;
      using MixtureBase::nSolvent_;
      using MixtureBase::nBlock_;
      using MixtureBase::vMonomer_;
 
   };
 
   // Inline member functions
 
   template <class PT, class ST>
   inline PT& MixtureTmpl<PT,ST>::polymer(int id)
   {  
      UTIL_CHECK(id < nPolymer_);
      return polymers_[id];
   }
 
   template <class PT, class ST>
   inline PT const & MixtureTmpl<PT,ST>::polymer(int id) const
   {  
      UTIL_CHECK(id < nPolymer_);
      return polymers_[id];
   }
 
   template <class PT, class ST>
   inline ST& MixtureTmpl<PT,ST>::solvent(int id)
   {  
      UTIL_CHECK(id < nSolvent_);
      return solvents_[id]; 
   }
 
   template <class PT, class ST>
   inline ST const & MixtureTmpl<PT,ST>::solvent(int id) const
   {  
      UTIL_CHECK(id < nSolvent_);
      return solvents_[id]; 
   }
 
}
#endif