GeomObj.h

// -*-c++-*-
//----------------------------------------------------------------------------
//  xc utils library; general purpose classes and functions.
//
//  Copyright (C)  Luis C. Pérez Tato
//
//  XC utils 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 3 of the License, or 
//  (at your option) any later version.
//
//  This software 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, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------------
//GeomObj.h

#ifndef GEOMOBJ_H
#define GEOMOBJ_H

#include <iostream>
#include "ProtoGeom.h"
#include "lists/PolyPos.h"

class Pos2d;
class Pos3d;
class Vector2d;
class Trf2d;
class Plotter;


//
class GeomObj: public ProtoGeom
  {
  public:
    typedef PolyPos<Pos3d> list_Pos3d;
    typedef PolyPos<Pos2d> list_Pos2d;

    GeomObj(void);
    virtual ~GeomObj(void) {}
    virtual GeomObj *getCopy(void) const=0;
    virtual unsigned short int Dimension(void) const= 0;
    virtual GEOM_FT getLength(void) const= 0;
    virtual GEOM_FT getArea(void) const= 0;
    virtual GEOM_FT getVolume(void) const= 0;
    inline virtual bool hasCenterOfMass(void) const
      { return false;}
    GEOM_FT getCenterOfMassArea(void) const;

    virtual void Print(std::ostream &os) const= 0;
    virtual void Plot(Plotter &) const {}
  };

std::ostream &operator<<(std::ostream &, const GeomObj &);


template <class input_iterator>
GEOM_FT area(input_iterator begin,input_iterator end)
  {
    GEOM_FT retval= 0.0;
    for(input_iterator i=begin;i!=end;i++)
      retval+= i->getArea();
    return retval;
  }

template <class input_iterator>
GEOM_FT volume(input_iterator begin,input_iterator end)
  {
    GEOM_FT retval= 0.0;
    for(input_iterator i=begin;i!=end;i++)
      retval+= i->getVolume();
    return retval;
  }

template <class input_iterator>
GEOM_FT length(input_iterator begin,input_iterator end)
  {
    GEOM_FT retval= 0.0;
    for(input_iterator i=begin;i!=end;i++)
      retval+= i->getLength();
    return retval;
  }

template <class input_iterator>
void print(std::ostream &os,input_iterator begin,input_iterator end)
  {
    if(begin!=end)
      {
        input_iterator i= begin;
        os << *i; i++;
        for(;i!=end;i++)
          os << ", " << *i;
      }
  }
#endif