rpg/field/WFields.h

#ifndef RPG_W_FIELDS_H
#define RPG_W_FIELDS_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 <prdc/field/WFieldsTmpl.h>     // base class template
#include <prdc/cpu/RField.h>            // template parameter
#include <rpg/field/FieldIo.h>          // template parameter
 
namespace Pscf {
namespace Rpg {
 
   using namespace Util;
   using namespace Pscf::Prdc;
   //using namespace Pscf::Prdc::Cuda;

   template <int D>
   class WFields : public WFieldsTmpl<D, Cuda::RField<D>, FieldIo<D> >
   {
 
   public:

      typedef WFieldsTmpl<D, RField<D>, FieldIo<D> > Base;
 
      // Inherited public member functions
      using Base::setFieldIo;
      using Base::setNMonomer;
      using Base::allocateRGrid;
      using Base::allocateBasis;
      using Base::allocate;
      using Base::setBasis;
      using Base::setRGrid;
      using Base::readBasis;
      using Base::readRGrid;
      using Base::symmetrize;
      using Base::clear;
      using Base::writeBasis;
      using Base::writeRGrid;
      using Base::basis;
      using Base::rgrid;
      using Base::isAllocatedRGrid;
      using Base::isAllocatedBasis;
      using Base::hasData;
      using Base::isSymmetric;

      void setRGrid(DeviceArray<cudaReal>& fields);
 
   protected:
 
      using Base::meshDimensions;
      using Base::meshSize;
      using Base::nBasis;
      using Base::nMonomer;
      using Base::fieldIo;
 
   private:

      void assignRField(RField<D>& lhs, RField<D> const & rhs) const 
      override;
 
   };
 
   // Explicit instantiation declarations
   extern template class WFields<1>;
   extern template class WFields<2>;
   extern template class WFields<3>;
 
} // namespace Rpg
 
namespace Prdc {
   // Explicit instantiation declarations for base class
   extern template class WFieldsTmpl<1, Cuda::RField<1>, Rpg::FieldIo<1> >;
   extern template class WFieldsTmpl<2, Cuda::RField<2>, Rpg::FieldIo<2> >;
   extern template class WFieldsTmpl<3, Cuda::RField<3>, Rpg::FieldIo<3> >;
} // namespace Prdc
 
} // namespace Pscf
#endif