pspg/solvers/Solvent.h

#ifndef PSPG_SOLVENT_H
#define PSPG_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 <pspg/solvers/Propagator.h>       // typedefs
#include <pspg/field/RDField.h>
 
namespace Pscf {
   template <int D> class Mesh;
}
 
namespace Pscf { 
namespace Pspg { 
 
   using namespace Util;
 
   template <int D>
   class Solvent : public SolventDescriptor
   {
   public:
 
      Solvent();
   
      ~Solvent();
 
      void setDiscretization(Mesh<D> const & mesh);
   
      void compute(RDField<D> const & wField );
 
      RDField<D> const & concField() const;
  
      // Inherited public accessor functions 
      using Pscf::Species::phi;
      using Pscf::Species::mu;
      using Pscf::Species::q;
      using Pscf::Species::ensemble;
      using Pscf::SolventDescriptor::monomerId;
      using Pscf::SolventDescriptor::size;
 
   protected:
 
      // Inherited protected 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:
 
      RDField<D> concField_;
 
      Mesh<D> const *  meshPtr_;
 
   };
 
   /*
   * Get monomer concentration field for this solvent.
   */
   template <int D>
   inline RDField<D> const & Solvent<D>::concField() const
   {  return concField_;  }
 
}
} 
#endif