|| ||           || Doxygen's  

---

HxExportExtraIdentMaskStdev.h

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

#ifndef HxExportExtraIdentMaskStdev_h
#define HxExportExtraIdentMaskStdev_h

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


template<class ResultT, class ImValT, class ExtraValT>
class HxExportExtraIdentMaskStdev
{
public:
    typedef HxTagTransInVar     TransVarianceCategory;

    typedef HxTag1Phase         PhaseCategory;


                        HxExportExtraIdentMaskStdev(HxTagList& tags);

                        ~HxExportExtraIdentMaskStdev();

    void                doIt(const ImValT& imV, const ExtraValT& extraV);

    static HxString     className();

private:
    HxTagList&          _tags;
    ExtraValT           _maskVal;
    ResultT             _sum;
    ResultT             _sumSqr;
    HxScalarInt         _num;
    bool                _exportVariance;
};

template<class ResultT, class ImValT, class ExtraValT>

HxExportExtraIdentMaskStdev<ResultT,ImValT,ExtraValT>::HxExportExtraIdentMaskStdev(HxTagList& tags) : _tags(tags)
{
    _maskVal = HxGetTag<int>(tags, "maskVal");
    _exportVariance = HxGetTag(tags, "exportVariance",false);
    _sum = HxScalarInt(0);
    _sumSqr = HxScalarInt(0);
    _num = 0;
}

template<class ResultT, class ImValT, class ExtraValT>
HxExportExtraIdentMaskStdev<ResultT,ImValT,ExtraValT>::~HxExportExtraIdentMaskStdev()
{
    ResultT num = _num;
    ResultT num1 = _num - 1;
    ResultT result = (_num == 1) ? ResultT(HxScalarInt(0))
                                 : (_sumSqr - ((_sum * _sum) / num)) / (num1);
    if(!_exportVariance)
        result = result.sqrt();
    HxAddTag(_tags, "result", HxValue(result));
}

template<class ResultT, class ImValT, class ExtraValT>
inline void

HxExportExtraIdentMaskStdev<ResultT,ImValT,ExtraValT>    ::doIt(const ImValT& imV, const ExtraValT& extraV)
{
    if (extraV == _maskVal) {
        _sum += imV;
        _sumSqr += imV * imV;
        _num += 1;
    }
}

template<class ResultT, class ImValT, class ExtraValT>
inline HxString
HxExportExtraIdentMaskStdev<ResultT,ImValT,ExtraValT>::className()
{
    return HxString("identMaskStdev");
}


#endif

---