Amesos_Support.cpp

Go to the documentation of this file.

// @HEADER
// ***********************************************************************
//
//                Amesos: Direct Sparse Solver Package
//                 Copyright (2004) 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., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301
// USA
// Questions? Contact Michael A. Heroux (maherou@sandia.gov)
//
// ***********************************************************************
// @HEADER

#include "Amesos_Klu.h"
#include "Amesos_Support.h"

Amesos_StandardIndex::Amesos_StandardIndex( const Epetra_Map& OriginalMap ) 
{

#ifdef HAVE_AMESOS_EPETRAEXT
  int NumMyElements = OriginalMap.NumMyElements();
#ifndef EPETRA_NO_32BIT_GLOBAL_INDICES
  if(OriginalMap.GlobalIndicesInt()) {
    int NumGlobalElements = OriginalMap.NumGlobalElements();
    StdIndexMap_ = rcp( new Epetra_Map( NumGlobalElements, NumMyElements, 0, OriginalMap.Comm() ) );
  }
  else
#endif
#ifndef EPETRA_NO_64BIT_GLOBAL_INDICES
  if(OriginalMap.GlobalIndicesLongLong()) {
    long long NumGlobalElements = OriginalMap.NumGlobalElements64();
    StdIndexMap_ = rcp( new Epetra_Map( NumGlobalElements, NumMyElements, 0LL, OriginalMap.Comm() ) );
  }
  else
#endif
    throw "Amesos_StandardIndex::Amesos_StandardIndex: ERROR, GlobalIndices type unknown.";

  MatTrans_ = rcp( new EpetraExt::CrsMatrix_Reindex( *StdIndexMap_ ) );
  VecTrans_ = rcp( new EpetraExt::MultiVector_Reindex( *StdIndexMap_ ) );
#endif

}

#ifdef HAVE_AMESOS_EPETRAEXT
Epetra_MultiVector* Amesos_StandardIndex::StandardizeIndex( Epetra_MultiVector* OriginalMultiVector ) {

  return ( &((*VecTrans_)( *OriginalMultiVector )) );
  
}

Teuchos::RCP<Epetra_MultiVector> Amesos_StandardIndex::StandardizeIndex( Epetra_MultiVector & OriginalMultiVector ) {

  return VecTrans_->transform(OriginalMultiVector);
}
  

Epetra_CrsMatrix*  Amesos_StandardIndex::StandardizeIndex( Epetra_CrsMatrix* OriginalCrsMatrix ) {

    return &((*MatTrans_)( *OriginalCrsMatrix ));
}
#endif