I tried running in the debug mode (should have done that earlier) and the code aborted with the following error.
No index 128 in ghosted vector. Vector contains [1352,2604) And empty ghost array. [1] ./include/libmesh/petsc_vector.h, line 1176, compiled May 22 2013 at 22:27:31 It also was trying to get value for an off-processor element. So, this perhaps explains the hans for the optimized executable. I was assuming that the project solution will work with only the local elements, but that does not seem to be the case. Will try with serialization of the vector. Manav On May 24, 2013, at 10:08 AM, Derek Gaston <[email protected]> wrote: > What are you trying to do? There are many scenarios where this could > fail... and some where it might succeed. > > Derek > > Sent from my iPad > > On May 24, 2013, at 8:02 AM, Manav Bhatia <[email protected]> wrote: > >> >> >> On May 24, 2013, at 9:58 AM, Derek Gaston <[email protected]> wrote: >> >>> Yes. >>> >>> Ok, well, MeshFunction itself doesn't care (as you found out)... but >>> if you want to evaluate your solution in elements containing >>> off-processor degrees of freedom then you need to serialize that >>> vector. >> >> >> But if the element is on a given processor, then wouldn't the local portion >> of system.solution contain all the relevant dofs to calculate the solution >> on that element? >> >> If so, then passing it the system.solution should be alright. No? >> >> >>> Think of it this way: MeshFunction doesn't do parallel >>> communication... you are responsible for feeding it the values you >>> need. >>> >>> Derek >>> >>> Sent from my iPad >> ------------------------------------------------------------------------------ Try New Relic Now & We'll Send You this Cool Shirt New Relic is the only SaaS-based application performance monitoring service that delivers powerful full stack analytics. Optimize and monitor your browser, app, & servers with just a few lines of code. Try New Relic and get this awesome Nerd Life shirt! http://p.sf.net/sfu/newrelic_d2d_may _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
