/********************************************************************
**                Image Component Library (ICL)                    **
**                                                                 **
** Copyright (C) 2006-2013 CITEC, University of Bielefeld          **
**                         Neuroinformatics Group                  **
** Website: www.iclcv.org and                                      **
**          http://opensource.cit-ec.de/projects/icl               **
**                                                                 **
** File   : ICLFilter/src/ICLFilter/WienerOp.cpp                   **
** Module : ICLFilter                                              **
** Authors: Christof Elbrechter, Robert Haschke                    **
**                                                                 **
**                                                                 **
** GNU LESSER GENERAL PUBLIC LICENSE                               **
** This file may be used under the terms of the GNU Lesser General **
** Public License version 3.0 as published by the                  **
**                                                                 **
** Free Software Foundation and appearing in the file LICENSE.LGPL **
** included in the packaging of this file.  Please review the      **
** following information to ensure the license requirements will   **
** be met: http://www.gnu.org/licenses/lgpl-3.0.txt                **
**                                                                 **
** The development of this software was supported by the           **
** Excellence Cluster EXC 277 Cognitive Interaction Technology.    **
** The Excellence Cluster EXC 277 is a grant of the Deutsche       **
** Forschungsgemeinschaft (DFG) in the context of the German       **
** Excellence Initiative.                                          **
**                                                                 **
********************************************************************/

#include <ICLFilter/WienerOp.h>
#include <ICLCore/Img.h>

using namespace icl::utils;
using namespace icl::core;

namespace icl {
  namespace filter{
  
    
  #ifdef ICL_HAVE_IPP
    
    namespace{
      template<typename T, IppStatus (IPP_DECL *ippiFunc) (const T*, int, T*, int, IppiSize, IppiSize, IppiPoint, float[], icl8u*)>
      void ippiWienerCall (const Img<T> *src, 
                           Img<T> *dst,
                           const Size &oMaskSize, 
                           const Point &oAnchor,
                           const Point &roiOffset,
                           std::vector<icl8u> &vBuffer, 
                           icl32f fNoise) {
        // {{{ open
        int iBufferSize;
        ippiFilterWienerGetBufferSize(dst->getROISize(), oMaskSize, 1, &iBufferSize);
        vBuffer.reserve (iBufferSize);
        
        for(int c=src->getChannels()-1; c>=0; --c) {
          ippiFunc(src->getROIData (c, roiOffset), src->getLineStep(),
                   dst->getROIData (c), dst->getLineStep(),
                   dst->getROISize(), oMaskSize, oAnchor, &fNoise, &vBuffer[0]);
        };
      }
    
  // }}}
    } // end of anonymous namespace   
  
  
    void WienerOp::apply (const ImgBase *poSrc, ImgBase **ppoDst) {
      // {{{ open
      FUNCTION_LOG("");
      if (!prepare (ppoDst, poSrc)) return;
      
      switch(poSrc->getDepth()){
  #define ARGLIST(D) poSrc->asImg<icl##D>(),(*ppoDst)->asImg<icl##D>(),getMaskSize(),getAnchor(),getROIOffset(),m_vecBuffer,m_fNoise
        case depth8u: ippiWienerCall<icl8u,ippiFilterWiener_8u_C1R>(ARGLIST(8u)); break;
        case depth16s: ippiWienerCall<icl16s,ippiFilterWiener_16s_C1R>(ARGLIST(16s)); break;
        case depth32f: ippiWienerCall<icl32f,ippiFilterWiener_32f_C1R>(ARGLIST(32f)); break;
        default: ICL_INVALID_DEPTH;
  #undef ARGLIST
      }
    }
    
    // }}}
  
  
  #endif     
  } // namespace filter
}