hi, dear all,
i am using trilino's direct solver interface to solve a time-dependent
problem. I realized that
when the number of iterations increases, the memory usage also increases.
This increase
does not arise when turning the solver off, so just assembling the matrices.
At the beginning, i used the following subroutine at every time step,
{
deallog.push("DirectKLU");
TrilinosWrappers::SolverDirect::AdditionalData data;
data.solver_type = "Amesos_Klu";
SolverControl solver_control (1000, 1e-10);
TrilinosWrappers::SolverDirect solver(solver_control, data);
solver.solve (matrix, solution, rhs);
thick_constraints.distribute (thick_solution);
}
then i realized that the memory might be allocated every time when
defining TrilinosWrappers::SolverDirect solver(.....), which is not
released.
Then I predefine the solver in a namespace, which is the used by calling
solver.solve (matrix, solution, rhs) at each step. Still, the memory
leakage
occurs.
I also tried to use "Amesos_Mumps", the problem still exists. I tried to
see the documentation of Amesos, i realized that when amesos solver
finished solving, there is a procedure 'delete solver' after. May I ask
whether
such a member function exists in dealII wrappers? Thanks in advance,
best,
lailai
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
