On Fri, Jan 9, 2015 at 8:47 AM, Serge Van Criekingen < [email protected]> wrote:
> Hi, > > Using petsc4py, I wish to control the distribution of matrix rows/columns > among the processors. > I found out how to do this for vectors in the archives: > https://lists.mcs.anl.gov/mailman/htdig/petsc-users/ > 2013-February/016358.html > But how can I do this with matrices? In other words, how is the > "MatSetSizes" PETSc routine interfaced in petsc4py? > > To take an example, how can I specify the local number of rows/columns > when running this minimal code in parallel: > > mpi_rank = PETSc.COMM_WORLD.getRank() > M = PETSc.Mat().create(comm=PETSc.COMM_WORLD) > M.setType('aij') > M.setSizes((8,8),bsize=1) > For example M.setSizes(((locSize, PETSC_DETERMINE), (locSize, PETSC_DETERMINE)), bsize=1) Thanks, Matt > M.setUp() > PETSc.Sys.syncPrint("rank: ",mpi_rank," M.sizes: ",M.sizes) > PETSc.Sys.syncFlush() > > Thanks a lot, > Serge > > > > -- 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
