[Trilinos-Users] Query STK entities for a global index?
kevin071586 at yahoo.com
Tue Dec 27 22:28:06 EST 2016
Hi, I'm just beginning to explore Trilinos and in particular the STK package. The STK package is fantastic and I've been able to get it to do a lot of things with the provided examples like read a decomposed mesh, write transient data and fields, etc., but there is a very simple query that I'm having trouble with:
I would like to obtain a global contiguous index for a given entity such as NODE_RANK or ELEMENT_RANK. I've created my exodus file so that there is an element number map so that the identifier() function returns IDs which are not sequential. For example, I may have 5 elements distributed on two MPI processes so that:
All IDs = [200, 400, 600, 800, 900]Processor 1 ID's = [200, 400, 600]Processor 2 ID's = [800, 900]
I can query stkBulkData.local_id( entity ) to get [0,1,2] and [0,1].I can query stkBulkData.identifier( entity ) to get [200, 400, 600] and [800, 900]
Is there a way to query an entity for global index so that I get [0, 1, 2] and [3, 4]?
I've searched all the documentation, code, and examples and can't find anything like this so I have to imagine I'm either (a) overlooking something very obvious, or (b) trying to do something that is not an intended usage of the STK package. For what it's worth, I'm trying to set up the global indices for a global stiffness matrix where each NODE_RANK entity may have multiple degrees of freedom and using a 'global index' seemed like a natural way to reference the correct rows and columns in the global stiffness matrix.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Trilinos-Users