SystemAccess.h

#ifndef R1D_SYSTEM_ACCESS_H
#define R1D_SYSTEM_ACCESS_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 "System.h"                        // member
 
namespace Pscf {
namespace R1d {
 
   using namespace Util;

   class SystemAccess 
   {
 
   public:

      SystemAccess();

      SystemAccess(System& system);

      ~SystemAccess();

      virtual void setSystem(System& system);


      const System& system() const;

      System& system();

      const Mixture& mixture() const;

      Mixture& mixture();

      const Domain& domain() const;

      Domain& domain();

      const Interaction& interaction() const;

      Interaction& interaction();

      FileMaster& fileMaster();


      DArray<System::WField>& wFields();

      System::WField& wField(int monomerId);

      DArray<System::CField>& cFields();

      System::CField& cField(int monomerId);

 
   private:
 
      System* systemPtr_;
 
   };
 
   // Inline member functions
 
   /*
   * Get the parent System object.
   */
   inline const System& SystemAccess::system() const
   { 
      UTIL_ASSERT(systemPtr_);
      return *systemPtr_; 
   }
 
   /*
   * Get the parent System object.
   */
   inline System& SystemAccess::system()
   { 
      UTIL_ASSERT(systemPtr_);
      return *systemPtr_; 
   }
 
   /*
   * Get the associated Mixture object.
   */
   inline const Mixture& SystemAccess::mixture() const
   { 
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->mixture(); 
   }
 
   /*
   * Get the associated Mixture object.
   */
   inline Mixture& SystemAccess::mixture()
   { 
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->mixture(); 
   }
 
   /*
   * Get the spatial Domain.
   */
   inline const Domain& SystemAccess::domain() const
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->domain(); 
   }
 
   /*
   * Get the spatial Domain.
   */
   inline Domain& SystemAccess::domain()
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->domain(); 
   }
 
   /*
   * Get the Interaction (excess free energy model).
   */
   inline const Interaction& SystemAccess::interaction() const
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->interaction(); 
   }
 
   /*
   * Get the Interaction (excess free energy model).
   */
   inline Interaction& SystemAccess::interaction()
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->interaction(); 
   }
 
   /*
   * Get the FileMaster.
   */
   inline FileMaster& SystemAccess::fileMaster()
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->fileMaster(); 
   }
 
   /*
   * Get an array of all monomer excess chemical potential fields.
   */
   inline 
   DArray< System::WField >& SystemAccess::wFields()
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->wFields(); 
   }
 
   /*
   * Get a single monomer excess chemical potential field.
   */
   inline 
   System::WField& SystemAccess::wField(int id)
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->wField(id); 
   }
 
   /*
   * Get array of all monomer concentration fields.
   */
   inline
   DArray< System::CField >& SystemAccess::cFields()
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->cFields(); 
   }
 
   /*
   * Get a single monomer concentration field.
   */
   inline System::CField& SystemAccess::cField(int id)
   {  
      UTIL_ASSERT(systemPtr_);
      return systemPtr_->cField(id); 
   }
 
} // namespace R1d
} // namespace Pscf
#endif