FieldComparison.h

#ifndef PSCF_FIELD_COMPARISON_H
#define PSCF_FIELD_COMPARISON_H
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <util/containers/DArray.h>
 
namespace Pscf {
 
   using namespace Util;
 
   template <class FT>
   class FieldComparison {
 
   public:
 
      FieldComparison(int begin = 0);
 
      // Use compiler defined destructor and assignment operator.
 
      double compare(FT const& a, FT const& b);
 
      double compare(DArray<FT> const& a, DArray<FT> const& b);
 
      double maxDiff() const
      {  return maxDiff_; }
   
      double rmsDiff() const
      {  return rmsDiff_; }
   
   protected:
  
      // Maximum element-by-element difference. 
      double maxDiff_;
 
      // Room-mean-squared element-by-element difference. 
      double rmsDiff_;
   
      // Index of first element (0 or 1)
      int begin_;
 
   };
 
} // namespace Pscf
#include "FieldComparison.tpp"
#endif