[petsc-users] Create DM matrix

[Matthew Knepley]

On Wed, Feb 10, 2016 at 8:59 PM, Sang pham van <pvsang002 at gmail.com> wrote:

> 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?
>
> Just put the declaration right into your source.

   Matt

> 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
>>
>


-- 
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/b96c9c55/attachment.html>