o
    EbBg                     @   sx  d gZ ddlZddlZddlZddlZddlZddddddd	d	d
d	d	dddddZdddddddddddddddZi Z	dd  Z
d!d" Zd#d$ Zd%d& Zd'd( Zd)d* Zd+d, Zd-d. Zd/d0 Zd1d2 Zd3d4 Zd5d6 Zd7d8 ZG d9d: d:ZG d;d< d<eZd=d> Zd?d@ ZdAdB ZdCdD ZdEdF ZdGdH ZdIdJ ZdKdL Z dMdN Z!dOdP Z"dQdR Z#G dSdT dTe$Z%	U	UdWdVd Z&dS )Xreadsav    Nz>u1z>i2z>i4z>f4z>f8z>c8z|Oz>c16z>u2z>u4z>i8z>u8)                        	   
                  ZSTART_MARKERZCOMMON_VARIABLEVARIABLESYSTEM_VARIABLE
END_MARKER	TIMESTAMPZCOMPILEDZIDENTIFICATIONVERSIONHEAP_HEADER	HEAP_DATAZ	PROMOTE64NOTICEDESCRIPTION)r   r   r   r   r   r   r   r   r   r               c                 C   s.   |   }|d dkr| |d |d   dS )z+Align to the next 32-bit position in a filer   r   N)tellseek)fpos r#   //usr/lib/python3/dist-packages/scipy/io/_idl.py	_align_32K   s   r%   c                 C   s   |  | dS )zSkip `n` bytesNreadr!   nr#   r#   r$   _skip_bytesT   s   
r*   c                 C   s
   |  |S )zRead the next `n` bytesr&   r(   r#   r#   r$   _read_bytesZ   s   
r+   c                 C   s$   t td| ddd d S )zRead a single bytez>Br   Nr   r   )npZuint8structunpackr'   r!   r#   r#   r$   
_read_byte_      $r0   c                 C      t td| dd S )Read a signed 32-bit integer>lr   r   r,   Zint32r-   r.   r'   r/   r#   r#   r$   
_read_longd      r6   c                 C   $   t td| ddd d S )zRead a signed 16-bit integerz>hr   r   r   )r,   Zint16r-   r.   r'   r/   r#   r#   r$   _read_int16i   r1   r9   c                 C   r2   )r3   z>ir   r   r5   r/   r#   r#   r$   _read_int32n   r7   r:   c                 C   r2   )zRead a signed 64-bit integerz>qr
   r   )r,   Zint64r-   r.   r'   r/   r#   r#   r$   _read_int64s   r7   r;   c                 C   r8   )zRead an unsigned 16-bit integerz>Hr   r   r   )r,   Zuint16r-   r.   r'   r/   r#   r#   r$   _read_uint16x   r1   r<   c                 C   r2   )zRead an unsigned 32-bit integer>Ir   r   )r,   Zuint32r-   r.   r'   r/   r#   r#   r$   _read_uint32}   r7   r>   c                 C   r2   )zRead an unsigned 64-bit integerz>Qr
   r   )r,   Zuint64r-   r.   r'   r/   r#   r#   r$   _read_uint64   r7   r?   c                 C   r2   )zRead a 32-bit floatz>fr   r   )r,   Zfloat32r-   r.   r'   r/   r#   r#   r$   _read_float32   r7   r@   c                 C   r2   )zRead a 64-bit floatz>dr
   r   )r,   Zfloat64r-   r.   r'   r/   r#   r#   r$   _read_float64   r7   rA   c                   @   s   e Zd ZdZdd ZdS )PointerzClass used to define pointersc                 C   s
   || _ d S N)index)selfrD   r#   r#   r$   __init__   s   zPointer.__init__N)__name__
__module____qualname____doc__rF   r#   r#   r#   r$   rB      s    rB   c                   @   s   e Zd ZdZdS )ObjectPointerz$Class used to define object pointersN)rG   rH   rI   rJ   r#   r#   r#   r$   rK      s    rK   c                 C   s4   t | }|dkrt| |d}t|  |S d}|S )zRead a stringr   latin1 )r6   r+   decoder%   )r!   lengthcharsr#   r#   r$   _read_string   s   rQ   c                 C   s6   t | }|dkrt | }t| |}t|  |S d}|S )z.Read a data string (length is specified twice)r   rM   )r6   r+   r%   )r!   rO   Zstring_datar#   r#   r$   _read_string_data   s   
rR   c                 C   sL  |dkrt | dkrtdt| S |dkrt| S |dkr"t | S |dkr*t| S |dkr2t| S |dkrGt| }t| }t||d  S |d	krOt| S |d
krWtd|dkrlt| }t| }t	||d  S |dkrvt
t | S |dkrtt | S |dkrt| S |dkrt| S |dkrt| S |dkrt| S td| )z*Read a variable with a specified data typer   z*Error occurred while reading byte variabler   r   r   r   r   y              ?r	   r
   z'Should not be here - please report thisr   r   r   r   r   r   r   z)Unknown IDL type: %i - please report this)r:   	Exceptionr0   r9   r@   rA   r,   Z	complex64rR   Z
