shells/utilities/MinMaxSensor.h

Go to the documentation of this file.

// -*- C++ -*-
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// 
//                                   Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2004 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//

#ifndef MINMAXSENSOR_H
#define MINMAXSENSOR_H

#include "../driverCC/MShell.h"
#include "../driverCC/SVertexFunctors.h"

#include <functional>
#include <iterator>
#include <ostream>
#include <cassert>


namespace shells {
    class MShell;
}
 
namespace utilities {
    template<typename EXTR>
    class MinMaxSensor;
}


namespace utilities {

template<typename EXTR>
class MinMaxSensor {
public:
    MinMaxSensor(std::ostream& os, shells::MShell * const mShell): 
        _os(os), _mShell(mShell)
        {
            _os << std::scientific;
        }
    ~MinMaxSensor(){_os.flush();}
    
    void printData(double timeStamp);

private:
    MinMaxSensor(const MinMaxSensor &);
    const MinMaxSensor & operator=(const MinMaxSensor &);
    
private:
    typedef std::vector<typename EXTR::DataType>     _DataCont;

    _DataCont                      _variable;
    std::ostream&                  _os;
    shells::MShell     * const     _mShell;
};


// implementation    
template<typename EXTR>
void MinMaxSensor<EXTR> ::printData(double timeStamp) 
{
    const size_t numVar = EXTR::numVar;
    size_t numVertices = _mShell->numberOfVertices();
    
    _variable.clear();
    _variable.reserve(numVertices*numVar);
    shells::SVertexCollector<std::back_insert_iterator<_DataCont>, EXTR> getData;
    _mShell->iterateOverVertices(std::bind2nd(getData, 
                                              std::back_inserter(_variable)));  

    assert((numVertices*numVar)==_variable.size());
    
    // find the minum and maximum
    typename EXTR::DataType min[numVar];
    typename EXTR::DataType max[numVar];
    
    for (size_t j=0; j<numVar; ++j) {
        min[j] = _variable[j];
        max[j] = min[j];
    }
    
    for (size_t i=0; i<numVertices; ++i ) {
        for (size_t j=0; j<numVar; ++j) {
            min[j] = std::min(min[j],_variable[numVar*i+j]);
                max[j] = std::max(max[j],_variable[numVar*i+j]);
        }
    }
        
    _os << timeStamp << " ";
    for (size_t i=0; i<numVar; ++i) {
        _os << min[i] << " " << max[i] << " ";
    }
    _os << std::endl;
}

} // namespace utilities

#endif

---