[petsc-users] Create DM matrix

[Matthew Knepley]

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/20160210/28ff9f07/attachment.html>