column map in Mat MPIAdj

[Lisandro Dalcin]

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