[petsc-users] distribute and cells mapping.

[Olivier Bonnefon]

Hello,

Thanks for your answers, I'm now able to import and distribute a mesh:
if (!rank){
        ierr = 
DMPlexCreateFromCellList(comm,dim,obNbCells,obNbVertex,3,0,obCells,2,obVertex,dm);CHKERRQ(ierr);
         for (i=0;i<obNbBound;i++){
           ierr =DMPlexSetLabelValue(*dm, "marker", 
obBoundary[i]+obNbCells, 1);CHKERRQ(ierr);
         }
}else {
          ierr = 
DMPlexCreateFromCellList(comm,dim,0,0,3,0,obCells,2,obVertex,dm);CHKERRQ(ierr);
}

ierr = DMPlexDistribute(*dm, partitioner, 0, 
&distributedMesh);CHKERRQ(ierr);
if (distributedMesh) {
      ierr = DMDestroy(dm);CHKERRQ(ierr);
      *dm  = distributedMesh;
}

Is it possible to known the resulting partition ? ie, What is the 
mapping between the initial cell number and the local cell (used in 
DMPlexComputeResidualFEM)?
I need this to write an efficient implementation of the FEM struct 
functions f0 and g0, space depending.

Regards,

Olivier B

--
Olivier Bonnefon
INRA PACA-Avignon, Unité BioSP
Tel: +33 (0)4 32 72 21 58