Dear Roy, On Wed, 25 Feb 2009, Tim Kroeger wrote:
> I'm currently trying to test the ghosted vectors on my application, > but the queue on the cluster is full, so I have to wait. Update: My application crashes with the ghosted vectors, but I'm not yet sure whether it's a bug in the ghosted code or a bug in my application that just got triggered by the ghost dofs. I'll try to track this down, but that will be next week. >> Besides, this can't be the right fix, can it? Why should any parallel >> communication be required to copy one consistent ghosted vector to >> another? It seems as if VecCopy is simply failing to copy the whole >> vec, and I don't like not understanding why. > > Good question. Perhaps one should apply VecCopy to the local form > (obtained via VecGhostGetLocalForm())? Also, I wonder whether other > methods, such as VecScale(), might possibly suffer from similar > problems? > > I'll ask in the petsc-users list and let you know what they say. The PETSc people say that the local form has to be used (or VecGhostUpdate*(), as we used in PetscVector::close()). See also Jed's reply. I'll go through the PetscVector code once more and change this in all methods that are expected to work without parallel communication. This will most probably also happen next week. Best Regards, Tim -- Dr. Tim Kroeger tim.kroe...@mevis.fraunhofer.de Phone +49-421-218-7710 tim.kroe...@cevis.uni-bremen.de Fax +49-421-218-4236 Fraunhofer MEVIS, Institute for Medical Image Computing Universitaetsallee 29, 28359 Bremen, Germany ------------------------------------------------------------------------------ Open Source Business Conference (OSBC), March 24-25, 2009, San Francisco, CA -OSBC tackles the biggest issue in open source: Open Sourcing the Enterprise -Strategies to boost innovation and cut costs with open source participation -Receive a $600 discount off the registration fee with the source code: SFAD http://p.sf.net/sfu/XcvMzF8H _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel