osiris-utils 1.0.5__py3-none-any.whl → 1.0.9__py3-none-any.whl

osiris_utils/__init__.py

@@ -1,5 +1,6 @@
 from .utils import (time_estimation, filesize_estimation, transverse_average, integrate, animate_2D,
-                    save_data, read_data, mft_decomposition, courant2D)
+                    save_data, read_data, courant2D)
 from .data_readers import open1D, open2D, open3D, read_osiris_file
 from .gui import LAVA_Qt, LAVA
-from .data import OsirisGridFile, OsirisRawFile
+from .data import OsirisGridFile, OsirisRawFile, OsirisData
+from .mean_field_theory import MeanFieldTheory

osiris_utils/data.py

@@ -1,158 +1,477 @@
 import numpy as np
 import h5py
 
-class OsirisGridFile():
+class OsirisData():
+    '''
+    Parent class for Osiris data files.
+
+    Input:
+    ------
+        - filename: the path to the HDF5 file
+    Attributes:
+    ----------
+        - filename - the path to the HDF5 file
+            str
+        - file - the HDF5 file
+            h5py.File
+        - verbose - if True, the class will print messages
+            bool
+        - dt - the time step
+            float
+        - dim - the number of dimensions
+            int
+        - time - the time and its units
+            list [time, units]
+            list [float, str]
+        - iter - the iteration number
+            int
+        - name - the name of the data
+            str
+        - type - the type of data
+            str
+    '''
+    def __init__(self, filename):
+        self._filename = filename
+        self._file = None
+
+        self._verbose = False
+
+        self._open_file(self._filename)
+        self._load_basic_attributes(self._file)
+        
+    def _load_basic_attributes(self, f: h5py.File) -> None:
+        '''Load common attributes from HDF5 file'''
+        self._dt = float(f['SIMULATION'].attrs['DT'][0])
+        self._dim = int(f['SIMULATION'].attrs['NDIMS'][0])
+        self._time = [float(f.attrs['TIME'][0]), f.attrs['TIME UNITS'][0].decode('utf-8')]
+        self._iter = int(f.attrs['ITER'][0])
+        self._name = f.attrs['NAME'][0].decode('utf-8')
+        self._type = f.attrs['TYPE'][0].decode('utf-8')
+    
+    def verbose(self, verbose: bool = True):
+        '''
+        Set the verbosity of the class
+
+        Parameters
+        ----------
+        verbose : bool, optional
+            If True, the class will print messages, by default True when calling (False when not calling)
+        '''
+        self._verbose = verbose
+
+    def _open_file(self, filename):
+        '''
+        Open the HDF5 file.
+
+        Parameters
+        ----------
+        filename : str
+            The path to the HDF5 file.
+        '''
+        if self._verbose: print(f'Opening file > {filename}')
+        self._file = h5py.File(self._filename, 'r')
+
+    def _close_file(self):
+        '''
+        Close the HDF5 file.
+        '''
+        if self._verbose: print('Closing file')
+        if self._file:
+            self._file.close()
+        
+    @property
+    def dt(self):
+        '''
+        Returns
+        -------
+        float
+            The time step
+        '''
+        return self._dt
+    @property
+    def dim(self):
+        '''
+        Returns
+        -------
+        int
+            The number of dimensions
+        '''
+        return self._dim
+    @property
+    def time(self):
+        '''
+        Returns
+        -------
+        list
+            The time and its units
+        '''
+        return self._time
+    @property
+    def iter(self):
+        '''
+        Returns
+        -------
+        int
+            The iteration number
+        '''
+        return self._iter
+    @property
+    def name(self):
+        '''
+        Returns
+        -------
+        str
+            The name of the data
+        '''
+        return self._name
+    @property
+    def type(self):
+        '''
+        Returns
+        -------
+        str
+            The type of data
+        '''
+        return self._type
+
+class OsirisGridFile(OsirisData):
     '''
     Class to read the grid data from an OSIRIS HDF5 file.
     
     Input:
+    -----
         - filename: the path to the HDF5 file
         
     Attributes:
-        - grid: the grid data ((x1.min, x1.max), (x2.min, x2.max), (x3.min, x3.max))
-            numpy.ndarray
-        - nx: the number of grid points (nx1, nx2, nx3)
-            numpy.ndarray
-        - dx: the grid spacing (dx1, dx2, dx3)
-            numpy.ndarray
-        - axis: the axis data [(name_x1, units_x1, long_name_x1, type_x1), ...]
-            list of dictionaries
-        - data: the data (numpy array) with shape (nx1, nx2, nx3) (Transpose to use `plt.imshow`)
-            numpy.ndarray
-        - dt: the time step
+    ----------
+        - filename - the path to the HDF5 file
+            str
+        - file - the HDF5 file
+            h5py.File
+        - verbose - if True, the class will print messages
+            bool
+        - dt - the time step
             float
-        - dim: the number of dimensions
+        - dim - the number of dimensions
             int
-        - time: the time and its units
+        - time - the time and its units
             list [time, units]
             list [float, str]
-        - iter: the iteration number
+        - iter - the iteration number
             int
-        - name: the name of the data
+        - name - the name of the data
+            str
+        - type - the type of data
             str
-        - units: the units of the data
+        - grid - the grid data ((x1.min, x1.max), (x2.min, x2.max), (x3.min, x3.max))
+            numpy.ndarray
+        - nx - the number of grid points (nx1, nx2, nx3)
+            numpy.ndarray
+        - dx - the grid spacing (dx1, dx2, dx3)
+            numpy.ndarray
+        - axis - the axis data [(name_x1, units_x1, long_name_x1, type_x1), ...]
+            list of dictionaries
+        - data: the data (numpy array) with shape (nx1, nx2, nx3) (Transpose to use `plt.imshow`)
+            numpy.ndarray
+        - units - the units of the data
             str
-        - label: the label of the data (LaTeX formatted)
+        - label - the label of the data (LaTeX formatted)
         
-    Example:
-        example = utils.OsirisGridFile("grid/file/osiris/output-000000.h5")
-        plt.imshow(example.data.T, aspect='auto', origin='lower', extent=[example.grid[0][0], example.grid[0][-1], example.grid[1][0], example.grid[1][-1]])
-        plt.xlabel(rf"${example.axis[0]["long_name"]}$ [${example.axis[0]["units"]}$]")
-        plt.ylabel(rf"${example.axis[1]["long_name"]}$ [${example.axis[1]["units"]}$]")
-        plt.title(rf"${example.label}$, t={example.time[0]} ${example.time[1]}$")
     '''
     def __init__(self, filename):
