/* Copyright (c) 1997-2005
   Ewgenij Gawrilow, Michael Joswig (Technische Universitaet Berlin, Germany)
   http://www.math.tu-berlin.de/polymake,  mailto:polymake@math.tu-berlin.de

   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, or (at your option) any
   later version: http://www.gnu.org/licenses/gpl.txt.

   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.
*/

#ifndef _POLYMAKE_CONNECTED_SUM_H
#define _POLYMAKE_CONNECTED_SUM_H "$Project: polymake $$Id: connected_sum.h,v 1.5 2005/09/22 17:12:10 gawrilow Exp $"

#if defined(__GNUC__) && __GNUC__==3 && __GNUC_MINOR__<3 && !defined(__APPLE__) && !defined(__COMO__)
#pragma interface
#endif

#include <Array.h>
#include <Set.h>

#include <string>
#include <list>
#include <ext/hash_map>


namespace polymake { namespace topaz {

   typedef std_ext::hash_map<int,int> map;

   /// Computes the connected sum of two complexes and adjusts the labels.
   /// The client writes the new lables into L1.
   /// If the labels are empty, the client just ignores them.
   /// The automorphism P determines to which of the vertices of f1 the
   /// vertices of f2 get glued.
   /// If the map is empty, the trivial map is used.
   template <typename Complex_1, typename Complex_2>
   std::list< Set<int> > connected_sum(const Complex_1& C1,
				       const Complex_2& C2,
				       const int f1, const int f2,
				       Array<std::string>& L1,
				       const Array<std::string>& L2,
				       map& P);
   
   /// Computes the connected sum of two complexes,
   /// using the first facet of each complex and the
   /// trivial map.
   template <typename Complex_1, typename Complex_2>
   std::list< Set<int> > connected_sum(const Complex_1& C1,
				       const Complex_2& C2)
   {
      map P;
      Array<std::string> L;
      
      return connected_sum(C1,C2,0,0,L,L,P);
   }
   
   
} }

#include <connected_sum.tcc>

#endif // _POLYMAKE_CONNECTED_SUM_H

// Local Variables:
// mode:C++
// c-basic-offset:3
// End:
