[petsc-users] Matrix Decomposition

Matthew Knepley knepley at gmail.com
Wed May 15 06:51:23 CDT 2019

On Wed, May 15, 2019 at 7:35 AM Eda Oktay via petsc-users <
petsc-users at mcs.anl.gov> wrote:

> Hello,
> I am trying to divide a matrix into unequal sized parts into different
> processors (for example I want to divide 10*10 matrix into 4*4 and 6*6
> submatrix in two processors). When my program reads a matrix from file, it
> automatically divides it into equal parts and then I can't change local
> sizes.
> How can I decompose a matrix that is read from a file?

MatLoad() takes a matrix argument. I believe you can use MatSetSizes()
before loading to get the distribution you want.


> Thanks,
> Eda

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.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190515/1193f82b/attachment.html>

More information about the petsc-users mailing list