On 12/07/2017 03:12 AM, Lucas Campos wrote:
Currently I am using a direct LU solver via PETSc/MUMPS to solve my matrix.
However, I have noticed that I spend a lot of time in allocation, at every
step. Is it possible (or useful) to preallocate the internal structures
necessary to solve the matrix? According to [1], it is possible to do it if I
use UMFPack, but it seems I would need to change a bit more code to still work
with MPI, so would be simpler to do it while using PETSc/MUMPS.
Lucas -- I don't think I entirely understand the question. Where is the time
allocating memory spent? When you build the matrix, or when MUMPS works on it
to compute the LU decomposition?
You have no control over what MUMPS does -- it's a black box from our
perspective. Or do you know whether PETSc allows setting parameters that are
then passed to MUMPS?
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.