-        with h5py.File(filename, 'r+') as f:
-            # Get the data 
-            known_keys = {"AXIS", "SIMULATION"}
-            all_keys = set(f.keys())
-            variable_key = (all_keys - known_keys).pop()
+        super().__init__(filename)
+            
+        variable_key = self._get_variable_key(self._file)
+        
+        self._units = self._file.attrs['UNITS'][0].decode('utf-8')
+        self._label = self._file.attrs['LABEL'][0].decode('utf-8')
+        self._FFTdata = None
+        
+        data = np.array(self._file[variable_key][:])
+
+        axis = list(self._file['AXIS'].keys())
+        if len(axis) == 1:
+            self._grid = self._file['AXIS/' + axis[0]][()]
+            self._nx = len(data)
+            self._dx = (self.grid[1] - self.grid[0] ) / self.nx
+            self._x = np.arange(self.grid[0], self.grid[1], self.dx)
+        else: 
+            grid = []
+            for ax in axis: grid.append(self._file['AXIS/' + ax][()])
+            self._grid = np.array(grid)
+            self._nx = self._file[variable_key][()].transpose().shape
+            self._dx = (self.grid[:, 1] - self.grid[:, 0])/self.nx
+            self._x = [np.arange(self.grid[i, 0], self.grid[i, 1], self.dx[i]) for i in range(self.dim)]
+
+        self._axis = []
+        for ax in axis:
+            axis_data = {
+                'name': self._file['AXIS/'+ax].attrs['NAME'][0].decode('utf-8'),
+                'units': self._file['AXIS/'+ax].attrs['UNITS'][0].decode('utf-8'),
+                'long_name': self._file['AXIS/'+ax].attrs['LONG_NAME'][0].decode('utf-8'),
+                'type': self._file['AXIS/'+ax].attrs['TYPE'][0].decode('utf-8'),
+                'plot_label': rf'${self._file["AXIS/"+ax].attrs["LONG_NAME"][0].decode("utf-8")}$ $[{self._file["AXIS/"+ax].attrs["UNITS"][0].decode("utf-8")}]$',
+            }
+            self._axis.append(axis_data)
+        
+        self._data = np.ascontiguousarray(data.T)
+
+        self._close_file()
+
+    def _load_basic_attributes(self, f: h5py.File) -> None:
+        '''Load common attributes from HDF5 file'''
+        self._dt = float(f['SIMULATION'].attrs['DT'][0])
+        self._dim = int(f['SIMULATION'].attrs['NDIMS'][0])
+        self._time = [float(f.attrs['TIME'][0]), f.attrs['TIME UNITS'][0].decode('utf-8')]
+        self._iter = int(f.attrs['ITER'][0])
+        self._name = f.attrs['NAME'][0].decode('utf-8')
+        self._type = f.attrs['TYPE'][0].decode('utf-8')
             
