On Thu, 14 Jun 2012, Dmitry Karpeev wrote:
What are the node iterators I need for this?
MeshBase::local_nodes_begin(),local_nodes_end() should give you the
range (on this processor) for nodes;
active_local_elements_begin)()/end() for elements.
I assume I'll need to call dof_map().get_indices() hereand something
else on the elements?
Nope; call DofObject::n_comp() and DofObject::dof_number() for both
elements and nodes.
---
Roy
On Wed, Jun 13, 2012 at 1:06 PM, Roy Stogner <[email protected]> wrote:
Suggestion: Do one loop over nodes adding all indices for each
variable at each node, then a loop over elements adding all indices
for each variable at each element (not at the nodes of each element,
but at the element-considered-as-DofObject itself). That way you
don't get any duplicates, don't need std::set, and you're
O(N) rather than O(N log N).
------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and
threat landscape has changed and how IT managers can respond. Discussions
will include endpoint security, mobile security and the latest in malware
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
