On Tue, 7 Feb 2012, Derek Gaston wrote:
The problem is that that version of the function explicitly builds full index
sets forcing complete communication. If we would
have just said "current_local_solution = solution" it would actually do the
correct ghosted assignment.
Interesting. So the quick fix is to pass the send_list in or use
operator= in project_solution; but really we ought to be updating that
localize() so that it works better on ghosted vectors too.
Yep. Dmitry was actually trying to compile libMesh against PETSc
built with 64bit integers just today! Apparently it didn't work out
because of petscint to int conversions in libMesh... oops!
The solution is clearly to make our DoF indexes all be a separate
typedef'ed variable, default to petscint when that's available,
default to int but configure overrideable otherwise. Crank up the
compiler warnings to catch local variables etc. and change their
declarations... all straightforward, but sounds like lots of work; let
me know how it goes! ;-)
---
Roy
------------------------------------------------------------------------------
Keep Your Developer Skills Current with LearnDevNow!
The most comprehensive online learning library for Microsoft developers
is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5, CSS3, MVC3,
Metro Style Apps, more. Free future releases when you subscribe now!
http://p.sf.net/sfu/learndevnow-d2d
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
