rpg/fts/perturbation/EinsteinCrystalPerturbation.h

#ifndef RPG_EINSTEIN_CRYSTAL_PERTURBATION_H
#define RPG_EINSTEIN_CRYSTAL_PERTURBATION_H
 
#include "Perturbation.h"      // base class
#include <prdc/cuda/RField.h>
#include <util/containers/DArray.h> 
 
namespace Pscf {
namespace Rpg {
 
   using namespace Util;
   using namespace Prdc;
   using namespace Prdc::Cuda;
 
   template <int D> class Simulator;
 
   template <int D>
   class EinsteinCrystalPerturbation : public Perturbation<D>
   {
 
   public:
 
      EinsteinCrystalPerturbation(Simulator<D>& simulator);
 
      virtual ~EinsteinCrystalPerturbation();
 
      virtual void readParameters(std::istream& in);
 
      virtual void setup();
 
      virtual double hamiltonian(double unperturbedHamiltonian);
      
      virtual void incrementDc(DArray< RField<D> >& dc);
      
      virtual double df();
      
      virtual void saveState();
 
      virtual void restoreState();
      
      // Inherited functions
      using ParamComposite::setClassName;
      using ParamComposite::read;
      using ParamComposite::readOptionalDArray;
      using Perturbation<D>::simulator;
      using Perturbation<D>::system;
      using Perturbation<D>::lambda;
      
   protected:
      
      // Inherited protected data members
      using Perturbation<D>::lambda_;
      
   private:
      
      // Parameters used in Einstein crystal integration
      DArray<double> epsilon_;
      
      // Reference w field
      DArray< RField<D> > w0_;
      
      // Eigenvector components of the reference w fields
      DArray< RField<D> > wc0_;
      
      // Current Einstein crystal Hamiltonian 
      double ecHamiltonian_;
      
      // Current unperturbed Hamiltonian 
      double unperturbedHamiltonian_;
      
      // Saved Einstein crystal Hamiltonian  
      double stateEcHamiltonian_;
      
      // Saved unperturbed Hamiltonian 
      double stateUnperturbedHamiltonian_;
      
      // Have epsilon values been set?
      bool hasEpsilon_;
      
      // Reference FieldFileName
      std::string referenceFieldFileName_;
      
      // Compute eigenvector components of the reference field
      void computeWcReference();
      
   };
 
   #ifndef RPG_EINSTEIN_CRYSTAL_PERTURBATION_TPP
   // Suppress implicit instantiation
   extern template class EinsteinCrystalPerturbation<1>;
   extern template class EinsteinCrystalPerturbation<2>;
   extern template class EinsteinCrystalPerturbation<3>;
   #endif
 
}
}
#endif