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 "Polymer.h"
#include "Solvent.h"
#include <pscf/solvers/MixtureTmpl.h>
#include <pscf/inter/Interaction.h>
#include <util/containers/DArray.h>
 
namespace Pscf {
namespace R1d
{
 
   class Domain;

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

      typedef Propagator::FieldT FieldT;
 
      // 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 MixtureTmpl< Polymer, Solvent >::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