Teko_mlutils.hpp

#ifndef __Teko_mlutils_hpp__
#define __Teko_mlutils_hpp__
 
#include "ml_operator.h"
#include "Epetra_Map.h"
#include "Epetra_CrsMatrix.h"
 
#include "Teuchos_RCP.hpp"
 
#include "Teko_Utilities.hpp"
 
namespace Teko {
 
class InverseLibrary;
 
namespace mlutils {

Teuchos::RCP<Epetra_Map> buildRowMap(ML_Operator * mlOp);

Teuchos::RCP<Epetra_CrsMatrix> convertToCrsMatrix(ML_Operator * mlOp,
                                                  const Teuchos::RCP<Epetra_Map> & rowMap=Teuchos::null);
 
Teko::LinearOp buildTekoBlockOp(ML_Operator * mlOp,int level);
 

struct SmootherData {
   Teuchos::RCP<Epetra_Operator> Amat;
   Teuchos::RCP<Epetra_Operator> smootherOperator;
};

int smoother(ML_Smoother *mydata, int leng1, double x[], int leng2, 
             double rhs[]);
 
extern "C" 
int ML_Gen_Smoother_Teko(ML *ml, int level, int pre_or_post, int ntimes, const Teuchos::RCP<const Teuchos::ParameterList> & tekoPL,
                         const Teuchos::RCP<const Teko::InverseLibrary> & invLib, const std::string & inverse, bool isBlocked);
 
}
}
 
#endif