-            # The data
-            self.data = np.array(f[variable_key][:])
-
-            keys = list(f.keys())
-            # Now get the infos
-            axis = list(f["AXIS"].keys())
-            if len(axis) == 1:
-                self.grid = f["AXIS/"+axis[0]][()]
-                self.nx = len(self.data)
-                self.dx = (self.grid[1] - self.grid[0] ) / self.nx
-            else: 
-                grid = []
-                for ax in axis: grid.append(f["AXIS/"+ax][()])
-                self.grid = np.array(grid)
-                
-                self.nx = f[variable_key][()].transpose().shape
-                self.dx = (self.grid[:, 1] - self.grid[:, 0])/self.nx
-            self.axis = []
-            for ax in axis:
-                axis_data = {
-                    "name": f["AXIS/"+ax].attrs["NAME"][0].decode('utf-8'),
-                    "units": f["AXIS/"+ax].attrs["UNITS"][0].decode('utf-8'),
-                    "long_name": f["AXIS/"+ax].attrs["LONG_NAME"][0].decode('utf-8'),
-                    "type": f["AXIS/"+ax].attrs["TYPE"][0].decode('utf-8'),
-                }
-                self.axis.append(axis_data)
-                    
-            # NOW WORK ON THE SIMULATION DATA
-            self.dt = float(f["SIMULATION"].attrs["DT"][0])
-            self.dims_simulation = int(f["SIMULATION"].attrs["NDIMS"][0])
-            self.time = [float(f.attrs["TIME"][0]), f.attrs["TIME UNITS"][0].decode('utf-8')]
-            self.iter = int(f.attrs["ITER"][0])
-            self.name = f.attrs["NAME"][0].decode('utf-8')
-            self.units = f.attrs["UNITS"][0].decode('utf-8')
-            self.label = f.attrs["LABEL"][0].decode('utf-8')
-            self.type = f.attrs["TYPE"][0].decode('utf-8')
-            self.data = self.data.T
-
-
-class OsirisRawFile():
+    def _get_variable_key(self, f: h5py.File) -> str:
+        return next(k for k in f.keys() if k not in {'AXIS', 'SIMULATION'})
+    
+    # Getters
+    @property
+    def grid(self):
+        '''
+        Returns
+        -------
+        numpy.ndarray
+            The grid data ((x1.min, x1.max), (x2.min, x2.max), (x3.min, x3.max
+        '''
+        return self._grid
+    @property
+    def nx(self):
+        '''
+        Returns
+        -------
+        numpy.ndarray
+            The number of grid points (nx1, nx2, nx3)
+        '''
+        return self._nx
+    @property
+    def dx(self):
+        '''
+        Returns
+        -------
+        numpy.ndarray
+            The grid spacing (dx1, dx2, dx3)
+        '''
+        return self._dx
+    @property
+    def x(self):
+        '''
+        Returns
+        -------
+        numpy.ndarray
+            The grid points in each axis
+        '''
+        return self._x
+    @property
+    def axis(self):
+        '''
+        Returns
+        -------
+        list of dictionaries
+            The axis data [(name_x1, units_x1, long_name_x1, type_x1), ...]
+        '''
+        return self._axis   
+    @property
+    def data(self):
+        '''
+        Returns
+        -------
+        numpy.ndarray
+            The data (numpy array) with shape (nx1, nx2, nx3) (Transpose to use `plt.imshow`)
+        '''
+        return self._data
+    @property
+    def units(self):
+        '''
+        Returns
+        -------
+        str
+            The units of the data (LaTeX formatted)
+        '''
+        return self._units
+    @property
+    def label(self):
+        '''
+        Returns
+        -------
+        str
+            The label of the data (LaTeX formatted)
+        '''
+        return self._label
+    @property
+    def FFTdata(self):
+        '''
+        Returns
+        -------
+        numpy.ndarray
+            The FFT of the data
+        '''
+        if self._FFTdata is None:
+            raise ValueError('The FFT of the data has not been computed yet. Compute it using the FFT method.')
+        return self._FFTdata
+    # Setters
+    @data.setter
+    def data(self, data):
+        '''
+        Set the data attribute
+        '''
+        self._data = data
+
+    def __str__(self):
+        # write me a template to print with the name, label, units, time, iter, grid, nx, dx, axis, dt, dim in a logical way
+        return rf'{self.name}' + f'\n' + rf'Time: [{self.time[0]} {self.time[1]}], dt = {self.dt}' + f'\n' + f'Iteration: {self.iter}' + f'\n' + f'Grid: {self.grid}' + f'\n' + f'dx: {self.dx}' + f'\n' + f'Dimensions: {self.dim}D'
+    
+
+    def __array__(self):
+        return np.asarray(self.data)
+    
+
+    def _yeeToCellCorner1d(self, boundary):
+        '''
+        Converts 1d EM fields from a staggered Yee mesh to a grid with field values centered on the corner of the cell (the corner of the cell [1] has coordinates [1])
+        '''
+
+        if self.name.lower() in ['b2', 'b3', 'e1']:
+            if boundary == 'periodic': return 0.5 * (np.roll(self.data, shift=1) + self.data) 
+            else: return 0.5 * (self.data[1:] + self.data[:-1])
+        elif self.name.lower() in ['b1', 'e2', 'e3']:
+            if boundary == 'periodic': return self.data 
+            else: return  self.data[1:]
+        else: 
+            raise TypeError(f'This method expects magnetic or electric field grid data but received \'{self.name}\' instead')
+    
+
+    def _yeeToCellCorner2d(self, boundary):
+        '''
+        Converts 2d EM fields from a staggered Yee mesh to a grid with field values centered on the corner of the cell (the corner of the cell [1,1] has coordinates [1,1])
+        '''
+
+        if self.name.lower() in ['e1', 'b2']:
+            if boundary == 'periodic': return 0.5 * (np.roll(self.data, shift=1, axis=0) + self.data)
+            else: return 0.5 * (self.data[1:, 1:] + self.data[:-1, 1:])
+        elif self.name.lower() in ['e2', 'b1']:
+            if boundary == 'periodic': return 0.5 * (np.roll(self.data, shift=1, axis=1) + self.data)
+            else: return 0.5 * (self.data[1:, 1:] + self.data[1:, :-1])
+        elif self.name.lower() in ['b3']:
+            if boundary == 'periodic': 
+               return 0.5 * (np.roll((0.5 * (np.roll(self.data, shift=1, axis=0) + self.data)), shift=1, axis=1) + (0.5 * (np.roll(self.data, shift=1, axis=0) + self.data)))
+            else:
+                return 0.25 * (self.data[1:, 1:] + self.data[:-1, 1:] + self.data[1:, :-1] + self.data[:-1, :-1])
+        elif self.name.lower() in ['e3']:
+            if boundary == 'periodic': return self.data
+            else: return self.data[1:, 1:]
+        else:
+            raise TypeError(f'This method expects magnetic or electric field grid data but received \'{self.name}\' instead')
+        
+
+    def _yeeToCellCorner3d(self, boundary):
+        '''
+        Converts 3d EM fields from a staggered Yee mesh to a grid with field values centered on the corner of the cell (the corner of the cell [1,1,1] has coordinates [1,1,1])
+        '''
+        if boundary == 'periodic':
+            raise ValueError('Centering field from 3D simulations considering periodic boundary conditions is not implemented yet')
+        if self.name.lower() == 'b1':
+            return 0.25 * (self.data[1:, 1:, 1:] + self.data[1:, :-1, 1:] + self.data[1:, 1:, :-1] + self.data[1:, :-1, :-1])
+        elif self.name.lower() == 'b2':
+            return 0.25 * (self.data[1:, 1:, 1:] + self.data[:-1, 1:, 1:] + self.data[1:, 1:, :-1] + self.data[:-1, 1:, :-1])
+        elif self.name.lower() == 'b3':
+            return 0.25 * (self.data[1:, 1:, 1:] + self.data[:-1, 1:, 1:] + self.data[1:, :-1, 1:] + self.data[:-1, :-1, 1:])
+        elif self.name.lower() == 'e1':
+            return 0.5 * (self.data[1:, 1:, 1:] + self.data[:-1, 1:, 1:])
+        elif self.name.lower() == 'e2':
+            return 0.5 * (self.data[1:, 1:, 1:] + self.data[1:, :-1, 1:])
+        elif self.name.lower() == 'e3':
+            return 0.5 * (self.data[1:, 1:, 1:] + self.data[1:, 1:, :-1])
+        else:
+            raise TypeError(f'This method expects magnetic or electric field grid data but received \'{self.name}\' instead')
+        
+    def yeeToCellCorner(self, boundary=None):
+        ''''
+        Converts EM fields from a staggered Yee mesh to a grid with field values centered on the corner of the cell.'
+        Can be used for 1D, 2D and 3D simulations.'
+        Creates a new attribute `data_centered` with the centered data.'
+        '''
+        
+        cases = {'b1', 'b2', 'b3', 'e1', 'e2', 'e3'}
+        if self.name not in cases:
+            raise TypeError(f'This method expects magnetic or electric field grid data but received \'{self.name}\' instead')
+        
+        if self.dim == 1:
+            self.data_centered = self._yeeToCellCorner1d(boundary)
+            return self.data_centered
+        elif self.dim == 2:
+            self.data_centered = self._yeeToCellCorner2d(boundary)
+            return self.data_centered
+        elif self.dim == 3:
+            self.data_centered = self._yeeToCellCorner3d(boundary)
+            return self.data_centered
+        else:
+            raise ValueError(f'Dimension {self.dim} is not supported')
+        
+    def FFT(self, axis=(0, )):
+        '''
+        Computes the Fast Fourier Transform of the data along the specified axis and shifts the zero frequency to the center.
+        Transforms the data to the frequency domain. A(x, y, z) -> A(kx, ky, kz)
+        '''
+        datafft = np.fft.fftn(self.data, axes=axis)
+        self._FFTdata = np.fft.fftshift(datafft, axes=axis)
+
+
+
+class OsirisRawFile(OsirisData):
     '''
     Class to read the raw data from an OSIRIS HDF5 file.
     
     Input:
         - filename: the path to the HDF5 file
-        
+    
     Attributes:
-
-        - axis: a dictionary where each key is a dataset name, and each value is another dictionary containing:
-                name (str): The name of the quantity (e.g., 'x1', 'ene').
-                units (str): The units associated with that dataset in LaTeX (e.g., 'c/\\omega_p', 'm_e c^2').
-                long_name (str): The name of the quantity in LaTeX (e.g., 'x_1', 'En2').
+        - axis - a dictionary where each key is a dataset name, and each value is another dictionary containing
+            name (str): The name of the quantity (e.g., r'x1', r'ene').
+            units (str): The units associated with that dataset in LaTeX (e.g., r'c/\\omega_p', r'm_e c^2').
+            long_name (str): The name of the quantity in LaTeX (e.g., r'x_1', r'En2').
             dictionary of dictionaries
-        - data: a dictionary where each key is a dataset name, and each value is the data
+        - data - a dictionary where each key is a dataset name, and each value is the data
             dictionary of np.arrays
-        - dim: the number of dimensions
+        - dim - the number of dimensions
             int
-        - dt: the time step
+        - dt - the time step
             float
-        - grid: maximum and minimum coordinates of the box, for each axis 
+        - grid - maximum and minimum coordinates of the box, for each axis 
             numpy.ndarray(dim,2)
-        - iter: the iteration number
+        - iter - the iteration number
             int
-        - name: the name of the species
+        - name - the name of the species
             str
-        - time: the time and its units
+        - time - the time and its units
             list [time, units]
             list [float, str]
-        -type: type of data (particles in the case of raw files) 
+        - type - type of data (particles in the case of raw files)
             str
 
-    Example:
-
-        electrons = utils.OsirisRawFile("MS/RAW/electrons/RAW-electrons-000002.h5")
-        plt.hist2d(electrons.data["x1"], electrons.data["x2"], bins=100, cmap='plasma', norm=mpl.colors.LogNorm())
-        plt.xlabel(r"${}$".format(electrons.axis["x1"]["long_name"]) + r"$\quad[{}]$".format(electrons.axis["x1"]["units"]))
-        plt.ylabel(r"${}$".format(electrons.axis["x2"]["long_name"]) + r"$\quad[{}]$".format(electrons.axis["x2"]["units"]))
-        plt.title(r"$t={}$".format(electrons.time[0]) + r"$\quad[{}]$".format(electrons.time[1]))
-        plt.show()
     '''
     
     def __init__(self, filename):
