a
    <bC)                      @   s   d Z ddlmZ ddlmZmZ ddlmZ ddlm	Z	m
Z
 ddlZdZdd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'ZeD ]Zeee< qejd(ejd)ejd*ejd+ejd,iZd-d.d/ZejgZG d0d1 d1eZd6d2d3Zd4d5 ZdS )7a	  
Maple code printer

The MapleCodePrinter converts single SymPy expressions into single
Maple expressions, using the functions defined in the Maple objects where possible.


FIXME: This module is still under actively developed. Some functions may be not completed.
    )S)IntegerIntegerConstant)CodePrinter)
precedence
PRECEDENCEN) sincostansecZcscZcotsinhcoshtanhZsechZcschZcothexpfloor	factorialZ	bernoulliZeulerZ	fibonaccigcdlcm	conjugateZCiChiZEiZLiZSiZShierferfcZharmonicZLambertWsqrtabslnZarcsinZarccosZarctanZarcsecZarccscZarccotZarcsinhZarccoshZarctanhZarcsechZarccschZarccothceilmaxminZdoublefactorialZ
pochhammerZBesselIZBesselJZBesselKZBesselYZHankelH1ZHankelH2ZAiryAiZAiryBiZAppellF1ZFresnelCZFresnelSZLerchPhi)ZAbslogasinacosatanZasecZacscZacotasinhacoshatanhZasechZacschZacothZceilingZMaxZMinZ
factorial2ZRisingFactorialZbesseliZbesseljZbesselkZbesselyZhankelh1Zhankelh2ZairyaiZairybiZappellf1ZfresnelcZfresnelsZlerchphiPizexp(1)Catalangammaz(1/2 + (1/2)*sqrt(5))=z<>)z==z!=c                       s  e Zd ZdZdZdZddddddZdC f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d 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*d+ Zd,d- Zd.d/ ZdDd1d2Zd3d4 Zd5d6 Zd7d8 Zd9d: Zd;d< Z d=d> Z!d?d@ Z"dAdB Z#  Z$S )EMapleCodePrinterzF
    Printer which converts a SymPy expression into a maple code.
    Z_mapleZmapleNautoT)orderZ	full_precZhumaninlineZallow_unknown_functionsc                    s@   |d u rt  }t | t t| _|di }| j| d S )NZuser_functions)dictsuper__init__known_functionsgetupdate)selfsettingsZ	userfuncs	__class__ d/Users/vegardjervell/Documents/master/model/venv/lib/python3.9/site-packages/sympy/printing/maple.pyr/   c   s    
zMapleCodePrinter.__init__c                 C   s   d| S )Nz%s;r7   )r3   Z
codestringr7   r7   r8   _get_statementk   s    zMapleCodePrinter._get_statementc                 C   s
   d |S )Nz# {})format)r3   textr7   r7   r8   _get_commentn   s    zMapleCodePrinter._get_commentc                 C   s   d ||| jd S )Nz	{} := {};Z	precision)r:   ZevalfZ	_settings)r3   namevaluer7   r7   r8   _declare_number_constq   s    z&MapleCodePrinter._declare_number_constc                 C   s   |S Nr7   )r3   linesr7   r7   r8   _format_codeu   s    zMapleCodePrinter._format_codec                 C   s   |  t|S r@   _printlistr3   exprr7   r7   r8   _print_tuplex   s    zMapleCodePrinter._print_tuplec                 C   s   |  t|S r@   rC   rF   r7   r7   r8   _print_Tuple{   s    zMapleCodePrinter._print_Tuplec                 C   s&   |  |j}|  |j}dj||dS )Nz{lhs} := {rhs})lhsrhs)rD   rJ   rK   r:   )r3   rG   rJ   rK   r7   r7   r8   _print_Assignment~   s    z"MapleCodePrinter._print_Assignmentc                 K   s   t |}|jdkr$d| |j| S |jdtjfv rDd| |j S |jdtj fv rfd| |j S dj| |j|| |j|dS d S )	Nz1/%sg      ?zsqrt(%s)g      z
