[Trilinos-Users] ILU IFPACK preconditioner

[Lucia Mirabella]

Hi,

thanks to a Simone Deparis' suggestion, I discovered that the result of 
the dynamic cast is a null pointer... so maybe that's the problem, but I 
still don't know why it happens...

thanks for your help,

Lucia



lucia at mathcs.emory.edu wrote:
> 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
>
>
>
>