[Trilinos-Users] ILU IFPACK preconditioner

[lucia at mathcs.emory.edu]

Hi Mike,

> If you haven't called the Compute() method, then you could see this kind
> of error.  Might this be the problem?

actually I did call the Compute() method. Please find attached the whole
piece of code I'm trying to compile & run

"
std::string PrecType = "ILU";

int OverlapLevel = 1;

Prec = Teuchos::rcp( Factory.Create(PrecType, &*B_rcp, OverlapLevel) );

ifpackList.set("fact: drop tolerance", drop_tol_IFPACK);

ifpackList.set("fact: level-of-fill", 1);

Prec->SetParameters(ifpackList);

Prec->Initialize();

Prec->Compute();

Ifpack_ILU * Prec_ILU = dynamic_cast<Ifpack_ILU *>(&*Prec);

Epetra_CrsMatrix L_factor= Prec_ILU->L();

//^^^ here I get the segfault
"

many thanks again for your help

Lucia