rpc/fts/analyzer/HamiltonianAnalyzer.tpp

#ifndef RPC_HAMILTONIAN_ANALYZER_TPP
#define RPC_HAMILTONIAN_ANALYZER_TPP
 
/*
* 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 "HamiltonianAnalyzer.h"
 
#include <rpc/System.h>
#include <rpc/fts/simulator/Simulator.h>
//#include <util/accumulators/Average.h>
 
namespace Pscf {
namespace Rpc {
 
   using namespace Util;
 
   /*
   * Constructor.
   */
   template <int D>
   HamiltonianAnalyzer<D>::HamiltonianAnalyzer(Simulator<D>& simulator, 
                                               System<D>& system)
    : AverageListAnalyzer<D>(system),
      simulatorPtr_(&simulator),
      systemPtr_(&(simulator.system())),
      hasAnalyzeChi_(false),
      idealId_(-1),
      fieldId_(-1),
      totalId_(-1)
   {  setClassName("HamiltonianAnalyzer"); }
 
   /*
   * Read interval and outputFileName. 
   */
   template <int D>
   void HamiltonianAnalyzer<D>::readParameters(std::istream& in) 
   {
      AverageListAnalyzer<D>::readParameters(in);
 
      idealId_ = 0;
      fieldId_ = 1;
      totalId_ = 2;
      AverageListAnalyzer<D>::initializeAccumulators(3);
 
      setName(idealId_, "ideal");
      setName(fieldId_, "field");
      setName(totalId_, "total");
   }
 
   /*
   * Output energy to file
   */
   template <int D>
   void HamiltonianAnalyzer<D>::compute() 
   {
      UTIL_CHECK(system().w().hasData());
      if (!system().hasCFields()) {
         system().compute();
      }
      if (!simulator().hasWc()){
         simulator().computeWc();
      }
      UTIL_CHECK(simulator().hasWc());
 
      #if 0
      if (!simulator().hasWc()){
         if (!hasAnalyzeChi_){
            simulator().analyzeChi();
            hasAnalyzeChi_ = true;
         }
         system().compute();
         simulator().computeWc();
      }
      #endif
 
      if (!simulator().hasHamiltonian()) {
         simulator().computeHamiltonian();
      }
 
      double ideal = simulator().idealHamiltonian();
      setValue(idealId_, ideal);
   
      double field = simulator().fieldHamiltonian();
      setValue(fieldId_, field);
   
      double total = simulator().hamiltonian();
      setValue(totalId_, total);
   }
   
}
}
#endif