cpc/solvers/Polymer.h

#ifndef CPC_POLYMER_H
#define CPC_POLYMER_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/PolymerTmpl.h>  // base class template
 
// Forward declarations
namespace Util {
   template <typename T> class DArray;
}
namespace Pscf { 
   namespace Cpc {
      template <int D> class Block;
      template <int D> class Propagator;
   }
   namespace Prdc { 
      namespace Cpu { 
         template <int D> class CField;
      }
   }
}
 
// Explicit instantiation declarations for base class
namespace Pscf { 
   extern template 
   class PolymerTmpl< Cpc::Block<1>, Cpc::Propagator<1>, std::complex<double> >;
   extern template 
   class PolymerTmpl< Cpc::Block<2>, Cpc::Propagator<2>, std::complex<double> >;
   extern template 
   class PolymerTmpl< Cpc::Block<3>, Cpc::Propagator<3>, std::complex<double> >;
}
 
namespace Pscf {
namespace Cpc {
 
   using namespace Prdc;
   using namespace Pscf::Prdc::Cpu;

   template <int D>
   class Polymer 
    : public PolymerTmpl< Block<D>, Propagator<D>, std::complex<double> >
   {
 
   public:
 
      // Public type name aliases

      using Base = PolymerTmpl< Block<D>, Propagator<D>, std::complex<double> >;

      using PolymerSpeciesT = PolymerSpecies< std::complex<double> >;

      using SpeciesT = Species< std::complex<double> >;

      using BlockT = Block<D>;

      using PropagatorT = Propagator<D>;
 
      // Public member functions

      Polymer();

      ~Polymer();
 
      void clearUnitCellData();
 
      void compute(DArray< CField<D> > const & wFields);
 
      // Inherited public member functions
      
      using Base::edge;
      using Base::block;
      using Base::propagator;
      using PolymerSpeciesT::vertex;
      using PolymerSpeciesT::propagatorId;
      using PolymerSpeciesT::path;
      using PolymerSpeciesT::nBlock;
      using PolymerSpeciesT::nVertex;
      using PolymerSpeciesT::nPropagator;
      using PolymerSpeciesT::length;
      using PolymerSpeciesT::nBead;
      using PolymerSpeciesT::type;
      using SpeciesT::phi;
      using SpeciesT::mu;
      using SpeciesT::q;
      using SpeciesT::ensemble;
      using SpeciesT::setPhi;
      using SpeciesT::setMu;
 
   private: 
 
      // Restricting access to inherited functions
      using Base::solve;
      using SpeciesT::setQ;
 
   };
 
   // Explicit instantiation declarations
   extern template class Polymer<1>;
   extern template class Polymer<2>;
   extern template class Polymer<3>;
 
}
}
#endif