Derek Gaston wrote: > Wouldn't the above loop over a _lot_ (like millions in some cases) of > unnecessary nodes when using Serial mesh?
At the moment, I think so, but that just makes it a wash efficiency-wise right now, and it wouldn't be hard to adjust our iterators to be more efficient on a partitioned SerialMesh later. > 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? You loop over every component from 0 to n_comp(s,var). Usually that'll mean "you use comp=0", but non-isoparametric Lagrange elements can have no components on some nodes, and non-Lagrange elements can have more than one on some nodes. > 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.... These don't sound like good options. --- Roy ------------------------------------------------------------------------- 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
