Re: [petsc-users] MatConvert changes distribution of local rows

[Marius Buerkle]

Matrix types is from MATMPIDENSE to MATSCALAPACK, but I think it happens also for other matrix types IIRC.

 

 

Gesendet: Freitag, 13. Januar 2023 um 16:58 Uhr
Von: "Pierre Jolivet" <pie...@joliv.et>
An: "Marius Buerkle" <mbuer...@web.de>
Cc: petsc-users@mcs.anl.gov
Betreff: Re: [petsc-users] MatConvert changes distribution of local rows

 

On 13 Jan 2023, at 8:49 AM, Marius Buerkle <mbuer...@web.de> wrote:

 

Hi,

 

I have a matrix A for which I defined the number of local rows per process manually using MatSetSizes. When I use MatConvert to change the matrix type it changes the number of local rows (to what one would get if MatSetSize is called with PETSC_DECIDE for number of local rows), which causes problems when doing MatVec producs and stuff like that. Is there any way to preserve the the number of local rows when using MatConvert?

 

This is most likely a bug, it’s not handled properly in some MatConvert() implementations.

Could you please share either the matrix types or a minimal working example?

 

Thanks,

Pierre

 

 

Best,

Marius