Molecule.h

#ifndef PSCF_FLORY_HUGGINS_MOLECULE_H
#define PSCF_FLORY_HUGGINS_MOLECULE_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 <pscf/floryHuggins/Clump.h>      // member template argument
#include <util/containers/Pair.h>        // member template
#include <util/containers/DArray.h>      // member template
 
#include <cmath>
 
namespace Pscf { 
   namespace FloryHuggins { 
   
      using namespace Util;
   
      class Molecule : public Util::ParamComposite
      {
   
      public:
   
         Molecule();
    
         ~Molecule();
   
         virtual void readParameters(std::istream& in);

         void setNClump(int nClump);
   
         void computeSize();
   

   
         Clump& clump(int id);
   
         Clump const & clump(int id) const;
   
         int nClump() const; 
   
         double size() const;
   
   
      private:
   
         DArray<Clump> clumps_;
   
         int nClump_;
   
         double size_;

         bool hasSize_;
   
      };
   
   }
 
   /*
   * Number of clumps.
   */
   inline int FloryHuggins::Molecule::nClump() const
   {  return nClump_; }
 
   /*
   * Total size of all clumps = volume / reference volume
   */
   inline double FloryHuggins::Molecule::size() const
   {
      UTIL_CHECK(hasSize_);
      return size_; 
   }
 
   /*
   * Get a specified Clump (non-constant reference)
   */
   inline FloryHuggins::Clump& FloryHuggins::Molecule::clump(int id)
   {  return clumps_[id]; }
 
   /*
   * Get a specified Clump (constant reference)
   */
   inline const FloryHuggins::Clump& FloryHuggins::Molecule::clump(int id) const
   {  return clumps_[id]; }
 
}
#endif