Re: [petsc-users] MatCompositeMerge + MatCreateRedundantMatrix

ok, so it seems there is no straight forward way to transfer data between PETSc matrices on different subcomms. Probably doing it by "hand" extracting the matricies on the subcomms create a MPI_INTERCOMM transfering the data to PETSC_COMM_WORLD and assembling them in a new PETSc matrix would be pos

Re: [petsc-users] MatCompositeMerge + MatCreateRedundantMatrix

I see. This would work if the matrices are on different subcommumicators ? Is it possible to add this functionality ?   marius     You basically need the inverse of MatCreateSubmatrices(). I do not think we have that right now, but it could probably be done without too much trouble by looking

Re: [petsc-users] Identifying matching points in differently distributed DMPlexs

On Tue, Feb 19, 2019 at 11:00 AM Lawrence Mitchell via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear petsc-users, > > I have two (different) distributions of the (topologically) same DMPlex > object (DM_a and DM_b). > > I would like to identify the map from points(DM_a) to points(DM_b) such

[petsc-users] Identifying matching points in differently distributed DMPlexs

Dear petsc-users, I have two (different) distributions of the (topologically) same DMPlex object (DM_a and DM_b). I would like to identify the map from points(DM_a) to points(DM_b) such that I can transfer fields between the two. Does such a facility exist? Lawrence

Re: [petsc-users] MatCompositeMerge + MatCreateRedundantMatrix

You basically need the inverse of MatCreateSubmatrices(). I do not think we have that right now, but it could probably be done without too much trouble by looking at that code. Thanks, Matt On Tue, Feb 19, 2019 at 6:15 AM Marius Buerkle via petsc-users < petsc-users@mcs.anl.gov> wrote: >

Re: [petsc-users] PCMG: level dependent smoothing steps

Hi Patrick, It worked, thank you! Cheers from Lugano, Pietro Il giorno 19-feb-2019, alle ore 16:19, Patrick Sanan mailto:patrick.sa...@gmail.com>> ha scritto: Will it work for you to use https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMGGetSmoother.html to pull out the KSP

Re: [petsc-users] PCMG: level dependent smoothing steps

Will it work for you to use https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMGGetSmoother.html to pull out the KSP for the level you're interested in (Note link to PCMGGetSmootherDown() if you want to control up and down smoothers separately) and then use https://www.mcs.anl.go

[petsc-users] PCMG: level dependent smoothing steps

Dear PETSc team, Is there a way to set the number of smoothing steps done inside a specific multigrid level l ? Something like: PCMGSetNumberSmoothDown(PC

[petsc-users] MatCompositeMerge + MatCreateRedundantMatrix

Hi !   Is there some way to combine MatCompositeMerge with MatCreateRedundantMatrix? I basically want to create copies of a matrix from PETSC_COMM_WORLD to subcommunicators, do some work on each subcommunicator and than gather the results back to PETSC_COMM_WORLD, namely  I want to sum the  the i