-        self.filename = filename
-
-        with h5py.File(self.filename, "r") as f: 
-
-            self.dt = float(f["SIMULATION"].attrs["DT"][0])
-            self.dim = int(f["SIMULATION"].attrs["NDIMS"][0])
-            self.time = [float(f.attrs["TIME"][0]), f.attrs["TIME UNITS"][0].decode('utf-8')]
-            self.iter = int(f.attrs["ITER"][0])
-            self.name = f.attrs["NAME"][0].decode('utf-8')
-            self.type = f.attrs["TYPE"][0].decode('utf-8')
-            self.grid = np.array([f["SIMULATION"].attrs["XMIN"], f["SIMULATION"].attrs["XMAX"]]).T
-
-            self.data = {}
-            self.axis = {}
-            for key in f.keys():
-                if key == "SIMULATION": continue
-
-                self.data[key] = np.array(f[key][()])
-
-                idx = np.where(f.attrs["QUANTS"] == str(key).encode('utf-8'))
-                axis_data = {
-                    "name": f.attrs["QUANTS"][idx][0].decode('utf-8'),
-                    "units": f.attrs["UNITS"][idx][0].decode('utf-8'),
-                    "long_name": f.attrs["LABELS"][idx][0].decode('utf-8'),
-                }
-                self.axis[key] = axis_data
+        super().__init__(filename)
+
+        self.grid = np.array([self._file['SIMULATION'].attrs['XMIN'], self._file['SIMULATION'].attrs['XMAX']]).T
+
+        self.data = {}
+        self.axis = {}
+        for key in self._file.keys():
+            if key == 'SIMULATION': continue
+
+            self.data[key] = np.array(self._file[key][()])
+
+            idx = np.where(self._file.attrs['QUANTS'] == str(key).encode('utf-8'))
+            axis_data = {
+                'name': self._file.attrs['QUANTS'][idx][0].decode('utf-8'),
+                'units': self._file.attrs['UNITS'][idx][0].decode('utf-8'),
+                'long_name': self._file.attrs['LABELS'][idx][0].decode('utf-8'),
+            }
+            self.axis[key] = axis_data

osiris_utils/data_readers.py

@@ -11,9 +11,9 @@ def read_osiris_file(filename, pressure = False):
     f = h5py.File(filename, 'r+')
     atr = f.attrs
     k = [key for key in f.keys()]
-    if "SIMULATION" in k:
+    if 'SIMULATION' in k:
         attrs1 = atr
-        attrs2 = f["SIMULATION"].attrs
+        attrs2 = f['SIMULATION'].attrs
         attrs = {}
         for i in range(len(attrs1)):
             attrs[[key for key in attrs1][i]] = [value for value in attrs1.values()][i]
@@ -24,21 +24,21 @@ def read_osiris_file(filename, pressure = False):
     axis = []
     for i in range(leanx):
         axis.append(ax.get([key for key in ax.keys()][i]))
-    if "SIMULATION" in k and pressure == False:
+    if 'SIMULATION' in k and pressure == False:
         data = f.get([key for key in f.keys()][2])
-        data.attrs["UNITS"] = attrs1["UNITS"]
-        data.attrs["LONG_NAME"] = attrs1["LABEL"]
-    elif "SIMULATION" in k and pressure == True:
+        data.attrs['UNITS'] = attrs1['UNITS']
+        data.attrs['LONG_NAME'] = attrs1['LABEL']
+    elif 'SIMULATION' in k and pressure == True:
         data = f.get([key for key in f.keys()][1])