complex128rB   rK   r<   r>   r;   r?   )r!   dtyperealimagr#   r#   r$   
_read_data   sJ   rW   c           
      C   sz  |d }|d }g }|D ]<}|d s|d r'| |d  |d ftjf q|d tv rA| |d  |d ft|d  f qtd|d  tj|f|d}t|D ]J}|D ]E}|d }|d r|t| |d	 |d  |d
 |d  ||d  |< qY|d rt	| ||d	 |d  ||d  |< qYt
| |||d  |< qYqU|d dkr|d dt|d  }	|	  ||	}|S )z
    Read a structure, with the array and structure descriptors given as
    `array_desc` and `structure_desc` respectively.
    	nelementstagtable	structurearraynametypecodez Variable type %i not implementedrT   arrtablestructtablendimsr   dimsN)appendlowerr,   object_
DTYPE_DICTrS   recarrayrange_read_structure_read_arrayrW   intreversereshape)
r!   
array_descstruct_descZnrowscolumnsrT   colrZ   irb   r#   r#   r$   ri      sB   "

ri   c           	      C   s   |dv r&|dkrt | }||d krtd tj| |d t| d}n9|dv r@tj| |d d t| dddd }ng }t|d	 D ]}|}t| |}|	| qHtj
|tjd}|d
 dkrz|d dt|d
  }|  ||}t|  |S )z`
    Read an array of type `typecode`, with the array descriptor given as
    `array_desc`.
    )	r   r   r   r   r   r   r   r   r   r   nbytesz.Not able to verify number of bytes from headerr^   )r   r   r   NrX   ra   rb   )r:   warningswarnr,   Z
frombufferr'   rf   rh   rW   rc   r[   re   rk   rl   rm   r%   )	r!   r]   rn   rs   r[   rr   rT   datarb   r#   r#   r$   rj     s4   


rj   c                    s  dt  i}t }|t d 7 }t d |d tvr%td|d  t|d  |d< |d dv r|d dkr@t |d< nt  |d< t d t }|d	 d
kre|  krad|d< ntdt  }|dkrqtd|d rt	 |d |d |d< n|d rt
 |d	 |d |d< n|d	 }t ||d< n|d dkrt d t |d< t |d< t |d< n|d dkrt  |d< t |d< t |d< t |d< n|d dkrt |d< t |d < t |d!< n|d d"krt |d#< n||d d$krt |d%< nn|d d&kr.t  |d'<  fd(d)t|d' D |d*< nQ|d d+krQt  |d,< t |d-<  fd.d)t|d, D |d/< n.|d d0kr]d1|d2< n"|d d3krjtd4 n|d d5krwtd6 ntd7|d   | |S )8z!Function to read in a full recordrectype        r   zUnknown RECTYPE: %i)r   r   r   varname
heap_indexr]   r   Nrv   zUnexpected type code: 0r	   zVARSTART is not 7rZ   rn   ro   r[   r   i   dateuserhostr   formatarchosreleaseIDENTIFICATONauthortitleidcoder   Znoticer   descriptionr   Znvaluesc                       g | ]}t  qS r#   r6   .0_r/   r#   r$   
<listcomp>      z _read_record.<locals>.<listcomp>indicesZCOMMONBLOCKZnvarsr\   c                    r   r#   rQ   r   r/   r#   r$   r     r   varnamesr   TendZUNKNOWNzSkipping UNKNOWN recordr   zSkipping SYSTEM_VARIABLE recordz$record['rectype']=%s not implemented)r6   r>   r*   RECTYPE_DICTrS   rQ   _read_typedescr   
ValueErrorri   rj   rW   rR   rh   rt   ru   r    )r!   recordnextrecZrectypedescZvarstartrT   r#   r/   r$   _read_record;  s   






  

r   c                 C   s   t | t | d}|d d@ dkrtd|d d@ dk|d< |d d@ dk|d< |d r;t| |d	< t| |d
< |S |d rEt| |d	< |S )z%Function to read in a type descriptor)r]   varflagsr   r   z System variables not implementedr   r[       rZ   rn   ro   )r6   rS   _read_arraydesc_read_structdesc)r!   Ztypedescr#   r#   r$   r     s   r   c                    s,  dt  i}|d dkr?t d t  |d< t  |d< t  |d< t d t  |d<  fdd	t|d D |d
< |S |d dkrtd t d t |d< t |d< t  |d< t d d|d< g |d
< t|d D ]}t  }|dkrtd|d
 t   qt|S td|d  )z'Function to read in an array descriptorZarrstartr
   r   rs   rX   ra   Znmaxc                    r   r#   r   r   r/   r#   r$   r     r   z#_read_arraydesc.<locals>.<listcomp>rb      z$Using experimental 64-bit array readr   zExpected a zero in ARRAY_DESCzUnknown ARRSTART: %i)r6   r*   rh   rt   ru   r?   rS   rc   )r!   Z	arraydescdvr#   r/   r$   r     s4   




r   c                    sx  i }t  }|dkrtdt |d< t  }t  |d< t  |d< |d@ |d< |d@ |d	< |d
@ |d< |d s fddt|d D |d< |d D ]}t |d< qM fdd|d D |d<  fdd|d D |d< |d	 sx|d rt |d< t  |d<  fddt|d D |d<  fddt|d D |d< |t|d < |S |d tvrtdt|d  }|S )z*Function to read in a structure descriptorr   zSTRUCTSTART should be 9r\   Zntagsrs   r   predefr   Zinheritsr   Zis_superc                    r   r#   )_read_tagdescr   r/   r#   r$   r     s    z$_read_structdesc.<locals>.<listcomp>rY   c                    "   i | ]}|d  r|d t  qS )r[   r\   )r   r   tagr/   r#   r$   
<dictcomp>      z$_read_structdesc.<locals>.<dictcomp>r_   c                    r   )rZ   r\   r   r   r/   r#   r$   r   	  r   r`   Z	classnameZnsupclassesc                    r   r#   r   r   r/   r#   r$   r         Zsupclassnamesc                    r   r#   r   r   r/   r#   r$   r     r   Zsupclasstablez"PREDEF=1 but can't find definition)r6   rS   rQ   rh   STRUCT_DICT)r!   Z
structdescZstructstartr   r   r#   r/   r$   r     sL   












	r   c                 C   sl   dt | i}|d dkrt| |d< t | |d< t | }|d@ dk|d< |d@ dk|d< |d tv |d< |S )	z$Function to read in a tag descriptoroffsetr]   r   r[   r   rZ   Zscalar)r6   r?   rf   )r!   ZtagdescZtagflagsr#   r#   r$   r   !  s   r   c                 C   sV  t | tr8t | tr)| jdkrd } n| j|v r|| j } ntd d } t | ts
