FieldBasisConverter.h

#ifndef PRDC_CPU_FIELD_BASIS_CONVERTER_H
#define PRDC_CPU_FIELD_BASIS_CONVERTER_H
 
/*
* PSCF Package 
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <util/containers/DMatrix.h>     // member
 
namespace Util {
   template <typename T> class DArray;
}
 
namespace Pscf {
namespace Prdc {
namespace Cpu {
 
   template <int D> class RField;
   template <int D> class RFieldDft;
 
   using namespace Util;
 
   template <int D>
   class FieldBasisConverter 
   {
 
   public:
 
      FieldBasisConverter();
 
      FieldBasisConverter(DMatrix<double> basis);
 
      virtual ~FieldBasisConverter();
 
      void setBasis(DMatrix<double> basis);
 
      double maxBasisError(double normSq = 1.0) const;
 
      void convertToBasis(DArray< RField<D> > const & in, 
                          DArray< RField<D> > & out,
                          double prefactor = 1.0) const;
 
      void convertFromBasis(DArray< RField<D> > const & in, 
                            DArray< RField<D> > & out,
                            double prefactor = 1.0) const;
 
      int nMonomer() const
      {  return nMonomer_; }
 
      bool hasBasis() const
      {  return (nMonomer_ > 0); }
 
   private:
 
      DMatrix<double> basis_;
 
      int nMonomer_;
 
   };
 
   #ifndef PRDC_CPU_FIELD_BASIS_CONVERTER_TPP
   // Suppress implicit instantiation
   extern template class FieldBasisConverter<1>;
   extern template class FieldBasisConverter<2>;
   extern template class FieldBasisConverter<3>;
   #endif
 
} // namespace Pscf::Prdc::Cpu
} // namespace Pscf::Prdc
} // namespace Pscf
#endif