-        data.attrs["UNITS"] = attrs1["UNITS"]
-        data.attrs["LONG_NAME"] = attrs1["LABEL"]
+        data.attrs['UNITS'] = attrs1['UNITS']
+        data.attrs['LONG_NAME'] = attrs1['LABEL']
     else:
         data = f.get([key for key in f.keys()][1])
     
     return attrs, axis, data
 
 def open1D(filename, pressure = False):
-    """ 
+    ''' 
     Open a 1D OSIRIS file and return the x axis and the data array.
 
     Parameters
@@ -52,7 +52,7 @@ def open1D(filename, pressure = False):
         The x axis.
     data_array : numpy.ndarray
         The data array.
-    """
+    '''
     attrs, axes, data = read_osiris_file(filename, pressure)
     datash = data.shape
     ax1 = axes[0]
@@ -61,7 +61,7 @@ def open1D(filename, pressure = False):
     return x, data_array, [attrs, axes, data]
 
 def open2D(filename, pressure = False):
-    """
+    '''
     Open a 2D OSIRIS file and return the x and y axes and the data array.
 
     Parameters
@@ -77,7 +77,7 @@ def open2D(filename, pressure = False):
         The y axis.
     data_array : numpy.ndarray
         The data array.
-    """
+    '''
     attrs, axes, data = read_osiris_file(filename, pressure)
     datash = data.shape
     ax1 = axes[0]
@@ -88,7 +88,7 @@ def open2D(filename, pressure = False):
     return x, y, data_array, [attrs, axes, data]
 
 def open3D(filename):
-    """
+    '''
     Open a 3D OSIRIS file and return the x, y and z axes and the data array.
 
     Parameters
@@ -106,7 +106,7 @@ def open3D(filename):
         The z axis.
     data_array : numpy.ndarray
         The data array.
-    """
+    '''
     attrs, axes, data = read_osiris_file(filename)
     datash = data.shape
     ax1 = axes[0]

osiris_utils/gui.py

@@ -14,7 +14,7 @@ from .utils import integrate, transverse_average  # Update import as needed
 class LAVA_Qt(QMainWindow):
     def __init__(self):
         super().__init__()
-        self.setWindowTitle("LAVA (LabAstro Visualization Assistant) - OSIRIS Data Grid Viewer")
+        self.setWindowTitle('LAVA (LabAstro Visualization Assistant) - OSIRIS Data Grid Viewer')
         self.setGeometry(100, 100, 1000, 600)
         
         # Initialize data
@@ -39,28 +39,28 @@ class LAVA_Qt(QMainWindow):
         control_layout = QHBoxLayout(control_frame)
         
         # Buttons
-        self.browse_btn = QPushButton("Browse Folder")
+        self.browse_btn = QPushButton('Browse Folder')
         self.browse_btn.clicked.connect(self.load_folder)
-        self.save_btn = QPushButton("Save Plot")
+        self.save_btn = QPushButton('Save Plot')
         self.save_btn.clicked.connect(self.save_plot)
         
         # File selector
         self.file_selector = QComboBox()
-        self.file_selector.setPlaceholderText("Select file...")
+        self.file_selector.setPlaceholderText('Select file...')
         self.file_selector.currentIndexChanged.connect(self.file_selection_changed)
         self.file_selector.view().setVerticalScrollBarPolicy(Qt.ScrollBarAlwaysOn)
         self.file_selector.setSizeAdjustPolicy(QComboBox.AdjustToContents)
         
         # Plot type combo box
         self.plot_combo = QComboBox()
-        self.plot_combo.addItem("Select Plot Type")
+        self.plot_combo.addItem('Select Plot Type')
         self.plot_combo.currentTextChanged.connect(self.plot_data)
         
         control_layout.addWidget(self.browse_btn)
         control_layout.addWidget(self.save_btn)
-        control_layout.addWidget(QLabel("Files:"))
+        control_layout.addWidget(QLabel('Files:'))
         control_layout.addWidget(self.file_selector)
-        control_layout.addWidget(QLabel("Plot Type:"))
+        control_layout.addWidget(QLabel('Plot Type:'))
         control_layout.addWidget(self.plot_combo)
         self.main_layout.addWidget(control_frame)
 
@@ -80,11 +80,11 @@ class LAVA_Qt(QMainWindow):
         self.ylabel_edit.textChanged.connect(self.update_plot_labels)
         
         
-        labels_layout.addWidget(QLabel("Title:"))
+        labels_layout.addWidget(QLabel('Title:'))
         labels_layout.addWidget(self.title_edit)
-        labels_layout.addWidget(QLabel("X Label:"))
+        labels_layout.addWidget(QLabel('X Label:'))
         labels_layout.addWidget(self.xlabel_edit)
-        labels_layout.addWidget(QLabel("Y Label:"))
+        labels_layout.addWidget(QLabel('Y Label:'))
         labels_layout.addWidget(self.ylabel_edit)
         
         # define the size of the labels frame
@@ -94,7 +94,7 @@ class LAVA_Qt(QMainWindow):
     def load_folder(self):
         folder_dialog = QFileDialog()
         folderpath = folder_dialog.getExistingDirectory(
-            self, "Select Folder with HDF5 Files"
+            self, 'Select Folder with HDF5 Files'
         )
         
         if not folderpath:
@@ -119,16 +119,16 @@ class LAVA_Qt(QMainWindow):
             h5_files.sort(key=sort_key)
             
             if not h5_files:
-                raise ValueError("No HDF5 files found in selected folder")
+                raise ValueError('No HDF5 files found in selected folder')
             
             self.file_selector.addItems(h5_files)
             self.file_selector.setCurrentIndex(0)
             
         except Exception as e:
-            QMessageBox.critical(self, "Error", str(e))
+            QMessageBox.critical(self, 'Error', str(e))
 
     def file_selection_changed(self, index):
-        """Handle file selection change in the combo box"""
+        '''Handle file selection change in the combo box'''
         if index >= 0 and self.current_folder:
             filename = self.file_selector.itemText(index)
             self.process_file(filename)
@@ -140,32 +140,32 @@ class LAVA_Qt(QMainWindow):
             self.dims = len(gridfile.axis)
             self.type = gridfile.type
             
-            if self.type == "grid":
+            if self.type == 'grid':
                 if self.dims == 1:
