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