Sorry, my comment about PETSC_COMM_SELF was a nonsense :-(
Let me write the rule in C-Python-pseudocode
rank = mat->comm->rank
size = mat->comm->size
cmap->N = # total_columns
for p in range(size):
if p == rank:
mat->cmap->n = cmap->N
else:
mat->cmap->n = 0
# and finally
PetscMapSetUp(mat->cmap)
That would be fine?
On Wed, Dec 3, 2008 at 1:19 PM, Lisandro Dalcin <dalcinl at gmail.com> wrote:
> Can someone help me to review this routine?
>
> MatCreate_MPIAdj(Mat B)
>
> I understand that MATMPIADJ is rather special, but the call above is never
> made
>
> PetscMapSetUp(B->cmap)
>
> then It seems that some stuff (ranges) in the column map will be
> uninitialized.
>
> What values should have the entries in "cmap" for mpiadj matrices?
> Should the "comm" in "cmap" be MPI_COMM_SELF, then we can call
> PetscMapSetUp?
>
>
>
>
> --
> Lisandro Dalc?n
> ---------------
> Centro Internacional de M?todos Computacionales en Ingenier?a (CIMEC)
> Instituto de Desarrollo Tecnol?gico para la Industria Qu?mica (INTEC)
> Consejo Nacional de Investigaciones Cient?ficas y T?cnicas (CONICET)
> PTLC - G?emes 3450, (3000) Santa Fe, Argentina
> Tel/Fax: +54-(0)342-451.1594
>
--
Lisandro Dalc?n
---------------
Centro Internacional de M?todos Computacionales en Ingenier?a (CIMEC)
Instituto de Desarrollo Tecnol?gico para la Industria Qu?mica (INTEC)
Consejo Nacional de Investigaciones Cient?ficas y T?cnicas (CONICET)
PTLC - G?emes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594
