[petsc-dev] Add natural-to-global and global-to-natural operations for DM?

[Lisandro Dalcin]

On 20 March 2013 06:38, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> On Tue, Mar 19, 2013 at 9:54 PM, Richard Tran Mills <rtm at eecs.utk.edu>
> wrote:
>>
>
> So I would be reluctant to add the Natural concept since it does not make
> sense with all DMs. What do you do with the natural vector coming out of
> here? Why is it easier to implement dispatch in DM rather than in your code?
>

IMHO, global<->natural mapping does make sense on sparse problems.
What about making global<->natural to be the identity map (i.e, just
use VecCopy or VecAXPY depending on InsertMode) for any DM where the
global<->natural does not make sense?



--
Lisandro Dalcin
---------------
CIMEC (INTEC/CONICET-UNL)
Predio CONICET-Santa Fe
Colectora RN 168 Km 472, Paraje El Pozo
3000 Santa Fe, Argentina
Tel: +54-342-4511594 (ext 1011)
Tel/Fax: +54-342-4511169