r1d/solvers/Solvent.h

#ifndef R1D_SOLVENT_H
#define R1D_SOLVENT_H
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <pscf/chem/SolventDescriptor.h>   // base class
//#include <pscf/chem/Species.h>             // base class
//#include <util/param/ParamComposite.h>     // base class
#include <r1d/solvers/Propagator.h>
#include <r1d/domain/Domain.h>
 
namespace Pscf { 
namespace R1d { 
 
   using namespace Util;
 
   //class Solvent : public Species, public ParamComposite
   class Solvent : public SolventDescriptor
   {
 
   public:
 
      typedef typename Propagator::CField CField;
 
      typedef typename Propagator::WField WField;
 
      Solvent();
   
      ~Solvent();
  
      void setDiscretization(Domain const & domain);
 
      void compute(WField const & wField );
 
 
      Domain const & domain() const;
 
      const CField& cField() const;
  
 
      // Inherited accessor functions 
      using Pscf::Species::phi;
      using Pscf::Species::mu;
      using Pscf::Species::q;
      using Pscf::Species::ensemble;
 
   protected:
 
      // Inherited data members
      using Pscf::Species::phi_;
      using Pscf::Species::mu_;
      using Pscf::Species::q_;
      using Pscf::Species::ensemble_;
      using Pscf::SolventDescriptor::monomerId_;
      using Pscf::SolventDescriptor::size_;
 
   private:
 
      // Concentration field for this solvent
      CField cField_;
 
      // Pointer to associated domain
      Domain const *  domainPtr_;
 
   };
   
   // Inline member functions
 
   inline const typename Solvent::CField& Solvent::cField() const
   {  return cField_;  }
 
   inline Domain const & Solvent::domain() const
   {   
      UTIL_ASSERT(domainPtr_);
      return *domainPtr_;
   }
 
}
}
#endif