1/sqrt(%s)z{base}^{exp})baser   )r   r   parenthesizerN   r   ZHalfrD   r:   )r3   rG   kwargsPRECr7   r7   r8   
_print_Pow   s    
zMapleCodePrinter._print_Powc                    sT   |j d jdur*|j d jtjkr*td fdd|j D }d|}dj|dS )	NrM   TzAll Piecewise expressions must contain an (expr, True) statement to be used as a default condition. Without one, the generated expression may not evaluate to anything under some condition.c                    sN   g | ]F\}}|d ur8|t jur8dj | |dndj |dqS )Tz{c}, {e})cez{e})rT   )r   BooleanTruer:   rD   ).0rT   rS   r3   r7   r8   
<listcomp>   s   z5MapleCodePrinter._print_Piecewise.<locals>.<listcomp>z, zpiecewise({_inbrace}))_inbrace)argsZcondr   rU   
ValueErrorjoinr:   )r3   rG   Z
_coup_listrY   r7   rW   r8   _print_Piecewise   s    "

z!MapleCodePrinter._print_Piecewisec                 C   s,   t |jt |j }}djt|t|dS )Nz{p}/{q})pq)intr^   r_   r:   str)r3   rG   r^   r_   r7   r7   r8   _print_Rational   s    z MapleCodePrinter._print_Rationalc                 C   sJ   t |}| |j|}| |j|}|j}|tv r:t| }dj|||dS )Nz{lhs} {rel_op} {rhs})rJ   rel_oprK   )r   rO   rJ   rK   rc   spec_relational_opsr:   )r3   rG   rQ   Zlhs_codeZrhs_codeopr7   r7   r8   _print_Relational   s    z"MapleCodePrinter._print_Relationalc                 C   s   t | S r@   )number_symbolsrF   r7   r7   r8   _print_NumberSymbol   s    z$MapleCodePrinter._print_NumberSymbolc                 C   s   dS )Nz	-infinityr7   rF   r7   r7   r8   _print_NegativeInfinity   s    z(MapleCodePrinter._print_NegativeInfinityc                 C   s   dS )NZinfinityr7   rF   r7   r7   r8   _print_Infinity   s    z MapleCodePrinter._print_Infinityc                 C   s   |  |jS r@   )rD   labelrF   r7   r7   r8   
_print_Idx   s    zMapleCodePrinter._print_Idxc                 C   s   dS )Ntruer7   rF   r7   r7   r8   _print_BooleanTrue   s    z#MapleCodePrinter._print_BooleanTruec                 C   s   dS )Nfalser7   rF   r7   r7   r8   _print_BooleanFalse   s    z$MapleCodePrinter._print_BooleanFalsec                 C   s   |rdS dS )Nrm   ro   r7   rF   r7   r7   r8   _print_bool   s    zMapleCodePrinter._print_boolc                 C   s   dS )N	undefinedr7   rF   r7   r7   r8   
_print_NaN   s    zMapleCodePrinter._print_NaNFc                 C   sF   t j|jv r"dj|rdndd}n dj| | |r:dndd}|S )NzMatrix([], storage = {storage})sparseZrectangular)storagez#Matrix({list}, storage = {storage}))rE   ru   )r   ZZeroshaper:   rD   tolist)r3   rG   rt   Z_strMr7   r7   r8   _get_matrix   s    

zMapleCodePrinter._get_matrixc                 C   s:   dj | j|jtd dd| |jd | |jd dS )Nz{parent}[{i_maple}, {j_maple}]ZAtomT)strict   )parentZi_mapleZj_maple)r:   rO   r{   r   rD   ijrF   r7   r7   r8   _print_MatrixElement   s
    z%MapleCodePrinter._print_MatrixElementc                 C   s   | j |ddS )NFrt   rx   rF   r7   r7   r8   _print_MatrixBase   s    z"MapleCodePrinter._print_MatrixBasec                 C   s   | j |ddS )NTr   r   rF   r7   r7   r8   _print_SparseRepMatrix   s    z'MapleCodePrinter._print_SparseRepMatrixc                 C   s8   t |jttfr | t|S dj| |jdS d S )Nz$Matrix({var_size}, shape = identity))Zvar_size)
