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
