r1d/solvers/Polymer.cpp

/*
* 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 "Polymer.h"
#include "Block.h"
#include "Propagator.h"
#include <pscf/solvers/PolymerTmpl.tpp>
 
// Explicit instantiation of base class
namespace Pscf {
   template class PolymerTmpl<R1d::Block, R1d::Propagator>;
}
 
namespace Pscf {
namespace R1d { 
 
   /*
   * Constructor
   */
   Polymer::Polymer()
   {  setClassName("Polymer"); }
 
   /*
   * Destructor
   */
   Polymer::~Polymer()
   {}
 
   /*
   * Compute solution to MDE and concentrations.
   */ 
   void Polymer::compute(DArray< DArray<double> > const & wFields)
   {
 
      // Setup solvers for all blocks
      int monomerId;
      for (int j = 0; j < nBlock(); ++j) {
         monomerId = block(j).monomerId();
         block(j).setupSolver(wFields[monomerId]);
      }
 
      // Solve MDE for all propagators
      solve();
 
      // Compute block concentration fields (thread model)
      double prefactor;
      prefactor = phi() / ( q() * length() );
      for (int i = 0; i < nBlock(); ++i) {
         block(i).computeConcentration(prefactor);
      }
 
   }
 
}
}