Commits

Stefano Zampini  committed 4843afd

PCIS: consider myself as a neighbor when classifying interface nodes

Local periodic nodes should be considered has interface nodes.

  • Participants
  • Parent commits 2a7da44

Comments (0)

Files changed (1)

File src/ksp/pc/impls/is/pcis.c

     ierr = VecGetSize(matis->x,&pcis->n);CHKERRQ(ierr);
     ierr = PetscMalloc(pcis->n*sizeof(PetscInt),&array);CHKERRQ(ierr);
     ierr = PetscMemzero(array,pcis->n*sizeof(PetscInt));CHKERRQ(ierr);
-    for (i=1;i<pcis->n_neigh;i++)
+    for (i=0;i<pcis->n_neigh;i++)
       for (j=0;j<pcis->n_shared[i];j++)
           array[pcis->shared[i][j]] += 1;