Author Topic: parfeap output  (Read 5370 times)

ace21

  • Jr. Member
  • **
  • Posts: 27
parfeap output
« on: February 06, 2018, 05:05:03 AM »
Dear feap users

Is it possible to get the global id of the elements, Actually I have a user element which  computes some element level quantity,

In order to extract it from the partitioned input files,I need the global element id corresponding to the  local one in each  partition

Or is it possible to write the output in a single output file

regards

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2647
Re: parfeap output
« Reply #1 on: February 06, 2018, 09:45:49 AM »
If you want the global node numbers they are in and array with name 'GN   ', located at mr(np(244)).

I do not think the global element numbers are available in each partition.

Prof. S. Govindjee

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 1160
Re: parfeap output
« Reply #2 on: February 06, 2018, 02:44:28 PM »
Prof. Taylor is correct.  Only the local node to global node numbering mapping is available within the elements, in the GN array.
If you want to have the global element numbers available you will have to edit parfeap/pmacr7.F to output this information to the parallel input files (and then also read it in).  To do so, you can output the information in the array IXINV( ) which is mr(np(112)) and maps global node numbers to element node numbers (see parpfeap/pmacr7.F(create_ixinv)  and parfeap/pmacr7.F(make_node_graph) for help on the structure of the data).

-sg

ace21

  • Jr. Member
  • **
  • Posts: 27
Re: parfeap output
« Reply #3 on: February 08, 2018, 02:39:42 AM »
Dear Prof Taylor and  Prof S.Govindjee

I will study the files that you have mentioned for extracting the global id of the elements, I have just started using parfeap so will take some time to understand the code,your advises are very helpful in this respect

If I am not wrong the partitioned input files contain the information block "LOCAl to GLOBal node numbers", the second column has the global nodal ids

So mean while for my problem , I have written a matlab script that: extract the global ids of nodes, replace the local one in the connectivity array of element present in each partition input file, compare the modified connectivity list with the one in a non-partitioned mesh/input file to get the global ids

I still have to take care of the duplicate/ghost  element ids, hope my approach is correct

regards