rp/solvers/Solvent.h

#ifndef RP_SOLVENT_H
#define RP_SOLVENT_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/SolventSpecies.h>  // base class
 
// Forward declaration
namespace Pscf {
   template <int D> class Mesh;
}
 
namespace Pscf {
namespace Rp {

   template <int D, class T>
   class Solvent : public SolventSpecies<double>
   {
 
   public:

      using SolventSpeciesT = SolventSpecies<double>;

      using SpeciesT = Species<double>;
 
      // Public member functions

      void associate(Mesh<D> const & mesh);

      void allocate();

      void compute(typename T::RField const & wField, double phiTot = 1.0);

      typename T::RField const & cField() const;
 
   protected:

      Solvent();

      ~Solvent();
 
   private:

      typename T::RField cField_;

      Mesh<D> const *  meshPtr_;
 
   };
 
   /*
   * Get monomer concentration field for this solvent.
   */
   template <int D, class T> inline
   typename T::RField const & Solvent<D,T>::cField() const
   {  return cField_;  }
 
}
}
#endif