56MLPreconditionerState::constructMLPreconditioner(
const Teuchos::ParameterList & mainList,
57 const std::vector<Teuchos::RCP<const Teuchos::ParameterList> > & coarseningParams)
59 TEUCHOS_ASSERT(isFilled());
60 std::vector<Teuchos::ParameterList> cpls(coarseningParams.size());
61 for(std::size_t i=0;i<coarseningParams.size();i++)
62 cpls[i] = *coarseningParams[i];
64 mlPreconditioner_ = rcp(
new ML_Epetra::MultiLevelPreconditioner(&*mlOp_, mainList, &diagonalOps_[0], &cpls[0], diagonalOps_.size()));
66 return mlPreconditioner_;
81 if(not mlState.isFilled())
82 fillMLPreconditionerState(blo,mlState);
83 TEUCHOS_ASSERT(mlState.isFilled());
85 Teuchos::RCP<ML_Epetra::MultiLevelPreconditioner> mlPrecOp =
86 mlState.constructMLPreconditioner(mainParams_,coarseningParams_);
89 return Teuchos::rcp(
new MLLinearOp(mlPrecOp));
97void MLPreconditionerFactory::fillMLPreconditionerState(
const BlockedLinearOp & blo,
MLPreconditionerState & mlState)
const
99 TEUCHOS_ASSERT(not mlState.isFilled());
101 EpetraExt::CrsMatrix_SolverMap RowMatrixColMapTrans;
108 ML_Comm_Create(&mlComm);
109 mlState.setMLComm(mlComm);
111 ML_Operator * mlBlkMat = ML_Operator_Create(mlComm);
112 ML_Operator_BlkMatInit(mlBlkMat,mlComm,rowCnt,colCnt,ML_DESTROY_EVERYTHING);
114 std::vector<Epetra_RowMatrix *> aggMats;
115 for(
int r=0;r<rowCnt;r++) {
116 for(
int c=0;c<colCnt;c++) {
117 ML_Operator * tmp = 0;
118 Epetra_RowMatrix * EMat = 0;
119 std::stringstream ss;
121 ss <<
"fine_"<< r <<
"_" << c;
124 Teuchos::RCP<Epetra_CrsMatrix> crsMat
125 = Teuchos::rcp_dynamic_cast<Epetra_CrsMatrix>(Thyra::get_Epetra_Operator(*blo->getNonconstBlock(r,c)));
126 EMat = ML_Epetra::ModifyEpetraMatrixColMap(*crsMat,RowMatrixColMapTrans, ss.str().c_str());
128 aggMats.push_back(EMat);
131 tmp = ML_Operator_Create(mlComm);
132 ML_Operator_WrapEpetraMatrix(EMat, tmp);
135 ML_Operator_BlkMatInsert(mlBlkMat, tmp, r,c);
138 ML_Operator_BlkMatFinalize(mlBlkMat);
141 mlState.setMLOperator(mlBlkMat);
144 mlState.setAggregationMatrices(aggMats);
146 mlState.setIsFilled(
true);
152void MLPreconditionerFactory::initializeFromParameterList(
const Teuchos::ParameterList & settings)
157 Teko_DEBUG_SCOPE(
"MLPreconditionerFactory::initializeFromParameterList",10);
160 coarseningParams_.clear();
163 blockRowCount_ = settings.get<
int>(
"Block Row Count");
167 mainParams_ = settings.sublist(
"Smoothing Parameters");
168 mainParams_.set<Teuchos::RCP<const Teko::InverseLibrary> >(
"smoother: teko inverse library",getInverseLibrary());
172 const Teuchos::ParameterList & aggSublist = settings.sublist(
"Block Aggregation");
174 for(
int block=0;block<blockRowCount_;block++) {
176 std::stringstream ss;
177 ss <<
"Block " << block;
178 std::string sublistName = ss.str();
181 RCP<Teuchos::ParameterList> userSpec = rcp(
new Teuchos::ParameterList(aggSublist.sublist(sublistName)));
182 coarseningParams_.push_back(userSpec);