fd1d/sweep/Sweep.h

#ifndef FD1D_SWEEP_H
#define FD1D_SWEEP_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 <pscf/sweep/SweepTmpl.h>             // base class
#include <fd1d/SystemAccess.h>                // base class
#include <fd1d/misc/HomogeneousComparison.h>  // member
#include <fd1d/misc/FieldIo.h>                // member
#include <util/containers/DArray.h>           // member
 
#include <util/global.h>
 
namespace Pscf {
namespace Fd1d
{
 
   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 Fd1d
} // namespace Pscf
#endif