[Trilinos-Users] segmentation fault in Ifpack or ML

[Jonas Thies]

Hi (again),

I'm encountering an issue I didn't have about a year ago (haven't been
running this solver in the mean time):

I use ML with Ifpack_AdditiveSchwarz as smoother (just 2 levels).
The first setup + solve go fine, but then before building the next
solver (in a Newton process) i get a segmentation fault.
I have narrowed it down to the function ML_Ifpack_Destroy in file
ml_ifpack_wrap.cpp,
line 316:

  const Epetra_MpiComm *comm = dynamic_cast<const
Epetra_MpiComm*>(&(Prec->Matrix().Comm()));

It seems that the call Prec->Matrix() causes the segfault,
though Prec is not NULL (Prec->Label() works)

I don't think it has to do with the machine/compiler I use, but is
an IBM P6 with IBM's XL compilers, version 10.1. It happens
also with a single process, and independent of the Ifpack preconditioner
used inside
the AdditiveSchwarz. It does not happen with the Aztec ILU as smoother.

I would be very grateful if someone can figure this out,
Jonas