t| |\}}|r4|} d| fS t | tjjj	rZt
| D ]\}}t||\}}|rU|| |< qDd| fS t | tjjjr|t
| D ]\}}t||\}}|rw|| |< qfd| fS t | tjr| jjtju rt| jD ]}t| ||\}}|r| || qd| fS d| fS )Nr   zNVariable referenced by pointer not found in heap: variable will be set to NoneTF)
isinstancerB   rD   rt   ru   _replace_heapr,   Zcorerecordsrg   	enumerater   ZndarrayrT   typere   rh   sizeitemZitemset)variableheapreplacenewZirr   Zivvaluer#   r#   r$   r   4  sH   





r   c                       sD   e Zd ZdZi fddZ fddZ fddZeZeZeZ	  Z
S )AttrDicta.  
    A case-insensitive dictionary with access via item, attribute, and call
    notations:

        >>> d = AttrDict()
        >>> d['Variable'] = 123
        >>> d['Variable']
        123
        >>> d.Variable
        123
        >>> d.variable
        123
        >>> d('VARIABLE')
        123
    c                 C   s   t | | d S rC   )dictrF   )rE   Zinitr#   r#   r$   rF        zAttrDict.__init__c                    s   t  | S rC   )super__getitem__rd   )rE   r\   	__class__r#   r$   r     r   zAttrDict.__getitem__c                    s   t  | |S rC   )r   __setitem__rd   )rE   keyr   r   r#   r$   r     s   zAttrDict.__setitem__)rG   rH   rI   rJ   rF   r   r   __getattr____setattr____call____classcell__r#   r#   r   r$   r   v  s    r   Fc              	   C   s  g }|s|r	i }nt  }t| d}t|d}|dkr td| t|d}	|	dkr*n|	dkr|r4td |r<t|d}
