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
