Gauss Separation method

Python implementation of the Gauss Separation method in Cartesian coordinates.

About

Fourier implementation of the Gauss Separation method. The Fourier implementation of this method is described in ¹. The method assumes data are on a Cartesian mesh with uniform spacing.

Usage

The module is called gauss_separation, the primary method is just called fourier.

Example/test

The file ut_gauss_separation_lfff.py applies the method to a simple Linear Force-Free magnetic Field (LFFFF) test case. The LFFF is simple enough that the decomposition can be performed analytically and provides a test case to demonstrate the correct functioning of the code. The code should reproduce the test case to approximately machine precision.

See the test case for an example of how to use the code.

Publications

A version of this code is used in Iyer et al.: 2026, ApJ, 1004, 112

• Paper DOI: 10.3847/1538-4357/ae6a8b
• Zenodo DOI: https://doi.org/10.5281/zenodo.17706429

References

Footnotes

1. Welsch,B.T.: 2022, Photospheric Imprints of Coronal Electric Currents, I. Magnetic Structure Near Polarity Inversion Lines,Arxiv,DOI:10.48550/arXiv.2211.01911 ↩