r1d/domain/Domain.h

#ifndef R1D_DOMAIN_H
#define R1D_DOMAIN_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 <util/param/ParamComposite.h>     // base class
#include "GeometryMode.h"                  // member
#include <util/containers/DArray.h>        // function parameter
 
namespace Pscf {
namespace R1d {
 
   using namespace Util;
 
   class Domain : public ParamComposite
   {
 
   public:
 
      typedef DArray<double> Field;
 
      Domain();
 
      ~Domain();
 
 
      void readParameters(std::istream& in);
 
      void setPlanarParameters(double xMin, double xMax, int nx);
 
      void setShellParameters(GeometryMode mode, 
                              double xMin, double xMax, 
                              int nx);
 
      void setCylinderParameters(double xMax, int nx);
 
      void setSphereParameters(double xMax, int nx);
 
 
      double xMin() const;
 
      double xMax() const;
 
      double volume() const;
 
      double dx() const;
 
      int nx() const;
 
      GeometryMode const & mode() const;
 
      bool isShell() const;
 
 
      double spatialAverage(Field const & f) const;
 
      double innerProduct(Field const & f, Field const & g) const;
 
 
   private:
 
      double xMin_;
 
      double xMax_;
 
      double dx_;
 
      double volume_;
 
      int nx_;
 
      GeometryMode mode_;
 
      bool isShell_;
 
      mutable DArray<double> work_;
 
      void computeVolume();
 
   };
 
   // Inline member functions
 
   inline int Domain::nx() const
   {  return nx_; }
 
   inline double Domain::dx() const
   {  return dx_; }
 
   inline double Domain::xMin() const
   {  return xMin_; }
 
   inline double Domain::xMax() const
   {  return xMax_; }
 
   inline double Domain::volume() const
   {  return volume_; }
 
   inline GeometryMode const & Domain::mode() const
   {  return mode_; }
 
   inline bool Domain::isShell() const
   {  return isShell_; }
 
} // namespace R1d
} // namespace Pscf
#endif