Python implementation of the Gauss Separation method in Cartesian coordinates.
Fourier implementation of the Gauss Separation method. The Fourier implementation of this method is described in 1. The method assumes data are on a Cartesian mesh with uniform spacing.
The module is called gauss_separation, the primary method is just called fourier.
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.
Footnotes
-
Welsch,B.T.: 2022, Photospheric Imprints of Coronal Electric Currents, I. Magnetic Structure Near Polarity Inversion Lines,Arxiv,DOI:10.48550/arXiv.2211.01911 ↩