Dashboard
PHP:
8.3.33
OS:
Linux
User:
proje374
/
/
usr
/
local
/
lib64
/
python3.9
/
site-packages
/
numpy
/
f2py
/
__pycache__
📤 Upload
📝 New File
📁 New Folder
Close
Editing: crackfortran.cpython-39.pyc
a C!�i7B � @ s d Z ddlZddlZddlZddlZddlZddlZddlZddlZddl m Z z ddl Z W n e yr dZ Y n0 ddl mZ ddlT ddl mZ ejadadadadad ad ad adadada g a d a d a da i a!d a"da#da$da%t%g ia&d a'g a(d a)g a*da+d a,g a-g a.g a/d d � Z0�ddd�Z1de_2i Z3dD ]Z4ddie3e4< �q4dD ]Z4ddie3e4< �qLi Z5i Z6dD ] Z7e7d e5e7< e7e6e7d < �qldd� Z8dd� Z9d d � Z:d d � Z;e�<d ej=�j>Z?e�<d!ej=�j>Z@e�<d"ej=�j>ZAe�<d#ej=�jBZCg d$�ZDg d%�ZEd&d'� ZFd(d)� ZGeHdfd*d+�ZId,ZJd-ZKe�<eJd eKeKd.f ej=�d/fZLe�<eJd eKd0 eKd0 d.f ej=�ZMe�<eJd1 ej=�d2fZNe�<eJd3 ej=�d2fZOd4ZPe�<eJd ePePd.f ej=�d2fZQePd5 d6 ZRe�<eJd eReRd.f ej=�d2fZSd7ZTe�<eJd eTeTd.f ej=�d8fZUd9ZVe�<eJd:eVeVd.f ej=�d;fZWd<ZXe�<eJd eXeXd.f ej=�d=fZYe�<eJd> ej=�d?fZZe�<eJd@ ej=�dAfZ[e�<eJdB ej=�dCfZ\e�<eJdD ej=�dEfZ]e�<eJdF ej=�dGfZ^e�<eJdH ej=�dIfZ_e�<eJdJ ej=�dKfZ`e�<eJdL ej=�dMfZae�<eJdN ej=�dOfZbe�<eJdP ej=�dQfZce�<eJdR ej=�dSfZde�<eJdT ej=�dUfZee�<eJdV ej=�dWfZfe�<eJdX ej=�dYfZge�<eJdZ ej=�d[fZhe�<eJd\ ej=�d]fZie�<eJd^ ej=�d_fZje�<eJd` ej=�dafZke�<eJdb ej=ejlB �dcfZme�<ddejl�defZndfdg� Zodhdi� Zpe�<djej=�Zqe�<dkej=�Zre�<dlej=�Zs�ddmdn�Ztdodp� Zu�ddrds�Zvdtdu� Zw�ddvdw�Zxe�<dxej=�Zye�<dyej=�Zze�<dzej=�Z{e�<d{ej=�Z|e�<d|ej=�Z}e�<d}�Z~e�<d~�Ze�<dej=�Z�d�d�� Z�d�d�� Z�d�d�� Z�d�d�� Z�d�d�� Z�d�d�� Z�d�d�� Z�e�<d�ej=�Z�e�<d�ej=�Z�e�<d�ej=�Z�e�<d�ej=�Z�e�<d�ej=�Z�d�d�� Z�d�d�� Z�d�d�� Z�d�d�� Z��dd�d��Z��dd�d��Z��dd�d��Z��dd�d��Z�d�d�� Z�d�d�� Z��dd�d��Z��dd�d��Z��dd�d��Z�d�d�� Z�d�d�� Z��dd�d��Z�d�d�� Z��dd�d��Z�e�<d�ej=�Z�d�d�� Z�e�<d�ej=�Z�d�d�� Z�d�d�� Z�d�d�� Z�d�dÄ Z�d�dń Z��dd�dDŽZ�i fd�dɄZ�d�d˄ Z�e�<d̡jBZ�d�d΄ Z�d�dЄ Z�e�<d�ej=�Z��d d�dӄZ�d�dՄ Z�g fd�dׄZ�d�dل Z�e�<d�ej=�Z�e�<d�ej=�Z�e�<d�ej=�Z�e�<d�ej=�Z�e�<d�ej=�Z�d�d� Z�i fd�d�Z��d d�d�Z��d d�d�Z��d d�d�Z�d�d� Z��d d�d�Z�g a�d�d� Z�d�d� Z�d�d� Z�g dfd�d��Z�d�d�� Z�t���eá e�d�k� rg Z�g Z�dZ�dZ�dZ�dZ�ej�dd� D �]�Z�e�d k�rҐq�e�d d�k�r�dZȐq�e�d�k�r�dada�q�e�d�k� rd�ada�q�e�d�k� r:t� r.e1d�d� dada�q�e͐d k� rNda�q�e͐dk� rbda�q�e͐dk� rxdadan�e͐dk� r�da�dadan�e͐dk� r�dZ�n�e͐dk� r�dZ�n�e͐dk� r�dZ�n�e�d �dk� r�eΐd e�e̓ � n�eɐ r�dZ�e�anveʐ rdZ�e�anfeȐ rjze�e̓�ѡ eƠ�e͡ W n8 eҐ yf Z� z eΐd eӛd�� W Y dZ�[�n dZ�[�0 0 n eǠ�e͡ �q�t� s�t� r�t� s�e1�d d� e�eƃZ�t� r�e1�d e�t� d� e�eԃZ�e�t�d ��Z�eȠ�eա W d � n1 � s�0 Y eː reHeԃ dS (! a crackfortran --- read fortran (77,90) code and extract declaration information. Copyright 1999 -- 2011 Pearu Peterson all rights reserved. Copyright 2011 -- present NumPy Developers. Permission to use, modify, and distribute this software is given under the terms of the NumPy License. NO WARRANTY IS EXPRESSED OR IMPLIED. USE AT YOUR OWN RISK. Usage of crackfortran: ====================== Command line keys: -quiet,-verbose,-fix,-f77,-f90,-show,-h <pyffilename> -m <module name for f77 routines>,--ignore-contains Functions: crackfortran, crack2fortran The following Fortran statements/constructions are supported (or will be if needed): block data,byte,call,character,common,complex,contains,data, dimension,double complex,double precision,end,external,function, implicit,integer,intent,interface,intrinsic, logical,module,optional,parameter,private,public, program,real,(sequence?),subroutine,type,use,virtual, include,pythonmodule Note: 'virtual' is mapped to 'dimension'. Note: 'implicit integer (z) static (z)' is 'implicit static (z)' (this is minor bug). Note: code after 'contains' will be ignored until its scope ends. Note: 'common' statement is extended: dimensions are moved to variable definitions Note: f2py directive: <commentchar>f2py<line> is read as <line> Note: pythonmodule is introduced to represent Python module Usage: `postlist=crackfortran(files)` `postlist` contains declaration information read from the list of files `files`. `crack2fortran(postlist)` returns a fortran code to be saved to pyf-file `postlist` has the following structure: *** it is a list of dictionaries containing `blocks': B = {'block','body','vars','parent_block'[,'name','prefix','args','result', 'implicit','externals','interfaced','common','sortvars', 'commonvars','note']} B['block'] = 'interface' | 'function' | 'subroutine' | 'module' | 'program' | 'block data' | 'type' | 'pythonmodule' | 'abstract interface' B['body'] --- list containing `subblocks' with the same structure as `blocks' B['parent_block'] --- dictionary of a parent block: C['body'][<index>]['parent_block'] is C B['vars'] --- dictionary of variable definitions B['sortvars'] --- dictionary of variable definitions sorted by dependence (independent first) B['name'] --- name of the block (not if B['block']=='interface') B['prefix'] --- prefix string (only if B['block']=='function') B['args'] --- list of argument names if B['block']== 'function' | 'subroutine' B['result'] --- name of the return value (only if B['block']=='function') B['implicit'] --- dictionary {'a':<variable definition>,'b':...} | None B['externals'] --- list of variables being external B['interfaced'] --- list of variables being external and defined B['common'] --- dictionary of common blocks (list of objects) B['commonvars'] --- list of variables used in common blocks (dimensions are moved to variable definitions) B['from'] --- string showing the 'parents' of the current block B['use'] --- dictionary of modules used in current block: {<modulename>:{['only':<0|1>],['map':{<local_name1>:<use_name1>,...}]}} B['note'] --- list of LaTeX comments on the block B['f2pyenhancements'] --- optional dictionary {'threadsafe':'','fortranname':<name>, 'callstatement':<C-expr>|<multi-line block>, 'callprotoargument':<C-expr-list>, 'usercode':<multi-line block>|<list of multi-line blocks>, 'pymethoddef:<multi-line block>' } B['entry'] --- dictionary {entryname:argslist,..} B['varnames'] --- list of variable names given in the order of reading the Fortran code, useful for derived types. B['saved_interface'] --- a string of scanned routine signature, defines explicit interface *** Variable definition is a dictionary D = B['vars'][<variable name>] = {'typespec'[,'attrspec','kindselector','charselector','=','typename']} D['typespec'] = 'byte' | 'character' | 'complex' | 'double complex' | 'double precision' | 'integer' | 'logical' | 'real' | 'type' D['attrspec'] --- list of attributes (e.g. 'dimension(<arrayspec>)', 'external','intent(in|out|inout|hide|c|callback|cache|aligned4|aligned8|aligned16)', 'optional','required', etc) K = D['kindselector'] = {['*','kind']} (only if D['typespec'] = 'complex' | 'integer' | 'logical' | 'real' ) C = D['charselector'] = {['*','len','kind','f2py_len']} (only if D['typespec']=='character') D['='] --- initialization expression string D['typename'] --- name of the type if D['typespec']=='type' D['dimension'] --- list of dimension bounds D['intent'] --- list of intent specifications D['depend'] --- list of variable names on which current variable depends on D['check'] --- list of C-expressions; if C-expr returns zero, exception is raised D['note'] --- list of LaTeX comments on the variable *** Meaning of kind/char selectors (few examples): D['typespec>']*K['*'] D['typespec'](kind=K['kind']) character*C['*'] character(len=C['len'],kind=C['kind'], f2py_len=C['f2py_len']) (see also fortran type declaration statement formats below) Fortran 90 type declaration statement format (F77 is subset of F90) ==================================================================== (Main source: IBM XL Fortran 5.1 Language Reference Manual) type declaration = <typespec> [[<attrspec>]::] <entitydecl> <typespec> = byte | character[<charselector>] | complex[<kindselector>] | double complex | double precision | integer[<kindselector>] | logical[<kindselector>] | real[<kindselector>] | type(<typename>) <charselector> = * <charlen> | ([len=]<len>[,[kind=]<kind>]) | (kind=<kind>[,len=<len>]) <kindselector> = * <intlen> | ([kind=]<kind>) <attrspec> = comma separated list of attributes. Only the following attributes are used in building up the interface: external (parameter --- affects '=' key) optional intent Other attributes are ignored. <intentspec> = in | out | inout <arrayspec> = comma separated list of dimension bounds. <entitydecl> = <name> [[*<charlen>][(<arrayspec>)] | [(<arrayspec>)]*<charlen>] [/<init_expr>/ | =<init_expr>] [,<entitydecl>] In addition, the following attributes are used: check,depend,note TODO: * Apply 'parameter' attribute (e.g. 'integer parameter :: i=2' 'real x(i)' -> 'real x(2)') The above may be solved by creating appropriate preprocessor program, for example. � N)�Path� )� __version__)�*)�symbolic�fix� � ���c C sp da dadadadadadadadada g a da t g ia da dadag ai adadadag ag ag ag ad ad S )Nr r r r r r )� strictf77�sourcecodeform�quiet�verbose�tabchar� pyffilename� f77modulename� skipemptyends�ignorecontains� dolowercase�debug� groupcounter� grouplist� neededmodule� expectbegin�skipblocksuntil� usermodules� f90modulevars� gotnextfile�filepositiontext�currentfilename� skipfunctions� skipfuncs� onlyfuncs� include_paths�previous_context� r% r% �C/usr/local/lib64/python3.9/site-packages/numpy/f2py/crackfortran.py�reset_global_f2py_vars� s4 r' c C s, t sd S ts(|r tj�t� tj�| � d S �N)r r �sys�stdout�writer )�line�flagr% r% r& �outmess� s r. �2 zabcdefghopqrstuvwxyz$_�typespec�realZijklmn�integer)3�int�double�float�char�short�long�void�case�while�return�signedZunsigned�if�forZtypedef�sizeof�union�struct�static�register�new�breakZdo�goto�switch�continue�else�inlineZextern�delete�const�auto�lenZrank�shape�index�slen�sizeZ_i�max�minZflenZfshape�stringZcomplex_doubleZ float_double�stdin�stderrr* �type�defaultZ_bnc C s( | t v r$td| t | f � t | S | S )Nz&rmbadname1: Replacing "%s" with "%s". )�badnames�errmess��namer% r% r& � rmbadname1 s �r_ c C s dd� | D �S )Nc S s g | ] }t |��qS r% )r_ ��.0�_mr% r% r& � <listcomp> � z rmbadname.<locals>.<listcomp>r% ��namesr% r% r& � rmbadname s rg c C s( | t v r$td| t | f � t | S | S )Nz+undo_rmbadname1: Replacing "%s" with "%s". )� invbadnamesr\ r] r% r% r& �undo_rmbadname1 s �ri c C s dd� | D �S )Nc S s g | ] }t |��qS r% )ri r` r% r% r& rc rd z"undo_rmbadname.<locals>.<listcomp>r% re r% r% r& �undo_rmbadname s rj z-\*-\s*fortran\s*-\*-z-\*-\s*f90\s*-\*-z-\*-\s*fix\s*-\*-z[^c*]\s*[^\s\d\t])z.f90z.f95z.f03z.f08)z.forz.ftnz.f77z.fc C s� t durt �| ��� j}n�tdtj�| ��}t| d��`}|� |�}|� t j �rTd}n4|� t j t jf�rld}n |� t jt jf�r�d}nd}W d � n1 s�0 Y t| ||d�S ) aN Ensures that filename is opened with correct encoding parameter. This function uses charset_normalizer package, when available, for determining the encoding of the file to be opened. When charset_normalizer is not available, the function detects only UTF encodings, otherwise, ASCII encoding is used as fallback. N� �rbz UTF-8-SIGzUTF-32zUTF-16�ascii)�encoding)�charset_normalizer� from_path�bestrn rU �os�path�getsize�open�read� startswith�codecs�BOM_UTF8� BOM_UTF32_LE� BOM_UTF32_BE�BOM_LE�BOM_BE)�filename�modern �nbytes�fhandle�rawr% r% r& �openhook- s "r� c C s� d}t | �j�� tv rd}t| d���}|�� }d}t|�r@d}nt|�rPd}d}|dkr�|r�|d dkr�|�� r�|d8 }|d dkr�t |d d � �s�|d d � d kr�d}q�|�� }qPW d � n1 s�0 Y |S )z(Check if file is in free format Fortran.FT�r� r �!r � N� ���r �&) r �suffix�lower�COMMON_FREE_EXTENSIONSr� �readline� _has_f_header�_has_f90_header�strip�_free_f90_start)�fname�resultr� r, �nr% r% r& �is_free_formatL s& ,(r� c C s` |st ttttttttf }| g kr&dS t}d}d}d}t � dt j �}t � d�} t � d�} t � d�} |rp|dd� d \}} d gd d � t d �D � } dat j| td�}z |�� }W nH ty� } z0td|�� � d|�� � d|� d���W Y d}~n d}~0 0 |�s�q@|�� �r�da|�� ada |} dadatj�t�d }tt�j�� t v �rft |��sft |��sfdant t��r~t |��s~dat�r�t ant!at"dt#t�tt�r�d�p�df � |�$� �%dd �}|dk�s�|d dv�rܐq�|dd� }�q�t&|d�\}}|d 7 }|dd � �� d k�r6t&|d |d d� d�\}}|�'� dk�rVtdk�rPq�d}q�tdk�r�|d d v �r�|dd � �� d k�r�d!|d d� }nd}q�n t�r�t(|�d"k�r�|dd"� }|d | v�r�td#t#|� ��|�r�t�r t(|�d k�r |d d k�s ||d$d� }d}d}nj| �)|�}|�r8|�*d%�}|�r\|| �)|��*d%� }d}d}n,d!|d d� }|�r||�� }n|}|}|}n�tdk�rZ|�s�|d&k�r�| �)|��r�|d' }|�� }|�s�t+d(� �q�|| }| �)|��r��q�q�|�,� }| �)|�}|�r|�*d%�}|�r4|| �)|��*d%� }d}d}n |�rD|�� }n|}|}|}|du}nt-d)t#t� ��d*|�.� d t| f a|�)|�}|�r2|�*d+�}tj�/|��r�t0||dd,� nxtj�1t�gt2 }d}|D ]:}tj�3||�} tj�/| ��r�d}t0| |dd,� �q�q�|�s:t"d-t#|�tj4�3|�f � n||� |} q�|�rP|�� }n|}|}d*|�.� d t| f a|�)|�}|�r |�*d+�}tj�/|��r�t0||dd,� nxtj�1t�gt2 }d}|D ]:}tj�3||�} tj�/| ��r�d}t0| |dd,� �q��q�|�s(t"d-t#|�tj4�3|�f � n||� da|�5� |�rF|dd� n|\ a aaaaaaaadS ).z� Read fortran codes from files and 1) Get rid of comments, line continuations, and empty lines; lower cases. 2) Call dowithline(line) on every line. 3) Recursively call itself when statement "include '<filename>'" is met. NFr z*\s*include\s*(\'|")(?P<name>[^\'"]*)(\'|")z(?P<line>.*)&\s*\Zz(\s*&|)(?P<line>.*)z.*?'''r )r r � c S s g | ] }t |��qS r% )�strr` r% r% r& rc � rd z#readfortrancode.<locals>.<listcomp>� )r� zreadfortrancode: reading �#z failed with zu. It is likely that installing charset_normalizer package will help f2py determine the input file encoding correctly.r r r �freez Reading file %s (format:%s%s) z,strict� z r� r� z!f2pyr )r �cr� �Cr� Zf2pyz �H zxreadfortrancode: Found non-(space,digit) char in the first column. Are you sure that this code is in fix form? line=%s� r, �.pyf� z.Unexpected end of file when reading multiline z6Flag sourcecodeform must be either 'fix' or 'free': %szLine #%d in %s:"%s" r^ )� dowithline�istopzAreadfortrancode: could not find include file %s in %s. Ignoring. )6r r r r r � beginpatternr r r �re�compile�I�range� fileinputZ FileInputr� r� �UnicodeDecodeError� Exceptionr~ �linenoZ isfirstlinerr rs �splitextr r� r� �COMMON_FIXED_EXTENSIONSr� �_has_fix_headerr� �beginpattern77�beginpattern90r. �repr� expandtabs�replace�split_by_unquotedr� rO �match�groupr\ �rstrip� ValueErrorZ filelineno�isfile�readfortrancode�dirnamer# �join�pathsep�close) �ffiler� r� Z saveglobalsZlocaldolowercase�contZ finalline�llZ includelineZcont1Zcont2Z mline_mark�l1Z spacedigitsZfin�l�msg�ext�rl�_Z origfinalliner� �lc�m�fn� include_dirsZ foundfile�inc_dirZfn1r% r% r& r� f sv �� �� ���� �( � �� �� � � � � �r� zN\s*(?P<before>%s(?=\s*(\b(%s)\b)))\s*(?P<this>(\b(%s)\b))\s*(?P<after>%s)\s*\Zzqcharacter|logical|integer|real|complex|double\s*(precision\s*(complex|)|complex)|type(?=\s*\([\w\s,=(*)]*\))|byte�.*rY z|static|automatic|undefined)z([a-z]+[\w\s(=*+-/)]*?|)�functionr� r� �begin)z [a-z\s]*?� subroutiner� r� zprogram|block\s*dataz@|module(?!\s*procedure)|python\s*module|(abstract|)\s*interface|z type(?!\s*\()z\end|endprogram|endblockdata|endmodule|endpythonmodule|endinterface|endsubroutine|endfunction�endzDend\s*(if|do|where|select|while|forall|associate|critical|enum|team)z[\w]*?�endifzmodule\s*procedure�moduleprocedure)r �implicitr� r� r� )r �dimension|virtualr� r� � dimension)r �externalr� r� r� )r �optionalr� r� r� )r �requiredr� r� r� )r �publicr� r� r� )r �privater� r� r� )r � intrinsicr� r� r� )r �intent|depend|note|checkr� z \s*\(.*?\).*�intent)r � parameterr� z\s*\(.*r� )r �datar� r� r� )r �callr� r� r� )r �entryr� r� r� )r �callfunr� r� r� )r �commonr� r� r� )r �user� r� r� )r �containsr� r r� )r �formatr� r� r� )r �Kthreadsafe|fortranname|callstatement|callprotoargument|usercode|pymethoddefr� r� �f2pyenhancementsz2\s*(?P<before>''')(?P<this>.*?)(?P<after>''')\s*\Z� multilinec C sv t d�t |�@ rJ d��t�djd�t�|��d�t�|��ddd��}|�| �}|rn|�� }|d |d fS | d fS ) z� Splits the line into (line[:i], line[i:]), where i is the index of first occurrence of one of the characters not within quotes, or len(line) if no such index exists z"'z cannot split by unquoted quoteszR\A(?P<before>({single_quoted}|{double_quoted}|{not_quoted})*)(?P<after>{char}.*)\Zz[^"'{}]z[{}]z('([^'\\]|(\\.))*')z("([^"\\]|(\\.))*"))Z not_quotedr6 � single_quotedZ double_quoted�before�afterr )�setr� r� r� �escaper� � groupdict)r, � charactersr� r� �dr% r% r& r� � s �� r� c C sB g }t | ��d�D ]$}dD ]}|�|d�}q|�|� qd�|�S )N�@,@z(),r� �,)�markoutercomma�splitr� �appendr� )Zargsline�ar� r� r% r% r& � _simplifyargs� s r� z"\s*(?P<result>\b[a-z]+\w*\b)\s*=.*z \s*(?P<bind>\b[a-z]+\w*\b)\s*=.*zH\s*bind\(\s*(?P<lang>[^,]+)\s*,\s*name\s*=\s*"(?P<lang_name>[^"]+)"\s*\)c C s� t | d�\}}|r�td �| �s�td �| �s�|dksBJ t|���t | d�\} }|rvt| |� t |dd� d�\} }qPt| |� dS |dk r�datdiati iatg ia g tt d< i tt d< dtt d< dtt d < d a d a dS |dk�r d}t �rt tk�rd }t|k�rpt d tt�tt�f � t d � t td �tt � t t t td d d< t t= td a�qt �r�t tk�r�t td �tt � t t t td d d< t t= td at td �tt � t t t td d d< t t= td ad a dS | dk�rdS d}tttttttttttttt t t t t t!t"t#t$t%ttt&fD ](}|d �| �}|�rh �qt|d }�qJ|�s�t'} dt k�r�tk�r�n ndS dtt v �r�tt d D ]�} | t(v �r�t(| } dtt v �r�| tt d v �r�q�t)�d| t*| �t)j+�} | �r�| �| �,d��} t-| �,d��} | �rJd| | | �,d�f } n d| | f } t.d �| �}|�s�t dt| � � dS t/|d| � dS �q�t0dk�s�t0dk�r�t1�2� �3d��r�da4t dt � dS |d dk�r�dt k�r�tk �rn ntd at tk�rdS tdk�r$t5dt ��t d �| �} | �rp| �,d �tt k�spt5d t| �,d ��ttt �t6f ��t tk�r~d a t td �tt � t t t td d d< t t= td at7�s�da8n�|d d k�rdt k�r�tk�r�n n td adS da9t/||d | � da8n�|d d k�r&n�|d d k�rFt/||d | � nn|d d!k�r�t:�r^dS dt k�rvtk�r~n ndS ta n0dt k�r�tk�r�n ndS t/||d | � dS )"z� reset=-1 --- initialize reset=0 --- crack the line reset=1 --- final check if mismatch of blocks occurred Cracked data is saved in grouplist[0]. �;r r Nr �body�vars�blockr^ r � z(crackline: groupcounter=%s groupname=%s zYcrackline: Mismatch of blocks encountered. Trying to fix it by assuming "end" statement. � externals� interfacedz5(?P<before>[^"]*)\b%s\b\s*@\(@(?P<args>[^@]*)@\)@.*\Zr� �argszcallfun %s(%s) result (%s)r� zcallfun %s(%s)z8crackline: could not resolve function call for line=%s. r� r� z"crackline:%d: No pattern for line r� z>crackline: groupcounter(=%s) is nonpositive. Check the blocks.�thiszGcrackline: End group %s does not match with previous Begin group %s %sr� r� r� r� );r� �f2pyenhancementspatternr� �multilinepatternr� � crackliner � groupname� groupcacher r r r r. r� �dimensionpattern�externalpattern� intentpattern�optionalpattern�requiredpattern�parameterpattern� datapattern� publicpattern�privatepattern�intrinsicpattern� endifpattern� endpattern� formatpatternr� �functionpattern�subroutinepattern�implicitpattern� typespattern� commonpattern� callpattern� usepattern�containspattern� entrypattern�moduleprocedurepattern�crackline_re_1rh r� �markouterparenr� r� r� �callfunpattern� analyzeliner r r� �endswithr$ r� r r r r r )r, �resetr� Z has_semicolonZsemicolon_line�flr- �patr� �re_1r^ �m1�m2r� r% r% r& r � s& � �� � � � � $ ��� r c C sd d}d}| D ]R}|dkr4|d }|dkrV|d }q n"|dkrV|d }|dkrV|d }q || }q |S )Nr r �(r �@(@�)�@)@r% )r, r� �fr� r% r% r& r" O s r" r� c C s� d}d}t | |d �\}}||7 }|r�|d |krL|dkrL|d| d 7 }n6||d 7 }|d dkrn|d7 }n|d dkr�|d8 }t |dd � |d �\}}||7 }q"|r�J t|| |f���|S )Nr r �()�@r, r r. )r� r� )r, �commar� r0 r� r� r% r% r& r� a s r� c C s | � dd�� dd�}|S )Nr- r, r/ r. )r� )r, r� r% r% r& �unmarkouterparent s r4 c C s� | si } |s| S | |u r | S t |�� �D ]�}|dkrN|s@|| vr�|| | |< q(|dkrr|| D ]}t| ||�} q^q(|dkr�t| || |�} q(|dkr�t| || |�} q(|dv r�|s�|| vr�|| | |< q(|dkr�q(|dv r�td| � q(td t|� ��q(| S ) Nr0 �attrspec� kindselector� charselector)�=�typename�note)r� �checkr� r� r� �dependz"appenddecl: "%s" not implemented. z-appenddecl: Unknown variable definition key: )�list�keys� setattrspec�setkindselector�setcharselectorr\ r� r� )�declZdecl2�force�kr� r% r% r& � appenddecly s8 �rE zD\s*(?P<this>(@\(@.*?@\)@|\*[\d*]+|\*\s*@\(@.*?@\)@|))(?P<after>.*)\Zz[(?:,(?P<attributes>[\w(),]+))?(::)?(?P<name>\b[a-z$_][\w$]*\b)(?:\((?P<params>[\w,]*)\))?\Zz�\s*(?P<name>\b[\w$]+\b)\s*(@\(@\s*(?P<args>[\w\s,]*)\s*@\)@|)\s*((result(\s*@\(@\s*(?P<result>\b[\w$]+\b)\s*@\)@|))|(bind\s*@\(@\s*(?P<bind>(?:(?!@\)@).)*)\s*@\)@))*\s*\ZzF\s*(?P<scheme>(operator|assignment))@\(@\s*(?P<name>[^)]+)\s*@\)@\s*\Zz9\s*(?P<name>\b[\w$]+\b)\s*@\(@\s*(?P<args>.*)\s*@\)@\s*\Zz4([-+]?(?:\d+(?:\.\d*)?|\d*\.\d+))[dD]((?:[-+]?\d+)?)zA([-+]?((?:\d+(?:\.\d*)?|\d*\.\d+))[eE]((?:[-+]?\d+)?)|(\d+\.\d*))zintent\s*\(.*?\bcallback\bc C s& | � dg �D ]}t�|�r dS q dS )Nr5 r r )�get�_intentcallbackpatternr� )Zvdeclr� r% r% r& �_is_intent_callback� s rH c C sl d� | �� �} t�| �}t| |� |rb|�d�}|rHdd� |�d�D �ng }|�d�||�d�fS d g d fS )Nr � attributesc S s g | ] }|� � �qS r% )r� )ra r� r% r% r& rc � rd z*_resolvetypedefpattern.<locals>.<listcomp>r� r^ �params)r� r� �typedefpatternr� �printr� )r, r* �attrsr% r% r& �_resolvetypedefpattern� s rN c C sR t �dt j�}|�| �}d }|rJ|�d�}| d |�� � | |�� d � } | |fS )NzObind\(\s*(?P<lang>[^,]+)(?:\s*,\s*name\s*=\s*["\'](?P<name>[^"\']+)["\']\s*)?\)r )r� r� r� �searchr� �startr� )r, �patternr� Zbind_statementr% r% r& �parse_name_for_bind� s rR c C s� t | �\} }t| �} t�| �}|r@|�d�|�d�|�d�|fS t�| �}|rv|�d�d |�d� d }|g d d fS t�| �}|r�|�d�|�d�d d fS d g d d fS )Nr^ r r� �schemer, r. )rR r" �nameargspatternr� r� �operatorpattern�callnameargspattern)r, Z bind_cnamer* r^ r% r% r& �_resolvenameargspattern� s rW c G C s� | � d�}|dkrdatr�|dvr�ts�tdk r�tj�t�� d�d }t d| � da td ad t t< i t t< g tt< g t t d < i t t d < d t t d < |t t d < dt t d< da|dv �r|�� }t�d|tj�r�d}n,t�d|tj�r�d}nt�d|tj��rd}|dk�rFt| � d��\}}}t|d�t t d |< g }d} nt| � d��\}}} } |du �r�|dk�rvd}nd}|d v�r�t d � ||tfa|�r�td d � t|�� d �D ��}ng }d|v �r�d|v �r�|�d� �q�t d!� d} d} |d"v �rnd} d#t t v�rdS |t t d# v�r&dS tt D ]} | d |k�r. dS �q.|t t d$ v �r`dS d%d&d"�| }t�r�td'k�r�tdk�r�td( ad} |d)v�r�d} td ai t t< g tt< | �rPtdk�r�t d*t t� d� d+t t< dt t d < tt t d < dt t d< g t t d < g t t d,< g t t d$< i t t d < td ai t t< g tt< | �rtdk�rnt d-t d� d.t t< d.t t d < d/t t d < d0t td d t td d f t t d< g t t d < g t t d,< g t t d$< i t t d < td ai t t< g tt< |t t< |t t d < |�s,d1|� d2d3� }| � d4�t t d5< t |�t t d < | t t d6< tdk�rrtt t d< n^t�r�td7k�r�d0t td d tf t t d< n,d0t td d t td d f t t d< t t t � � �D ] }t t | �s�t t |= �q�|t t d#< g t t d < g t t d,< g t t d$< i t t d < i t t d8< |dk�r^g t t d9< |d"v �r�|t td( d, v�r�t td( d, �!|� t"�#t td( d �t t d < z,t t d | t t d | d: �$d;�= W n t%�y� Y n0 |d<v �r<| �rft�&t'| �}|�rf|i it t d=< |� d>�t t d= | d>< |� d?��rf|� d?�t t d= | d < z6t(t t d | t td( d d �t t d |< W n t%�y� Y n0 |d@k�r | �r | t t d v �r || k�s t(t t d | t t d | �t t d |< zt td( d$ �!|� W n t%�y: Y n0 |d&k�r�t)d �| � d4�d2 | �}|�r�t*|� d�|� d��\}}}}t+||||� |d"v �r�ttd �!t t � tt ttd d' d < tt= td attd �!t t � tt ttd d' d < tt= td a�n�|d8k�r�t| � d��\}}} }|du�r�|�rhtdAd � t|�� d �D ��}ng }| du �s�J t | ���|t t d8 |< d8|tfa�n8|dk�r�t*|| � d��\}}}}t+||||�}|du�r�dB|tfa�n�|dCv � r�t t d }| � d��,� }|�-dD�}|dk � r�|dEk� r�t.|��-dF�d( }|d|d � dD ||d d� }|�-dD�}||d� dDk� r�d#t t v � r�t dG| � d�|d|� f � |dH�/t t d# � }|dk � r�d}d}n |d|� �,� }||d( d� }t|�� d �}t0|�dk� r0|d }t dIdH�/|dd� � � d}dJd � t|�� d �D �D �]}t1�|�}|� s�|dKv � rrd}n$t2| �3� � t dL|t |�f � � qLnt |� d ��}|dKv � r�|dMk� s�|dNk� r�||� d�7 }||v� r�i ||< |dOk� r�||� d� } |dEk� r�| � d�| } t4�| �� r�|t t d# v� r�tdk� r�dPt td( d v� rZt dQ� |t t d k� r�t dR|t t d f � t t d# �!|� n t5dS| � n t5dT| � |dUv � r�|} d:|| v � r�|| d: �!| � n| g|| d:< |d;k� rZt t d d k� rt dV� � qL|t t d# v� r.� qLd,t t v� rHg t t d,< t t d, �!|� |}� qL|t t d < |du�r�dB|tfa� nR|dWk� r�dXd � | � d�� dH�D �t t dY< � n"|dZk�rt t d }| � d��,� dd'� }d}t|�� d �D �]�}z d[d � |� d\�D �\}} W n* t%� y@ t d]||f � Y � q�Y n0 t6|�} t |�}||v� rdi ||< d\|| v � r�|| d\ | k� s�t d^||| d\ | f � t7| | �}|�rH|�8d_�d`k�rt | �} t9�:| �D ]:} t | | �;� | �<� � �� � dadb��| | �;� | �<� �< � q�d�/| �} n0|�8d_�dck�rH| dd� �� � dadb�� dHdd�} zt=| i | �} W nJ t>t?t@f�y� }! z*t5de| |!f � W Y d}!~!� q�W Y d}!~!n d}!~!0 0 t | �|| d\< d:|| v �r�|| d: �!dZ� ndZg|| d:< |}� q�|t t d < |du�r�dB|tfa� n�|dfk�rH| � d��,� �� dgk�rBdt t df< �q�| � d��r�dft t v �rjt t df }"ni }"|"du �r�t dh� i }"t| � d��� d �D �]�}i }#t�di|tj�}|�s�t dj| � �q�tA�|� d��}$|$�s�t dk| � �q�t*|$� d�|$� d��\}}}}tB||�\}%}&}'||#d_< |%|#dl< |&|#dm< |'|#dn< t |#� � �D ]}|#| �sH|#|= �qHt|� d��� d �D ]�}(do|(v �r�z dpd � |(� do�D �\})}*W n& t%�y� t dq|( � Y �qtY n0 n |(�,� })}*t0|)�t0|*� k�r�dk�sn t dr|( � �qttCtD|)�tD|*�d �D ]}+|#|"tE|+�< �q �qt�q�|"t t df< �n�|dsk�rg }d},d}-d}.d}/d}0| � d�D ]�}1|0�s�|1dtk�r�|/ }/|1duk�r�|/�r�|.d }.�qt|1dvk�r�|0d }0n|1dwk�r�|0d }0|.dk�r�|,|1 },n\|.dk�r�|-|1 }-nH|.d(k�rt|,�,� },|,�FdH��r&|,dd� �,� },|�!|,|-g� |1},d}-d}.�qt|.d(k�r�|,�,� },|,�FdH��rr|,dd� �,� },|�!|,|-g� t t �8d i �}2d}|D �]T}3|3d �,� |3d �,� |3d< |3d< |3d �FdH��r�|3d dd� |3d< |3d �Fdv��rt dx|3d � �q�tGtdyd � t|3d �� d �D ���D �]�\}4} | �Fdv��rRt dx| � �q,dz|3d v �rtt d{|3d � �q,|2�H| i � |2| �8d_�}5tI|2| �}6|5dck�r�t�Jd||3d �n |3d � dH�}7z"|6�r�d}�Kd~�/|7��n|7|4 }8W n� tL�y� tMdd�� |7D ���r�g }9|7D ]n}:d�|:v �rnz*|:� d��\};}<|9�N|<�,� gtO|;� � W n" tP�yj |9�!|:�,� � Y n0 n|9�!|:�,� � �q|9}7|6�r�d}�Kd~�/|7��n|7|4 }8Y n0 |2| �8d\�}=|=�r�|=|8k�r�t d�| |=|8f � |8|2| d\< | }�q,�q�|2t t d < |�r�dB|tfa�n�|d�k�r�| � d��,� }|d duk�s>d�| }g }>d}.d}?d}@|D ]r}1|1duk�rl|.d }.�qR|.d7k�r�|?�,� }?|?�s�d�}?|>�!|?|@g� |.d( }.d}?d}@|.d( �r�|?|1 }?n|@|1 }@�qR|?�,� }?|?�s�d�}?|>�!|?|@g� i }Ad�t t v �rt t d� }A|>D ]\}1|1d |Av�r&g |A|1d < d�d � t|1d �� d �D �D ] }|�rB|A|1d �!|� �qB�q|At t d�< d�|?tfa�nX|d�k�r�t�d�| � d�tj�}|�r�|�3� }Bd�t t v�r�i t t d�< |� d �}i t t d� |< d}Cd�|Bv �r�|Bd� du�r�d�|Bv �r|Bd� du �rd}C|Ct t d� | d�< d�d � |Bd� � dH�D �}i }D|D ]r}3d\|3v �r�t�d�|3tj�}$|$�r�|$� d���,� |D|$� d���,� < nt d�t |3� � n|3|D|3< |Dt t d� | d�< �qNn nt2| �3� � t d�� n�|d�v �rpd�t t v�r i t t d�< t t d� }E| � d�d�k�rZd�|Ev �rZtQ|Ed� tR��rD|Ed� g|Ed�< |Ed� �!| � d�� n| � d�|E| � d�< nh|dk�r�tdu �r�t�r�t d�� dS t}FtSt |F tdd(� | � d�� n tdk�r�t2| �3� � t d�� dS )�a� Reads each line in the input file in sequence and updates global vars. Effectively reads and collects information from the input file to the global variable groupcache, a dictionary containing info about each part of the fortran module. At the end of analyzeline, information is filtered into the correct dict keys, but parameter values and dimensions are not yet interpreted. r r� N)r� r� r� rY r �.r zBanalyzeline: no group yet. Creating program group with name "%s". �programr� r� r� r^ Zfromsky�from)r� r� r� z block\s*data� block datazpython\s*module� python modulezabstract\s*interface�abstract interfacerY r� )r5 � _BLOCK_DATA_r )� interfacer[ r] z2analyzeline: No name/args pattern found for line. c S s g | ] }|� � �qS r% �r� �ra �xr% r% r& rc s �z analyzeline.<locals>.<listcomp>r� z<analyzeline: argument list is malformed (missing argument). )r� r� r r r� r� r r �r_ r] z&analyzeline: Creating module block %s �moduler zDanalyzeline: Creating additional interface block (groupcounter=%s). r_ �unknown_interfacez%s:%sZunknown_r� r� r� �prefixr� � r� �varnamesr5 r� �r� r� Zbindlang�langZ lang_namer� c S s g | ] }|� � �qS r% r` ra r% r% r& rc � s ��variable)r� r� r� r� r� r� r� r� �::r� r/ z'All arguments will have attribute %s%s r� zXanalyzeline: cannot handle multiple attributes without type specification. Ignoring %r. c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd �r� r� zEanalyzeline: no name pattern found in %s statement for %s. Skipping. �operatorZ assignmentr� �__user__z8analyzeline: missing __user__ module (could be nothing) z;analyzeline: appending intent(callback) %s to %s arguments z,analyzeline: intent(callback) %s is ignored z=analyzeline: intent(callback) %s is already in argument list )r� r� r� r� r� r� z(analyzeline: ignoring program arguments r� c S s g | ] }|� � �qS r% r` ra r% r% r& rc rd � implementedbyr� c S s g | ] }|� � �qS r% r` ra r% r% r& rc rd r8 zManalyzeline: could not extract name,expr in parameter statement "%s" of "%s" zGanalyzeline: Overwriting the value of parameter "%s" ("%s") with "%s". r0 r1 r� �e�complexz+1j*(z1analyzeline: Failed to evaluate %r. Ignoring: %s r� �nonez<analyzeline: Overwriting earlier "implicit none" statement. z9\s*(?P<this>.*?)\s*(\(\s*(?P<after>[a-z-, ]+)\s*\)\s*|)\ZzDanalyzeline: could not extract info of implicit statement part "%s" zManalyzeline: could not extract types pattern of implicit statement part "%s" r6 r7 r9 �-c S s g | ] }|� � �qS r% r` ra r% r% r& rc _ rd zZanalyzeline: expected "<char>-<char>" instead of "%s" in range list of implicit statement z^analyzeline: expected "<char>-<char>" instead of "%s" in range list of implicit statement (2) r� �'�/r, r. z>analyzeline: implied-DO list "%s" is not supported. Skipping. c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd r� z=Comment line in declaration "%s" is not supported. Skipping. z\(.*?\)z(/{}/)�, c s s | ]}d |v V qdS )r Nr% )ra r� r% r% r& � <genexpr>� rd z analyzeline.<locals>.<genexpr>r z=analyzeline: changing init expression of "%s" ("%s") to "%s" r� z//�_BLNK_c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd r� zR\A\s*(?P<name>\b\w+\b)\s*((,(\s*\bonly\b\s*:|(?P<notonly>))\s*(?P<list>.*))|)\s*\Zr= Znotonly�onlyc S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd z7\A\s*(?P<local>\b\w+\b)\s*=\s*>\s*(?P<use>\b\w+\b)\s*\Z�localz0analyzeline: Not local=>use pattern found in %s �mapz0analyzeline: Could not crack the use statement. )r� r� Zusercodez-analyzeline: No context for multiline block. z+analyzeline: No code implemented for line. )Tr� r$ r r r rr rs �basenamer r� r. r r r r r� r� r� r� rN �dictrW rg r� �remover r r r� r� r_ r= r> r� �copy�deepcopyrQ r� rO �crackline_bindlangrE r �cracktypespec0� updatevarsr� �findr" r� rO � namepatternrL r� rG r\ �get_parameters�determineexprtyperF � real16pattern�finditerrP r� �eval� SyntaxError� NameError� TypeError�typespattern4implicit� cracktypespecr� �ord�chrrw � enumerate� setdefaultZ getdimension�findallr� � IndexError�any�extendr3 r� � isinstancer� �appendmultiline)Gr� r: r, r� Znewnamer^ rM r� r r� Z bindclineZ needmoduleZ needinterface�itrD Zbindcdat�tr0 �selector�attr�edecl� last_namer� �i�pl�chrq r* ZapZinitexprrJ �tt�vr� �implrB r+ � kindselect� charselectr9 r� ZbegcZendc�o�dlZilr0 Zfc�inpr� r� r� �idx�vtypeZvdim�matches�new_valZ expanded_listr� Z multiplier�valueZ current_valZclZbn�olZ commonkey�mmZisonlyr� r� �gcr% r% r& r$ � s� ��� � � �� �� � � � � � � � � � � $ � � � �� ���� � � � � ��$ ����� � � "� ". (" $ � �� � � r$ c C s: d| vri | d< | d }||vr(g ||< || � |� d S )NZf2pymultilines�r� )r� Z context_name�mlr� r% r% r& r� * s r� c C s� d }d }t �d| t j�r d} n"t �d| t j�r4d} n | �� �� } t�t|��}|s^td� d S |�� }t |� � �D ]}t || �||< qr| dv r�|d }|d }|� d �}|d kr�|d |� �� }||d d � }| |||fS ) Nzdouble\s*complex�double complexzdouble\s*precisionzdouble precisionz>cracktypespec0: no kind/char_selector pattern found for line. )rr r2 �logicalr1 � characterrY r r� rl r r ) r� r� r� r� r� � selectpatternr" r. r� r= r> r4 r� )r0 r� r� r� r* r� rD r� r% r% r& r� 4 s0 � r� z)\s*(?P<name>\b\w+\b)\s*(?P<after>.*)\s*\ZzB\s*(\(\s*(kind\s*=)?\s*(?P<kind>.*)\s*\)|\*\s*(?P<kind2>.*?))\s*\Zz2\s*(\((?P<lenkind>.*)\)|\*\s*(?P<charlen>.*))\s*\Zz�\s*(kind\s*=\s*(?P<kind>.*?)\s*(@,@\s*len\s*=\s*(?P<len>.*)|)|(len\s*=\s*|)(?P<len2>.*?)\s*(@,@\s*(kind\s*=\s*|)(?P<kind2>.*)|(f2py_len\s*=\s*(?P<f2py_len>.*))|))\s*\Zz�\s*(@\(@\s*(?!/)\s*(?P<array>.*?)\s*@\)@\s*\*\s*(?P<len>.*?)|(\*\s*(?P<len2>.*?)|)\s*(@\(@\s*(?!/)\s*(?P<array2>.*?)\s*@\)@|))\s*(=\s*(?P<init>.*?)|(@\(@|)/\s*(?P<init2>.*?)\s*/(@\)@|)|)\s*\Zc C s� | � � } t| �dkr| S | d }tdt| �d �D ]>}| | dkrd| |d dv s2| |d dv rdq2|| | }q2|| d }|S )Nr r r� z ()[]{}=+-/* r )r� rO r� )�exprZexpr2r� r% r% r& � removespaces[ s �� r� c C s~ d}d}d}d}| D ]d}|dkr6|dv r6||7 }|}q|sF|dv rF|}||krV| }n|dkrl|rl|d7 }q||7 }|}q|S ) a( The function replace all spaces in the input variable line which are surrounded with quotation marks, with the triplet "@_@". For instance, for the input "a 'b c'" the function returns "a 'b@_@c'" Parameters ---------- line : str Returns ------- str r FN�\)r� ru �")ru r� r� �@_@r% )r, �fragmentZinsideZ current_quote�escapedr� r% r% r& �markinnerspacesj s& r� c C sB d}t | |�\}}}|r�dd� t|��d�D �}g }t�d�} |D ]D} | sLqB| �| �} | r|| �d��� } | | t| �d� } |� | � qB|}dd� t|��d�D �} g }| D ]D}dd� tt t |��d d ��d �D �D ]}|r�|� |� d d �� q�q�|D �]D}t �|�} | �s td t|� � q�t| �d��}i }|tt d v �rXtt d | �� }d|v}|�rn| |d< n*| �r�| |d k�s�td||d | f � d|v�r�t�|�|d< nv|�r(t|�� �D ]b}||d v �r|| |d | k�std|||d | || f � nt�|| �|d |< �q�d|v�rZ|�rZ|�rH||d< ntd||f � nv|�r�t|�� �D ]b}||d v �r�|| |d | k�s�td|||d | || f � nt�|| �|d |< �qld|v�r�||d< n*|�r|d |k�std||d |f � d|v�r(t�|�|d< n.|�r�|D ]"} | |d v�r2|d � | � �q2n@t�| �|d< t�|�|d< t�|�|d< ||d< t�|�|d< d|�d��p�g v �r�|tt d v �r�d tt v�r�g tt d < tt d � |� | �d ��r�t�t| �d ���}|�r�|�� }d D ]2}||d du�r ||d ||< ||d = �q t|�� �D ],}|| du�r�t|| �||< n||= �q^d |v �rPd!|v �rP|d d"k�r�|d! |d < |d!= n�| d#k�rd|v�s�|d �s�i |d< d |d v �r |d d = |d |d d$< |d = n8|d! d% |d |d!< |d = td&| || ||d! f � d |v �r�| d'v �r�d|v�sx|d �s�i |d< |d |d d$< |d = nT| d#k�r�d|v�s�|d �s�i |d< d |d v �r�|d d = |d |d d$< |d = d(|v �r:d)|v �r.|d) |d( k�s.td*||d) |d( f � n |d( |d)< d!|v �r�d+|d! }d|v�sd|d �sp|g|d< n\|d � |� |d D ]D}|dd,� d-k�r�||k�r�|d d.= td/|||f � �q�q�ntd0|| �d � � t|�� �D ]}|| �s�||= �q�|tt d |< d1tt v �r8tt d1 � |� |}q�|S )2z� Returns last_name, the variable name without special chars, parenthesis or dimension specifiers. Alters groupcache to add the name, typespec, attrspec (and possibly value) of current variable. Nc S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd z updatevars.<locals>.<listcomp>r� z(?P<start>[a-zA-Z]+)rP c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd r� �r3 z@ @r� z;updatevars: no name pattern found for entity=%s. Skipping. r^ r� r0 zIupdatevars: attempt to change the type of "%s" ("%s") to "%s". Ignoring. r6 zVupdatevars: attempt to change the kindselector "%s" of "%s" ("%s") to "%s". Ignoring. r7 zEupdatevars:%s: attempt to change empty charselector to %r. Ignoring. zVupdatevars: attempt to change the charselector "%s" of "%s" ("%s") to "%s". Ignoring. r9 zMupdatevars: attempt to change the typename of "%s" ("%s") to "%s". Ignoring. r5 r� r r r� )rO �array�init�2rO r� r r� r r� z-updatevars: "%s %s" is mapped to "%s %s(%s)" �rr r2 r� r1 r� r8 zTupdatevars: attempt to change the init expression of "%s" ("%s") to "%s". Ignoring. � dimension(%s)� r� r z5updatevars:%s: attempt to change %r to %r. Ignoring. z?updatevars: could not crack entity declaration "%s". Ignoring. rh ) r� r� r� r� r� r� r� r� rO r� r� r� r� r� r. r� r_ r r r� r= r> r\ rF �lenarraypatternr" r� r4 )r0 r� r5 Z entitydeclr� r� r� r9 r� r� r� r� �s�elZel1rq �e1Zenamer� Znot_has_typespecrD r* �d1�lk�dmZdm1r% r% r& r� � s0 ( � � $� �$� � ( � � � � r� c C s d }d }d }|�r| dv r�t �|�}|s<tdt|� � d S |�� }|d |d< |d= t|�� �D ]}|| sb||= qbt|�� �D ]\}}t|�||< q��nf| dk�r�t �|�}|s�tdt|� � d S |�� }|d |d< |d= |d �r`t �t |d ��}|�� }d D ]:}||d �r,||d ||< || ||< ||d = �q |d d u�r`|d |d < |d= t|�� �D ]}|| �sr||= �qrt|�� �D ]\}}t|�||< �q�nR| d k�r�t �d |t j �}|�r�|�d�}ntdt| | � � ntdt|� � |||fS )Nr� z4cracktypespec: no kindselector pattern found for %s Zkind2r r� z4cracktypespec: no charselector pattern found for %s Zcharlen�lenkind)rO �kindr� Zf2py_lenrY z\s*\(\s*(?P<name>\w+)\s*\)r^ z'cracktypespec: no typename found in %s z'cracktypespec: no selector used for %s )r6 r� r. r� r� r= r> �itemsr_ r7 �lenkindpatternr� r� r� r� ) r0 r� r� r� r9 rD r� r� r� r% r% r& r� 8 st � � � ��r� c C s� | si } |s| S d| vr&|g| d< | S |r8| d � |� || d v rH| S |dkrld| d vrl| d � |� nz|dkr�d| d vr�| d � |� nV|dkr�d| d vr�| d � |� n2|dkr�d| d vr�| d � |� n| d � |� | S )Nr5 rC � automaticr� r� r� )rB r� rC r% r% r& r? t s. r? c C sZ | si } |s| S d| vr$|| d< | S t |�� �D ]$}|sD|| d vr0|| | d |< q0| S )Nr6 �r= r> �rB �selrC rD r% r% r& r@ � s r@ c C sZ | si } |s| S d| vr$|| d< | S t |�� �D ]$}|sD|| d vr0|| | d |< q0| S )Nr7 r� r� r% r% r& rA � s rA �unknownc C s d| v r| d S |S )Nr^ r% )r� r� r% r% r& � getblockname� s r� c C s0 zd| d | d f a W n ty* Y n0 d S )Nz In: %s:%s rZ r^ )r r� )r� r% r% r& � setmesstext� s r� c C s2 i }d| v rt | d �}d| v r.|�| d � |S )N� parent_blockr� )� get_usedict�update)r� �usedictr% r% r& r� � s r� c C s� |d u r i }t | �}|s |S t|�� �D ]�\}}|�� }|tvrXtd|| �d�f � q(t| }t|�}|snq(|r~td| � t|�� �D ]0\}}||v r�tdt |�t |�f � |||< q�q(|S )Nz0get_useparameters: no module %s info used by %s r^ z,get_useparameters: mapping for %s not impl. zEget_useparameters: overriding parameter %s with value from module %s ) r� r= r� r� r r. rF r� r\ r� ) r� � param_mapr� Zusename�mapping�mvarsrJ rD r� r% r% r& �get_useparameters� s2 � � r� c s� t s| S t| t�r*� �fdd�| D �}|S t| � td�| d f d� � d u rXt| �� � d ur�d| v r�| d }t|�� �D ]@}|| }d|v r||d }d|v r||d }|� v r|� | |d< q|� �fd d�| d D �} | | d < | S ) Nc s g | ]}t |�d � d��qS �r� )�tabr� �� postcrack2)ra �g�r� r� r% r& rc � s �z postcrack2.<locals>.<listcomp>� %sBlock: %s r^ r r� r6 r� c s g | ]}t |�d � d��qS r� r� �ra �br� r% r& rc � s �r� )r r� r= r� r. r� r> ) r� r� r� �retr� r� �varr� �valZnew_bodyr% r� r& r� � s4 � �r� c C sP t | t�rfg }g }| D ]F}t|� t||d d�}d|v rRd|d v rR|�|� q|�|� q|| S t| � t | t�s�d| vr�tdt| � ��d| v r�| d dks�td|| d f d � t | �} t | �} t | �| d < t | d �| d < d | v r�| d r�| d }t | ||d�| d < g }d| v �rP| d }t|�� �D ]}d|v �r2|�|� �q2ni }d} d| v �rj| d } d| v �r>| d �r>g } d| v �r�| d } t�| d �} | �r�| d } nd} | |v �r�d} d| | f |v �r�| d } �q�d| | f } dg i | d d�}| d D ]�}|| v �r�g }d}| d D ]�}|d }|d dk�r&d} |d D ]@}| d } d|v �rL|d |k�rLt�|�}|d | = �q��qL|�r&|d �s�| d |= | | �|�= �q��q&|d �|� n(|| v �rt| | ��s| | |d |< �q|d �s|d �r>| | d< d|gi | | d d�}i || < t�|� |�rL|| d< | S ) zi TODO: function return values determine expression types if in argument list r� �r� r^ ro r� z0postcrack: Expected block dictionary instead of re r� r r� Zsortvarsr r� r� r r r Z__user__routinesZunknown__user__routinesr z%s_%ir_ Z_user_interface)r� r� r� r^ r r\ )r� r� r� r^ r )r� r= r� � postcrackr� r~ r� r� r. � analyzeargs� analyzecommon� analyzevars� sortvarnames� analyzebodyr> r� rQ Z isexternalr )r� r r� ZgretZuretr� Z userisdefinedZuseblockrD r^ r r� Zmnamer� r_ rq Zedef�jr� ZbbZmblockr% r% r& r� s� � � � � r� c C s� g }g }t | �� �D ]2}d| | v r<| | d r<|�|� q|�|� qt|�}d}|r�|d }d}|dd � D ] }|| | d v rpd} q�qp|r�|dd � |g }|d }||kr�tdd�|� d � || }q�qT|�|� |dd � }t|�}d}qT|S )Nr< r r zTsortvarnames: failed to compute dependencies because of cyclic dependencies between rw r� )r= r> r� rO r\ r� )r� Zindep�depr� r� r� r' �wr% r% r& r� d s> �� r� c C s� t | �s | S g }t| d �� �D �]<}g }| d | D �]}t�d|tj�}|�r2g }|�d�r~dd� t|�d���d�D �}t |�d�� � �}|| d v r�d | d | v r�| d | d � d d � |� � n d d � |� g| d | d < n0|�rd d d � |� gi| d |< n i | d |< ||v�rF|� |� n|}t d ||f � |� |� q6|| d |< q d | v�rt|| d < n| d | | d < | S )Nr� z2\A\s*\b(?P<name>.*?)\b\s*(\((?P<dims>.*?)\)|)\s*\Z�dimsc S s g | ] }|� � �qS r% r` ra r% r% r& rc � s �z!analyzecommon.<locals>.<listcomp>r� r^ r� r5 r� r� zNanalyzecommon: failed to extract "<name>[(<dims>)]" from "%s" in common /%s/. � commonvars)Z hascommonr= r> r� r� r� r� r� r� r_ r� r� r� r\ )r� r� rD Zcomvarsrq r� r� r� r% r% r& r� � sN � � � �� � r� c C sD t | � dd� | d �� D �}g }| d D �]}| |d< |d dv r�|d ur\|d |vr\q*n|d }|d |�� v r�t�|d � |d tv r�q*tr�|d tvr�q*t|d d d �|d < n|}t|||d d�}|d dv r�|d s�|�d�s�d|vr�q*|d � dd�dk�rt �|� q*|d dk�r4|d t |d < |�|� q*|S )Nc S s* i | ]"\}}d |vs d|d vr||�qS )r5 r� r% )ra �keyr� r% r% r& � <dictcomp>� s �z analyzebody.<locals>.<dictcomp>r� r� r� r� ri r^ r z T�� as_interfaceZsaved_interfacer� r� rc rp r� r� r � pythonmodulerd ) r� r� r> r! r� r" �crack2fortrangenr� rF r� r r )r� r r� Z maybe_privater� r� Zas_r% r% r& r� � sF � � �� r� c C s� t | � t}i }d| v r�| d d u rFd }tdkr�tdt| d � � nRt| d �� �D ]@}| d | �d�dvr�| d | ||< qV| d | d ||< qV||fS )Nr� r z6buildimplicitrules: no implicit rules for routine %s. r^ r0 )rC r� )r� �defaultimplicitrulesr r. r� r= r> rF )r� � implicitrules� attrrulesrD r% r% r& �buildimplicitrules� s �r c C s0 t | ||�}t|�ttfv r |S td| ��dS )z2 Like `eval` but returns only integers and floats zr=%rN)r� rY r3 r5 r� )rq r� r� r� r% r% r& �myeval� s r z \A\b\w+\b\Zc C s z t t| i i ��}d|dfW S ty. Y n0 t�| �rDdd| fS t| �}|D �]�}t|�|krdqPt�d| d | �rzqPt�d| d tj �}|�| �}|rP�zF|�| �}|r�d|� d �d|� d �f }|�|�}q�t|i i �} |�| �}|�rd|� d �d|� d �f }|�|�}q�t|i i �| } |�| �}|�r`d|� d �d |� d �f }|�|�}�q2t|i i �}|�| �}|�r�d|� d �d |� d �f }|�|�}�qvt|i i �} | d | |k�r�| d | | k�r�| | |fW S W n t�y� Y n0 �qqPd S )a Obtain ``a`` and ``b`` when ``e == "a*x+b"``, where ``x`` is a symbol in xset. >>> getlincoef('2*x + 1', {'x'}) (2, 1, 'x') >>> getlincoef('3*x + x*2 + 2 + 1', {'x'}) (5, 3, 'x') >>> getlincoef('0', {'x'}) (0, 0, None) >>> getlincoef('0*x', {'x'}) (0, 0, 'x') >>> getlincoef('x*x', {'x'}) (None, None, None) This can be tricked by sufficiently complex expressions >>> getlincoef('(x - 0.5)*(x - 1.5)*(x - 1)*x + 2*x + 3', {'x'}) (2.0, 3.0, 'x') r Nr z\w\s*\([^)]*\b�\b�(?P<before>.*?)\b�\b(?P<after>.*)z%s(%s)%sr� r� g �?g �?)NNN) r3 r r� �getlincoef_re_1r� rO r� rO r� r� r� ) rq Zxsetr� Zlen_erb r) r� r* �eer� r� �c2r% r% r& � getlincoef� sf � � � � $r z\b[a-z][\w$]*\bc C s� | |v r�|| � dg �}d|| v rjt|| �sjt�|| d �D ]&}||vrB||v rB|| krB|�|� qB|d d � D ]4}|� |g �p�t|||�D ]}||vr�|�|� q�qvntdt| � � g }||| < |S )Nr< r8 z,_get_depend_dict: no dependence info for %s )rF �isstring� word_patternr� r� �_get_depend_dictr. r� )r^ r� �deps�words�wordr� r% r% r& r F s �r c C s* t | �� �}i }|D ]}t|| |� q|S r( )r= r> r )r� rf � depend_dictr� r% r% r& �_calc_depend_dict\ s r c s� t ��� g }t� �� �D ] }� | s|�|� � |= q� r�t� �� �D ]8\}}� fdd�|D �}|sv|�|� � |= qF|� |< qFq6�fdd�|D �S )Nc s g | ]}|� v r|�qS r% r% )ra r� )r r% r& rc m rd z$get_sorted_names.<locals>.<listcomp>c s g | ]}|� v r|�qS r% r% )ra r^ �r� r% r& rc s rd )r r= r> r� r� )r� rf r^ �lstZnew_lstr% )r r� r& �get_sorted_namesd s r c C s@ | d dv r| dd� } t �| �r&dS t�| �r4dS d| d S ) Nr �'"r r � � zkind(r. )r� r� � real8pattern)rV r% r% r& � _kind_funcv s r c C sL d| }|dkrdS |dkr dS |dkr,dS |dkr8d S |dd krHd S d S ) Nr� � r i r l r l r � � r r% )r� r� r% r% r& �_selected_int_kind_func� s r c C sX | dk r dS | dk rdS t �� �� }|�d�r<| dkrTdS n| dk rHdS | dkrTdS d S ) N� r r r ) �aarch64�alpha�arm64Z loongarchZmips�power�ppcZriscv�s390xZsparc�! � r� r )�platform�machiner� rw )�pr� �radixr+ r% r% r& �_selected_real_kind_func� s r. c C sF t � |�}t � |�}dtfdtfdtffD ]\}}||vr*|||< q*g }t| �D ]*}d| | v rPd| | d v rP|�|� qPt�dtj�}t�dtj�} t�dtj�} |D �]�}d | | v �r$| | d } t | | �r�| � � } d D ]} | j | � } q�|� d | �} | � d | �} d } d| | v �r^d| | d v �r^t | �}| � d| | d d d�} t | �|k } | �s�| �| ��s�| �d�}t |�dk�r�d�|d d� �� � � |d � � d�} t| | ��rt| �}t�| �D ]:}t| |�� |�� � � � � dd��||�� |�� �< �q�d�|�} n t| | ��r4td| � d�� dd� | | d D ��pPd gd }t�| ��rhd} nt�| ��rxd} zt| |||d �||< W nD t �y� } z*| ||< td |� d |�d�� W Y d }~n d }~0 0 t | | ��rt || t ��rt || �||< |� � }||k�r@|| ||< q�t!| | � td |�d �� q�|S )!Nr� Zselected_int_kindZselected_real_kindr5 r� z \bkind\s*\(\s*(?P<value>.*)\s*\)z-\bselected_int_kind\s*\(\s*(?P<value>.*)\s*\)z4\bselected_(int|real)_kind\s*\(\s*(?P<value>.*)\s*\)r8 ))z.false.�False)z.true.�Truez kind("\1")zselected_int_kind(\1)Fr6 r� r r r r� rq zAget_parameters[TODO]: implement evaluation of complex expression r� c S s$ g | ] }|� d �r|�d ��� �qS )r� )rw �lstripr� �ra r� r% r% r& rc � s �z"get_parameters.<locals>.<listcomp>r r r )�dimspeczget_parameters: got "�" on zget_parameters:parameter z does not have value?! )"r� r r r. r r� r� r� r� Z islogicalr� r� �subrO r� r� r� Zisdoubler= r� r� rP r� Z iscomplexr. rO r � param_evalr� r r� r3 r� rL )r� Z global_paramsrJ �g_paramsr^ �funcZ param_namesr� Zkind_reZselected_int_kind_reZselected_kind_rer� �replZ is_replacedZ orig_v_lenZv_r� r� r3 r� �nlr% r% r& r� � s� � �� & � � �� , r� c C s | dv r dS t | |�S )N)z(:)�(*)r r; )� _eval_scalar)�lengthrJ r% r% r& � _eval_length s r> z\d+_c C s� t | �r| �d�d } z&t| i |�} t| t�r2tnt| �} W nX tttfyX | Y S t y� } z&t d|| t |� � �f � W Y d }~n d }~0 0 | S )Nr� r z,"%s" in evaluating %r (available names: %s) ) �_is_kind_numberr� r� r� r� r� r� r� r� r� r\ r= r> )r� rJ r� r% r% r& r< s � r< c > s� t | � t| �\}}t�| d ��| d dkrF| d �vrFi �| d < d| d v rĈd= d| d d v r�| d d d }t��tdd� | d D ��B D ],}d D ]"}||v r�t��|i �|��|< q�q�g }| d }|D ]0}z�| |�|� W q� t�y Y q�0 q�t�� � �D ]}||v�r|�|� �qt �t | ��} i } t � d �j} t�� � �D ]b} | | �} | �rX| | �� | �� � }z | | W n, t�y� t � d | t j�j| |< Y n0 �qX|D � ]�}|d t|� � �v �r�t�| ||d ��|< d�| v�r�d�| v �r$d�| d v �s�|�r�|d �� }t|| � � �D ]z}|dk�rj|| | dk�rj�qF|�| v�r�|| | �| |< n0|dk�rF|| | D ]}t�| |��|< �q��qFn&|| d v �r�tdt|�| d f � d�| v �rTd�| d v �rT�| d d }ztt|i | ��}W n t�yB Y n0 |�| d d< d�| v �r�d�| d v �r��| d d }ztt|i | ��}W n t�y� Y n0 |�| d d< i }d�| v �r��| d }|�� g �| d< d\}}}}}|D �]�}|dd� dk�r4|dd� �� dd � }n�|dd � d k�r`|d d� �� dd � }n�|dd � d k�r�|d d� �� dd � }nj|dd � d!k�r�|d d� �� dd � }n>|dd"� d#k�r�|d"d� �� dd � }nt�| |��|< |�rjd �| v�rg �| d < d$d%� t|��d&�D �D ]6}|� d'd�}|�| d v�r.�| d �|� �q.d}|�r�|� d(d)�}|� d*d+�}d#�| v�r�|g�| d#< n�| d# �|� d}|du�r,d �| v�r�g �| d < t d,d%� t|��d&�D ��D ]*}|�| d v�r��| d �|� �q�d}|du�rd!�| v�rPg �| d!< d-d%� t|��d&�D �D ]*}|�| d! v�rh�| d! �|� �qhd}�q|�r�d�| v�r�g �| d< t d.d%� t|��d&�D ��D �]�}zt || �}W n* t t tf�y td/|�d+�� Y n0 |d0k�r&d0nd1}||k�r<|g}nt|d0��d2�}t!|�d3k�rnd1|v �rnd1g}d1}t!|�dk�r�|d |k�r�d4|d g}t!|�d3k�r�t"t#j$j%|�\} } | | d } | j&t#j'j(d5�}i } | d D ]�} t#�)| �} | �*| ��r�z:| �+| �\}}!||!fd6d7�}"t|�,� �}#|#�-|!�,� � W n6 t.�yh }$ z d}"t| �,� �}#W Y d}$~$n d}$~$0 0 t�fd8d�|#D ��}%|"t|%�f| | < �q�| ||< �| d �|� �q�d!�| v� rdd | v � rd|| d v � rd�| �d g �}&g }'t/t0t1t2��| �}(t3�| �� r�t4�| d �D �]x\})}|�|�� � du � r@�n�|(� r� �5� D �]�\} \}*}+� �fd9d:��t� },�| |,� | |&v � sNd;�| v � sNd �| v � r�� qN|*du� r6| |,v� r6d<}-|*t#�)d=|� d>|)� d?���}.|.j&t#j'j(d5�}.|.�| d;< |g|+ �| d < d!�| v� r�d=|� d>|)� d@|� �g�| d!< nhdA}-d �| v� rTg �| d < dB�| d v� rx�| d �dB� |&�| � |'�d=|� d>|)� d@|� �� �| �dg �}/dC|/v � s�dD|/v � s�|/�|-� r�dDndC� |/� rN|/�| d< � qN� du� r � �5� D ]�\} \}*}+�| �d g �}%�| �dg �D ]8}0|0�6d �� r*d�7|0�� �}0|%�8|0dEd � �dF�� � q*|%� r~tt|%���| d < ||%v� r�|&�| � � q�� q n�t9�| �� r8d�| v � r8d1�| d v � r�t:�| d d1 | �}1|1�| d d1< nFd�| d v � r8t:�| d d | �}1�| d d= |1�| d d1< |'� rJ|'�| d!< |&� rdtt|&���| d < d;�| v �r�d�| v� r�g �| d< dC�| d v� r�dD�| d v� r| d �dC� d �| v� r0g �| d < t| �5� �D ].\} } | �| d; �� r�| d �| � � q�| d � s0�| d = t;�| ��r�t<�| d; | ��| d;< �q�t�� � �D �]�}|| d k� rhd#�| v � r��| d# | d#< | d dk� rhdG| v � r�| dG �v � r�t=�| �| dG ��|< dH| v � rh| dH }2|2� dId�}3|2|3k }4|3� dJd�}2|2|3k }5t>d �|2�} | �r:t?| �@dK�| �@dL��\}6}7}}8tA|6|7�\}9}:};|6�| d< z | dG �r�|6�| dG d< W n t�y� Y n0 |9�r�d|9v �r�zt|9d i | �|9d< W n t�y� Y n0 |9�| d< |:�r�|:�| d< |;�r|;�| dM< |4�r t�| dI��|< |5�rNt�| dJ��|< ntdNt| dH � � � qh| d dOv�r�dP| v �r�t�| d | dP �}<nt�| d �}<t�� � �D ]$}t/tBtC��| ��r�|<�|� �q�dQ| v �r"|<�8t| dQ � � �� t| dQ � � �D ].}| dQ | D ]}||<v�r|<�|� �q�q�| d dk�rXdG| v �rJ|<�| dG � n|<�| d � | d dRv �r�| d }=|=�v �r�d �|= v �r��|= d | d < | d dSk�r�|<�8t�� � �� t�� � �D ]}||<v�r|= �qS )TzH Sets correct dimension information for each variable/parameter r� r� r� r^ r r5 c s s | ]}|d V qdS )r^ Nr% r� r% r% r& rx 2 rd z analyzevars.<locals>.<genexpr>r� rm r z[A-Za-z][\w$]*z.*\b%s\br r0 r� Z undefinedzCanalyzevars: typespec of variable %s is not defined in routine %s. r7 rO r6 r� )NNNNNNr� r� r r r� r� r< r� r; r r: c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd z analyzevars.<locals>.<listcomp>r� r� z\n\nz z\n r� c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd z4analyzevars: could not parse dimension for variable �:r z@:@r �1)�languagec S s | | | S r( r% )r� r� r� r% r% r& �solve_v� s z analyzevars.<locals>.solve_vc 3 s | ]}|j � v r|j V qd S r( )r� r2 r r% r& rx � s �c s: � � | d g g�d D ] }||vr|�|� �||� qd S )Nr )rF �add)r� r Zv1)�coeffs_and_deps� compute_depsr% r& rF � s z!analyzevars.<locals>.compute_depsr8 Fzshape(rw r. z) == T�inr� r� r! r� r� rf Zpure� recursiver r� r9 z'analyzevars: prefix (%s) were not used )rd r� r\ r[ r� r� )r� r� rY )Dr� r r� r� r? rF r� �KeyErrorr= r> r� r� r� r� r� rP r� r� r� r. r� r� r� r� �reverser� r� r� r� rg � param_parser� r� rO r| r �Expr�parse�tostringZLanguager� Z as_symbolr� Z linear_solve�symbolsr� � RuntimeErrorZl_orZ isintent_inZisintent_inoutZisintent_inplaceZisarrayr� r� rw r� r� r r> Zisscalarr< rE r r� r� r� �isintent_callbackZ isintent_aux)>r� r r �genr� rD Zsvarsr r� rJ Z dep_matchesZ name_matchr� r� Zln0r� Zdimension_exprsr� Zdimr� r< r; r: r� �tmpr� Zdim_charr� r� Zd2ZdsizeZsolver_and_depsr� r� rC Z all_symbolsr� Zv_depsZn_depsZn_checksZ n_is_inputr� Zsolverr Zall_depsZ is_requiredr� Zv_attrZaar= �prZpr1ZispureZisrecr0 r� r� r� r� r9 Z neededvarsr^ r% )rE rF r� r& r� " s� " " � � �� " �&�� � ��� � � � � �� � � � � � r� z\A[a-z]+[\w$]*\Zc C s� |du r\zt | ||�}W n> tyV } z&| }td|� d| �d�� W Y d}~n d}~0 0 |S t|�dk s�|ddt|�d � dkr�td|� d ���|dd � �d �}t|�dk�r"|d �d �}t|�dkr�t|d |�}tdt|�d �}n2t|d |�}t|d |�} tt|�t| �d �}ntd��| � d��rB| dd� n| �d �} g } | D ]^} zt | ||�} W n< t�y� } z"td|� d| �d�� W Y d}~n d}~0 0 | � | � �qTt t || ��}|S )a� Creates a dictionary of indices and values for each parameter in a parameter array to be evaluated later. WARNING: It is not possible to initialize multidimensional array parameters e.g. dimension(-3:1, 4, 3:5) at this point. This is because in Fortran initialization through array constructor requires the RESHAPE intrinsic function. Since the right-hand side of the parameter declaration is not executed in f2py, but rather at the compiled c/fortran extension, later, it is not possible to execute a reshape of a parameter array. One issue remains: if the user wants to access the array parameter from python, we should either 1) allow them to access the parameter array using python standard indexing (which is often incompatible with the original fortran indexing) 2) allow the parameter array to be accessed in python as a dictionary with fortran indices as keys We are choosing 2 for now. Nzparam_eval: got "r4 r� r r r1 zparam_eval: dimension z can't be parsedr r� r r@ zEparam_eval: multidimensional array parameters {dimspec} not supportedz(/r� ) r� r� r. rO r� r� rK r� r3 rw r� r~ �zip) r� r7 rJ r3 r, r� Zdimrange�boundZlboundZuboundZv_eval�itemr% r% r& r6 � s: ,& ",r6 c C s� d| v rT| d| � d�� }| | � d�d | �d�� }tt||��}t|| | �S | |v rht|| �S |D ]T}t�d| d tj�}|�| �}|rl|� d�t|| � |� d� } |�| �}q�ql| S dS ) a� Recursively parse array dimensions. Parses the declaration of an array variable or parameter `dimension` keyword, and is called recursively if the dimension for this array is a previously defined parameter (found in `params`). Parameters ---------- d : str Fortran expression describing the dimension of an array. params : dict Previously parsed parameters declared in the Fortran source file. Returns ------- out : str Parsed dimension expression. Examples -------- * If the line being analyzed is `integer, parameter, dimension(2) :: pa = (/ 3, 5 /)` then `d = 2` and we return immediately, with >>> d = '2' >>> param_parse(d, params) 2 * If the line being analyzed is `integer, parameter, dimension(pa) :: pb = (/1, 2, 3/)` then `d = 'pa'`; since `pa` is a previously parsed parameter, and `pa = 3`, we call `param_parse` recursively, to obtain >>> d = 'pa' >>> params = {'pa': 3} >>> param_parse(d, params) 3 * If the line being analyzed is `integer, parameter, dimension(pa(1)) :: pb = (/1, 2, 3/)` then `d = 'pa(1)'`; since `pa` is a previously parsed parameter, and `pa(1) = 3`, we call `param_parse` recursively, to obtain >>> d = 'pa(1)' >>> params = dict(pa={1: 3, 2: 5}) >>> param_parse(d, params) 3 r, Nr r. r r r� r� ) r� �rfindr3 rK r� r� r� r� r� r� )r� rJ ZdnameZddimsrQ r, r) r� r% r% r& rK � s( 9 � ��rK c C sf | }t �| � }|r�t|�\}}t| |d |�}d}| D ](} | �� } | tjtj vrXd} || }q8|d dkrx|d }n|d }|} | |d v s�| |d v r�| d } q�| |v r�d } | t| � |v r�| d } q�| t| � } |r�||d | < nx| |d v�r(||d v �r |d | |d | < n i |d | < d |v �rb||d |d v �rbt |d | d �|d | < | S ) Nr� Ze_r� r rq �_er r� r r r r� ) �analyzeargs_re_1r� r r� r� rV �ascii_lowercase�digitsr� r? ) r� r� r Zorig_aZ a_is_exprr r �at�nar� rD r% r% r& � expr2name: s> r_ c C s t | � t| �\}}d| vr$g | d< g }| d D ]}t|| |�}|�|� q0|| d< d| v r�t| d �� �D ]*\}}|D ] }|| d vrxi | d |< qxql| d D ]B}|d |v r�d| vr�g | d< |d | d vr�| d �|d � q�d| v �r| d | d v�ri | d | d < | S )Nr r� r� r� r^ r r� )r� r r_ r� r= r� )r� r r� r r� rD Zargs1r� r% r% r& r� _ s. r� z\A\(.+?,.+?\)\Zz \A[+-]?\d+(_(?P<name>\w+)|)\Zz*\A[+-]?[\d.]+[-\d+de.]*(_(?P<name>\w+)|)\Zz \A\(.*\)\Zz \A(?P<name>\w+)\s*\(.*?\)\s*\Zc C sT t | t�rddiS t | t�r$ddiS t | t�r6ddiS t | t�rD| S tt| ���d S )Nr0 r2 r1 rr )r� r3 r5 rr r~ �AssertionErrorr� )r� r% r% r& �_ensure_exprdict� s ra c C s� | |v rt || �S | �� } t�| �r.ddiS t�| �}|rjd|�� v rb|�d�rbtdt| � � ddiS t �| �}|r�d|�� v r�|�d�r�tdt| � � ddiS dD ]J}dd � t | |d �� d | d �D �D ] }||v r�t || � S q�q�i }t �| ��r t | d d � ||�}nft�| �}|�r�|�d�}t |�d�||�}|�r^d|v �r^|d= |�s�|d |v �r�t ||d �S | d dv �r�dddid�S |�s�tdt| � � |S )Nr0 rr r^ z:determineexprtype: selected kind types not supported (%s) r2 r1 )�+rt r rv c S s g | ] }|� � �qS r% r` ra r% r% r& rc � rd z%determineexprtype.<locals>.<listcomp>r� r2 r r r5 r r r� r )r0 r7 z>determineexprtype: could not determine expressions (%s) type. )ra r� �determineexprtype_re_1r� �determineexprtype_re_2r� r� r. r� �determineexprtype_re_3r� r� �determineexprtype_re_4r� �determineexprtype_re_5)r� r� �rulesr� �oprq r� Zrnr% r% r& r� � sT � �( �r� r� Fc C s� t | � d}t| t�rf| D ]F}|rN|d dv rN|d tv r<qtrN|d tvrNq|t|||d� }q|S d}d}d}| d }|dkr�dS g } d| v r�| d }d| v r�| d } | d D ]&} t| | | �} t| | �s�| �| � q�| d d ks�| r�d d � | � }d} d | v �r:t| d � � �D ]"} d | |t | | d | f } �q| � dg �d d � }|d k�rld|v �rl|� d� |�r�d| |t d � |�|f } d}d| v �r�t| d |t �}d}d| v �r�t| d |t �}|dk�r�d}d}d| v �rd| d }| d | v�r| �| d � t| d |t |d�}t| | d | |t |d�} d}d| v �r`|�s`d| d }d| v �r�d}t| d �� �D ]$\} }d||t | d � |�f }�q~|| }|dk�r�|d k�r�d}d |||||||| || |||||f }|S ) Nr r� ri r^ r� rY r r� r� z(%s)r� r� z %s%s%s %sr� �callbackz%s%sintent(%s) %sr� r� re r� z result (%s)r� rZ z! in %sr� z%s%sentry %s(%s)r[ r^ z%%s%s%s %s%s%s %s%s%s%s%s%s%send %s %s)r� r� r= r! r" r r_ rQ r� r� r> r rF r � use2fortran�common2fortran� vars2fortranr� )r� r� r� r� r� rf r^ r Z blocktypeZargslr� r� r� rD Z intent_lstr� r� r� r� ZmessZ entry_stmtsr� r% r% r& r � s� � �� � � �r c C sV d}t | �� �D ]@}|dkr6d||d�| | �f }qd|||d�| | �f }q|S )Nr ry z %s%scommon %sr� z%s%scommon /%s/ %s)r= r> r� )r� r� r� rD r% r% r& rl s rl c C s d}t | �� �D ]�}d|||f }| | i krL|r|d dkr|d d� }qd| | v rl| | d rld| }d| | v r�| | d r�d}t | | d �� �D ]L}|| | d | kr�d |||f }d}q�d |||| | d | f }d}q�|r|d dkr|d d� }q|S ) Nr z %s%suse %s,r r� rz z%s only:r| r� �%s%s%sz %s%s%s=>%sr� )r� r� r� r� r� rD r% r% r& rk s( rk c C sR | d }g }|D ]<}zt � d| }W n ty8 Y q0 || �r|�|� q|S )Nr� z isintent_%s)�globalsrI r� )r� r r� r� r0 r% r% r& �true_intent_list' s rp c C s, t | � d}g }|D ]}|| d v r|�|� qd| v rn| d D ],}||v r`||vrl|�|� q@td| � q@d| v r�|�| d � |s�t|�� �D ]}||vr�|�|� q�|D �]t}d|| v �r || d D ]8}||v r�d|| v r�||| d v r�td||f � q�d| v �r�|| d v �r�t|| ��r>d |||f }d |||f }t|| ��rhd |||f }||v �r�d || v�r�q�d } | d D ],} || d k�r�| d dk�r�d} �q��q�| �r�q�||v�r�t|� t d| � q�|| d k�r | d dks�| � d��r q�d || v�rld|| v �rRd|| d v �rR||v r�d |||f }q�t|| � t d| � q�|| d } | dk�r�d|| v �r�d| || d f } i } d|| v �r�|| d } nd || v �r�|| d } d | v �r| d d v �rd | | d f } nd | | d f } nZd!| v �rZd"| | d! f } d#| v �rPd$| | d# f } nd%| } nd#| v �rtd&| | d# f } d'} d|| v �r�d(d)� || d D �}|�r�d*|v �r�d+|v �r�|� d+� |�r�d,| d-� |�f } d-} d.|| v �r d/| | d-� || d. �f } d-} d0|| v �rBt || �}|�r>d1| | d-� |�f } d-} d2|| v �rpd3| | d-� || d2 �f } d-} d|| v �r�d4| | d-� || d �f } d-} d5|| v �r || d5 }|| d d6v �r�z t|�}d7|j|jf }W n t�y� Y n0 d8| ||f } n d9| |f } d:||| f }q�|S );Nr r� r� z7vars2fortran: Confused?!: "%s" is not defined in vars. rh r< zHvars2fortran: Warning: cross-dependence between variables "%s" and "%s" r z%s%sintent(callback) %sz%s%sexternal %sz%s%soptional %sr0 r r� r^ r� r� r z/vars2fortran: No definition for argument "%s". r� r5 r� z-vars2fortran: No typespec for argument "%s". rY r9 z%s(%s)r6 r7 r )r r@ z%s*(%s)z%s*%srO z %s(len=%sr� z %s,kind=%s)z%s)z %s(kind=%s)r� c S s g | ]}|d vr|�qS ))r� r% )ra r� r% r% r& rc � s �z vars2fortran.<locals>.<listcomp>z intent(in)z intent(out)z%s, %sr� r� z%s%sdimension(%s)r� z%s%sintent(%s)r; z %s%scheck(%s)z%s%sdepend(%s)r8 )rr r� z(%s,%s)z %s :: %s=%sz%s :: %srn )r� r� r\ r� r= r> rQ Z isoptional�showr. rF r r� rp r� r1 �imagr� )r� r� r r� r� r� Znoutr� r� r� r� Zvardefr� r� r� r r� r% r% r& rm 5 s� � $ � � � rm c C sz t dd� t| t� t dd� g attd �}t dd� tD ]"}t d|j� d�d� t||�}q<t dd� t |�}t| S )NzReading fortran codes... r zPost-processing... z"Applying post-processing hooks... z r� z Post-processing (stage 2)... ) r. r� r r r� r �post_processing_hooks�__name__�traverser� )�files�postlist�hookr% r% r& � crackfortran� s ry c C s$ t | �d }d}dt }|| | S )Nr� zE! -*- f90 -*- ! Note: the context of this file is case sensitive. z� ! This file was auto-generated with f2py (version:%s). ! See: ! https://web.archive.org/web/20140822061353/http://cens.ioc.ee/projects/f2py2e )r � f2py_version)r� �pyf�header�footerr% r% r& � crack2fortran� s �r~ c C s( t | t�o&t| �dko&t | d ttf�S )Nr r )r� �tuplerO r3 r� )�objr% r% r& �_is_visit_pair� s ��r� c O s@ t | �rX| d dkr| S || ||g|�R i |��}|durJt |�sFJ �|S | }| \}} n d| f}d}t| t�r�g }t| �D ]F\} } t| | f|g|�R ||g |d�|��\} } | durz|�| � qznft| t��r&t� }| �� D ]D\} } t| | f|g|�R ||g |d�|��\}}|dur�|||< q�n| }|du �r8|S ||fS )a� Traverse f2py data structure with the following visit function: def visit(item, parents, result, *args, **kwargs): """ parents is a list of key-"f2py data structure" pairs from which items are taken from. result is a f2py data structure that is filled with the return value of the visit function. item is 2-tuple (index, value) if parents[-1][1] is a list item is 2-tuple (key, value) if parents[-1][1] is a dict The return value of visit must be None, or of the same kind as item, that is, if parents[-1] is a list, the return value must be 2-tuple (new_index, new_value), or if parents[-1] is a dict, the return value must be 2-tuple (new_key, new_value). If new_index or new_value is None, the return value of visit is ignored, that is, it will not be added to the result. If the return value is None, the content of obj will be traversed, otherwise not. """ r r� N)�parentsr� )r� r� r= r� ru r� r~ r� )r� �visitr� r� r �kwargsZ new_result�parentZ result_keyrQ r� Z new_indexZnew_itemr� Znew_key� new_valuer% r% r& ru � sR ��� ��� ru c O s> |d \}}| \}}dd� } |dv rF|d d dks8J �|d d } n.|d krp|d d dksbJ �|d d } nd } d } | d ur�|} | � � D ]\} } t| �r�| | | �} q�nP|d kr�|d d d } |} | � � D ]*\} } t| �r�t�d | d d| | �} q�| d u�r:| |k�r2td|� d|� d|� d| � d� d� || fS d S )a� Previously, Fortran character was incorrectly treated as character*1. This hook fixes the usage of the corresponding variables in `check`, `dimension`, `=`, and `callstatement` expressions. The usage of `char*` in `callprotoargument` expression can be left unchanged because C `character` is C typedef of `char`, although, new implementations should use `character*` in the corresponding expressions. See https://github.com/numpy/numpy/pull/19388 for more information. r c S s0 t �d| d | |�}t �d| d | |�}|S )Nz[*]\s*\br z\b\s*[\[]\s*0\s*[\]])r� r5 )�varnamer� r% r% r& � fix_usage5 s �z8character_backward_compatibility_hook.<locals>.fix_usage)r� r; ���r r� r r8 r� NZ callstatementz (?<![&])\br r� zcharacter_bc_hook[rX z ]: replaced `z` -> `z` )r� Z ischaracterr� r5 r. )rW r� r� r r� Z parent_keyZ parent_valuer� r� r� Z vars_dictr� r� Zvdr% r% r& �%character_backward_compatibility_hook# sD � ���r� �__main__r@ z-quietz-verboser z-fixz?Use option -f90 before -fix if Fortran 90 code is in fix form. z-skipemptyendsz--ignore-containsz-f77z-f90r� z-hz-showz-mrt zUnknown option %s z OSError: a Warning: You have specified module name for non Fortran 77 code that should not need one (expect if you are scanning F90 code for non module blocks but then you should use flag -skipemptyends and also be sure that the files do not contain programs without program statement). z Writing fortran code to file %s r� )r )r )r� )r )r )r )r )r� )N)r N)Nr )r )NN)r r )N)r� F)r )r )r F)��__doc__r) rV r� r� rr r� r* rx �pathlibr ro � ImportErrorr r Zauxfuncsr �versionrz r r r r r r r r r r r r� r r r r r r r r r r# r r" r$ r r! r r r' r. � _MAXCACHEr r� r[ rh r� r_ rg ri rj r� r� rO r� r� r� r� r� r� r� r� r� rq r� ZbeforethisafterZ fortrantypesr r� r r Z groupbegins77r� Z groupbegins90r� Z groupendsr Zendifsr Zmoduleproceduresr r r r r r r r r r r r r r r# r r r r �Sr r r� r� r! Zcrackline_bind_1r� r r" r� r4 rE r� rK rT rU rV r� r rG rH rN rR rW r$ r� r� r� r6 r7 r� r� r� r� r� r� r? r@ rA r� r� r� r� r� r� r� r� r� r r r r r r r r r r r. r� r> r? r< r� rZ r6 rK r_ r� rc rd re rf rg ra r� r rl rk rp rm rs ry r~ r� ru r� r� rt rv �funcsr0 �f2Zf3Z showblocklist�argvr� r\ r� ru r� �OSErrorZdetailrw r{ r+ r% r% r% r& �<module> sL & Z���� �� �������������� �� �� �� �� �� �� �� �� �� �� �� �� �� �� �� �� �� �� ���� ������� W ����% *< `"+ - H e >P %� - K B< � *�*
Save
Cancel