[petsc-users] Matload with given parallel layout

Matthew Knepley knepley at gmail.com
Wed Oct 25 04:40:12 CDT 2017

On Wed, Oct 25, 2017 at 4:32 AM, Karin&NiKo <niko.karin at gmail.com> wrote:

> Dear PETSc team,
> I have a code that creates a parallel matrix based on domain
> decomposition. I serialize this matrix with MatView.
> Then, I would like to relaod it with MatLoad but not leaving PETSc decide
> the parallel layout but rather use the original distribution of the degrees
> of freedom.
> Does MatLoad help to do that? Shall  I use the IS of the local/global
> mapping ?

You can prescribe the layout of the Mat that is passed to MatLoad().
However, if you are asking that the original
layout be stored in the file, we do not do that. It could be made to work
by writing the layout to the .info file and
then having the option call MatSetSizes() on load. Should not be much code.



> I look forward to reading you,
> Nicolas

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171025/593803f2/attachment.html>

More information about the petsc-users mailing list