Import coordinates

fei · 6 June 2022 14:20

Hello, I want to import the coordinates of a bunch of points, and store the coordinates in an array as global variables. Is it possible? could anyone give me some suggestions?

Thaks a lot.

fdw · 6 June 2022 19:22

Not without modifying the code. From then it really depends what you want to do with the points. There is a big difference between using the points in say a Python plugin compared with code on the backend.

Regards, Freddie.

fei · 7 June 2022 15:04

Yes, I want to modify the code. Can you give me some ideas?

fdw · 7 June 2022 19:52

The first thing you need to do is figure out how exactly what you want to do fits within the high-order FR framework. Only once you’ve done this and come up with a formulation which does not impact the order of accuracy is it sensible to start thinking about how to code up the functionality.

Regards, Freddie.

fei · 12 June 2022 14:57

Thanks.
I want to use the coordinates of solution points in the physical space, how can I use the dictionary vars to get all solution points coordinates in other py file? Are there any examples for my reference?

Best regards

fdw · 13 June 2022 13:24

It depends what Python file and in what context.

However, a quick search shows several values which make use of these coordinates:

github.com

PyFR/PyFR/blob/c9d02b004e70a71536cb423a0e0f39f66c59462a/pyfr/plugins/integrate.py#L44


      
                        for k in self.cfg.items(cfgsect)
                        if k.startswith('int-')]
          
          
# Integration region pre-processing
          esetmask = self._prepare_esetmask(intg)
          
          
# Gradient pre-processing
          self._init_gradients(intg)
          
          
# Save a reference to the physical solution point locations
          self.plocs = system.ele_ploc_upts
          
          
# Integration parameters
          self.nsteps = self.cfg.getint(cfgsect, 'nsteps')
          
          
# The root rank needs to open the output file
          if rank == root:
              header = ['t'] + [k for k in self.cfg.items(cfgsect)
                                if k.startswith('int-')]
          
          
    # Open

github.com

PyFR/PyFR/blob/d561c878dc0b517409409f4371bec0d08c399cbe/pyfr/plugins/sampler.py#L194


      
                  'hex': 'gauss-legendre',
                  'pri': 'williams-shunn~gauss-legendre',
                  'pyr': 'gauss-legendre',
                  'tet': 'shunn-ham'
              }
          
          
    for etype, eles in intg.system.ele_map.items():
                  pts = get_quadrule(etype, qrule_map[etype], eles.basis.nupts).pts
          
          
        tlocs.append(pts)
                  plocs.append(eles.ploc_at_np(pts).swapaxes(1, 2))
          
          
    return tlocs, plocs
          
          
def _refine_pts(self, intg, elepts):
              elelist = intg.system.ele_map.values()
              ptsinfo = []
          
          
    # Loop over all the points for each element type
              for etype, (eles, epts) in enumerate(zip(elelist, elepts)):
                  if not epts:

for example.

Regards, Freddie.

fei · 22 May 2023 15:50

In this way, x = plocs[0], y = plocs[1], z = plocs[2], are these right?

And when I use the physical coor of solution point, is it necessary to enter the element type?

WillT · 23 May 2023 07:46

The element type is handled in the loop, eles is an element object for a single element type.

As to where x, y and z are in that array, I imagine you’re right but check the shape of the array.

Topic		Replies	Views
Accessing coordinates of solution points Just Starting	3	348	13 July 2021
Sampling pyfrs files at certain coordinates General	4	222	21 May 2022
Solutions at specific points General	1	171	2 June 2014
Data location in high order cells and calculation of dimensionless wall normal distance General	1	194	25 June 2018
Computation of time step for FR General	6	219	23 October 2014

Import coordinates

Related topics