Epetra_OskiPermutation.cpp

Go to the documentation of this file.
00001 
00002 //@HEADER
00003 // ************************************************************************
00004 // 
00005 //               Epetra: Linear Algebra Services Package 
00006 //                 Copyright (2001) Sandia Corporation
00007 // 
00008 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00009 // license for use of this work by or on behalf of the U.S. Government.
00010 // 
00011 // This library is free software; you can redistribute it and/or modify
00012 // it under the terms of the GNU Lesser General Public License as
00013 // published by the Free Software Foundation; either version 2.1 of the
00014 // License, or (at your option) any later version.
00015 //  
00016 // This library is distributed in the hope that it will be useful, but
00017 // WITHOUT ANY WARRANTY; without even the implied warranty of
00018 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00019 // Lesser General Public License for more details.
00020 //  
00021 // You should have received a copy of the GNU Lesser General Public
00022 // License along with this library; if not, write to the Free Software
00023 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00024 // USA
00025 // Questions? Contact Michael A. Heroux (maherou@sandia.gov) 
00026 // 
00027 // ************************************************************************
00028 //@HEADER
00029 
00030 #include "Epetra_ConfigDefs.h"
00031 
00032 #ifdef HAVE_OSKI
00033 #ifdef HAVE_EPETRA_TEUCHOS
00034 #include "Epetra_OskiPermutation.h"
00035 
00036 //=============================================================================
00037 
00038 
00039 Epetra_OskiPermutation::Epetra_OskiPermutation() 
00040   : Permutation_(NULL) {}
00041 
00042 Epetra_OskiPermutation::Epetra_OskiPermutation(const Epetra_OskiPermutation& Source) 
00043   : Permutation_(Source.Permutation_) {}
00044 
00045 Epetra_OskiPermutation::Epetra_OskiPermutation(bool RowPerm, const Epetra_OskiMatrix& Source) {
00046   if(RowPerm)
00047     *this = Source.ViewRowPermutation();
00048   else
00049     *this = Source.ViewColumnPermutation();
00050 }
00051 
00052 Epetra_OskiPermutation::~Epetra_OskiPermutation() {
00053   if(Permutation_ != NULL)
00054     delete this;
00055 }
00056 
00057 void Epetra_OskiPermutation::ReplacePermutation (const oski_perm_t& InPerm) {
00058   Permutation_ = &InPerm;
00059 }
00060 
00061 int Epetra_OskiPermutation::PermuteVector(const bool TransA, Epetra_OskiMultiVector& Vector) const {
00062   int ReturnVal;
00063   if(TransA)
00064     ReturnVal = oski_PermuteVecView(*Permutation_, OP_TRANS, Vector.Oski_View());
00065   else
00066     ReturnVal = oski_PermuteVecView(*Permutation_, OP_NORMAL, Vector.Oski_View());
00067   if(ReturnVal)
00068     std::cerr << "Error in PermuteVector.\n";
00069   return ReturnVal;
00070 }
00071 
00072 #endif
00073 #endif

Generated on Wed May 12 21:41:05 2010 for Epetra Package Browser (Single Doxygen Collection) by  doxygen 1.4.7