ntjd	d
}
|rKtd|
j  |
d 	 t	|}|
t
dt| t|}|t|d 7 }|d}t| dkr|
t
dt|d  |
t
dt||d  d  |
| nA| }t||| }|
 t| d }|
t
dt|d  |
t
dt||d  d  |
| |
| qQ|  |
}|d ntd|	 	 t|}|| d|v r|d rnq|  i }|D ]}|d dkr|d ||d < q|D ]$}|d dkrDt|d |\}}|r:||d< |d ||d  < q!|r;|D ](}|d dkrrtd td|d   td |d!   td"|d#    nqK|D ]0}|d d$krtd td%|d&   td'|d(   td)|d*   td+|d,    nqv|D ](}|d d-krtd td.|d/   td0|d1   td2|d3    nq|D ]}|d d4krtd td5|d6    nqtd td7t|  d8d9 |D }t|D ]}|dkrtd:|||f  qtd d|v r;td; |D ]}td<|t|| f  q'td |rL|D ]	}|| ||< q@|S |S )=a	  
    Read an IDL .sav file.

    Parameters
    ----------
    file_name : str
        Name of the IDL save file.
    idict : dict, optional
        Dictionary in which to insert .sav file variables.
    python_dict : bool, optional
        By default, the object return is not a Python dictionary, but a
        case-insensitive dictionary with item, attribute, and call access
        to variables. To get a standard Python dictionary, set this option
        to True.
    uncompressed_file_name : str, optional
        This option only has an effect for .sav files written with the
        /compress option. If a file name is specified, compressed .sav
        files are uncompressed to this file. Otherwise, readsav will use
        the `tempfile` module to determine a temporary filename
        automatically, and will remove the temporary file upon successfully
        reading it in.
    verbose : bool, optional
        Whether to print out information about the save file, including
        the records read, and available variables.

    Returns
    -------
    idl_dict : AttrDict or dict
        If `python_dict` is set to False (default), this function returns a
        case-insensitive dictionary with item, attribute, and call access
        to variables. If `python_dict` is set to True, this function
        returns a Python dictionary with all variable names in lowercase.
        If `idict` was specified, then variables are written to the
        dictionary specified, and the updated dictionary is returned.

    Examples
    --------
    >>> from os.path import dirname, join as pjoin
    >>> import scipy.io as sio
    >>> from scipy.io import readsav

    Get the filename for an example .sav file from the tests/data directory.

    >>> data_dir = pjoin(dirname(sio.__file__), 'tests', 'data')
    >>> sav_fname = pjoin(data_dir, 'array_float32_1d.sav')

    Load the .sav file contents.

    >>> sav_data = readsav(sav_fname)

    Get keys of the .sav file contents.

    >>> print(sav_data.keys())
    dict_keys(['array1d'])

    Access a content with a key.

    >>> print(sav_data['array1d'])
    [0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
     0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
     0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
     0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
     0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
     0. 0. 0.]

    rbr   s   SRzInvalid SIGNATURE: %ss    s    zIDL Save file is compressedzw+bz.sav)suffixz -> expanding to %ss   SR Tr4   rx   r   r   r=   r   zInvalid RECFMT: %sr   rw   r   rv   rz   r   ry   r   z2--------------------------------------------------zDate: %sr{   zUser: %sr|   zHost: %sr}   r   z
Format: %sr~   zArchitecture: %sr   zOperating System: %sr   zIDL Version: %sr   r   z
Author: %sr   z	Title: %sr   zID Code: %sr   r   zDescription: %sr   z&Successfully read %i records of which:c                 S   s   g | ]}|d  qS )rw   r#   )r   rr#   r#   r$   r   q  r   zreadsav.<locals>.<listcomp>z - %i are of type %szAvailable variables:z
 - %s [%s])r   openr+   rS   printtempfileZNamedTemporaryFiler\   writer6   r-   Zpackrk   r>   r'   r   r   zlib
decompresslencloser    r   rc   r   rd   setcountr   )	file_nameZidictZpython_dictZuncompressed_file_nameverboser   Z	variablesr!   Z	signatureZrecfmtZfoutrw   r   unknownr"   Z
rec_stringr   r   r   r   r   ZrectypesZrtvarr#   r#   r$   r     s   F




"
"

$


	

)NFNF)'__all__r-   Znumpyr,   r   r   rt   rf   r   r   r%   r*   r+   r0   r6   r9   r:   r;   r<   r>   r?   r@   rA   rB   rK   rQ   rR   rW   ri   rj   r   r   r   r   r   r   r   r   r   r#   r#   r#   r$   <module>   s   	*-/m/6B