Sent from my iPad On May 12, 2011, at 5:22 PM, Roy Stogner <[email protected]> wrote:
> operator= is implemented correctly for src.type == dst.type, and as of > mid-last year it's implemented correctly in PetscVector when one type > is PARALLEL and the other is GHOSTED. That's it. The new > System::update() probably doesn't even work correctly for Petsc+MPI > when --enable-ghosted-vectors is off. Why would that be? It would just run the old code with the localize().... if it worked before it would still work now. As for operator=, David and I shored it up for ghosted = parallel a while ago. We've been using that capability without issue for a while. If there is a bug in it, I would like to get to the bottom of it. We rely fairly heavily on it. All of our regression tests pass with the new System::update()... Can you pinpoint a problem that doesn't work properly? We might be able to take a look at it. I do agree that there might be better ways... and I'm committed to looking into them... but we need to understand why this wouldn't work because that would mean there is a serious implementation issue with operator=. Derek ------------------------------------------------------------------------------ Achieve unprecedented app performance and reliability What every C/C++ and Fortran developer should know. Learn how Intel has extended the reach of its next-generation tools to help boost performance applications - inlcuding clusters. http://p.sf.net/sfu/intel-dev2devmay _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
