ScftReal.h

#ifndef PRDC_SCFT_REAL_H
#define PRDC_SCFT_REAL_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 "SystemConstRefReal.h"
#include <util/global.h>
#include <iostream>
 
namespace Pscf {
namespace Prdc {

   template <int D, class ST>
   class ScftReal : protected SystemConstRefReal<ST>
   {
 
   public:

      using Base = SystemConstRefReal<ST>;

      using SystemT = typename Base::SystemT;
 
      // Public member functions

      ScftReal(SystemT const & system);

      virtual ~ScftReal();


      void compute();

      void clear();


      bool hasData() const;

      double fHelmholtz() const;

      double fIdeal() const;

      double fInter() const;

      double fExt() const;

      double pressure() const;

      void write(std::ostream& out);

 
   protected:
 
      // Protected inherited type name aliases.
      using MixtureT = typename Base::MixtureT;
      using InteractionT = typename Base::InteractionT;
      using DomainT = typename Base::DomainT;
      using WFieldContainerT = typename Base::WFieldContainerT;
      using CFieldContainerT = typename Base::CFieldContainerT;
      using MaskT = typename Base::MaskT;
      using FieldT = typename Base::FieldT;
 
      // Protected inherited member functions
      using Base::system;
      using Base::mixture;
      using Base::interaction;
      using Base::domain;
      using Base::c;
      using Base::w;
      using Base::h;
      using Base::mask;
      using Base::fileMaster;

      virtual
      double innerProduct(FieldT const & A, FieldT const & B) const
      {  return 0.0; };
 
   private:

      double fHelmholtz_;

      double fIdeal_;

      double fInter_;

      double fExt_;

      double pressure_;

      bool hasData_;
 
   };
 
   // Inline member functions
 
   // Get the Helmholtz free energy per monomer / kT.
   template <int D, class ST>
   inline double ScftReal<D,ST>::fHelmholtz() const
   {
      UTIL_CHECK(hasData_);
      return fHelmholtz_;
   }
 
   // Get the ideal gas contribution to fHelmholtz.
   template <int D, class ST>
   inline double ScftReal<D,ST>::fIdeal() const
   {
      UTIL_CHECK(hasData_);
      return fIdeal_;
   }
 
   // Get the interaction contribution to fHelmholtz.
   template <int D, class ST>
   inline double ScftReal<D,ST>::fInter() const
   {
      UTIL_CHECK(hasData_);
      return fInter_;
   }
 
   // Get the external field contribution to fHelmholtz.
   template <int D, class ST>
   inline double ScftReal<D,ST>::fExt() const
   {
      UTIL_CHECK(hasData_);
      return fExt_;
   }
 
   // Get the precomputed pressure (units of kT / monomer volume).
   template <int D, class ST>
   inline double ScftReal<D,ST>::pressure() const
   {
      UTIL_CHECK(hasData_);
      return pressure_;
   }
 
   // Have free energies and pressure been computed?
   template <int D, class ST>
   inline bool ScftReal<D,ST>::hasData() const
   {  return hasData_; }
 
} // namespace Prdc
} // namespace Pscf
#endif