mitkEqual.h

Go to the documentation of this file.

/*============================================================================
 
The Medical Imaging Interaction Toolkit (MITK)
 
Copyright (c) German Cancer Research Center (DKFZ)
All rights reserved.
 
Use of this source code is governed by a 3-clause BSD license that can be
found in the LICENSE file.
 
============================================================================*/
 
#ifndef mitkEqual_h
#define mitkEqual_h
 
#include <cmath>
#include <iomanip>
 
#include "mitkLog.h"
#include "mitkNumericConstants.h"
 
namespace mitk
{
  template <typename DifferenceType>
  inline bool DifferenceBiggerOrEqualEps(DifferenceType diff, mitk::ScalarType epsilon = mitk::eps)
  {
    return std::fabs(diff) >= epsilon;
  }
 
  template <typename ElementToOutput1, typename ElementToOutput2>
  inline void ConditionalOutputOfDifference(
    ElementToOutput1 elem1, ElementToOutput2 elem2, mitk::ScalarType eps, bool verbose, bool isEqual)
  {
    if (verbose && !isEqual)
    {
      MITK_INFO << typeid(ElementToOutput1).name() << " and " << typeid(ElementToOutput2).name()
                << " not equal. Lefthandside " << std::setprecision(12) << elem1 << " - Righthandside " << elem2
                << " - epsilon " << eps;
    }
  }
 
  inline bool Equal(ScalarType scalar1, ScalarType scalar2, ScalarType eps = mitk::eps, bool verbose = false)
  {
    bool isEqual(!DifferenceBiggerOrEqualEps(scalar1 - scalar2, eps));
 
    ConditionalOutputOfDifference(scalar1, scalar2, eps, verbose, isEqual);
 
    return isEqual;
  }
}
 
#endif