Skip to content

Commit

Permalink
MatSetLocalToGlobalMapping update for petsc-dev
Browse files Browse the repository at this point in the history
  • Loading branch information
jedbrown committed Dec 21, 2010
1 parent 93502df commit d150390
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
4 changes: 2 additions & 2 deletions src/fs/interface/fs.c
Original file line number Diff line number Diff line change
Expand Up @@ -659,8 +659,8 @@ dErr dFSGetMatrix(dFS fs,const MatType mtype,Mat *inJ)
err = MatMPIAIJSetPreallocation(J,bs*27,NULL,bs*25,NULL);dCHK(err);
}
err = MatSetBlockSize(J,bs);dCHK(err);
err = MatSetLocalToGlobalMappingBlock(J,fs->bmapping);dCHK(err);
err = MatSetLocalToGlobalMapping(J,fs->mapping);dCHK(err);
err = MatSetLocalToGlobalMappingBlock(J,fs->bmapping,fs->bmapping);dCHK(err);
err = MatSetLocalToGlobalMapping(J,fs->mapping,fs->mapping);dCHK(err);

/* We want the resulting matrices to be usable with matrix-free operations based on this FS */
err = PetscObjectCompose((dObject)J,"DohpFS",(dObject)fs);dCHK(err);
Expand Down
4 changes: 2 additions & 2 deletions src/util/util.c
Original file line number Diff line number Diff line change
Expand Up @@ -40,12 +40,12 @@ dErr dRealTableView(dInt m,dInt n,const dReal mat[],const char *name,dViewer vie
if (name) {
err = PetscViewerASCIIPrintf(viewer,"%10s[%2d][%2d:%2d] ",name,i,0,n);dCHK(err);
}
err = PetscViewerASCIIUseTabs(viewer,PETSC_NO);dCHK(err);
err = PetscViewerASCIIUseTabs(viewer,PETSC_FALSE);dCHK(err);
for (dInt j=0; j<n; j++) {
err = PetscViewerASCIIPrintf(viewer," % 9.5f",mat[i*n+j]);dCHK(err);
}
err = PetscViewerASCIIPrintf(viewer,"\n");dCHK(err);
err = PetscViewerASCIIUseTabs(viewer,PETSC_YES);dCHK(err);
err = PetscViewerASCIIUseTabs(viewer,PETSC_TRUE);dCHK(err);
}
dFunctionReturn(0);
}

0 comments on commit d150390

Please sign in to comment.