isinstancerowsr   r   rD   sympyZSparseMatrixr:   rF   r7   r7   r8   _print_Identity   s    z MapleCodePrinter._print_Identityc                    s   t | t|j}d }t|d tjtjtjtjfsJ|d |dd   }}|d u sZ|dkrtd	 fdd|D S dj
|d	 fdd|D dS d S )	Nr   rz   .c                 3   s   | ]} | V  qd S r@   rO   rV   _mrQ   r3   r7   r8   	<genexpr>       z1MapleCodePrinter._print_MatMul.<locals>.<genexpr>z{c}*{m}c                 3   s   | ]} | V  qd S r@   r   r   r   r7   r8   r      r   )rS   m)r   rE   rZ   r   r   Z
MatrixBaseZ
MatrixExprZMatrixSliceZMatrixSymbolr\   r:   )r3   rG   
_fact_listZ_constr7   r   r8   _print_MatMul   s    
zMapleCodePrinter._print_MatMulc                 C   s   dj | |j| |jdS )NzMatrixPower({A}, {n}))An)r:   rD   rN   r   rF   r7   r7   r8   _print_MatPow   s    zMapleCodePrinter._print_MatPowc                    s,   t | t|j}d fdd|D S )N*c                 3   s   | ]} | V  qd S r@   r   r   r   r7   r8   r      r   z:MapleCodePrinter._print_HadamardProduct.<locals>.<genexpr>)r   rE   rZ   r\   )r3   rG   r   r7   r   r8   _print_HadamardProduct   s    
z'MapleCodePrinter._print_HadamardProductc                 C   sX   |j \}\}}|dkr2dj| || |d}ndj| |d}dj| ||dS )Nrz   z{var}${order})varr+   z{var})r   zdiff({func_expr}, {sec_arg}))Z	func_exprZsec_arg)rZ   r:   rD   )r3   rG   Z_fZ_var_orderZ_second_argr7   r7   r8   _print_Derivative   s    z"MapleCodePrinter._print_Derivative)N)F)%__name__
__module____qualname____doc__ZprintmethodlanguageZ_default_settingsr/   r9   r<   r?   rB   rH   rI   rL   rR   r]   rb   rf   rh   ri   rj   rl   rn   rp   rq   rs   rx   r~   r   r   r   r   r   r   r   __classcell__r7   r7   r5   r8   r)   T   sL   	

r)   c                 K   s   t || |S )aU  Converts ``expr`` to a string of Maple code.

    Parameters
    ==========

    expr : Expr
        A SymPy expression to be converted.
    assign_to : optional
        When given, the argument is used as the name of the variable to which
        the expression is assigned.  Can be a string, ``Symbol``,
        ``MatrixSymbol``, or ``Indexed`` type.  This can be helpful for
        expressions that generate multi-line statements.
    precision : integer, optional
        The precision for numbers such as pi  [default=16].
    user_functions : dict, optional
        A dictionary where keys are ``FunctionClass`` instances and values are
        their string representations.  Alternatively, the dictionary value can
        be a list of tuples i.e. [(argument_test, cfunction_string)].  See
        below for examples.
    human : bool, optional
        If True, the result is a single string that may contain some constant
        declarations for the number symbols.  If False, the same information is
        returned in a tuple of (symbols_to_declare, not_supported_functions,
        code_text).  [default=True].
    contract: bool, optional
        If True, ``Indexed`` instances are assumed to obey tensor contraction
        rules and the corresponding nested loops over indices are generated.
        Setting contract=False will not generate loops, instead the user is
        responsible to provide values for the indices in the code.
        [default=True].
    inline: bool, optional
        If True, we try to create single-statement code instead of multiple
        statements.  [default=True].

    )r)   Zdoprint)rG   Z	assign_tor4   r7   r7   r8   
maple_code  s    $r   c                 K   s   t t| fi | dS )a&  Prints the Maple representation of the given expression.

    See :func:`maple_code` for the meaning of the optional arguments.

    Examples
    ========

    >>> from sympy import print_maple_code, symbols
    >>> x, y = symbols('x y')
    >>> print_maple_code(x, assign_to=y)
    y := x
    N)printr   )rG   r4   r7   r7   r8   print_maple_code*  s    r   )N)r   Z
sympy.corer   Zsympy.core.numbersr   r   Zsympy.printing.codeprinterr   Zsympy.printing.precedencer   r   r   Z_known_func_same_namer0   Z_funcr%   ZExp1r&   Z
EulerGammaZGoldenRatiorg   rd   ZComplexInfinityZnot_supported_symbolr)   r   r   r7   r7   r7   r8   <module>   sn   

$
 0
'