o *ήcAS @sBdZddlmZmZddlmZmZmZm Z ddl m Z idddddd d d d d dddfddddfddddfddddfddddfddddfddddfgdd dddfgd!d"d#d$d%d&d'd'd(d(d)d)d*d*d+d+d,d,d-d-d.d.d/d0d1d2d3dd4Z gd5Z Gd6d7d7e Zd=d9d:Zd;d<Zd8S)>a Rust code printer The `RustCodePrinter` converts SymPy expressions into Rust expressions. A complete code generator, which uses `rust_code` extensively, can be found in `sympy.utilities.codegen`. The `codegen` module can be used to generate complete source code files. )AnyDict)SRationalFloatLambda) CodePrinterfloorceilingceilAbsabssignsignumPowcC |tj kSNrOnebaseexpr:/tmp/pip-target-vg8gfxp4/lib/python/sympy/printing/rust.pyG rrecipcCs |tjkSrrHalfrrrrrHs sqrtcCrrrrrrrrIrz sqrt().recipcCs|tddkS)N)rrrrrrJcbrtcCs|tjdkS)NrrrrrrrKr#exp2r"cC|jSr is_integerrrrrrLspowir!cCs|j Srr'rrrrrMspowfrcCdS)NTr)rrrrrNsloglnMaxmaxMinminsincostanasinacosatanatan2sinhcoshtanhasinhacoshatanh)r:r;r<r=r>r )4abstractalignofasbecomeboxbreakconstcontinuecratedoelseenumexternfalsefinalfnforifimplinletloopmacromatchmodmovemutoffsetofoverrideprivprocpubpurerefreturnSelfselfsizeofstaticstructsupertraittruetypetypeofunsafeunsizedusevirtualwherewhileyieldc sleZdZdZdZdZdddiddedd dd Zifd d Zd dZ ddZ ddZ ddZ ddZ ddZddZddZddZdd ZdOfd!d" ZdOfd#d$ Zd%d&Zd'd(Zd)d*Zd+d,Zd-d.ZdOd/d0ZdOd1d2ZdOd3d4ZdOd5d6ZdOd7d8ZdOd9d:ZdOd;d<Z dOd=d>Z!d?d@Z"dAdBZ#dCdDZ$dEdFZ%dGdHZ&fdIdJZ'dKdLZ(dMdNZ)Z*S)PRustCodePrinterz>A printer to convert SymPy expressions to strings of Rust code _rust_codeRustNautoTF_) order full_prec precisionuser_functionshumancontract dereferenceerror_on_reservedreserved_word_suffixinlinecCsNt||tt|_|di}|j|t|dg|_tt|_dS)Nr|r) r__init__dictknown_functionsgetupdateset _dereferencereserved_words)rcsettings userfuncsrrrrs    zRustCodePrinter.__init__cCs|dS)Nr)rcprrr_rate_index_positionz$RustCodePrinter._rate_index_positioncCd|S)Nz%s;r)rc codestringrrr_get_statementrzRustCodePrinter._get_statementcCr)Nz// %sr)rctextrrr _get_commentrzRustCodePrinter._get_commentcCs d||fS)Nzconst %s: f64 = %s;r)rcnamevaluerrr_declare_number_const z%RustCodePrinter._declare_number_constcC ||Sr) indent_code)rclinesrrr _format_codes zRustCodePrinter._format_codecs |j\}fddt|DS)Nc3s&|]}tD]}||fVqqdSr)range).0ijcolsrr s$z;RustCodePrinter._traverse_matrix_indices..)shaper)rcmatrowsrrr_traverse_matrix_indicess z(RustCodePrinter._traverse_matrix_indicesc CsXg}g}d}|D]}||||||j||jdd|dq||fS)Nz#for %(var)s in %(start)s..%(end)s {r!)varstartend})append_printlowerupper)rcindices open_lines close_lines loopstartrrrr_get_loop_opening_endings  z(RustCodePrinter._get_loop_opening_endingcCs>t|jdkrd||dS|jr|j|ddS||S)Nr!()T)_type)lenargsr is_numberrcexprrrr_print_caller_vars  z!RustCodePrinter._print_caller_varc s|jjjvrj|jj}d}d}t|tr|}n|D] \}}}||jr(nq|dur|dkrSd|jd|t|jdkrL|jdddndd}|S|dkrfd |jd|d }|S|d kryd |jd|d }|Sd ||jdd }|SdSt |drt|j t r |j |jS|jjj vrɈj |jj\}}|rtfdd|DrLj ||SdSdS|S)a< basic function for printing `Function` Function Style : 1. args[0].func(args[1:]), method with arguments 2. args[0].func(), method without arguments 3. args[1].func(), method without arguments (e.g. (e, x) => x.exp()) 4. func(args), function with arguments Nr!z%(var)s.%(method)s(%(args)s)r, )rmethodrrz%(var)s.%(method)s())rrr"z%(func)s(%(args)s))funcr_imp_c3|]}|VqdSr) _can_print)rfrcrrrQz2RustCodePrinter._print_Function..)r__name__r isinstancestrrrr stringifyhasattrrrr_rewriteable_functionsrallrewrite_print_not_supported) rcr cond_funcrstylecondrettarget_f required_fsrrr_print_FunctionsX   &    zRustCodePrinter._print_FunctioncCs:|jjr|jjst|t|j|j}||S||Sr)rr(rrjrrrrrrr _print_PowVs  zRustCodePrinter._print_Powct|}|r |dS|S)N_f64)rg _print_Floatrcrrr __class__rrr\ zRustCodePrinter._print_Floatcr)N_i32)rg_print_IntegerrrrrrcrzRustCodePrinter._print_IntegercCs"t|jt|j}}d||fS)Nz %d_f64/%d.0)intrq)rcrrrrrr_print_Rationaljs zRustCodePrinter._print_RationalcCs,||j}||j}|j}d|||S)Nz{} {} {})rlhsrhsrel_opformat)rcrlhs_coderhs_codeoprrr_print_Relationalns  z!RustCodePrinter._print_RelationalcCs`|j}tj}tj}tt|jD]}||j||7}|||9}qd||j j ||fS)N%s[%s]) rrZerorreversedrrankrrrlabel)rcrdimselemoffsetrrrr_print_IndexedtszRustCodePrinter._print_IndexedcCs|jjSr)rrrrrr _print_Idx~rzRustCodePrinter._print_IdxcCr&r)rrrrr _print_DummyszRustCodePrinter._print_DummycCr+)NErrcrrrrr _print_Exp1zRustCodePrinter._print_Exp1cCr+)NPIrrrrr _print_PirzRustCodePrinter._print_PicCr+)NINFINITYrrrrr_print_InfinityrzRustCodePrinter._print_InfinitycCr+)N NEG_INFINITYrrrrr_print_NegativeInfinityrz'RustCodePrinter._print_NegativeInfinitycCr+)Nrirrrrr_print_BooleanTruerz"RustCodePrinter._print_BooleanTruecCr+)NrLrrrrr_print_BooleanFalserz#RustCodePrinter._print_BooleanFalsecCs t|Sr)rrrrrr _print_boolrzRustCodePrinter._print_boolcCr+)NNANrrrrr _print_NaNrzRustCodePrinter._print_NaNcCs|jdjdkr tdg}t|jD]G\}\}}|dkr(|d||n#|t|jdkr>|dkr>|dd7<n |dd||7<||}|||d q|jd red |Sd |S) NTzAll 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.rz if (%s) {r!z else {z else if (%s) {rr  ) rr ValueError enumeraterrr _settingsjoin)rcrrreccode0rrr_print_Piecewises      z RustCodePrinter._print_PiecewisecCs ddlm}||j|ddS)Nr) PiecewiseF)deep)sympy.functionsrrr)rcrrrrr _print_ITEs zRustCodePrinter._print_ITEcs.|jdkrddfdd|DStd)Nr!z[%s]rc3rr)r)rarrrrrz4RustCodePrinter._print_MatrixBase..zLFull Matrix Support in Rust need Crates (https://crates.io/keywords/matrix).)rrr )rcArrr_print_MatrixBases z!RustCodePrinter._print_MatrixBasecCrr)r)rcrrrr_print_SparseRepMatrixs z&RustCodePrinter._print_SparseRepMatrixcCs"d|j|j|j|jjdfS)Nrr!)parentrrrrrrr_print_MatrixElementsz$RustCodePrinter._print_MatrixElementcs"t|}||jvrd|S|S)Nz(*%s))rg _print_Symbolr)rcrrrrrrs  zRustCodePrinter._print_SymbolcCshddlm}|j}|j}|jdr!||s||r!|||S||}||}|d||fS)Nr) IndexedBaser~z%s = %s) sympy.tensor.indexedrrrr has_doprint_loopsrr)rcrrrrrrrrr_print_Assignments    z!RustCodePrinter._print_Assignmentc st|tr||d}d|Sd}dddd|D}fdd|D}fd d|D}g}d }t|D]%\}} | d vrG|| q9|||8}|d ||| f|||7}q9|S) z0Accepts a string of code or a list of code linesTrz ){rz{ z( )rrcSsg|]}|dqS)z )lstriprlinerrr sz/RustCodePrinter.indent_code..c g|] }ttt|jqSr)ranymapendswithr%) inc_tokenrrr's cr(r)rr)r* startswithr%) dec_tokenrrr'sr)rr z%s%s)rrr splitlinesrr r) rccode code_linestabincreasedecreaseprettylevelnr&r)r.r,rrs*     zRustCodePrinter.indent_code)F)+r __module__ __qualname____doc__ printmethodlanguager_default_settingsrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr"r __classcell__rrrrrssb   8           rsNcKst|||S)aConverts an expr to a string of Rust 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 is helpful in case of line-wrapping, or for expressions that generate multi-line statements. precision : integer, optional The precision for numbers such as pi [default=15]. user_functions : dict, optional A dictionary where the keys are string representations of either ``FunctionClass`` or ``UndefinedFunction`` instances and the values are their desired C string representations. Alternatively, the dictionary value can be a list of tuples i.e. [(argument_test, cfunction_string)]. See below for examples. dereference : iterable, optional An iterable of symbols that should be dereferenced in the printed code expression. These would be values passed by address to the function. For example, if ``dereference=[a]``, the resulting code would print ``(*a)`` instead of ``a``. 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]. Examples ======== >>> from sympy import rust_code, symbols, Rational, sin, ceiling, Abs, Function >>> x, tau = symbols("x, tau") >>> rust_code((2*tau)**Rational(7, 2)) '8*1.4142135623731*tau.powf(7_f64/2.0)' >>> rust_code(sin(x), assign_to="s") 's = x.sin();' Simple custom printing can be defined for certain types by passing a dictionary of {"type" : "function"} to the ``user_functions`` kwarg. Alternatively, the dictionary value can be a list of tuples i.e. [(argument_test, cfunction_string)]. >>> custom_functions = { ... "ceiling": "CEIL", ... "Abs": [(lambda x: not x.is_integer, "fabs", 4), ... (lambda x: x.is_integer, "ABS", 4)], ... "func": "f" ... } >>> func = Function('func') >>> rust_code(func(Abs(x) + ceiling(x)), user_functions=custom_functions) '(fabs(x) + x.CEIL()).f()' ``Piecewise`` expressions are converted into conditionals. If an ``assign_to`` variable is provided an if statement is created, otherwise the ternary operator is used. Note that if the ``Piecewise`` lacks a default term, represented by ``(expr, True)`` then an error will be thrown. This is to prevent generating an expression that may not evaluate to anything. >>> from sympy import Piecewise >>> expr = Piecewise((x + 1, x > 0), (x, True)) >>> print(rust_code(expr, tau)) tau = if (x > 0) { x + 1 } else { x }; Support for loops is provided through ``Indexed`` types. With ``contract=True`` these expressions will be turned into loops, whereas ``contract=False`` will just print the assignment expression that should be looped over: >>> from sympy import Eq, IndexedBase, Idx >>> len_y = 5 >>> y = IndexedBase('y', shape=(len_y,)) >>> t = IndexedBase('t', shape=(len_y,)) >>> Dy = IndexedBase('Dy', shape=(len_y-1,)) >>> i = Idx('i', len_y-1) >>> e=Eq(Dy[i], (y[i+1]-y[i])/(t[i+1]-t[i])) >>> rust_code(e.rhs, assign_to=e.lhs, contract=False) 'Dy[i] = (y[i + 1] - y[i])/(t[i + 1] - t[i]);' Matrices are also supported, but a ``MatrixSymbol`` of the same dimensions must be provided to ``assign_to``. Note that any expression that can be generated normally can also exist inside a Matrix: >>> from sympy import Matrix, MatrixSymbol >>> mat = Matrix([x**2, Piecewise((x + 1, x > 0), (x, True)), sin(x)]) >>> A = MatrixSymbol('A', 3, 1) >>> print(rust_code(mat, A)) A = [x.powi(2), if (x > 0) { x + 1 } else { x }, x.sin()]; )rsdoprint)r assign_torrrr rust_codesmrAcKstt|fi|dS)z3Prints Rust representation of the given expression.N)printrA)rrrrrprint_rust_codemsrCr)r:typingrrtDict sympy.corerrrrsympy.printing.codeprinterrrrrsrArCrrrrsv!         !"#$%&'+,l6 & p