[petsc-dev] -ksp_view_binary with DA vectors

[Barry Smith]

On Jun 22, 2010, at 10:53 AM, Lisandro Dalcin wrote:

> On 22 June 2010 11:17, Jed Brown <jed at 59a2.org> wrote:
>> The matrix is always written in PETSc ordering, but DA vectors are
>> written in natural order.
> 
> I rally hate the view/load overriding hackery in DACreateGlobalVector().

    Hack? That's no hack, it is called inheritance :-)

    Actually we should add the same support for matrices and then the inconsistency would disappear.

   Barry

> 
>>  This means that loading a problem with ex10
>> produces a genuinely different RHS than the user had initially, and the
>> solution is truly different (not just different by a permutation) when
>> they write the problem out on a different number of procs.  Should
>> -ksp_view_binary push PETSC_VIEWER_NATIVE so that the vector is always
>> consistent with the matrix ordering?
>> 
> 
> I would very much prefer to remove the hackery in
> DACreateGlobalVector() and perhaps add a helper routine (let say
> DAViewGlobalVector()). However, I understand the backward
> compatibility issues.
> 
> 
> 
> -- 
> Lisandro Dalcin
> ---------------
> CIMEC (INTEC/CONICET-UNL)
> Predio CONICET-Santa Fe
> Colectora RN 168 Km 472, Paraje El Pozo
> Tel: +54-342-4511594 (ext 1011)
> Tel/Fax: +54-342-4511169