Actually I am using fortran. But I also tried with C and yes I am casting A to PetscObject.
To confirm, are you casting A to PetscObject? e.g.
MPI_Comm comm;
/* ... */
PetscObjectGetComm((PetscObject)A,&comm);
Am Mi., 22. Apr. 2020 um 07:11 Uhr schrieb Marius Buerkle <mbuer...@web.de>:
Hi,
What is PetscObjectGetComm expected to return? I thought it would give the MPI communicator the object lives on. So if I create A matrix on PETSC_COMM_WORLD a call of PetscObjectGetComm for A it would return PETSC_COMM_WORLD? But it seems to return something else, and while most of the nodes return a similar communicator some are giving a different one. That said, is there a way to get the MPI communicator a matrix lives on?
Best,
Marius
