[Trilinos-Users] AztecOO and ~OperatorData()
Williams, Alan B
william at sandia.gov
Tue Dec 13 13:01:59 MST 2005
Simone,
No, AztecOO does not destroy your P operator. When you use
SetPrecOperator, AztecOO wraps your operator in an internal struct.
AztecOO later destroys that struct, but it should not destroy your
operator.
Alan
> -----Original Message-----
> From: trilinos-users-bounces at software.sandia.gov
> [mailto:trilinos-users-bounces at software.sandia.gov] On Behalf
> Of Simone Deparis
> Sent: Tuesday, December 13, 2005 12:26 PM
> To: trilinos-users at software.sandia.gov
> Subject: [Trilinos-Users] AztecOO and ~OperatorData()
>
> Hi,
>
> I have the following issue/understanding:
>
> I have a preconditioner P of type Prec derived from Epetra_Operator:
>
> Prec P;
>
> I set this as preconditioner to Aztec00:
>
> {
> AztecOO Solver(Problem); // Problem is of class Epetra_LinearProblem
> Solver.SetPrecOperator(&P);
> M_Solver->Iterate(numiter, tol);
> }
>
> Now the instance Solver of Aztec00 has been destroyed (since
> outside the brackets).
> If I look through Aztec00.cpp, the destruction of Solver has
> involved the distruction of my preconditioner P ( through
> AztecOO::DeleteMemory() and AztecOO::~AztecOO(void) ).
>
> Do I understand well what's going on? If I am right, is this
> the behavior that we wish?
> I don't! I would prefer that I decide when to destroy my
> preconditioner.
> Is there a reason for this?
>
> Thank you
>
> Simone
>
>
> --
> _____________________________________________________________________
> Simone Deparis
> .~. Mech Eng Dept - MIT
> /V\ L I N U X 77 Mass Ave Room 3-264
> // \\ ================= Cambridge MA 02139
> /( )\ USA
> ^^-^^
> phone : +1 617 452 3285 mailto:deparis at mit.edu
> fax : +1 617 258 8559 http://www.mit.edu/~deparis
> _____________________________________________________________________
>
>
>
> _______________________________________________
> Trilinos-Users mailing list
> Trilinos-Users at software.sandia.gov
> http://software.sandia.gov/mailman/listinfo/trilinos-users
>
More information about the Trilinos-Users
mailing list