MixtureBase.h

#ifndef PSCF_MIXTURE_BASE_H
#define PSCF_MIXTURE_BASE_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/Monomer.h>
#include <util/containers/DArray.h>
 
namespace Pscf
{
 
   class PolymerSpecies;
   class SolventSpecies;
 
   using namespace Util;

   class MixtureBase 
   {
 
   public:

      MixtureBase();

      ~MixtureBase();

      void setVmonomer(double vMonomer);


      int nMonomer() const;

      int nPolymer() const;

      int nSolvent() const;

      int nBlock() const;

      double vMonomer() const;

      bool isCanonical() const;


      Monomer const & monomer(int id) const;

      virtual
      PolymerSpecies const & polymerSpecies(int id) const = 0;

      virtual
      SolventSpecies const & solventSpecies(int id) const = 0;

 
   protected:

      Monomer& monomer(int id);

      DArray<Monomer> monomers_;

      int nMonomer_; 

      int nPolymer_;

      int nSolvent_;

      int nBlock_;

      double vMonomer_;
 
   };
 
   // Inline public member functions
 
   inline int MixtureBase::nMonomer() const
   {  return nMonomer_; }
 
   inline int MixtureBase::nPolymer() const
   {  return nPolymer_; }
 
   inline int MixtureBase::nSolvent() const
   {  return nSolvent_; }
 
   inline int MixtureBase::nBlock() const
   {  return nBlock_; }
 
   inline 
   double MixtureBase::vMonomer() const
   {  return vMonomer_; }
 
   inline 
   Monomer const & MixtureBase::monomer(int id) const
   {  
      UTIL_CHECK(id < nMonomer_);
      return monomers_[id]; 
   }
 
   // Inline protected member function
 
   inline 
   Monomer& MixtureBase::monomer(int id)
   {  
      UTIL_CHECK(id < nMonomer_);
      return monomers_[id]; 
   }
 
}
#endif