r1d/sweep/Sweep.h

#ifndef R1D_SWEEP_H
#define R1D_SWEEP_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 <pscf/sweep/SweepTmpl.h>             // base class
#include <r1d/SystemAccess.h>                // base class
#include <r1d/misc/HomogeneousComparison.h>  // member
#include <r1d/misc/FieldIo.h>                // member
#include <util/containers/DArray.h>           // member
 
#include <util/global.h>
 
namespace Pscf {
namespace R1d
{
 
   using namespace Util;

   class Sweep : public SweepTmpl< DArray<System::WField> >, 
                 public SystemAccess
   {
 
   public:
 
      typedef DArray<System::WField> State;
 
      typedef SweepTmpl<State>  Base;

      Sweep(System& system);

      ~Sweep();

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

      virtual void outputSolution(std::string const & stateFileName);

      virtual void outputSummary(std::ostream& outFile);
 
   protected:

      int homogeneousMode_;

      virtual void checkAllocation(State & state);

      virtual void setup();

      virtual void setParameters(double s);

      virtual void extrapolate(double s);

      virtual int solve(bool isContinuation);

      virtual void reset();

      virtual void getSolution();

      virtual void cleanup();
 
   private:

      HomogeneousComparison comparison_;

      FieldIo fieldIo_;

      std::ofstream logFile_;

      void assignFields(State& lhs, State const & rhs) const;
 
   };
 
} // namespace R1d
} // namespace Pscf
#endif