*h7 SrSSKJr SSKJr SSKJr SSKJr SSK J r J r SSK J r 0SS _S S _S S _S S _S S _SS_SS_SS_SS_SS_SS_SS_SS_SS_SS_SS_SS_SSSSSSS S!S"S#S$S%. Er/S&Qr"S'S(\5rS,S*jrS+rg))-z R code printer The RCodePrinter converts single SymPy expressions into single R expressions, using the functions defined in math.h where possible. ) annotations)Any) equal_valued) CodePrinter) precedence PRECEDENCE)RangeAbsabssincostanasinacosatanatan2explogerfsinhcoshtanhasinhacoshatanhfloorceilingsignmaxmin factorialgammadigammatrigammabetasqrt) rrrMaxMinr!r"r#r$r%r&)ifelserepeatwhilefunctionforinnextbreakTRUEFALSENULLInfNaNNA NA_integer_NA_real_ NA_complex_ NA_character_volatilec $^\rSrSr%SrSrSr\"\R40S0S\ "5S.D6r S\ S 'S S S S .r 0r S\ S'04SjrSrSrSrSrSrSrSrSrSrSrSrSrSrSrSrS rS!rS"r U4S#jr!S$r"S%r#S&r$S'r%S(r&U=r'$)) RCodePrinterNz;A printer to convert SymPy expressions to strings of R code_rcodeRT) precisionuser_functionscontract dereferencezdict[str, Any]_default_settings&|!)andornotzdict[str, str] _relationalsc[R"X5 [[5UlUR S05nURR U5 [ UR S/55Ul[ [5Ulg)NrDrF) r__init__dictknown_functionsgetupdateset _dereferencereserved_words)selfsettings userfuncss Q/home/mande/repo/quber/.venv/lib/python3.13/site-packages/sympy/printing/rcode.pyrPRCodePrinter.__init__asgT,#O4LL!126  ##I. ]B ?@!.1c US-$)N)rXps r[_rate_index_position!RCodePrinter._rate_index_positionis s r]c SU-$)Nz%s;r`)rX codestrings r[_get_statementRCodePrinter._get_statementls z!!r]c$SRU5$)Nz// {}format)rXtexts r[ _get_commentRCodePrinter._get_commentos~~d##r]c$SRX5$)Nz{} = {};ri)rXnamevalues r[_declare_number_const"RCodePrinter._declare_number_constrs  --r]c$URU5$N) indent_code)rXliness r[ _format_codeRCodePrinter._format_codeus&&r]cL^URunmU4Sj[U55$)Nc3P># UHn[T5Ho!U4v M M g7frt)range).0ijcolss r[ 8RCodePrinter._traverse_matrix_indices..zsA 1U4[A[ s#&)shaper{)rXmatrowsrs @r[_traverse_matrix_indices%RCodePrinter._traverse_matrix_indicesxsYY dAd AAr]c /n/nSnUH}nURUURUR5URURS-5URURS-5S.-5 URS5 M X#4$)zPReturns a tuple (open_lines, close_lines) containing lists of codelines z#for (%(var)s in %(start)s:%(end)s){)varstartend})append_printlabellowerupper)rXindices open_lines close_lines loopstartr}s r[_get_loop_opening_ending%RCodePrinter._get_loop_opening_ending|s  9 A   i{{177+QWWQY/{{177Q;/+11 2   s # &&r]cSUR;aURU5$[U5n[URS5(aSUR UR U5-$[URS5(aSURUR 5-$UR UR U5<SUR URU5<3$)NPowz1.0/%sg?zsqrt(%s)^)rR_print_Functionrrr parenthesizebaser)rXexprPRECs r[ _print_PowRCodePrinter._print_Pows D(( (''- -$ " % %t00DAB B $((C ( ( DII 66 6"// 4@!%!2!2488T!BD Dr]c`[UR5[UR5p2SX#4-$)Nz %d.0/%d.0)intraq)rXrrars r[_print_RationalRCodePrinter._print_Rationals&466{CK1aV##r]cURVs/sHo RU5PM nnURURR5<SSR U5<S3$s snf)N[z, ])rrrrjoin)rXrr}indss r[_print_IndexedRCodePrinter._print_IndexedsI)-7AQ7;;tyy74II8sA'cg)Nzexp(1)r`rXrs r[ _print_Exp1RCodePrinter._print_Exp1sr]cg)Npir`rs r[ _print_PiRCodePrinter._print_Pisr]cg)Nr5r`rs r[_print_InfinityRCodePrinter._print_Infinitysr]cg)Nz-Infr`rs r[_print_NegativeInfinity$RCodePrinter._print_NegativeInfinitysr]cLSSKJn SSKJn SSKJn UR nURn[XS5(a`/nURU5H8upU"XXU 4XhU 45n URU 5n URU 5 M: SRU5$URS(a=URU5(dURU5(aURXe5$URU5n URU5n UR!U <SU <35$)Nr) Assignment) MatrixSymbol) IndexedBase rEz = )sympy.codegen.astr"sympy.matrices.expressions.matexprrsympy.tensor.indexedrlhsrhs isinstancerrrr _settingshas_doprint_loopsrf)rXrrrrrrrvr}r~tempcode0lhs_coderhs_codes r[_print_AssignmentRCodePrinter._print_Assignments0C4hhhh c ( (E77<!#d)SAY7 D) U#=99U# # ^^J 'SWW[-A-A $$&&s0 0{{3'H{{3'H&&Hh'GH Hr]cURSRS:Xa,SURURSR5-nOWSURURSR5<SURURSR5<S3nUn[ URSS5H4upESURU5<SURU5<S3U-S-nM6 U$)NrTz%szifelse(,z,NA)))argscondrrreversed)rXr last_linecodeecs r[_print_PiecewiseRCodePrinter._print_Piecewises 99R=   %t{{499R=+=+=>>I &/3kk$))B-:L:L.Mt{{[_[d[deg[h[m[mOnoI TYYs^,DAq%)[[^T[[^DTI#MD- r]cNSSKJn URURU55$)Nr) Piecewise)sympy.functionsrrrewrite)rXrrs r[ _print_ITERCodePrinter._print_ITEs-{{4<< 233r]cSRURUR[SSS9URUR URR S--5$)Nz{}[{}]AtomT)strictr)rjrparentrr~r}rrs r[_print_MatrixElement!RCodePrinter._print_MatrixElements[t00j>P 1 &&466$++*;*;A*>#>>@ @r]cf>[TU]U5nXR;aSRU5$U$)Nz(*{}))super _print_SymbolrVrj)rXrro __class__s r[rRCodePrinter._print_Symbols2w$T* $$ $>>$' 'Kr]cURUR5nURUR5nURnSR X$U5$)Nz{} {} {})rrrrel_oprj)rXrrrops r[_print_RelationalRCodePrinter._print_RelationalsB;;txx(;;txx( [[  x88r]cURUR5nURnURUR5nSR X#U5$)Nz {} {} {};)rrrrrj)rXrrrrs r[_print_AugmentedAssignment'RCodePrinter._print_AugmentedAssignmentsB;;txx( WW;;txx(!!(99r]cURUR5n[UR[5(aURR up4nO [ S5eURUR5nSRX#US- XVS9$)Nz*Only iterable currently supported is RangezCfor({target} in seq(from={start}, to={stop}, by={step}){{ {body} }}r)targetrstopstepbody) rrriterabler rNotImplementedErrorrrj)rXrrrrrrs r[ _print_ForRCodePrinter._print_Fors|T[[) dmmU + + $ 2 2 E%&RS S{{499%V]]ek!V$^3 3r]c [U[5(a1URURS55nSR U5$SnSnSnUVs/sHofR S5PM nnUVs/sH+n[ [[URU555PM- nnUVs/sH+n[ [[URU555PM- nn/n Sn [U5HEupUS;aU RU5 MXU -n U RX:-<U<35 XU - n MG U $s snfs snfs snf) z0Accepts a string of code or a list of code linesTz ){(z{ z( )rrz r)rr) rstrru splitlinesrlstripranymapendswith startswith enumerater) rXr code_linestab inc_token dec_tokenlineincreasedecreaseprettylevelns r[ruRCodePrinter.indent_codes3 dC ))$//$*?@J77:& &,  046U#6IMOSS :;<O!%(!%S)<=>!% ( GAz! d# a[ E MMCIt4 5 a[ E ' 7O(sD712D<)2E)rVrRrW)(__name__ __module__ __qualname____firstlineno____doc__ printmethodlanguagerQrrGrU__annotations__ _operatorsrNrPrbrfrlrqrwrrrrrrrrrrrrrrrrrru__static_attributes__ __classcell__)rs@r[r>r>NsEKH(,[-J-J)u O)~J $&L.% "2"$.'B ' D$J%IN 4@9 : 3r]r>Nc 6[U5RX5$)aConverts an expr to a string of r 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 R string representations. Alternatively, the dictionary value can be a list of tuples i.e. [(argument_test, rfunction_string)] or [(argument_test, rfunction_formater)]. 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]. Examples ======== >>> from sympy import rcode, symbols, Rational, sin, ceiling, Abs, Function >>> x, tau = symbols("x, tau") >>> rcode((2*tau)**Rational(7, 2)) '8*sqrt(2)*tau^(7.0/2.0)' >>> rcode(sin(x), assign_to="s") 's = sin(x);' 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"), ... (lambda x: x.is_integer, "ABS")], ... "func": "f" ... } >>> func = Function('func') >>> rcode(func(Abs(x) + ceiling(x)), user_functions=custom_functions) 'f(fabs(x) + CEIL(x))' or if the R-function takes a subset of the original arguments: >>> rcode(2**x + 3**x, user_functions={'Pow': [ ... (lambda b, e: b == 2, lambda b, e: 'exp2(%s)' % e), ... (lambda b, e: b != 2, 'pow')]}) 'exp2(x) + pow(3, x)' ``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(rcode(expr, assign_to=tau)) tau = ifelse(x > 0,x + 1,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])) >>> rcode(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(rcode(mat, A)) A[0] = x^2; A[1] = ifelse(x > 0,x + 1,x); A[2] = sin(x); )r>doprint)r assign_torYs r[rcoder$"sV  ! ) )$ ::r]c .[[U40UD65 g)z0Prints R representation of the given expression.N)printr$)rrYs r[ print_rcoder's % ! !"r]rt)r __future__rtypingrsympy.core.numbersrsympy.printing.codeprinterrsympy.printing.precedencerrsympy.sets.fancysetsr rRrWr>r$r'r`r]r[r.s:#+2<& 5 5 5   5   F  F F W 5 5 5 F F F  W!" W#$ W%&      ;F,Q;Qhk;\#r]