a a5O@sLdZddlZddlZddlZddlZddlmZm Z m Z m Z m Z m Z mZddlmZmZddlmZejdkZdZdZdZd Zd Zd Zed ed ededededdgdddiZeeeeeeeeeeed ZdZdZdZ ddddddZ!edede diZ"Gdd d Z#Gd!d"d"Z$Gd#d$d$eZ%d,d&d'Z&Gd(d)d)Z'Gd*d+d+Z(dS)-z1 Classes for read / write of matlab (TM) 4 files N) MatFileReader docfillermatdims read_dtypeconvert_dtypes arr_to_charsarr_dtype_number)squeeze_elementchars_to_strings)reducelittlef8f4i4i2u2u1header))moptr)mrowsr)ncolsr)imagfr)namlenrU1) rZc32Zc24Zc16rZc8rrrrS1<>z VAX D-floatz VAX G-floatZCray)rrrrrdoublecharsparsec@seZdZdZdZddZdS) VarHeader4FcCs"||_||_||_||_||_dSN)namedtypemclassdims is_complex)selfr'r(r)r*r+r-d/Users/vegardjervell/Documents/master/model/venv/lib/python3.9/site-packages/scipy/io/matlab/mio4.py__init__Ps zVarHeader4.__init__N)__name__ __module__ __qualname__Z is_logical is_globalr/r-r-r-r.r%Ksr%c@sTeZdZdZddZddZdddZdd d Zd d Zd dZ ddZ ddZ dS) VarReader4z" Class to read matlab 4 variables cCs*||_|j|_|j|_|j|_|j|_dSr&) file_reader mat_streamdtypeschars_as_strings squeeze_me)r,r5r-r-r.r/`s zVarReader4.__init__c Cst|j|jd}|jt|dd}|ddksD|ddkrLtdt|dd\}}|d vrzt d t |t t|d \}}|dkrtd t|d \}}|}|d|df}|ddk} |j|} t || ||| S)z% Read and return header for variable rrrrz.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? rrrr) rr6r7readintstrip ValueErrordivmodwarningswarn order_codes UserWarningr%) r,datar'MrestOPTr*r+r(r-r-r. read_headergs4  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) mxFULL_CLASSread_full_array mxCHAR_CLASSread_char_arrayr8r mxSPARSE_CLASSread_sparse_array TypeErrorr9r )r,hdrprocessr)arrr-r-r.array_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*itemsizer6r?r@lenrBr'npndarraycopy) r,rWrcdtr* num_bytesdr]rYr-r-r.read_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)rcy?)r+rg)r,rWresZres_jr-r-r.rQs  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-1rr\r(r]) rgastyperaZuint8tobytesdecoderbr*r(arrayrc)r,rWrYSr-r-r.rSs 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)rpr)rprrZintcr(rrrfloatcomplex) rgr@raZascontiguousarrayr\imagscipyr$Z coo_matrix)r,rWrhtmpr*IJVr-r-r.rUs 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. Nrprrrr-rjrOcSsg|]}|dkr|qS)rr-).0xr-r-r. 'z0VarReader4.shape_from_header..)r)rPtuplemapr@r*rRr8rTr(r`r6seekr_rarbr?rVr9)r,rWr)r\rdr*rowscolsr-r-r.shape_from_headers2 $   zVarReader4.shape_from_headerN)T)T) r0r1r2__doc__r/rNrZrgrQrSrUrr-r-r-r.r4]s  $*r4csXeZdZdZefddZddZddZdd Zdd d Z dddZ ddZ Z S)MatFile4Readerz Reader for Mat4 files cs$tj|g|Ri|d|_dS)zR Initialize matlab 4 file reader %(matstream_arg)s %(load_args)s N)superr/_matrix_reader)r,r6argskwargs __class__r-r.r/-szMatFile4Reader.__init__cCs^|jdt|jtd}|jd|dkr6dS|dksF|dkrRtrNdpPdStrZdp\dS)Nrrr r;r!)r6rrrar(SYS_LITTLE_ENDIAN)r,rr-r-r.guess_byte_order7s   zMatFile4Reader.guess_byte_ordercCstt|j|_t||_dS)za Run when beginning read of variables Sets up readers from parameters in `self` N)rmdtypes_templateZ byte_orderr7r4r)r,r-r-r.initialize_readCszMatFile4Reader.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||Sr&r-)r{yr-r-r.[r}z0MatFile4Reader.read_var_header..rr) rrNr r*r(r_r+r)rTr6tell)r,rWnZremaining_bytes next_positionr-r-r.read_var_headerKs  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`. )rrZ)r,rrXr-r-r.read_var_arraybszMatFile4Reader.read_var_arrayNcCst|tr|g}n|dur"t|}|jd|i}|s|\}}|jdur\dn |j d}|dur||vr|j|q:| |||<|j||dur:| |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) isinstancestrlistr6rr end_of_streamrr'rmrremover`)r,Zvariable_namesmdictrWrr'r-r-r. get_variablests(       zMatFile4Reader.get_variablescCs|jd|g}|s~|\}}|jdur:dn |jd}|j|}t |j d}| |||f|j|q|S)z list variables from stream rNrrunknown) r6rrrrr'rmrr mclass_infogetr)append)r,varsrWrr'r\infor-r-r.list_variabless   zMatFile4Reader.list_variables)T)N) r0r1r2rrr/rrrrrr __classcell__r-r-rr.r+s    rrowcCs(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 rz=Matlab 4 files cannot save arrays with more than 2 dimensions)rr`rBZreshape)rYoned_asr*r-r-r. arr_to_2ds  rc@sTeZdZddZddZddZeedfdd Zd d Z d d Z ddZ ddZ dS) VarWriter4cCs|j|_|j|_dSr&) file_streamr)r,Z file_writerr-r-r.r/szVarWriter4.__init__cCs|j|jdddS)Nr[)r^)rwriterl)r,rYr-r-r. write_bytesszVarWriter4.write_bytescCs|j|dSr&)rr)r,sr-r-r. 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-rrr<r=r>rrrrrrrN)raemptyrrr`rrencode) r,r'r\rLrMrrrIrKrHr-r-r. write_headers"   zVarWriter4.write_headercCstj|r|||dSt|}|j}|jsB|| d}|j }|tj ur\t dn4|tj urpt dn |tjtjfvr|||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)rur$issparse write_sparseraZasarrayr(ZisnativerkZ newbyteordertypeZobject_rVvoidZunicode_Zstring_ write_char write_numeric)r,rYr'rdZdttr-r-r.rs        zVarWriter4.writecCst||j}|jjdk}zt|jjdd}Wn0ty`|rN|d}n |d}t}Yn0|j ||j |t |d|r| |j | |jn | |dS)NcrZc128r)rLrMr)rrr(kind np_to_mtypesrKeyErrorrkmiDOUBLErr\rPrrealrt)r,rYr'rrLr-r-r.rs&       zVarWriter4.write_numericcCst|}t||j}|j}|j||ttd|jjdkrvt |}t j dt |||d}| d}t j |d|d}||dS)NrLrMUr-rjrir)rrrr\rmiUINT8rRr(rraprodrbr itemrr)r,rYr'r*Zn_charsZst_arrstr-r-r.rs$   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 rrrrrqNrprrr)Ztocoor(rraZzerosZnnzrcolrHrrtr\rrrTr)r,rYr'ArZijvr-r-r.r.s$ zVarWriter4.write_sparseN) r0r1r2r/rrrrPrrrrrr-r-r-r.rs!rc@s$eZdZdZdddZdddZdS) MatFile4Writerz) Class for writing matlab 4 format files NcCs"||_|durd}||_d|_dS)Nr)rr_matrix_writer)r,rrr-r-r.r/Is 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-r. put_variablesPs zMatFile4Writer.put_variables)N)N)r0r1r2rr/rr-r-r-r.rGs r)r))rsysrDZnumpyraZ scipy.sparseruZmiobaserrrrrrr Z mio_utilsr r functoolsr byteorderrrZmiSINGLEZmiINT32ZmiINT16ZmiUINT16rrrrPrRrTrFrr%r4rrrrr-r-r-r.sp$   Oy