[petsc-users] How to construct Jacobian for multicomponent problem ?

[Christophe Ortiz]

>
>
>>>
>>>
>>
>> Ok. But what I still don't understand is the ordering of the residual in
>> the multicomponent case. Let's say dof=2, components u and v. Which is the
>> correct ordering in the following ?
>>
>> Fu0, Fu1, Fu2,...Fumx-1, Fv0, Fv1, Fv2, ...Fvmx-1
>>
>> or
>>
>> Fu0, Fv0, Fu1, Fv1, Fu2, Fv2, ...    ?
>>
>
> This one.
>

Ok.


>
>
>> It is this point that is not clear to me.
>>
>> Also, which subroutine do you advice to use to fill in the Matrix ?
>> MatSetValues ? MatSetValuesBlocked ?
>>
>
> MatSetValuesStencil.
>


Thanks a lot Matt for your prompt replies !

Christophe


>
>    Matt
>
>
>> Christophe
>>
>>
>>
>
>
>
> --
> 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/20130925/c19f5d57/attachment.html>