BinaryRelaxIterator.h

#ifndef R1D_BINARY_RELAX_ITERATOR_H
#define R1D_BINARY_RELAX_ITERATOR_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 "Iterator.h"
#include <r1d/solvers/Mixture.h>
#include <pscf/math/LuSolver.h>
#include <util/containers/Array.h>
#include <util/containers/DArray.h>
#include <util/containers/DMatrix.h>
 
 
namespace Pscf {
namespace R1d
{
 
   using namespace Util;

   class BinaryRelaxIterator : public Iterator
   {
    
   public:
    
      typedef Mixture::FieldT FieldT;
    
      BinaryRelaxIterator(System& system);
    
      virtual ~BinaryRelaxIterator();
    
      void readParameters(std::istream& in);
    
      int solve(bool isContinuation = false);
   
   private:
    
    
      DArray<FieldT> wFieldsNew_;
    
      DArray<FieldT> cFieldsNew_;
      
      DArray<double> cArray_;
    
      DArray<double> wArray_;
    
      DArray<double> residual_;
    
      DArray<double> residualNew_;
      
      DArray<FieldT> dW_;
      
      DArray<FieldT> dWNew_;
      
      double dWNorm_;
      
      double dWNormNew_;
      
      double epsilon_;
      
      double lambdaPlus_;
      
      double lambdaMinus_;
      
      int maxItr_;
    
      bool isAllocated_;
    
      bool isCanonical_;
    
      void allocate();
    
      void computeDW(Array<FieldT> const & wOld, 
                     Array<FieldT> const & cFields,
                     Array<FieldT> & dW,
                     double & dWNorm);
     
      void updateWFields(Array<FieldT> const & wOld,
                         Array<FieldT> const & dW_,
                         Array<FieldT> & wNew);
                         
   };
 
} // namespace R1d
} // namespace Pscf
#endif