On Wed, Feb 10, 2016 at 8:44 PM, Sang pham van <pvsang...@gmail.com> wrote:
> That is because my matrix has some rows which need more entries than usual. > If its only a few, you could just turn off the allocation error. > Where can i find source code of DMCreateMatrix()? > > https://bitbucket.org/petsc/petsc/src/827b69d6bb12709ff9b9a0dede31640477fc2b74/src/dm/impls/da/fdda.c?at=master&fileviewer=file-view-default#fdda.c-1024 Matt > Pham. > On Feb 11, 2016 8:35 AM, "Matthew Knepley" <knep...@gmail.com> wrote: > >> On Wed, Feb 10, 2016 at 6:14 PM, Sang pham van <pvsang...@gmail.com> >> wrote: >> >>> Hi, >>> >>> I am trying to create a DM matrix with DMCreateMatrix_DA_3d_MPIAIJ() >>> instead of using DMCreateMatrix(). >>> >> >> Why, that should be called automatically by DMCreateMatrix()? >> >> Matt >> >> >>> Which header file should I include to use that routine? also, what is >>> the source file containing the DMCreateMatrix() routine? >>> >>> Many thanks in advance. >>> Pham >>> >> >> >> >> -- >> What most experimenters take for granted before they begin their >> experiments is infinitely more interesting than any results to which their >> experiments lead. >> -- Norbert Wiener >> > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener