A full-fledged micromagnetic code in fewer than 70 lines of NumPy

Abstract

We present a complete micromagnetic finite-difference code in fewer than 70 lines of Python. The code makes a large use of the NumPy library and computes the exchange field by finite differences and the demagnetization field with a fast convolution algorithm. Since the magnetization in finite-difference micromagnetics is represented by a multi-dimensional array and the NumPy library features a rich interface for this data structure, the code we present is an ideal starting point for the development of novel algorithms.