|| ||           || Doxygen's  

---

HxNgbBernsen.h

/*
 *  Copyright (c) 2002, University of Amsterdam, The Netherlands.
 *  All rights reserved.
 *
 *  Author(s):
 *  Leon Todoran            (todoran@science.uva.nl)
 */

#ifndef HxNgbBernsen_h
#define HxNgbBernsen_h

#include "HxClassName.h"
#include "HxSizes.h"
#include "HxCategories.h"
#include "HxTagList.h"

#include <algorithm>

template <class ArgType, class ResType>
class HxNgbBernsen 
{
public:

    typedef HxTagLoop        IteratorCategory;

    typedef HxTag1Phase      PhaseCategory;


//                          HxNgbBernsen();
                            HxNgbBernsen(HxTagList& tags);
                            ~HxNgbBernsen();

    HxSizes                 size();
    void                    init(int, int, ArgType);

    void                    next(int x, int y, ArgType value);
    ResType                  result() const;

    static HxString         className();

private:
    int                     _pctIdx;
    int                     _size;  //this 
    int                     _conn;  //connectivity 4 or 8
    int                     _windowSz;
    int                     _uniformTh;
    bool                    _uniformLow;

    ArgType*                 _values;
    int                     _i;
};

template <class ArgType, class ResType>
ostream&
operator<<(ostream& os, const HxNgbBernsen<ArgType, ResType>& pct)
{
    return pct.put(os);
}

template <class ArgType, class ResType>
inline HxString
HxNgbBernsen<ArgType, ResType>::className()
{
    static HxString s("bernsen");
    return s;
}

template <class ArgType, class ResType>
inline void
HxNgbBernsen<ArgType, ResType>::init(int, int, ArgType)
{
    _i = 0;
}

template <class ArgType, class ResType>
inline void
HxNgbBernsen<ArgType, ResType>::next(int, int, ArgType value)
{
    _values[_i++] = value;
}

template <class ArgType, class ResType>
inline HxSizes
HxNgbBernsen<ArgType, ResType>::size()
{
    return HxSizes(_size, _size, 1);
}

template <class ArgType, class ResType>
inline ResType
HxNgbBernsen<ArgType, ResType>::result() const
{
    std::sort(&_values[0], &_values[_size*_size]);
    ArgType minval=_values[0];
    ArgType maxval=_values[_size*_size-1];
    ArgType t = 0.5*(minval*maxval);
    ArgType d = maxval-minval;
    if(d>_uniformTh)
    {
        return _values[_pctIdx]<t?0:1;
    }
    else
    {
        return _uniformLow?0:1;
    }

    return _values[_pctIdx];
}

#ifdef INC_TEMPLATE_SRC
#include "HxNgbBernsen.c"
#endif


#endif

---