r1d/solvers/Mixture.h

#ifndef R1D_MIXTURE_H
#define R1D_MIXTURE_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/solvers/MixtureTmpl.h>         // base class
 
// Forward declarations
namespace Util {
   template <typename T> class DArray;
}
namespace Pscf {
   namespace R1d {
      class Domain;
      class Polymer;
      class Solvent;
   }
   extern template class MixtureTmpl<R1d::Polymer, R1d::Solvent>;
}
 
namespace Pscf {
namespace R1d {

   class Mixture : public MixtureTmpl<Polymer, Solvent>
   {
 
   public:
 
      // Public typedefs

      using FieldT = DArray<double>;
 
      // Public member functions

      Mixture();

      ~Mixture();

      void readParameters(std::istream& in);

      void setDomain(Domain const & domain);

      void setKuhn(int monomerId, double kuhn);

      void 
      compute(DArray<FieldT> const & wFields, DArray<FieldT>& cFields);
 
 
      // Inherited public member functions with non-dependent names
      using MixtureTmpl< Polymer, Solvent >::nMonomer;
      using MixtureTmpl< Polymer, Solvent >::nPolymer;
      using MixtureTmpl< Polymer, Solvent >::nSolvent;
      using MixtureTmpl< Polymer, Solvent >::nBlock;
      using MixtureTmpl< Polymer, Solvent >::polymer;
      using MixtureTmpl< Polymer, Solvent >::monomer;
      using MixtureTmpl< Polymer, Solvent >::solvent;
      using MixtureTmpl< Polymer, Solvent >::vMonomer;
 
   protected:
 
      // Inherited protected member functions with non-dependent names
      using ParamComposite::setClassName;
      using ParamComposite::read;
      using ParamComposite::readOptional;
 
   private:

      double ds_;

      Domain const * domainPtr_;

      Domain const & domain() const;
 
   };
 
   // Inline member function
 
   /*
   * Get Domain by constant reference (private).
   */
   inline Domain const & Mixture::domain() const
   {   
      UTIL_ASSERT(domainPtr_);
      return *domainPtr_;
   }
 
} // namespace R1d
} // namespace Pscf
#endif