fd1d/solvers/Mixture.h

#ifndef FD1D_MIXTURE_H
#define FD1D_MIXTURE_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 "Polymer.h"
#include "Solvent.h"
#include <pscf/solvers/MixtureTmpl.h>
#include <pscf/inter/Interaction.h>
#include <util/containers/DArray.h>
 
namespace Pscf {
namespace Fd1d
{
 
   class Domain;
 
   class Mixture : public MixtureTmpl<Polymer, Solvent>
   {
 
   public:
 
      // Public typedefs
 
      typedef Propagator::WField WField;
 
      typedef Propagator::CField CField;
 
      // Public member functions
 
      Mixture();
 
      ~Mixture();
 
      void readParameters(std::istream& in);
 
      void setDomain(Domain const & domain);
 
      void setKuhn(int monomerId, double kuhn);
 
      void 
      compute(DArray<WField> const & wFields, DArray<CField>& 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 Fd1d
} // namespace Pscf
#endif