SallyLLG Data analysis provides a series of different function to compute the hysteresis cycle from the data obtained with the micromagnetic code SallyLLG in the .h5 format. It can be used both with the results provided by the 2D 2.5D and 3D code.
Copy the .py files from this repository in the working folder
Three main scripts can be found:
hdf5_test.py, which contains the routines and classes to compute the hysteresis cycle of from the 2.5d and 2d simulations;
hdf5_test_3Dv2.py, which contains the routines to compute the hysteresis cycle of from the 3d simulations;
distributions4_corrections.py, which allows to plot the data obtained from the previous scripts;
2d3d_Compute_Hyst_SpecificEnergy.py, Example script puts together hdf5_test class and integral class and uses them to compute the Hysteresis loop and specific energy from the loop.
""" Basic script Open a HDF5 file of my simulation and compute the Hysteresis loop, saves data in an opportune file with specific naming pattern at the end plots the calculated loop. Is a hdf5_test class, initialized with mainDir filename, saves values such as output file names and things like that. """
class hdf5_test:
def __init__(self,mainDir,filename,numObj=1,versoreu = np.array([[1],[0],[0]]),versorev = np.array([[1],[0],[0]]),versorew = np.array([[1],[0],[0]]))has now two methods compute and plot to calculate hysteresis loop and plot it
calcoloMagnMedia(time,file,Volumes,versoreu,versorev,versorew)Return a the average magnetization for each time step / equilibrium point of the simulation.
calcoloMagnMediaDisks(time,file,Volumes,numObj,versoreu,versorev,versorew)Return a the average magnetization for each time step / equilibrium point of the simulation, for each object in a 2.5D distribution of N objects.
calcoloEnergia(time,file,Volumes):Computes the Zeeman and magnetostatic energy of the system for each time step / equilibrium point.
""" Basic script Open a HDF5 file of my simulation and compute the Hysteresis loop, saves data in an opportune file with specific naming pattern at the end plots the calculated loop. """
MagnetizationCalc(mainDir,filename):Sets up the class and loads the filename.h5 file from the folder in mainDir.
computeData(self, hysteresis=True, dt=1, versoreu = np.array([[1],[0],[0]]), versorev = np.array([[0],[1],[0]]), versorew = np.array([[0],[0],[1]])):Evaluate the Hysteresis loop or the average magnetization temporal behaviour along the direction specified by the versors versoreu, versorev, verosew.
calcoloMagnMediaVsappField(time, file, versoreu, versorev, versorew):Used in computeData to evaluate the Hysteresis loop.
calcoloMagnMedia(time, file, versoreu, versorev, versorew):Used in computeData to evaluate the temporal behaviour of the magnetization