h o m e

d o c u m e n t a t i o n

c l a s s h i e r a r c h y

---

Interface1D.h

Go to the documentation of this file.

//
//  Filename         : Interface1D.h
//  Author(s)        : Emmanuel Turquin
//  Purpose          : Interface to 1D elts
//  Date of creation : 01/07/2003
//


//
//  Copyright (C) : Please refer to the COPYRIGHT file distributed 
//   with this source distribution. 
//
//  This program is free software; you can redistribute it and/or
//  modify it under the terms of the GNU General Public License
//  as published by the Free Software Foundation; either version 2
//  of the License, or (at your option) any later version.
//
//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software
//  Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
//

#ifndef  INTERFACE1D_H
# define INTERFACE1D_H

# include <string>
# include <iostream>
# include <float.h>
# include "../system/Id.h"
# include "../system/Precision.h"
# include "../winged_edge/Nature.h"
# include "Functions0D.h"

using namespace std;
// Integration method
typedef enum {
  MEAN,
  MIN,
  MAX,
  FIRST,
  LAST
} IntegrationType;

template <class T>
T integrate(UnaryFunction0D<T>& fun,
            Interface0DIterator it,
            Interface0DIterator it_end,
            IntegrationType integration_type = MEAN) {
  T res;
  T res_tmp;
  unsigned size;
  switch (integration_type) {
  case MIN:
    res = fun(it);++it;
    for (; !it.isEnd(); ++it) {
      res_tmp = fun(it);
      if (res_tmp < res)
        res = res_tmp;
    }
    break;
  case MAX:
    res = fun(it);++it;
    for (; !it.isEnd(); ++it) {
      res_tmp = fun(it);
      if (res_tmp > res)
        res = res_tmp;
    }
    break;
  case FIRST:
    res = fun(it);
    break;
  case LAST:
    res = fun(--it_end);
    break;
  case MEAN:
  default:
    res = fun(it);++it;
    for (size = 1; !it.isEnd(); ++it, ++size)
      res += fun(it);
    res /= (size ? size : 1);
    break;
  }
  return res;
}

//
// Interface1D
//

class Interface1D
{
public:

  Interface1D() {_timeStamp=0;}

  virtual string getExactTypeName() const {
    return "Interface1D";
  }

  // Iterator access

  virtual Interface0DIterator verticesBegin() = 0;
  virtual Interface0DIterator verticesEnd() = 0;
  virtual Interface0DIterator pointsBegin(float t=0.f) = 0;
  virtual Interface0DIterator pointsEnd(float t=0.f) = 0;

  // Data access methods

  virtual real getLength2D() const {
    cerr << "Warning: method getLength2D() not implemented" << endl;
    return 0;
  }

  virtual Id getId() const {
    cerr << "Warning: method getId() not implemented" << endl;
    return Id(0, 0);
  }

  // FIXME: ce truc n'a rien a faire la...(c une requete complexe qui doit etre ds les Function1D)
  virtual Nature::EdgeNature getNature() const {
    cerr << "Warning: method getNature() not implemented" << endl;
    return Nature::NO_FEATURE;
  }
  
  virtual unsigned getTimeStamp() const {
    return _timeStamp;
  }
  
  inline void setTimeStamp(unsigned iTimeStamp){
    _timeStamp = iTimeStamp;
  }

protected:
  unsigned _timeStamp;
};

#endif // INTERFACE1D_H

---

Last modified 3 Mar 2008

© Stephane Grabli