r1d/solvers/Solvent.h

#ifndef R1D_SOLVENT_H
#define R1D_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
#include <r1d/solvers/Propagator.h>     // typedefs
 
namespace Pscf { 
namespace R1d { 
 
   class Domain;
   using namespace Util;

   class Solvent : public SolventSpecies
   {
 
   public:
 
      // Public typename alias

      using FieldT = typename Propagator::FieldT;
 
      // Public member functions

      Solvent();
   
      ~Solvent();
  
      void setDiscretization(Domain const & domain);

      void compute(FieldT const & wField );


      Domain const & domain() const;

      const FieldT& cField() const;
  
 
   private:
 
      // Concentration field for this solvent
      FieldT cField_;
 
      // Pointer to associated domain
      Domain const *  domainPtr_;
 
   };
   
   // Inline member functions

   inline const typename Solvent::FieldT& Solvent::cField() const
   {  return cField_;  }

   inline Domain const & Solvent::domain() const
   {   
      UTIL_ASSERT(domainPtr_);
      return *domainPtr_;
   }
 
}
}
#endif