U mdP@s~dZddlZddlZddlZddlZddlmZm Z m Z m Z m Z m Z mZddlmZmZddlmZddd d d d d dddddddddddddgZejdkZdZdZdZdZdZdZed ed!ed"ed#ed$ed%d&d'd(d)d*d+gd,d,iZeeeeeeeeeeed- ZdZdZ dZ!d.d/d0d1d2d3Z"ed4e d5e!d6iZ#Gd7d d Z$Gd8d d Z%Gd9ddeZ&d>d;d Z'Gdz VAX D-floatz VAX G-floatZCray)rrr"r#r$doublecharsparsec@seZdZdZdZddZdS)rFcCs"||_||_||_||_||_dSN)namedtypemclassdims is_complex)selfr:r;r<r=r>r@N/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/scipy/io/matlab/_mio4.py__init__Ys zVarHeader4.__init__N)__name__ __module__ __qualname__Z is_logical is_globalrBr@r@r@rArTsc@sTeZdZdZddZddZdddZdd d Zd d Zd dZ ddZ ddZ dS)rz" Class to read matlab 4 variables cCs*||_|j|_|j|_|j|_|j|_dSr9) file_reader mat_streamdtypeschars_as_strings squeeze_me)r?rGr@r@rArBis zVarReader4.__init__c Cst|j|jd}|jt|dd}|ddksD|ddkrLtdt|dd\}}|d krzt d t |t t|d \}}|dkrtd t|d \}}|}|d|df}|ddk} |j|} t || ||| S)z% Read and return header for variable r,r1r-rz.Mat 4 mopt wrong format, byteswapping problem?)rrzBWe do not support byte ordering '%s'; returned data may be corruptdz,O in MOPT integer should be 0, wrong format? r.r/r0r) rrHrIreadintstrip ValueErrordivmodwarningswarnr UserWarningr) r?datar:MrestOPTr=r>r;r@r@rA read_headerps4  zVarReader4.read_headerTcCst|j}|tkr||}nD|tkr@||}|r^|jr^t|}n|tkrR||St d||rp|j rpt |S|S)NNo reader for class code %s) r<rread_full_arrayrread_char_arrayrJr rread_sparse_array TypeErrorrKr )r?hdrprocessr<arrr@r@rAarray_from_headers       zVarReader4.array_from_headerc Csp|j}|j}|j}|D] }||9}q|jt|}t||krNtd|jt j |||dd}|rl| }|S)a Mat4 read using header `hdr` dtype and dims Parameters ---------- hdr : object object with attributes ``dtype``, ``dims``. dtype is assumed to be the correct endianness copy : bool, optional copies array before return if True (default True) (buffer is usually read only) Returns ------- arr : ndarray of dtype given by `hdr` ``dtype`` and shape given by `hdr` ``dims`` zNot enough bytes to read matrix '%s'; is this a badly-formed file? Consider listing matrices with `whosmat` and loading named matrices with `variable_names` kwarg to `loadmat`F)shaper;bufferorder) r;r=itemsizerHrQrRlenrTr:npndarraycopy) r?rerqdtr= num_bytesdrkrgr@r@rAread_sub_arrays$  zVarReader4.read_sub_arraycCs8|jr.|j|dd}|j|dd}||dS||S)aM Full (rather than sparse) matrix getter Read matrix (array) can be real or complex Parameters ---------- hdr : ``VarHeader4`` instance Returns ------- arr : ndarray complex array if ``hdr.is_complex`` is True, otherwise a real numeric array F)rqy?)r>ru)r?reresZres_jr@r@rAras  zVarReader4.read_full_arraycCsB||tj}|d}tj|jtdt |d S)z latin-1 text matrix (char matrix) reader Parameters ---------- hdr : ``VarHeader4`` instance Returns ------- arr : ndarray with dtype 'U1', shape given by `hdr` ``dims`` latin-1r2rjr;rk) ruastyperoZuint8tobytesdecoderpr=r;arrayrq)r?rergSr@r@rArbs zVarReader4.read_char_arraycCs||}|ddddf}t|dt|df}tj|dddfdd}tj|dddfdd}|d8}|d8}|jdd krtj|ddd fd d}n,tj|ddd fd d}|ddd f|_tj|||ff|S) a Read and return sparse matrix type Parameters ---------- hdr : ``VarHeader4`` instance Returns ------- arr : ``scipy.sparse.coo_matrix`` with dtype ``float`` and shape read from the sparse matrix data Notes ----- MATLAB 4 real sparse arrays are saved in a N+1 by 3 array format, where N is the number of non-zero values. Column 1 values [0:N] are the (1-based) row indices of the each non-zero value, column 2 [0:N] are the column indices, column 3 [0:N] are the (real) values. The last values [-1,0:2] of the rows, column indices are shape[0] and shape[1] respectively of the output matrix. The last value for the values column is a padding 0. mrows and ncols values from the header give the shape of the stored matrix, here [N+1, 3]. Complex data are saved as a 4 column matrix, where the fourth column contains the imaginary component; the last value is again 0. Complex sparse data do *not* have the header ``imagf`` field set to True; the fact that the data are complex is only detectable because there are 4 storage columns. N)r~r)r~rrZintcr;rr#r"floatcomplex) rurRroZascontiguousarrayrjimagscipyr8Z coo_matrix)r?rervtmpr=IJVr@r@rArcs zVarReader4.read_sparse_arraycCs2|j}|tkr ttt|j}n|tkrLttt|j}|jrJ|dd}n|tkr|j }|j}t |dkr|ddkr|ddksdS|j |j |dddtjd||j |j d}|j |j |dddtjd||j |j d}t|t|f}n td||jr.td d |D}|S) zwRead the shape of the array described by the header. The file position after this call is unspecified. Nr~r"rrr@rxr`cSsg|]}|dkr|qS)rr@).0xr@r@rA 0sz0VarReader4.shape_from_header..)r<rtuplemaprRr=rrJrr;rnrHseekrmrorprQrdrK)r?rer<rjrrr=rowscolsr@r@rAshape_from_headers2 $   zVarReader4.shape_from_headerN)T)T) rCrDrE__doc__rBr_rhrurarbrcrr@r@r@rArfs  $*csXeZdZdZefddZddZddZdd Zdd d Z dddZ ddZ Z S)r z Reader for Mat4 files cstj|f||d|_dS)zR Initialize matlab 4 file reader %(matstream_arg)s %(load_args)s N)superrB_matrix_reader)r?rHargskwargs __class__r@rArB6szMatFile4Reader.__init__cCs^|jdt|jtd}|jd|dkr6dS|dksF|dkrRtrNdpPdStrZdp\dS)Nrr(r4rMr5)rHrrror;r)r?r-r@r@rAguess_byte_order@s   zMatFile4Reader.guess_byte_ordercCstt|j|_t||_dS)za Run when beginning read of variables Sets up readers from parameters in `self` N)rrZ byte_orderrIrr)r?r@r@rAinitialize_readLszMatFile4Reader.initialize_readcCsV|j}tdd|jd}|jj|}|jr@|jtks@|d9}|j |}||fS)ab Read and return header, next position Parameters ---------- None Returns ------- header : object object that can be passed to self.read_var_array, and that has attributes ``name`` and ``is_global`` next_position : int position in stream of next variable cSs||Sr9r@)ryr@r@rAdz0MatFile4Reader.read_var_header..rr") rr_r r=r;rmr>r<rrHtell)r?renZremaining_bytes next_positionr@r@rAread_var_headerTs  zMatFile4Reader.read_var_headerTcCs|j||S)a Read array, given `header` Parameters ---------- header : header object object with fields defining variable header process : {True, False}, optional If True, apply recursive post-processing during loading of array. Returns ------- arr : array array with post-processing applied or not according to `process`. )rrh)r?r,rfr@r@rAread_var_arraykszMatFile4Reader.read_var_arrayNcCst|tr|g}n|dk r"t|}|jd|i}|s|\}}|jdkr\dn |j d}|dk r||kr|j|q:| |||<|j||dk r:| |t |dkr:qq:|S)a, get variables from stream as dictionary Parameters ---------- variable_names : None or str or sequence of str, optional variable name, or sequence of variable names to get from Mat file / file stream. If None, then get all variables in file. NrNonelatin1) isinstancestrlistrHrr end_of_streamrr:r{rremovern)r?Zvariable_namesmdictrerr:r@r@rA get_variables}s(       zMatFile4Reader.get_variablescCs|jd|g}|s~|\}}|jdkr:dn |jd}|j|}t |j d}| |||f|j|q|S)z list variables from stream rNrrunknown) rHrrrrr:r{rrrgetr<append)r?varsrerr:rjinfor@r@rAlist_variabless   zMatFile4Reader.list_variables)T)N) rCrDrErrrBrrrrrr __classcell__r@r@rrAr 4s    rowcCs(t||}t|dkrtd||S)a Make ``arr`` exactly two dimensional If `arr` has more than 2 dimensions, raise a ValueError Parameters ---------- arr : array oned_as : {'row', 'column'}, optional Whether to reshape 1-D vectors as row vectors or column vectors. See documentation for ``matdims`` for more detail Returns ------- arr2d : array 2-D version of the array r"z=Matlab 4 files cannot save arrays with more than 2 dimensions)rrnrTZreshape)rgoned_asr=r@r@rArs  c@sTeZdZddZddZddZeedfdd Zd d Z d d Z ddZ ddZ dS)rcCs|j|_|j|_dSr9) file_streamr)r?Z file_writerr@r@rArBszVarWriter4.__init__cCs|j|jdddS)Nri)rl)rwriterz)r?rgr@r@rA write_bytesszVarWriter4.write_bytescCs|j|dSr9)rr)r?sr@r@rA write_stringszVarWriter4.write_stringrc Cstdtd}t }d}|d|d|d||d<|d|d<|d |d <||d <t|d |d <|||d } || ddS)a) Write header for given data options Parameters ---------- name : str name of variable shape : sequence Shape of array as it will be read in matlab P : int, optional code for mat4 data type, one of ``miDOUBLE, miSINGLE, miINT32, miINT16, miUINT16, miUINT8`` T : int, optional code for mat4 matrix class, one of ``mxFULL_CLASS, mxCHAR_CLASS, mxSPARSE_CLASS`` imagf : int, optional flag indicating complex r@r,rrNrOrPr-r.rr/r0r1rN)roemptyrrrnrrencode) r?r:rjr]r^r0r,rZr\rYr@r@rA write_headers"   zVarWriter4.write_headercCstj|r|||dSt|}|j}|jsB|| d}|j }|tj kr\t dn4|tj krpt dn |tjtjfkr|||dS|||dS)z Write matrix `arr`, with name `name` Parameters ---------- arr : array_like array to write name : str name in matlab workspace N=z!Cannot save object arrays in Mat4zCannot save void type arrays)rr8issparse write_sparseroZasarrayr;ZisnativeryZ newbyteordertypeZobject_rdvoidZunicode_Zstring_ write_char write_numeric)r?rgr:rrZdttr@r@rArs        zVarWriter4.writecCst||j}|jjdk}zt|jjdd}Wn2tk rb|rP|d}n |d}t}YnX|j ||j |t |d|r| |j | |jn | |dS)NcrZc128r&)r]r^r0)rrr;kindrrKeyErrorryrrrjrrrealr)r?rgr:r0r]r@r@rArs&      zVarWriter4.write_numericcCst|}t||j}|j}|j||ttd|jjdkrvt |}t j dt |||d}| d}t j |d|d}||dS)Nr]r^Ur@rxrwr3)rrrrjrrrr;rroprodrpr itemrr)r?rgr:r=Zn_charsZst_arrstr@r@rAr$s$   zVarWriter4.write_charcCs|}|jjdk}tj|jdd|fdd}|j|dddf<|j|dddf<|dddd fd7<|r|jj |ddd f<|jj |dddf<n|j|ddd f<|j |ddd f<|j ||j t td ||dS) zY Sparse matrices are 2-D See docstring for VarReader4.read_sparse_array rrr#r&rNr~rr"r)Ztocoor;rroZzerosZnnzrcolrYrrrjrrrr)r?rgr:Ar0Zijvr@r@rAr7s$ zVarWriter4.write_sparseN) rCrDrErBrrrrrrrrrr@r@r@rArs!c@s$eZdZdZdddZdddZdS) rz) Class for writing matlab 4 format files NcCs"||_|dkrd}||_d|_dS)Nr)rr_matrix_writer)r?rrr@r@rArBRs zMatFile4Writer.__init__cCs.t||_|D]\}}|j||qdS)a Write variables in `mdict` to stream Parameters ---------- mdict : mapping mapping with method ``items`` return name, contents pairs where ``name`` which will appeak in the matlab workspace in file load, and ``contents`` is something writeable to a matlab file, such as a NumPy array. write_header : {None, True, False} If True, then write the matlab file header before writing the variables. If None (the default) then write the file header if we are at position 0 in the stream. By setting False here, and setting the stream position to the end of the file, you can append variables to a matlab file N)rritemsr)r?rrr:varr@r@rA put_variablesYs zMatFile4Writer.put_variables)N)N)rCrDrErrBrr@r@r@rArPs )r)*rsysrVnumpyroZ scipy.sparserZ_miobaserrrrrrr Z _mio_utilsr r functoolsr __all__ byteorderrrrrrrrrrrrrr rrrr rrrr@r@r@rAs$   Oy