On Tue, 24 Feb 2009, Roy Stogner wrote:

> On Wed, 18 Feb 2009, Roy Stogner wrote:
>
>> I'm swamped right now.  The next step is to reduce ex10 to a simpler
>> case (maybe solving du/dt = 0 in nD with some non-zero ICs?) where I
>> can step through the whole projection in gdb, but it's likely to be
>> weeks before I can get to that.
>
> I apparently underestimated the productivity benefits of trying to
> avoid twiddling my thumbs while the supercomputer I'm using undergoes
> maintenance.

Sounds great, although that sentence seems to have a wrong parity of 
negations.  Anyway, great that you found the time to work on ghosted 
vectors.

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.

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

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

Reply via email to