Teko_mlutils.hpp

// @HEADER
// *****************************************************************************
//      Teko: A package for block and physics based preconditioning
//
// Copyright 2010 NTESS and the Teko contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#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);

}  // namespace mlutils
}  // namespace Teko

#endif