00001 // $Id: Anasazi_LOCA_Sort.C,v 1.6 2008/06/18 21:51:25 cgbaker Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-loca/src/Anasazi_LOCA_Sort.C,v $ 00003 00004 //@HEADER 00005 // ************************************************************************ 00006 // 00007 // NOX: An Object-Oriented Nonlinear Solver Package 00008 // Copyright (2002) Sandia Corporation 00009 // 00010 // LOCA: Library of Continuation Algorithms Package 00011 // Copyright (2005) Sandia Corporation 00012 // 00013 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive 00014 // license for use of this work by or on behalf of the U.S. Government. 00015 // 00016 // This library is free software; you can redistribute it and/or modify 00017 // it under the terms of the GNU Lesser General Public License as 00018 // published by the Free Software Foundation; either version 2.1 of the 00019 // License, or (at your option) any later version. 00020 // 00021 // This library is distributed in the hope that it will be useful, but 00022 // WITHOUT ANY WARRANTY; without even the implied warranty of 00023 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00024 // Lesser General Public License for more details. 00025 // 00026 // You should have received a copy of the GNU Lesser General Public 00027 // License along with this library; if not, write to the Free Software 00028 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00029 // USA 00030 // 00031 // Questions? Contact Roger Pawlowski (rppawlo@sandia.gov) or 00032 // Eric Phipps (etphipp@sandia.gov), Sandia National Laboratories. 00033 // ************************************************************************ 00034 // CVS Information 00035 // $Source: /space/CVS/Trilinos/packages/nox/src-loca/src/Anasazi_LOCA_Sort.C,v $ 00036 // $Author: cgbaker $ 00037 // $Date: 2008/06/18 21:51:25 $ 00038 // $Revision: 1.6 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #include "Anasazi_LOCA_Sort.H" 00043 #include "LOCA_EigenvalueSort_Strategies.H" 00044 #include "LOCA_GlobalData.H" 00045 #include "LOCA_ErrorCheck.H" 00046 00047 Anasazi::LOCASort::LOCASort( 00048 const Teuchos::RCP<LOCA::GlobalData>& global_data, 00049 const Teuchos::RCP<LOCA::EigenvalueSort::AbstractStrategy>& strategy_) 00050 : globalData(global_data), 00051 strategy(strategy_) 00052 { 00053 } 00054 00055 Anasazi::LOCASort::~LOCASort() 00056 { 00057 } 00058 00059 void 00060 Anasazi::LOCASort::sort(std::vector<double>& evals, 00061 Teuchos::RCP<std::vector<int> > perm, 00062 int n) const 00063 { 00064 if (n == -1) { 00065 n = evals.size(); 00066 } 00067 NOX::Abstract::Group::ReturnType res = strategy->sort(n, &evals[0], perm.get()); 00068 globalData->locaErrorCheck->checkReturnType(res, "Anasazi::LOCASort::sort()"); 00069 } 00070 00071 void 00072 Anasazi::LOCASort::sort(std::vector<double>& r_evals, 00073 std::vector<double>& i_evals, 00074 Teuchos::RCP<std::vector<int> > perm, 00075 int n) const 00076 { 00077 if (n == -1) { 00078 n = r_evals.size(); 00079 } 00080 NOX::Abstract::Group::ReturnType res = 00081 strategy->sort(n, &r_evals[0], &i_evals[0], perm.get()); 00082 globalData->locaErrorCheck->checkReturnType(res, 00083 "Anasazi::LOCASort::sort()"); 00084 }