cpu/CField.h

#ifndef PRDC_CPU_C_FIELD_H
#define PRDC_CPU_C_FIELD_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/cpu/FftwDArray.h>    // base class
#include <pscf/math/IntVec.h>       // member
#include <fftw3.h>
 
namespace Pscf {
namespace Prdc {
namespace Cpu {
 
   using namespace Util;
   using namespace Pscf;

   template <int D>
   class CField : public FftwDArray<fftw_complex>
   {
 
   public:
 
      // Type aliases
 
      using FftwDArray<fftw_complex>::ValueType;

      using RealType = double;
 
      // Member functions

      CField();

      CField(const CField& other);

      virtual ~CField();

      CField& operator = (const CField& other);

      void allocate(const IntVec<D>& meshDimensions);

      virtual void deallocate();

      const IntVec<D>& meshDimensions() const;

      template <class Archive>
      void serialize(Archive& ar, const unsigned int version);
 
   private:
 
      // Vector containing number of grid points in each direction.
      IntVec<D> meshDimensions_;
 
      using FftwDArray<fftw_complex>::allocate;
 
   };
 
   /*
   * Return mesh dimensions by constant reference.
   */
   template <int D>
   inline const IntVec<D>& CField<D>::meshDimensions() const
   {  return meshDimensions_; }
 
   /*
   * Serialize a Field to/from an Archive.
   */
   template <int D>
   template <class Archive>
   void CField<D>::serialize(Archive& ar, const unsigned int version)
   {
      FftwDArray<fftw_complex>::serialize(ar, version);
      ar & meshDimensions_;
   }
 
   // Explicit instantiation declarations
   extern template class CField<1>;
   extern template class CField<2>;
   extern template class CField<3>;
 
}
}
}
#endif