# Development

**URL:** https://pyfr.discourse.group/c/dev/14.md?page=1

[Latest](https://pyfr.discourse.group/latest.md) · [Categories](https://pyfr.discourse.group/categories.md) · [Tags](https://pyfr.discourse.group/tags.md)

**Page:** 2

---

## [How does pyfr organize its data](https://pyfr.discourse.group/t/how-does-pyfr-organize-its-data/903)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 1\
**Last updated:** [29 August 2023 16:33 UTC](https://pyfr.discourse.group/t/how-does-pyfr-organize-its-data/903 "2023-08-29T16:33:02Z")

</div>

Hi, everyone . I really want to know how Pyfr organizes its data? Where does it appear in the code?

---

## [Import and update a matrix during runtime](https://pyfr.discourse.group/t/import-and-update-a-matrix-during-runtime/887)

<div class="topic-metadata">

**Author:** [@Edmund](https://pyfr.discourse.group/u/Edmund)\
**Replies:** 5\
**Last updated:** [3 August 2023 12:57 UTC](https://pyfr.discourse.group/t/import-and-update-a-matrix-during-runtime/887 "2023-08-03T12:57:29Z")

</div>

Hi, everyone I wonder if I can update some data in the kernel run like time. Because I’m trying to use the actuator line method in PyFR, which requires updating some parameters in every step.I can only transfer individu…

---

## [Source term for some particular solution points](https://pyfr.discourse.group/t/source-term-for-some-particular-solution-points/821)

<div class="topic-metadata">

**Author:** [@fei](https://pyfr.discourse.group/u/fei)\
**Replies:** 6\
**Last updated:** [5 June 2023 11:32 UTC](https://pyfr.discourse.group/t/source-term-for-some-particular-solution-points/821 "2023-06-05T11:32:38Z")

</div>

Hello, The source-terms-plugin adds the source term to all solution points, but I want to add source term for some particular solution points, could someone give me some tips?

---

## [Zonal P-order in pyFR](https://pyfr.discourse.group/t/zonal-p-order-in-pyfr/809)

<div class="topic-metadata">

**Author:** [@Frankx9](https://pyfr.discourse.group/u/Frankx9)\
**Replies:** 1\
**Last updated:** [3 May 2023 07:56 UTC](https://pyfr.discourse.group/t/zonal-p-order-in-pyfr/809 "2023-05-03T07:56:44Z")

</div>

Hi, in pyFR is there any possibility to use different polynomial order for different zones in the domain? Regards

---

## [Import coordinates](https://pyfr.discourse.group/t/import-coordinates/619)

<div class="topic-metadata">

**Author:** [@fei](https://pyfr.discourse.group/u/fei)\
**Replies:** 7\
**Last updated:** [23 May 2023 07:46 UTC](https://pyfr.discourse.group/t/import-coordinates/619 "2023-05-23T07:46:45Z")

</div>

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.

---

## [Mathlib in pyfr question](https://pyfr.discourse.group/t/mathlib-in-pyfr-question/803)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 1\
**Last updated:** [20 April 2023 09:34 UTC](https://pyfr.discourse.group/t/mathlib-in-pyfr-question/803 "2023-04-20T09:34:00Z")

</div>

As far as I know, the math library currently in use by pyfr is gimmik and blas. I noticed that nvidia has a matrix multiplying tool that uses the sparse feature of the Ampere architecture called cuSPARSELt. Can it be use…

---

## [Cuda mpi aware question](https://pyfr.discourse.group/t/cuda-mpi-aware-question/797)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 3\
**Last updated:** [19 April 2023 09:39 UTC](https://pyfr.discourse.group/t/cuda-mpi-aware-question/797 "2023-04-19T09:39:46Z")

</div>

As shown, there is a code in pyfr to determine whether to use cuda aware and then determine whether to send the data to the host, is this necessary? When cuda aware can’t send data directly between GPUs, won’t it auto…

---

## [Restart simulation with different resources](https://pyfr.discourse.group/t/restart-simulation-with-different-resources/798)

<div class="topic-metadata">

**Author:** [@Frankx9](https://pyfr.discourse.group/u/Frankx9)\
**Replies:** 4\
**Last updated:** [17 April 2023 10:23 UTC](https://pyfr.discourse.group/t/restart-simulation-with-different-resources/798 "2023-04-17T10:23:48Z")

</div>

Hi, is there any plan on allowing a restart of a given simulation using different resources (i.e. number of partitions)? Say, for example, that I’m running with cuda with X GPUs and want to switch to openmp with Y CPUs …

---

## [What's the k in AoSoA?](https://pyfr.discourse.group/t/whats-the-k-in-aosoa/749)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 4\
**Last updated:** [7 March 2023 04:31 UTC](https://pyfr.discourse.group/t/whats-the-k-in-aosoa/749 "2023-03-07T04:31:00Z")

</div>

Hello, I have some trouble in understanding the format of AoSoA. From the introduction of pyfr data structure, I know that the memory arrangement of AoSoA is as shown in the figure.I guess the k in the picture correspond…

---

## [Correcting discontinuous field](https://pyfr.discourse.group/t/correcting-discontinuous-field/759)

<div class="topic-metadata">

**Author:** [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Replies:** 1\
**Last updated:** [2 March 2023 14:54 UTC](https://pyfr.discourse.group/t/correcting-discontinuous-field/759 "2023-03-02T14:54:49Z")

</div>

Hi, I wonder if there’s anyway to corrected discontinued field especially for unstructured mesh? No matter in visualization or simply doing exact interpolation, results from unstructured mesh always give strong disconti…

---

## [Plugins: source terms and performance](https://pyfr.discourse.group/t/plugins-source-terms-and-performance/753)

<div class="topic-metadata">

**Author:** [@MohamedKamra](https://pyfr.discourse.group/u/MohamedKamra)\
**Replies:** 5\
**Last updated:** [21 February 2023 21:57 UTC](https://pyfr.discourse.group/t/plugins-source-terms-and-performance/753 "2023-02-21T21:57:24Z")

</div>

Greetings everyone, I have two questions: The first question is: Is it possible for a plug-in to specifically evaluate the solver’s source term? My second question: by default, plug-ins are executed as Python code; is…

---

## [Where does pyfr use sympy?](https://pyfr.discourse.group/t/where-does-pyfr-use-sympy/748)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 2\
**Last updated:** [21 February 2023 07:52 UTC](https://pyfr.discourse.group/t/where-does-pyfr-use-sympy/748 "2023-02-21T07:52:38Z")

</div>

The paper “PyFR: An open source framework for solving advection–diffusion type problems on streaming architectures using the flux” says that： To keep the codebase compact PyFR makes extensive use of symbolic manipulat…

---

## [Writing gradients to disk](https://pyfr.discourse.group/t/writing-gradients-to-disk/739)

<div class="topic-metadata">

**Author:** [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Replies:** 1\
**Last updated:** [11 January 2023 15:35 UTC](https://pyfr.discourse.group/t/writing-gradients-to-disk/739 "2023-01-11T15:35:14Z")

</div>

Hi guys, I wonder if there’s a stable way to write gradients to disk? For my case I need collecting a lot of snapshots as well as gradients for post processing. For now, I modified tavg plugin to output gradient without…

---

## [Gradient export flag](https://pyfr.discourse.group/t/gradient-export-flag/527)

<div class="topic-metadata">

**Author:** [@fdw](https://pyfr.discourse.group/u/fdw)\
**Replies:** 2\
**Last updated:** [20 January 2022 00:01 UTC](https://pyfr.discourse.group/t/gradient-export-flag/527 "2022-01-20T00:01:10Z")

</div>

Is anyone actively using the gradient export functionality of PyFR? That is to say export -g …? In my experience it tends to produce worse results than current versions of ParaView’s gradient of unstructured data set f…

---

## [Matrix dimensions and cublas kernels](https://pyfr.discourse.group/t/matrix-dimensions-and-cublas-kernels/721)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 2\
**Last updated:** [13 December 2022 09:13 UTC](https://pyfr.discourse.group/t/matrix-dimensions-and-cublas-kernels/721 "2022-12-13T09:13:40Z")

</div>

Dear PyFR team, I am trying to use cublas in PyFR and I am using the Euler vortex case for testing. However, I am confused because when I print the matrix that is input to cublas, I found that it is three-dimensional (6…

---

## [View matrices in interface](https://pyfr.discourse.group/t/view-matrices-in-interface/717)

<div class="topic-metadata">

**Author:** [@Zhenyang](https://pyfr.discourse.group/u/Zhenyang)\
**Replies:** 2\
**Last updated:** [1 December 2022 08:47 UTC](https://pyfr.discourse.group/t/view-matrices-in-interface/717 "2022-12-01T08:47:45Z")

</div>

Hi guys, I have a weird question regarding implementing linear solver. After some discussion we decide to use (rho’, rho\*u’, p’) as our variables ( ’ here means perturbation and without ’ is baseflow). So we need to dea…

---

## [Meaning of variable names in PyFR source code](https://pyfr.discourse.group/t/meaning-of-variable-names-in-pyfr-source-code/713)

<div class="topic-metadata">

**Author:** [@Yubo](https://pyfr.discourse.group/u/Yubo)\
**Replies:** 2\
**Last updated:** [4 November 2022 15:42 UTC](https://pyfr.discourse.group/t/meaning-of-variable-names-in-pyfr-source-code/713 "2022-11-04T15:42:05Z")

</div>

Hi Developers, I am reading the source code of PyFR. But I am confused by the abbreviations of some variables. For example, in PyFR/pyfr/shapes.py: 1. self.upts 2. self.nupts 3. self.qpts 4. self.nqpts 5. self.fpts 6. …

---

## [Calculating cell averages in Mako macros](https://pyfr.discourse.group/t/calculating-cell-averages-in-mako-macros/710)

<div class="topic-metadata">

**Author:** [@fei](https://pyfr.discourse.group/u/fei)\
**Replies:** 16\
**Last updated:** [31 October 2022 16:47 UTC](https://pyfr.discourse.group/t/calculating-cell-averages-in-mako-macros/710 "2022-10-31T16:47:15Z")

</div>

Hi, Assuming cell i, I want to compute cell average of solution value u\_{i}, u\_{i+1} in macro environment. How to implement this? And are there any examples I can refer to? Second, Is it possible to install two differe…

---

## [Persson shock sensor, implementation and config](https://pyfr.discourse.group/t/persson-shock-sensor-implementation-and-config/709)

<div class="topic-metadata">

**Author:** [@randommc](https://pyfr.discourse.group/u/randommc)\
**Replies:** 3\
**Last updated:** [23 October 2022 07:17 UTC](https://pyfr.discourse.group/t/persson-shock-sensor-implementation-and-config/709 "2022-10-23T07:17:41Z")

</div>

Hi I have a question about the shock detector in FR, I know the method is from “P. O. Persson and J. Peraire. Sub-cell shock capturing for discontinuous Galerkin methods. 44th AIAA Aerospace Sciences Meeting and Exhibit…

---

## [Gradients without writing to disk](https://pyfr.discourse.group/t/gradients-without-writing-to-disk/436)

<div class="topic-metadata">

**Author:** [@tsas](https://pyfr.discourse.group/u/tsas)\
**Replies:** 13\
**Last updated:** [21 October 2022 07:47 UTC](https://pyfr.discourse.group/t/gradients-without-writing-to-disk/436 "2022-10-21T07:47:01Z")

</div>

Is there a way to access the gradients of a PyFR simulation from within Python, similar to integ.soln, without the need to write the solution file to disk and using pyfr export -g? I intend to run a large number of simu…

---

## [Installing mulitple pyfr packages](https://pyfr.discourse.group/t/installing-mulitple-pyfr-packages/694)

<div class="topic-metadata">

**Author:** [@fei](https://pyfr.discourse.group/u/fei)\
**Replies:** 3\
**Last updated:** [27 September 2022 07:56 UTC](https://pyfr.discourse.group/t/installing-mulitple-pyfr-packages/694 "2022-09-27T07:56:30Z")

</div>

My computer has already installed pyfr package, and I want to install another pyfr package. Should I delete the original pyfr installation package, then install the new one. If so, I copy the new pyfr package to the site…

---

## [Element-wise operations](https://pyfr.discourse.group/t/element-wise-operations/683)

<div class="topic-metadata">

**Author:** [@EitanA](https://pyfr.discourse.group/u/EitanA)\
**Replies:** 2\
**Last updated:** [29 August 2022 14:39 UTC](https://pyfr.discourse.group/t/element-wise-operations/683 "2022-08-29T14:39:44Z")

</div>

Hello all I am wondering how element wise operations can be achieved in pyfr (beyond matrix multiplication)? In particular, I am interested in extracting the element average/0th order polynomial, and the minimum/maximu…

---

## [Question about Gimmik flops and Arithmetic intensity](https://pyfr.discourse.group/t/question-about-gimmik-flops-and-arithmetic-intensity/673)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 1\
**Last updated:** [12 August 2022 11:03 UTC](https://pyfr.discourse.group/t/question-about-gimmik-flops-and-arithmetic-intensity/673 "2022-08-12T11:03:28Z")

</div>

I read paper about gimmik package GiMMiK—Generating bespoke matrix multiplication kernels for accelerators: Application to high-order Computational Fluid Dynamics - ScienceDirect . When we use high-order scheme, the kern…

---

## [Where are the "atomic" operations used in PyFR?](https://pyfr.discourse.group/t/where-are-the-atomic-operations-used-in-pyfr/672)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 1\
**Last updated:** [11 August 2022 11:01 UTC](https://pyfr.discourse.group/t/where-are-the-atomic-operations-used-in-pyfr/672 "2022-08-11T11:01:50Z")

</div>

I notice in pyfr/banckend/hip/driver.py file in class “HIPDevProps” , some atomic functions are listed, like this , where are these operations used in PyFR ?

---

## [Does pyfr support hybrid precision？](https://pyfr.discourse.group/t/does-pyfr-support-hybrid-precision/666)

<div class="topic-metadata">

**Author:** [@luli](https://pyfr.discourse.group/u/luli)\
**Replies:** 1\
**Last updated:** [9 August 2022 09:10 UTC](https://pyfr.discourse.group/t/does-pyfr-support-hybrid-precision/666 "2022-08-09T09:10:05Z")

</div>

Hybrid precision is a scheme to ensure both accuracy and efficiency ，does pyfr support hybrid precision? Such like hybrid fp32/fp64?

[Previous page](https://pyfr.discourse.group/c/dev/14.md)
