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, typename WT = double>
   class MixtureTmpl : public MixtureBase<WT>, public ParamComposite
   {
   public:
 
      // Public type aliases

      using SolventT = ST;

      using PolymerT = PT;
 
      // Public member functions
 
      /*
      * 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);
      */

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


      PolymerT& polymer(int id);

      PolymerT const & polymer(int id) const;

      PolymerSpecies<WT> const & polymerSpecies(int id) const final;

      SolventT& solvent(int id);

      SolventT const & solvent(int id) const;

      SolventSpecies<WT> const & solventSpecies(int id) const final;

 
   protected:

      MixtureTmpl();

      ~MixtureTmpl() override = default;

      using  MixtureBaseT = MixtureBase<WT>;
 
   private:

      DArray<PolymerT> polymers_;

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