PyMHD 0.1.0__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- pymhd-0.1.0/.gitignore +221 -0
- pymhd-0.1.0/LICENSE +21 -0
- pymhd-0.1.0/PKG-INFO +100 -0
- pymhd-0.1.0/README.md +79 -0
- pymhd-0.1.0/pymhd/__init__.py +31 -0
- pymhd-0.1.0/pymhd/derivatives/TENO.py +278 -0
- pymhd-0.1.0/pymhd/derivatives/WENO.py +323 -0
- pymhd-0.1.0/pymhd/derivatives/__init__.py +24 -0
- pymhd-0.1.0/pymhd/derivatives/compact.py +365 -0
- pymhd-0.1.0/pymhd/derivatives/derivative.py +926 -0
- pymhd-0.1.0/pymhd/numdiss.py +598 -0
- pymhd-0.1.0/pymhd/plot/__init__.py +48 -0
- pymhd-0.1.0/pymhd/plot/nd.py +1519 -0
- pymhd-0.1.0/pymhd/plot/slc.py +648 -0
- pymhd-0.1.0/pymhd/plot/spc.py +249 -0
- pymhd-0.1.0/pymhd/preprocess/Athena.py +847 -0
- pymhd-0.1.0/pymhd/preprocess/__init__.py +69 -0
- pymhd-0.1.0/pymhd/preprocess/helper/NOTICE +42 -0
- pymhd-0.1.0/pymhd/preprocess/helper/bin_convert.py +2000 -0
- pymhd-0.1.0/pymhd/preprocess/helper/make_athdf.py +45 -0
- pymhd-0.1.0/pymhd/spectra.py +376 -0
- pymhd-0.1.0/pymhd/turbulence.py +917 -0
- pymhd-0.1.0/pyproject.toml +44 -0
- pymhd-0.1.0/tests/numdiss.test.py +107 -0
- pymhd-0.1.0/tests/preprocess.test.py +108 -0
- pymhd-0.1.0/tests/slice.test.py +65 -0
- pymhd-0.1.0/tests/spectra.test.py +56 -0
pymhd-0.1.0/.gitignore
ADDED
|
@@ -0,0 +1,221 @@
|
|
|
1
|
+
# Byte-compiled / optimized / DLL files
|
|
2
|
+
__pycache__/
|
|
3
|
+
*.py[codz]
|
|
4
|
+
*$py.class
|
|
5
|
+
|
|
6
|
+
# C extensions
|
|
7
|
+
*.so
|
|
8
|
+
|
|
9
|
+
# Distribution / packaging
|
|
10
|
+
.Python
|
|
11
|
+
build/
|
|
12
|
+
develop-eggs/
|
|
13
|
+
dist/
|
|
14
|
+
downloads/
|
|
15
|
+
eggs/
|
|
16
|
+
.eggs/
|
|
17
|
+
lib/
|
|
18
|
+
lib64/
|
|
19
|
+
parts/
|
|
20
|
+
sdist/
|
|
21
|
+
var/
|
|
22
|
+
wheels/
|
|
23
|
+
share/python-wheels/
|
|
24
|
+
*.egg-info/
|
|
25
|
+
.installed.cfg
|
|
26
|
+
*.egg
|
|
27
|
+
MANIFEST
|
|
28
|
+
|
|
29
|
+
# PyInstaller
|
|
30
|
+
# Usually these files are written by a python script from a template
|
|
31
|
+
# before PyInstaller builds the exe, so as to inject date/other infos into it.
|
|
32
|
+
*.manifest
|
|
33
|
+
*.spec
|
|
34
|
+
|
|
35
|
+
# Installer logs
|
|
36
|
+
pip-log.txt
|
|
37
|
+
pip-delete-this-directory.txt
|
|
38
|
+
|
|
39
|
+
# Unit test / coverage reports
|
|
40
|
+
htmlcov/
|
|
41
|
+
.tox/
|
|
42
|
+
.nox/
|
|
43
|
+
.coverage
|
|
44
|
+
.coverage.*
|
|
45
|
+
.cache
|
|
46
|
+
nosetests.xml
|
|
47
|
+
coverage.xml
|
|
48
|
+
*.cover
|
|
49
|
+
*.py.cover
|
|
50
|
+
.hypothesis/
|
|
51
|
+
.pytest_cache/
|
|
52
|
+
cover/
|
|
53
|
+
|
|
54
|
+
# Translations
|
|
55
|
+
*.mo
|
|
56
|
+
*.pot
|
|
57
|
+
|
|
58
|
+
# Django stuff:
|
|
59
|
+
*.log
|
|
60
|
+
local_settings.py
|
|
61
|
+
db.sqlite3
|
|
62
|
+
db.sqlite3-journal
|
|
63
|
+
|
|
64
|
+
# Flask stuff:
|
|
65
|
+
instance/
|
|
66
|
+
.webassets-cache
|
|
67
|
+
|
|
68
|
+
# Scrapy stuff:
|
|
69
|
+
.scrapy
|
|
70
|
+
|
|
71
|
+
# Sphinx documentation
|
|
72
|
+
docs/_build/
|
|
73
|
+
|
|
74
|
+
# PyBuilder
|
|
75
|
+
.pybuilder/
|
|
76
|
+
target/
|
|
77
|
+
|
|
78
|
+
# Jupyter Notebook
|
|
79
|
+
.ipynb_checkpoints
|
|
80
|
+
|
|
81
|
+
# IPython
|
|
82
|
+
profile_default/
|
|
83
|
+
ipython_config.py
|
|
84
|
+
|
|
85
|
+
# pyenv
|
|
86
|
+
# For a library or package, you might want to ignore these files since the code is
|
|
87
|
+
# intended to run in multiple environments; otherwise, check them in:
|
|
88
|
+
# .python-version
|
|
89
|
+
|
|
90
|
+
# pipenv
|
|
91
|
+
# According to pypa/pipenv#598, it is recommended to include Pipfile.lock in version control.
|
|
92
|
+
# However, in case of collaboration, if having platform-specific dependencies or dependencies
|
|
93
|
+
# having no cross-platform support, pipenv may install dependencies that don't work, or not
|
|
94
|
+
# install all needed dependencies.
|
|
95
|
+
# Pipfile.lock
|
|
96
|
+
|
|
97
|
+
# UV
|
|
98
|
+
# Similar to Pipfile.lock, it is generally recommended to include uv.lock in version control.
|
|
99
|
+
# This is especially recommended for binary packages to ensure reproducibility, and is more
|
|
100
|
+
# commonly ignored for libraries.
|
|
101
|
+
# uv.lock
|
|
102
|
+
|
|
103
|
+
# poetry
|
|
104
|
+
# Similar to Pipfile.lock, it is generally recommended to include poetry.lock in version control.
|
|
105
|
+
# This is especially recommended for binary packages to ensure reproducibility, and is more
|
|
106
|
+
# commonly ignored for libraries.
|
|
107
|
+
# https://python-poetry.org/docs/basic-usage/#commit-your-poetrylock-file-to-version-control
|
|
108
|
+
# poetry.lock
|
|
109
|
+
# poetry.toml
|
|
110
|
+
|
|
111
|
+
# pdm
|
|
112
|
+
# Similar to Pipfile.lock, it is generally recommended to include pdm.lock in version control.
|
|
113
|
+
# pdm recommends including project-wide configuration in pdm.toml, but excluding .pdm-python.
|
|
114
|
+
# https://pdm-project.org/en/latest/usage/project/#working-with-version-control
|
|
115
|
+
# pdm.lock
|
|
116
|
+
# pdm.toml
|
|
117
|
+
.pdm-python
|
|
118
|
+
.pdm-build/
|
|
119
|
+
|
|
120
|
+
# pixi
|
|
121
|
+
# Similar to Pipfile.lock, it is generally recommended to include pixi.lock in version control.
|
|
122
|
+
# pixi.lock
|
|
123
|
+
# Pixi creates a virtual environment in the .pixi directory, just like venv module creates one
|
|
124
|
+
# in the .venv directory. It is recommended not to include this directory in version control.
|
|
125
|
+
.pixi
|
|
126
|
+
|
|
127
|
+
# PEP 582; used by e.g. github.com/David-OConnor/pyflow and github.com/pdm-project/pdm
|
|
128
|
+
__pypackages__/
|
|
129
|
+
|
|
130
|
+
# Celery stuff
|
|
131
|
+
celerybeat-schedule
|
|
132
|
+
celerybeat.pid
|
|
133
|
+
|
|
134
|
+
# Redis
|
|
135
|
+
*.rdb
|
|
136
|
+
*.aof
|
|
137
|
+
*.pid
|
|
138
|
+
|
|
139
|
+
# RabbitMQ
|
|
140
|
+
mnesia/
|
|
141
|
+
rabbitmq/
|
|
142
|
+
rabbitmq-data/
|
|
143
|
+
|
|
144
|
+
# ActiveMQ
|
|
145
|
+
activemq-data/
|
|
146
|
+
|
|
147
|
+
# SageMath parsed files
|
|
148
|
+
*.sage.py
|
|
149
|
+
|
|
150
|
+
# Environments
|
|
151
|
+
.env
|
|
152
|
+
.envrc
|
|
153
|
+
.venv
|
|
154
|
+
env/
|
|
155
|
+
venv/
|
|
156
|
+
ENV/
|
|
157
|
+
env.bak/
|
|
158
|
+
venv.bak/
|
|
159
|
+
|
|
160
|
+
# Spyder project settings
|
|
161
|
+
.spyderproject
|
|
162
|
+
.spyproject
|
|
163
|
+
|
|
164
|
+
# Rope project settings
|
|
165
|
+
.ropeproject
|
|
166
|
+
|
|
167
|
+
# mkdocs documentation
|
|
168
|
+
/site
|
|
169
|
+
|
|
170
|
+
# mypy
|
|
171
|
+
.mypy_cache/
|
|
172
|
+
.dmypy.json
|
|
173
|
+
dmypy.json
|
|
174
|
+
|
|
175
|
+
# Pyre type checker
|
|
176
|
+
.pyre/
|
|
177
|
+
|
|
178
|
+
# pytype static type analyzer
|
|
179
|
+
.pytype/
|
|
180
|
+
|
|
181
|
+
# Cython debug symbols
|
|
182
|
+
cython_debug/
|
|
183
|
+
|
|
184
|
+
# PyCharm
|
|
185
|
+
# JetBrains specific template is maintained in a separate JetBrains.gitignore that can
|
|
186
|
+
# be found at https://github.com/github/gitignore/blob/main/Global/JetBrains.gitignore
|
|
187
|
+
# and can be added to the global gitignore or merged into this file. For a more nuclear
|
|
188
|
+
# option (not recommended) you can uncomment the following to ignore the entire idea folder.
|
|
189
|
+
# .idea/
|
|
190
|
+
|
|
191
|
+
# Abstra
|
|
192
|
+
# Abstra is an AI-powered process automation framework.
|
|
193
|
+
# Ignore directories containing user credentials, local state, and settings.
|
|
194
|
+
# Learn more at https://abstra.io/docs
|
|
195
|
+
.abstra/
|
|
196
|
+
|
|
197
|
+
# Visual Studio Code
|
|
198
|
+
# Visual Studio Code specific template is maintained in a separate VisualStudioCode.gitignore
|
|
199
|
+
# that can be found at https://github.com/github/gitignore/blob/main/Global/VisualStudioCode.gitignore
|
|
200
|
+
# and can be added to the global gitignore or merged into this file. However, if you prefer,
|
|
201
|
+
# you could uncomment the following to ignore the entire vscode folder
|
|
202
|
+
# .vscode/
|
|
203
|
+
# Temporary file for partial code execution
|
|
204
|
+
tempCodeRunnerFile.py
|
|
205
|
+
|
|
206
|
+
# Ruff stuff:
|
|
207
|
+
.ruff_cache/
|
|
208
|
+
|
|
209
|
+
# PyPI configuration file
|
|
210
|
+
.pypirc
|
|
211
|
+
|
|
212
|
+
# Marimo
|
|
213
|
+
marimo/_static/
|
|
214
|
+
marimo/_lsp/
|
|
215
|
+
__marimo__/
|
|
216
|
+
|
|
217
|
+
# Streamlit
|
|
218
|
+
.streamlit/secrets.toml
|
|
219
|
+
|
|
220
|
+
# macOS
|
|
221
|
+
.DS_Store
|
pymhd-0.1.0/LICENSE
ADDED
|
@@ -0,0 +1,21 @@
|
|
|
1
|
+
MIT License
|
|
2
|
+
|
|
3
|
+
Copyright (c) 2026 Yuyang HUA
|
|
4
|
+
|
|
5
|
+
Permission is hereby granted, free of charge, to any person obtaining a copy
|
|
6
|
+
of this software and associated documentation files (the "Software"), to deal
|
|
7
|
+
in the Software without restriction, including without limitation the rights
|
|
8
|
+
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
|
9
|
+
copies of the Software, and to permit persons to whom the Software is
|
|
10
|
+
furnished to do so, subject to the following conditions:
|
|
11
|
+
|
|
12
|
+
The above copyright notice and this permission notice shall be included in all
|
|
13
|
+
copies or substantial portions of the Software.
|
|
14
|
+
|
|
15
|
+
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
|
16
|
+
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
17
|
+
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
|
18
|
+
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
|
19
|
+
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
|
20
|
+
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
|
|
21
|
+
SOFTWARE.
|
pymhd-0.1.0/PKG-INFO
ADDED
|
@@ -0,0 +1,100 @@
|
|
|
1
|
+
Metadata-Version: 2.4
|
|
2
|
+
Name: PyMHD
|
|
3
|
+
Version: 0.1.0
|
|
4
|
+
Summary: A Python package for post-processing MHD turbulence simulations
|
|
5
|
+
Project-URL: Homepage, https://github.com/PlasmaHua/PyMHD
|
|
6
|
+
Project-URL: Repository, https://github.com/PlasmaHua/PyMHD.git
|
|
7
|
+
Project-URL: Issues, https://github.com/PlasmaHua/PyMHD/issues
|
|
8
|
+
Author-email: "Yuyang Hua (华宇阳)" <yuyanghua@stu.pku.edu.cn>
|
|
9
|
+
License-Expression: MIT
|
|
10
|
+
License-File: LICENSE
|
|
11
|
+
Classifier: Operating System :: OS Independent
|
|
12
|
+
Classifier: Programming Language :: Python :: 3
|
|
13
|
+
Requires-Python: >=3.10
|
|
14
|
+
Requires-Dist: h5py>=3.10
|
|
15
|
+
Requires-Dist: jax>=0.4.26
|
|
16
|
+
Requires-Dist: kdepy>=1.1
|
|
17
|
+
Requires-Dist: matplotlib>=3.7
|
|
18
|
+
Requires-Dist: numpy>=1.26
|
|
19
|
+
Requires-Dist: yt>=4.2
|
|
20
|
+
Description-Content-Type: text/markdown
|
|
21
|
+
|
|
22
|
+
# PyMHD
|
|
23
|
+
|
|
24
|
+
[](https://doi.org/10.5281/zenodo.20749062)
|
|
25
|
+
[](https://github.com/PlasmaHua/PyMHD/blob/main/LICENSE)
|
|
26
|
+
|
|
27
|
+
|
|
28
|
+
A Python package for post-processing magnetohydrodynamic (MHD) turbulence simulations.
|
|
29
|
+
|
|
30
|
+
## Features
|
|
31
|
+
|
|
32
|
+
- [x] High-fidelity ***a posteriori* estimation of numerical dissipation** in MHD turbulence simulations;
|
|
33
|
+
- [x] Computation and visualization of turbulent energy spectra;
|
|
34
|
+
- [x] Visualization of 2D slices of MHD variables;
|
|
35
|
+
- [x] Built-in support for simulations from the Athena family, including [Athena++](https://github.com/PrincetonUniversity/athena), [AthenaK](https://github.com/IAS-Astrophysics/athenak), and [AthenaPK](https://github.com/parthenon-hpc-lab/athenapk);
|
|
36
|
+
- [x] Support for **driven MHD turbulence** (e.g., Alfvénic turbulence and small-scale dynamos), driven hydrodynamic turbulence, and **MRI-driven turbulence** (from shearing-box simulations) with triply periodic boundary conditions;
|
|
37
|
+
- [ ] Support for adiabatic and isothermal equations of state (EoS); **TODO:** add support for incompressible EoS;
|
|
38
|
+
- [ ] **TODO:** Computation and visualization of correlation functions;
|
|
39
|
+
- [ ] **TODO:** Spectral energy transfer analysis following [Grete, Philipp, et al. Physics of Plasmas 24.9 (2017)](https://doi.org/10.1063/1.4990613).
|
|
40
|
+
- [ ] **TODO:** Support JAX acceleration on GPUs and multi-node CPUs (currently, certain algorithms support parallelism on single-node, multi-core CPUs).
|
|
41
|
+
|
|
42
|
+
|
|
43
|
+
## Quick Start
|
|
44
|
+
|
|
45
|
+
### Installation
|
|
46
|
+
|
|
47
|
+
Install **PyMHD** from [PyPI](https://pypi.org/) via:
|
|
48
|
+
|
|
49
|
+
```bash
|
|
50
|
+
pip install pymhd
|
|
51
|
+
```
|
|
52
|
+
|
|
53
|
+
This will automatically install PyMHD and its required dependencies, including:
|
|
54
|
+
- [NumPy](https://github.com/numpy/numpy) for core computation;
|
|
55
|
+
- [JAX](https://github.com/jax-ml/jax) for parallel acceleration on single-node, multi-core CPUs; therefore, [the CPU version of JAX](https://docs.jax.dev/en/latest/installation.html#pip-installation-cpu) is required;
|
|
56
|
+
- [yt](https://github.com/yt-project/yt) and [h5py](https://github.com/h5py/h5py) for extracting data from HDF5 output files, specifically Athena++ `.athdf` files, AthenaK `.bin` files, and AthenaPK `.phdf` files;
|
|
57
|
+
- [Matplotlib](https://github.com/matplotlib/matplotlib) for plotting functionality in the `pymhd.plot` modules;
|
|
58
|
+
- [KDEpy](https://github.com/tommyod/KDEpy) for plotting smoothed histograms of numerical dissipation data in the `pymhd.plot.nd` module.
|
|
59
|
+
|
|
60
|
+
### Built-in Examples
|
|
61
|
+
|
|
62
|
+
For a quick start, clone the source repository, for example:
|
|
63
|
+
|
|
64
|
+
```bash
|
|
65
|
+
git clone https://github.com/PlasmaHua/PyMHD.git
|
|
66
|
+
```
|
|
67
|
+
|
|
68
|
+
Then navigate to the `./examples` subdirectory using `cd PyMHD/examples`, where two AthenaK simulation examples are provided.
|
|
69
|
+
|
|
70
|
+
Under construction.
|
|
71
|
+
|
|
72
|
+
### Analyze Your Own Athena Simulations
|
|
73
|
+
|
|
74
|
+
PyMHD natively supports simulations from Athena++, AthenaK, and AthenaPK. Use the `pymhd.output2turbulence` function to extract output data and input-file metadata, then construct a PyMHD `Turbulence` object for further analysis.
|
|
75
|
+
|
|
76
|
+
## API Documentation
|
|
77
|
+
|
|
78
|
+
Under construction.
|
|
79
|
+
|
|
80
|
+
## Citation
|
|
81
|
+
|
|
82
|
+
If PyMHD contributes to your work in any way, please consider citing the package via:
|
|
83
|
+
|
|
84
|
+
```bibtex
|
|
85
|
+
@misc{Hua2026,
|
|
86
|
+
doi = {10.5281/zenodo.20749062},
|
|
87
|
+
url = {https://zenodo.org/doi/10.5281/zenodo.20749062},
|
|
88
|
+
author = {Hua, Yuyang},
|
|
89
|
+
title = {PyMHD: Python package for post-processing MHD turbulence simulations},
|
|
90
|
+
publisher = {Zenodo},
|
|
91
|
+
year = {2026},
|
|
92
|
+
copyright = {MIT License}
|
|
93
|
+
}
|
|
94
|
+
```
|
|
95
|
+
|
|
96
|
+
If you specifically use the framework for numerical dissipation estimation, please cite our method paper via:
|
|
97
|
+
|
|
98
|
+
```bibtex
|
|
99
|
+
|
|
100
|
+
```
|
pymhd-0.1.0/README.md
ADDED
|
@@ -0,0 +1,79 @@
|
|
|
1
|
+
# PyMHD
|
|
2
|
+
|
|
3
|
+
[](https://doi.org/10.5281/zenodo.20749062)
|
|
4
|
+
[](https://github.com/PlasmaHua/PyMHD/blob/main/LICENSE)
|
|
5
|
+
|
|
6
|
+
|
|
7
|
+
A Python package for post-processing magnetohydrodynamic (MHD) turbulence simulations.
|
|
8
|
+
|
|
9
|
+
## Features
|
|
10
|
+
|
|
11
|
+
- [x] High-fidelity ***a posteriori* estimation of numerical dissipation** in MHD turbulence simulations;
|
|
12
|
+
- [x] Computation and visualization of turbulent energy spectra;
|
|
13
|
+
- [x] Visualization of 2D slices of MHD variables;
|
|
14
|
+
- [x] Built-in support for simulations from the Athena family, including [Athena++](https://github.com/PrincetonUniversity/athena), [AthenaK](https://github.com/IAS-Astrophysics/athenak), and [AthenaPK](https://github.com/parthenon-hpc-lab/athenapk);
|
|
15
|
+
- [x] Support for **driven MHD turbulence** (e.g., Alfvénic turbulence and small-scale dynamos), driven hydrodynamic turbulence, and **MRI-driven turbulence** (from shearing-box simulations) with triply periodic boundary conditions;
|
|
16
|
+
- [ ] Support for adiabatic and isothermal equations of state (EoS); **TODO:** add support for incompressible EoS;
|
|
17
|
+
- [ ] **TODO:** Computation and visualization of correlation functions;
|
|
18
|
+
- [ ] **TODO:** Spectral energy transfer analysis following [Grete, Philipp, et al. Physics of Plasmas 24.9 (2017)](https://doi.org/10.1063/1.4990613).
|
|
19
|
+
- [ ] **TODO:** Support JAX acceleration on GPUs and multi-node CPUs (currently, certain algorithms support parallelism on single-node, multi-core CPUs).
|
|
20
|
+
|
|
21
|
+
|
|
22
|
+
## Quick Start
|
|
23
|
+
|
|
24
|
+
### Installation
|
|
25
|
+
|
|
26
|
+
Install **PyMHD** from [PyPI](https://pypi.org/) via:
|
|
27
|
+
|
|
28
|
+
```bash
|
|
29
|
+
pip install pymhd
|
|
30
|
+
```
|
|
31
|
+
|
|
32
|
+
This will automatically install PyMHD and its required dependencies, including:
|
|
33
|
+
- [NumPy](https://github.com/numpy/numpy) for core computation;
|
|
34
|
+
- [JAX](https://github.com/jax-ml/jax) for parallel acceleration on single-node, multi-core CPUs; therefore, [the CPU version of JAX](https://docs.jax.dev/en/latest/installation.html#pip-installation-cpu) is required;
|
|
35
|
+
- [yt](https://github.com/yt-project/yt) and [h5py](https://github.com/h5py/h5py) for extracting data from HDF5 output files, specifically Athena++ `.athdf` files, AthenaK `.bin` files, and AthenaPK `.phdf` files;
|
|
36
|
+
- [Matplotlib](https://github.com/matplotlib/matplotlib) for plotting functionality in the `pymhd.plot` modules;
|
|
37
|
+
- [KDEpy](https://github.com/tommyod/KDEpy) for plotting smoothed histograms of numerical dissipation data in the `pymhd.plot.nd` module.
|
|
38
|
+
|
|
39
|
+
### Built-in Examples
|
|
40
|
+
|
|
41
|
+
For a quick start, clone the source repository, for example:
|
|
42
|
+
|
|
43
|
+
```bash
|
|
44
|
+
git clone https://github.com/PlasmaHua/PyMHD.git
|
|
45
|
+
```
|
|
46
|
+
|
|
47
|
+
Then navigate to the `./examples` subdirectory using `cd PyMHD/examples`, where two AthenaK simulation examples are provided.
|
|
48
|
+
|
|
49
|
+
Under construction.
|
|
50
|
+
|
|
51
|
+
### Analyze Your Own Athena Simulations
|
|
52
|
+
|
|
53
|
+
PyMHD natively supports simulations from Athena++, AthenaK, and AthenaPK. Use the `pymhd.output2turbulence` function to extract output data and input-file metadata, then construct a PyMHD `Turbulence` object for further analysis.
|
|
54
|
+
|
|
55
|
+
## API Documentation
|
|
56
|
+
|
|
57
|
+
Under construction.
|
|
58
|
+
|
|
59
|
+
## Citation
|
|
60
|
+
|
|
61
|
+
If PyMHD contributes to your work in any way, please consider citing the package via:
|
|
62
|
+
|
|
63
|
+
```bibtex
|
|
64
|
+
@misc{Hua2026,
|
|
65
|
+
doi = {10.5281/zenodo.20749062},
|
|
66
|
+
url = {https://zenodo.org/doi/10.5281/zenodo.20749062},
|
|
67
|
+
author = {Hua, Yuyang},
|
|
68
|
+
title = {PyMHD: Python package for post-processing MHD turbulence simulations},
|
|
69
|
+
publisher = {Zenodo},
|
|
70
|
+
year = {2026},
|
|
71
|
+
copyright = {MIT License}
|
|
72
|
+
}
|
|
73
|
+
```
|
|
74
|
+
|
|
75
|
+
If you specifically use the framework for numerical dissipation estimation, please cite our method paper via:
|
|
76
|
+
|
|
77
|
+
```bibtex
|
|
78
|
+
|
|
79
|
+
```
|
|
@@ -0,0 +1,31 @@
|
|
|
1
|
+
# PyMHD: Python for Magnetohydrodynamic Turbulence.
|
|
2
|
+
# Copyright (c) 2026 Yuyang Hua (华宇阳)
|
|
3
|
+
# License: MIT
|
|
4
|
+
|
|
5
|
+
"""
|
|
6
|
+
pymhd/__init__.py
|
|
7
|
+
-----------------
|
|
8
|
+
|
|
9
|
+
Unified interface for PyMHD.
|
|
10
|
+
"""
|
|
11
|
+
|
|
12
|
+
from .turbulence import Scalar, Vector, ScalarField, VectorField, Turbulence, avg, std
|
|
13
|
+
from .derivatives.derivative import Dx, Dy, Dz, grad, div, curl, laplacian, Algorithm
|
|
14
|
+
|
|
15
|
+
from .spectra import Spectrum, Spectrum1D, EnergySpectra
|
|
16
|
+
from .numdiss import NumericalDissipation
|
|
17
|
+
from .preprocess import output2turbulence
|
|
18
|
+
|
|
19
|
+
from .plot import plot
|
|
20
|
+
|
|
21
|
+
__all__ = [
|
|
22
|
+
"Scalar", "Vector",
|
|
23
|
+
"ScalarField", "VectorField", "avg", "std", "Turbulence",
|
|
24
|
+
"Algorithm",
|
|
25
|
+
"Dx", "Dy", "Dz", "grad", "div", "curl", "laplacian",
|
|
26
|
+
"Spectrum", "Spectrum1D",
|
|
27
|
+
"EnergySpectra",
|
|
28
|
+
"NumericalDissipation",
|
|
29
|
+
"output2turbulence",
|
|
30
|
+
"plot",
|
|
31
|
+
]
|