def METHOD_NAME(self):
    """
    Property to get parsed matrices converted to numpy arrays using atomic units.
    Returns:
        numpy.array: Array of the parsed GLE matrix with the shape normal_modes x s+1 x s+1, where normal_modes is 1
                     except in the case of the PIGLET thermostat and s is the number of degrees of freedom added via
                     GLE. If no matrix is found, None is returned.
    """
    # Write out last buffer
    if len(self._tmp_matrix) > 0:
        self._matrix.append(self._tmp_matrix)
    # Convert to numpy array
    _matrix = np.array(self._matrix)
    # Perform unit conversion
    if self.units is not None:
        return _matrix * self.units
    else:
        return None