fixed in master and next with the commit a9de710 Barry
On May 29, 2015, at 1:08 PM, Barry Smith <[email protected]> wrote:
I will fix this.
Barry
On May 29, 2015, at 10:40 AM, Mark Adams <[email protected]> wrote:
I have a 2 x 2 MatComposite and PCSetup is hanging because this branch below does "NEED" to be checked. This is in matrix.c around line 7500. Barry added this a year ago (and the code has been working so I am puzzled). I commented this code out and it works.
The 2nd block here is very small and is all on one processor, so all but one proc have (n == rend-rstart).
This code seems to work in maint, so I am puzzled. Perhaps this code was merged recently? Mark
if (!iscol || isrow == iscol) {
PetscBool stride;
ierr = PetscObjectTypeCompare((PetscObject)isrow,ISSTRIDE,&stride);CHKERRQ(ierr);
if (stride) {
PetscInt first,step,n,rstart,rend;
ierr = ISStrideGetInfo(isrow,&first,&step);CHKERRQ(ierr);
if (step == 1) {
ierr = MatGetOwnershipRange(mat,&rstart,&rend);CHKERRQ(ierr);
if (rstart == first) {
ierr = ISGetLocalSize(isrow,&n);CHKERRQ(ierr);
if (n == rend-rstart) {
/* special case grabbing all rows; NEED to do a global reduction to make sure all processes are doing this */
if (cll == MAT_INITIAL_MATRIX) {
*newmat = mat;
ierr = PetscObjectReference((PetscObject)mat);CHKERRQ(ierr);
}
PetscFunctionReturn(0);
}
}
}
}
}