[petsc-dev] -ksp_view_binary with DA vectors

[Lisandro Dalcin]

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().

> 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