-                    x = np.linspace(gridfile.grid[0], gridfile.grid[1], gridfile.nx)
-                    self.xlabel_edit.setText(r"$%s$ [$%s$]" % (gridfile.axis[0]["long_name"], gridfile.axis[0]["units"]))
-                    self.ylabel_edit.setText(r"$%s$ [$%s$]" % (gridfile.label, gridfile.units))
+                    x = np.arange(gridfile.grid[0], gridfile.grid[1], gridfile.dx)
+                    self.xlabel_edit.setText(r'$%s$ [$%s$]' % (gridfile.axis[0]['long_name'], gridfile.axis[0]['units']))
+                    self.ylabel_edit.setText(r'$%s$ [$%s$]' % (gridfile.label, gridfile.units))
                     self.data_info = (x, gridfile.data)
                 elif self.dims == 2:
-                    x = np.linspace(gridfile.grid[0][0], gridfile.grid[0][1], gridfile.nx[0])
-                    y = np.linspace(gridfile.grid[1][0], gridfile.grid[1][1], gridfile.nx[1])
-                    self.xlabel_edit.setText(r"$%s$ [$%s$]" % (gridfile.axis[0]["long_name"], gridfile.axis[0]["units"]))
-                    self.ylabel_edit.setText(r"$%s$ [$%s$]" % (gridfile.axis[1]["long_name"], gridfile.axis[1]["units"]))
+                    x = np.arange(gridfile.grid[0][0], gridfile.grid[0][1], gridfile.dx[0])
+                    y = np.arange(gridfile.grid[1][0], gridfile.grid[1][1], gridfile.dx[1])
+                    self.xlabel_edit.setText(r'$%s$ [$%s$]' % (gridfile.axis[0]['long_name'], gridfile.axis[0]['units']))
+                    self.ylabel_edit.setText(r'$%s$ [$%s$]' % (gridfile.axis[1]['long_name'], gridfile.axis[1]['units']))
                     self.data_info = (x, y, gridfile.data)
                 elif self.dims == 3:
-                    raise ValueError("3D not supported yet")
+                    raise ValueError('3D not supported yet')
                 else:
-                    raise ValueError("Unsupported dimensionality")
+                    raise ValueError('Unsupported dimensionality')
                 
-                self.title_edit.setText(r"$%s$ [$%s$]" %( gridfile.label, gridfile.units))
+                self.title_edit.setText(r'$%s$ [$%s$]' %( gridfile.label, gridfile.units))
                 self.update_plot_menu()
                 self.plot_data()
                 
             else:
-                QMessageBox.information(self, "Info", f"{self.type} data not supported yet")
+                QMessageBox.information(self, 'Info', f'{self.type} data not supported yet')
                 
         except Exception as e:
-            QMessageBox.critical(self, "Error", str(e))
+            QMessageBox.critical(self, 'Error', str(e))
 
     def create_plot_area(self):
         # Matplotlib figure and canvas
@@ -194,9 +194,9 @@ class LAVA_Qt(QMainWindow):
         self.current_ax = self.figure.add_subplot(111)
         plot_type = self.plot_combo.currentText()
         
-        if "Line" in plot_type:
+        if 'Line' in plot_type:
             self.current_ax.plot(x, data)
-        elif "Scatter" in plot_type:
+        elif 'Scatter' in plot_type:
             self.current_ax.scatter(x, data)
             
         self.current_ax.set_xlabel(self.xlabel_edit.text())
@@ -208,17 +208,17 @@ class LAVA_Qt(QMainWindow):
         self.current_ax = self.figure.add_subplot(111)
         plot_type = self.plot_combo.currentText()
         
-        if "Quantity" in plot_type:
-            img = self.current_ax.imshow(data, extent=(x[0], x[-1], y[0], y[-1]), origin='lower', aspect='auto')
+        if 'Quantity' in plot_type:
+            img = self.current_ax.imshow(data.T, extent=(x[0], x[-1], y[0], y[-1]), origin='lower', aspect='auto')
             self.figure.colorbar(img)
-        elif "Integral" in plot_type:
+        elif 'Integral' in plot_type:
             avg = integrate(transverse_average(data), x[-1]/len(x))
             self.current_ax.plot(x, avg)
-        elif "Transverse" in plot_type:
+        elif 'Transverse' in plot_type:
             avg = transverse_average(data)
             self.current_ax.plot(x, avg)
-        elif "Phase" in plot_type:
-            img = self.current_ax.imshow(np.abs(-data), extent=(x[0], x[-1], y[0], y[-1]), origin='lower', aspect='auto', norm=LogNorm())
+        elif 'Phase' in plot_type:
+            img = self.current_ax.imshow(np.abs(-data.T), extent=(x[0], x[-1], y[0], y[-1]), origin='lower', aspect='auto', norm=LogNorm())
             self.figure.colorbar(img)
         
         self.current_ax.set_xlabel(self.xlabel_edit.text())
@@ -233,9 +233,9 @@ class LAVA_Qt(QMainWindow):
         
         # Determine items based on dimensions
         if self.dims == 1:
-            items = ["Line Plot", "Scatter Plot"]
+            items = ['Line Plot', 'Scatter Plot']
         elif self.dims == 2:
-            items = ["Quantity Plot", "T. Average Integral", "Transverse Average", "Phase Space"]
+            items = ['Quantity Plot', 'T. Average Integral', 'Transverse Average', 'Phase Space']
         else:
             items = []
         
@@ -250,11 +250,11 @@ class LAVA_Qt(QMainWindow):
     def save_plot(self):
         file_dialog = QFileDialog()
         filepath, _ = file_dialog.getSaveFileName(
-            self, "Save Plot", "", "PNG Files (*.png);;PDF Files (*.pdf)"
+            self, 'Save Plot', '', 'PNG Files (*.png);;PDF Files (*.pdf)'
         )
         
         if filepath:
-            self.figure.savefig(filepath, dpi=800, bbox_inches="tight")
+            self.figure.savefig(filepath, dpi=800, bbox_inches='tight')
 
 def LAVA():
     app = QApplication(sys.argv)
@@ -262,5 +262,5 @@ def LAVA():
     window.show()
     sys.exit(app.exec())
 
-if __name__ == "__main__":
+if __name__ == '__main__':
     LAVA()

osiris_utils/mean_field_theory.py

