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.

Reply via email to