cp/field/Domain.h

#ifndef PRDC_CL_DOMAIN_H
#define PRDC_CL_DOMAIN_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 <util/param/ParamComposite.h>    // base class
#include <prdc/crystal/UnitCell.h>        // member
#include <pscf/mesh/Mesh.h>               // member
 
// Forward declarations
namespace Util {
   class FileMaster;
   template <typename T> class Signal;
   template <> class Signal<void>;
}
 
namespace Pscf {
namespace Cp {
 
   using namespace Util;
   using namespace Prdc;

   template <int D, class FFT, class WLT, class FIT>
   class Domain : public ParamComposite
   {
 
   public:

      Domain();

      ~Domain();


      void setFileMaster(FileMaster& fileMaster);

      virtual void readParameters(std::istream& in);

      void readFieldHeader(std::istream& in, int& nMonomer);


      Mesh<D>& mesh();

      Mesh<D> const & mesh() const;

      UnitCell<D>& unitCell();

      UnitCell<D> const & unitCell() const;

      FFT& fft();

      FFT const & fft() const;

      WLT& waveList();

      WLT const & waveList() const;

      FIT& fieldIo();

      FIT const & fieldIo() const;

      typename UnitCell<D>::LatticeSystem lattice() const;

 
   private:
 
      // Private member variables

      Mesh<D> mesh_;

      UnitCell<D> unitCell_;

      typename UnitCell<D>::LatticeSystem lattice_;

      FFT* fftPtr_;

      WLT* waveListPtr_;

      FIT* fieldIoPtr_;

      Signal<void>* signalPtr_;

      FileMaster* fileMasterPtr_;
 
      // Boolean flags

      bool isInitialized_;
 
      // Private member function
 
      /*
      * Get FileMaster as const reference.
      */
      FileMaster const & fileMaster() const
      {
         UTIL_CHECK(fileMasterPtr_);
         return * fileMasterPtr_;
      }
 
      // Members of parent class with non-dependent names
      using ParamComposite::read;
      using ParamComposite::readOptional;
 
   };
 
   // Inline member functions
 
   // Get the UnitCell by non-const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline UnitCell<D>& Domain<D,FFT,WLT,FIT>::unitCell()
   {  return unitCell_; }
 
   // Get the UnitCell by const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline UnitCell<D> const & Domain<D,FFT,WLT,FIT>::unitCell() const
   {  return unitCell_; }
 
   // Get the Mesh by non-const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline Mesh<D>& Domain<D,FFT,WLT,FIT>::mesh()
   {  return mesh_; }
 
   // Get the Mesh by const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline Mesh<D> const & Domain<D,FFT,WLT,FIT>::mesh() const
   {  return mesh_; }
 
   // Get the FFT by non-const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline FFT& Domain<D,FFT,WLT,FIT>::fft()
   {  return *fftPtr_; }
 
   // Get the FFT by const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline FFT const & Domain<D,FFT,WLT,FIT>::fft() const
   {  return *fftPtr_; }
 
   // Get the WaveList by non-const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline WLT& Domain<D,FFT,WLT,FIT>::waveList()
   {  return *waveListPtr_; }
 
   // Get the WaveList by const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline WLT const & Domain<D,FFT,WLT,FIT>::waveList() const
   {  return *waveListPtr_; }
 
   // Get the FieldIo by const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline FIT& Domain<D,FFT,WLT,FIT>::fieldIo()
   {  return *fieldIoPtr_; }
 
   // Get the FieldIo by const reference.
   template <int D, class FFT, class WLT, class FIT>
   inline FIT const & Domain<D,FFT,WLT,FIT>::fieldIo() const
   {  return *fieldIoPtr_; }
 
   // Get the lattice system enumeration value.
   template <int D, class FFT, class WLT, class FIT>
   inline
   typename UnitCell<D>::LatticeSystem Domain<D,FFT,WLT,FIT>::lattice()
   const
   {  return lattice_; }
 
} // namespace Cp
} // namespace Pscf
#endif