rpc/fts/perturbation/EinsteinCrystalPerturbation.h

#ifndef RPC_EINSTEIN_CRYSTAL_PERTURBATION_H
#define RPC_EINSTEIN_CRYSTAL_PERTURBATION_H
 
#include "Perturbation.h"      // base class
 
namespace Util {
   template <typename T> class DArray;
}
 
namespace Pscf {
namespace Prdc{
namespace Cpu{
   template <int D> class RField;
}
}
}
 
namespace Pscf {
namespace Rpc {
 
   using namespace Util;
   using namespace Prdc;
   using namespace Prdc::Cpu;
 
   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_;
 
      // Reference FieldFileName
      std::string referenceFieldFileName_;
      
      // Have epsilon values been set?
      bool hasEpsilon_;
      
      // Compute eigenvector components of the reference field
      void computeWcReference();
      
   };
 
   #ifndef RPC_EINSTEIN_CRYSTAL_PERTURBATION_TPP
   // Suppress implicit instantiation
   extern template class EinsteinCrystalPerturbation<1>;
   extern template class EinsteinCrystalPerturbation<2>;
   extern template class EinsteinCrystalPerturbation<3>;
   #endif
 
}
}
#endif