Isorropia_EpetraPartitioner.hpp

//@HEADER
/*
************************************************************************

              Isorropia: Partitioning and Load Balancing Package
                Copyright (2006) Sandia Corporation

Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
license for use of this work by or on behalf of the U.S. Government.

This library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as
published by the Free Software Foundation; either version 2.1 of the
License, or (at your option) any later version.

This library 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
Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public
License along with this library; if not, write to the Free Software
Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
USA
Questions? Contact Alan Williams (william@sandia.gov)
                or Erik Boman    (egboman@sandia.gov)

************************************************************************
*/
//@HEADER

#ifndef _Isorropia_EpetraPartitioner_hpp_
#define _Isorropia_EpetraPartitioner_hpp_

#include <Isorropia_configdefs.hpp>
#include <Teuchos_RefCountPtr.hpp>
#include <Teuchos_ParameterList.hpp>

#include <Isorropia_Partitioner.hpp>

#ifdef HAVE_EPETRA
class Epetra_Map;
class Epetra_BlockMap;
class Epetra_Import;
class Epetra_Vector;
class Epetra_MultiVector;
class Epetra_CrsGraph;
class Epetra_CrsMatrix;
class Epetra_RowMatrix;
class Epetra_LinearProblem;

namespace Isorropia {
  class CostDescriber;

namespace Epetra {

class Partitioner : public Isorropia::Partitioner {
public:
  Partitioner(Teuchos::RefCountPtr<const Epetra_CrsGraph> input_graph,
              const Teuchos::ParameterList& paramlist,
              bool compute_partitioning_now=true);

  Partitioner(Teuchos::RefCountPtr<const Epetra_CrsGraph> input_graph,
              Teuchos::RefCountPtr<const Isorropia::CostDescriber> costs,
              const Teuchos::ParameterList& paramlist,
              bool compute_partitioning_now=true);

  Partitioner(Teuchos::RefCountPtr<const Epetra_RowMatrix> input_matrix,
              const Teuchos::ParameterList& paramlist,
              bool compute_partitioning_now=true);

  Partitioner(Teuchos::RefCountPtr<const Epetra_RowMatrix> input_matrix,
              Teuchos::RefCountPtr<const Isorropia::CostDescriber> costs,
              const Teuchos::ParameterList& paramlist,
              bool compute_partitioning_now=true);

  virtual ~Partitioner();

  void setParameters(const Teuchos::ParameterList& paramlist);

  void compute_partitioning(bool force_repartitioning=false);

  bool partitioning_already_computed() const;

  int newPartitionNumber(int myElem) const;

  int numElemsInPartition(int partition) const;

  void elemsInPartition(int partition, int* elementList, int len) const;

private:

  Teuchos::RefCountPtr<const Epetra_BlockMap> input_map_;
  Teuchos::RefCountPtr<const Epetra_CrsGraph> input_graph_;
  Teuchos::RefCountPtr<const Epetra_RowMatrix> input_matrix_;
  Teuchos::RefCountPtr<const Isorropia::CostDescriber> costs_;
  Teuchos::ParameterList paramlist_;
  Teuchos::RefCountPtr<Epetra_Vector> weights_;

  std::map<int,int> exports_, imports_;
  std::vector<int> myNewElements_;

  bool partitioning_already_computed_;

};//class Partitioner

}//namespace Epetra
}//namespace Isorropia

#endif //HAVE_EPETRA

#endif