On Nov 6, 2008, at 8:57 AM, Kirk, Benjamin (JSC-EG) wrote:
Wouldn't the above loop over a _lot_ (like millions in some cases) of unnecessary nodes when using Serial mesh?Right...The fastest thing to do may be to loop over all local elements, and only look at their internal degreees of freedom, which are guaranteed to be owned by the processor. Then you can loop over all local nodes. Should be faster than the typical loop over elements / loop over element's nodes with-embedded-if-test...
Also I don't quite see how you get the dof_indices for a variable from a node.... if you use the dof_number()... how do you know the component?
My current algorithm loops over local_active_elements... calls dof_indices for that variable... but then I'm kinda stuck. I have all of the dof indices but I don't know if I own them. One option is just to pass them into the NumericVector to get added to the norm calculation... and just have NumericVector throw away the ones that it doesn't own....
What do you think? Derek
------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