@@ -0,0 +1,52 @@
+import numpy as np
+from .data import OsirisGridFile
+
+        
+class MeanFieldTheory(OsirisGridFile):
+    '''
+    Class to handle the mean field theory on data. Inherits from OsirisGridFile.
+    
+    Parameters
+    ----------
+    source : str or OsirisGridFile
+        The filename or an OsirisGridFile object.
+    axis : int
+        The axis to average over.
+    '''
+    def __init__(self, source, axis=1):
+        if isinstance(source, OsirisGridFile):
+            self.__dict__.update(source.__dict__)
+        else:
+            super().__init__(source)
+        self._compute_mean_field(axis=axis)
+
+    def _compute_mean_field(self, axis=1):
+        self._average = np.expand_dims(np.mean(self.data, axis=axis), axis=axis)
+        self._fluctuations = self.data - self._average
+
+    def __array__(self):
+        return self.data
+
+    @property
+    def average(self):
+        return self._average
+    
+    @property
+    def delta(self):
+        return self._fluctuations
+    
+    def __str__(self):
+        return super().__str__() + f'\nAverage: {self.average.shape}\nDelta: {self.delta.shape}'
+    
+    def derivative(self, field, axis=0):
+        '''
+        Compute the derivative of the average or the fluctuations.
+        
+        Parameters
+        ----------
+        field : MeanFieldTheory.average or MeanFieldTheory.delta
+            The field to compute the derivative.
+        axis : int
+            The axis to compute the derivative.
+        '''
+        return np.gradient(field, self.dx[axis], axis=0)

osiris_utils/utils.py

@@ -7,7 +7,7 @@ import scipy
 import pandas as pd
 
 def courant2D(dx, dy):
-    """
+    '''
     Compute the Courant number for a 2D simulation.
 
     Parameters
@@ -21,12 +21,12 @@ def courant2D(dx, dy):
     -------
     float
         The limit for dt.
-    """
+    '''
     dt = 1 / (np.sqrt(1/dx**2 + 1/dy**2))
     return dt
 
 def time_estimation(n_cells, ppc, push_time, t_steps, n_cpu, hours = False):
-    """
+    '''
     Estimate the simulation time.
 
     Parameters
@@ -48,7 +48,7 @@ def time_estimation(n_cells, ppc, push_time, t_steps, n_cpu, hours = False):
     -------
     float
         The estimated time in seconds or hours.
-    """
+    '''
     time = (n_cells*ppc*push_time*t_steps)/n_cpu
     if hours:
         return time/3600
@@ -59,7 +59,7 @@ def filesize_estimation(n_gridpoints):
     return n_gridpoints*4/(1024**2)
 
 def transverse_average(data):
-    """
+    '''
     Computes the transverse average of a 2D array.
     
     Parameters
@@ -74,14 +74,14 @@ def transverse_average(data):
         Dim: 1D.
         The transverse average.
 
-    """
+    '''
 
     if len(data.shape) != 2:
-        raise ValueError("The input data must be a 2D array.")
+        raise ValueError('The input data must be a 2D array.')
     return np.mean(data, axis = 1)
 
 def integrate(array, dx):
-    """
+    '''
     Integrate a 1D from the left to the right. This may be changed in the future to allow 
     for integration in both directions or for other more general cases.
 
@@ -98,24 +98,24 @@ def integrate(array, dx):
     numpy.ndarray
         Dim: 1D.
         The integrated array.
-    """
+    '''
 
     if len(array.shape) != 1:
-        raise ValueError(f"Array must be 1D\n Array shape: {array.shape}")
+        raise ValueError(f'Array must be 1D\n Array shape: {array.shape}')
     flip_array = np.flip(array)
     # int = -scipy.integrate.cumulative_trapezoid(flip_array, dx = dx, initial = flip_array[0])
     int = -scipy.integrate.cumulative_simpson(flip_array, dx = dx, initial = 0)
     return np.flip(int)
 
 def animate_2D(datafiles, frames, interval, fps, savename, **kwargs):
-    """
+    '''
     Animate 2D OSIRIS files.
 
     Parameters
     ----------
     datafiles : str
         The path to the files.
-        Must be of the type "path/to/file_%06d.h5".
+        Must be of the type 'path/to/file_%06d.h5'.
     kwargs : dict
         Additional keyword arguments for plotting.
 
@@ -123,7 +123,7 @@ def animate_2D(datafiles, frames, interval, fps, savename, **kwargs):
     -------
     matplotlib.animation.FuncAnimation
         The animation.
-    """
+    '''
     fig, ax = plt.subplots(figsize=(12, 6), tight_layout=True)
     im = 0
 
@@ -133,8 +133,8 @@ def animate_2D(datafiles, frames, interval, fps, savename, **kwargs):
         ax.clear()
         # Display image data, make sure data shape is valid for imshow
         im = ax.imshow(-data, extent=[x.min(), x.max(), y.min(), y.max()], aspect='auto', origin='lower', **kwargs)
-        plt.xlabel(r"x [c/$\omega_p$]")
-        plt.ylabel(r"y [c/$\omega_p$]")
+        plt.xlabel(r'x [c/$\omega_p$]')
+        plt.ylabel(r'y [c/$\omega_p$]')
         
 
     # Creating the animation, and frames should be updated accordingly
@@ -146,7 +146,7 @@ def animate_2D(datafiles, frames, interval, fps, savename, **kwargs):
     # Display the animation
     return ani
 
-def save_data(data, savename, option="numpy"):
+def save_data(data, savename, option='numpy'):
     """
     Save the data to a .txt (with Numpy) or .csv (with Pandas) file.
 
@@ -157,17 +157,17 @@ def save_data(data, savename, option="numpy"):
     savename : str
         The path to the file.
     option : str, optional
