Is there a smarter way to select and *get?

[francesco123212]

Hi all,

I have a question concerning speed. I have implemented a procedure that does the following:

1. select nodes on a line
2. for each node do the following:
   • select the adjacent node on the line (the one immediately before and after for computing moving average and creating a reference system later on in the code)
   • get all elements attached to that node; for each of these element:
   • get nodes attached to the element

To speed up things I interact with the database only to get element and node IDs, then I ask for the coordinates using mapdl.mesh.nodes which I noted is much much faster. So I kinda have the following pseudo-code:

for line in lines
    for node in nodes_on_line
       # select the node
       # select the adjacent nodes IDs
       get_IDsList("NODE")
       # select the elements attached to the node
       get_IDsList("ELEM")
       for element in elements_attached2node
           # select the nodes of the element
           get_IDsList("NODE")
      end
   end
end

where the function is defined as follows get_IDsList()

function get_IDs(string::String)
   count = get_count(string, check=false)
   IDsList = zeros(Int,count)
   ID = 0
   for i in 1:count
      ID = int(mapdl.get("tmp",string,ID,"NXTH"))
      IDsList[i] = ID
   end
   return count, IDsList
end

so I am basically interacting with the database a lot, and it seems to be very timeconsuming (#757). Is there a smarter way to do this? Maybe I can speedup the get_IDs bit by doing:

function get_IDs_nodes()
   IDsList = mapdl.mesh.nnum
   count = int(length(IDsList))
end

But still, I have all the slow nodes/elements selections.
Thanks all for ur support :)

[francesco123212]

I have modified my get_IDs(string) function as follows and it is already much faster, but the mapdl.nsel(), mapdl.lsel(), mapdl.nsle(),... are still very slow.

function get_IDs(string::String; method::String="database")
   if string == "NODE" && method == "apimesh"
      IDsList = mapdl.mesh.nnum
      count = int(length(IDsList))
      return count, IDsList
   elseif string == "ELEM" && method == "apimesh"
      IDsList = mapdl.mesh.enum
      count = int(length(IDsList))
      return count, IDsList
   elseif string == "KP" && method == "geometrymesh"
      IDsList = mapdl.geometry.knum
      count = int(length(IDsList))
      return count, IDsList
   elseif string == "LINE" && method == "geometrymesh"
      IDsList = mapdl.geometry.lnum
      count = int(length(IDsList))
      return count, IDsList
   elseif string == "AREA" && method == "geometrymesh"
      IDsList = mapdl.geometry.anum
      count = int(length(IDsList))
      return count, IDsList
   elseif string == "VOLU" && method == "geometrymesh"
      IDsList = mapdl.geometry.vnum
      count = int(length(IDsList))
      return count, IDsList
   end
   count = get_count(string, check=false)
   IDsList = zeros(Int,count)
   ID = 0
   for i in 1:count
      ID = int(mapdl.get("tmp",string,ID,"NXTH"))
      IDsList[i] = ID
   end
   return count, IDsList
end