cpu/RFieldDftComparison.h

#ifndef PRDC_CPU_R_FIELD_DFT_COMPARISON_H
#define PRDC_CPU_R_FIELD_DFT_COMPARISON_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.
*/
 
// Forward declarations
namespace Util {
   template <typename T> class DArray;
}
namespace Pscf {
   namespace Prdc {
      namespace Cpu {
         template <int D> class RFieldDft;
      }   
   }   
}
 
namespace Pscf {
namespace Prdc {
namespace Cpu {
 
   using namespace Util;

   template <int D>
   class RFieldDftComparison {
 
   public:

      RFieldDftComparison();
 
      // Use compiler defined destructor and assignment operator.
 
      double compare(RFieldDft<D> const& a, RFieldDft<D> const& b);
 
      double compare(DArray<RFieldDft<D> > const& a, 
                     DArray<RFieldDft<D> > const& b);

      double maxDiff() const
      {  return maxDiff_; }
   
      double rmsDiff() const
      {  return rmsDiff_; }
   
   private:
  
      // Maximum element-by-element difference. 
      double maxDiff_;
 
      // Room-mean-squared element-by-element difference. 
      double rmsDiff_;
   
   };
 
   // Explicit instantiation declarations
   extern template class RFieldDftComparison<1>;
   extern template class RFieldDftComparison<2>;
   extern template class RFieldDftComparison<3>;
 
} // namespace Cpu
} // namespace Prdc
} // namespace Pscf
#endif