-        The option for saving the data. The default is "numpy". Can be "numpy" or "pandas".
+        The option for saving the data. The default is 'numpy'. Can be 'numpy' or 'pandas'.
     """
-    if option == "numpy":
+    if option == 'numpy':
         np.savetxt(savename, data)
-    elif option == "pandas":
+    elif option == 'pandas':
         pd.DataFrame(data).to_csv(savename, index=False)
     else:
         raise ValueError("Option must be 'numpy' or 'pandas'.")
 
-def read_data(filename, option="numpy"):
-    """
+def read_data(filename, option='numpy'):
+    '''
     Read the data from a .txt file.
 
     Parameters
@@ -180,54 +180,5 @@ def read_data(filename, option="numpy"):
     numpy.ndarray
         Dim: 2D.
         The data.
-    """
-    return np.loadtxt(filename) if option == "numpy" else pd.read_csv(filename).values
-
-def mft_decomposition(filename, pressure =  False, xy = False, data = False):
-    """
-    Mean Field Theory decomposition of the data.
-    Considering that A = ⟨A⟩ + δA with ⟨δA⟩ = 0
-    This function returns ⟨A⟩ and δA from A. 
-    
-    Parameters
-    ----------
-    filename : str
-        The path to the file.
-        The data is 2D.
-    pressure : bool, optional
-        If True, the file is a pressure file. The default is False.
-    xy : bool, optional
-        If True, the function returns x and y axes. The default is False.
-    data : bool, optional
-        If True, the function returns the data (2D with no transformation). The default is False.
-        
-    Returns
-    -------
-    mean : numpy.ndarray
-        Dim: 1D.
-        The mean field ⟨A⟩.
-    fluctuation : numpy.ndarray
-        Dim: 2D.
-        The fluctuation δA.
-    x : numpy.ndarray
-        Dim: 1D.
-        The x axis.
-    y : numpy.ndarray
-        Dim: 1D.
-        The y axis.
-    data : numpy.ndarray
-        Dim: 2D.
-        The data.
-    """ 
-    x, y, d, _ = open2D(filename, pressure)
-    mean = transverse_average(d)
-    fluctuation = d - mean
-   
-    if xy and data:
-        return mean, fluctuation, x, y, d
-    elif xy:
-        return mean, fluctuation, x, y
-    elif data:
-        return mean, fluctuation, d
-    else:
-        return mean, fluctuation
+    '''
+    return np.loadtxt(filename) if option == 'numpy' else pd.read_csv(filename).values

{osiris_utils-1.0.5.dist-info → osiris_utils-1.0.9.dist-info}/METADATA

@@ -1,8 +1,8 @@
 Metadata-Version: 2.2
 Name: osiris_utils
-Version: 1.0.5
+Version: 1.0.9
 Summary: Utilities to manipulate and visualize OSIRIS framework output data
-Author: ['João Pedro Ferreira Biu', 'João Cândido', 'Rafael Russo de Almeida']
+Author: ['João Pedro Ferreira Biu', 'João Cândido']
 Author-email: ['joaopedrofbiu@tecnico.ulisboa.pt']
 License: MIT
 Project-URL: Issues Tracker, https://github.com/joaopedrobiu6/osiris_utils/issues
@@ -17,8 +17,9 @@ Classifier: Programming Language :: Python :: 3.8
 Classifier: Programming Language :: Python :: 3.9
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
 Classifier: Topic :: Scientific/Engineering
-Requires-Python: >=3.12
+Requires-Python: >=3.11
 Description-Content-Type: text/markdown
 License-File: LICENSE.txt
 Requires-Dist: numpy
@@ -27,6 +28,8 @@ Requires-Dist: pandas
 Requires-Dist: scipy
 Requires-Dist: h5py
 Requires-Dist: pyside6
+Requires-Dist: pyarrow
+Requires-Dist: tqdm
 Dynamic: author
 Dynamic: author-email
 Dynamic: classifier

osiris_utils-1.0.9.dist-info/RECORD

@@ -0,0 +1,11 @@
+osiris_utils/__init__.py,sha256=5qsicKjKzheUmdI89cIw7Egvjv6OZcGAQ5eEAiGZIcY,358
+osiris_utils/data.py,sha256=vj5WV7SubCqR3Ia3ylFuwR5hyotk31d1xZLfYyjC5ow,16491
+osiris_utils/data_readers.py,sha256=ZbzYQ-MK5NSv4Hfgtdc1YqXHs9yL_ZEKFFCiRn7GN5o,3264
+osiris_utils/gui.py,sha256=YyG4MBmUFDJ141-NeorVRnRU3fcym4cfDMSGsJUFbYA,10531
+osiris_utils/mean_field_theory.py,sha256=sJGvoJ5fTgW6fQZjNNsPswzXkh7dhMgO96qB8R5tl0M,1517
+osiris_utils/utils.py,sha256=dlW2BvYpVSsW48kA_A6NNQ77VvKfyvs8hVwnOfJu2qk,4736
+osiris_utils-1.0.9.dist-info/LICENSE.txt,sha256=Cawy2v7wKc7n8yL8guFu-cH9sQw9r1gll1pEFPFAB-Q,1084
+osiris_utils-1.0.9.dist-info/METADATA,sha256=YqRPvxzvsA5TlTHP0p961tkaDrUbP7n9hqh-jJJ7_5E,1463
+osiris_utils-1.0.9.dist-info/WHEEL,sha256=52BFRY2Up02UkjOa29eZOS2VxUrpPORXg1pkohGGUS8,91
+osiris_utils-1.0.9.dist-info/top_level.txt,sha256=mM-_dX5fjzIKB7te655PhZOrPACVY-bJmiASCqW1eOA,13
+osiris_utils-1.0.9.dist-info/RECORD,,

{osiris_utils-1.0.5.dist-info → osiris_utils-1.0.9.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.8.0)
+Generator: setuptools (76.0.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

osiris_utils-1.0.5.dist-info/RECORD

@@ -1,10 +0,0 @@
-osiris_utils/__init__.py,sha256=mXRRxiNNJ6pxmRGO23OCNR3Yn3MvrBZFiNI9Rui_un8,318
-osiris_utils/data.py,sha256=awf6v_zWsOTPjdbLQ_8E0SsaGGvp_UQZGUBh2UuCIyk,6678
-osiris_utils/data_readers.py,sha256=ZeqbYHsFmdKRxqUNkkauurjz-1VABW4GZ-WaPQREdfc,3264
-osiris_utils/gui.py,sha256=DPRqqFIYfgNZEMnt4rEo5TiHe5fogkWySrOllQvj0ak,10533
-osiris_utils/utils.py,sha256=SrHsMaOq1zw9Qx9YJV7jJvS1u6hr5RLV7IvikygM4pg,6097
-osiris_utils-1.0.5.dist-info/LICENSE.txt,sha256=Cawy2v7wKc7n8yL8guFu-cH9sQw9r1gll1pEFPFAB-Q,1084
-osiris_utils-1.0.5.dist-info/METADATA,sha256=6kPFf27IPybWbwPEyf36VxPgRyD6-Du81tfVj04DeSM,1396
-osiris_utils-1.0.5.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-osiris_utils-1.0.5.dist-info/top_level.txt,sha256=mM-_dX5fjzIKB7te655PhZOrPACVY-bJmiASCqW1eOA,13
-osiris_utils-1.0.5.dist-info/RECORD,,