[petsc-users] Create DM matrix

[Sang pham van]

The irregular rows is quite many. The matrix really needs to be
preallocated.
Could you show me how to use DMCreateMatrix_DA_3d_MPIAIJ() directly?

Pham
On Feb 11, 2016 9:52 AM, "Matthew Knepley" <knepley at gmail.com> wrote:

> On Wed, Feb 10, 2016 at 8:44 PM, Sang pham van <pvsang002 at 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" <knepley at gmail.com> wrote:
>>
>>> On Wed, Feb 10, 2016 at 6:14 PM, Sang pham van <pvsang002 at 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160211/cf5cf559/attachment.html>