FieldBasisConverter.h

#ifndef PRDC_CPU_FIELD_BASIS_CONVERTER_H
#define PRDC_CPU_FIELD_BASIS_CONVERTER_H
 
/*
* 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 <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