o o h6+@s<ddlmZddlmZddlmZmZddlmZddl m Z ddl m Z ddl mZddlmZdd lmZdd lmZmZmZmZmZmZmZdd lmZmZmZdd lm Z dd l!m"Z"m#Z#ddl$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z*m+Z+ddl,m-Z-ddl.m/Z/m0Z0m1Z1m2Z2ddl3m4Z4m5Z5m6Z6ddl7m8Z8ddl9m:Z:m;Z;ddle>re>j?@e>j?AgjBZCeCjDZEddZFddZGddZHddZIddZJd d!ZKd"d#ZLd$d%ZMd&d'ZNd(d)ZOd*d+ZPd,d-ZQd.d/ZRd0d1ZSd2d3ZTd4d5ZUd6d7ZVd8d9ZWd:d;ZXdd?ZZd@dAZ[dBdCZ\dDdEZ]dFdGZ^dHS)I)Sum)Mod)Equality Unequality)sqrt) Piecewise) BlockMatrix) MatrixSymbol)Identity)lambdify)xijabcd)FunctionPowSymbol) MatrixSolve) logaddexp logaddexp2)log1pexpm1hypotlog10exp2log2Sqrt)Array)ArrayTensorProductArrayAdd PermuteDims ArrayDiagonal) JaxPrinter_jax_known_constants_jax_known_functions)convert_matrix_to_array)skipraises) import_modulejaxcCsDt}tdtdkfd}||dksJ|jdhdiks JdS)z NumPyPrinter needs to print Piecewise()'s choicelist as a list to avoid breaking compatibility with numpy 1.8. This is not necessary in numpy 1.9+. See gh-9747 and gh-9749 for details. r)rTzKjax.numpy.select([jax.numpy.less(x, 0),True], [1,0], default=jax.numpy.nan)z jax.numpy>nanlessselectN)r%rr doprintmodule_imports)printerpr5q/var/www/html/construction_image-detection-poc/venv/lib/python3.10/site-packages/sympy/printing/tests/test_jax.pytest_jax_piecewise_regression(s r7cCs@ttt}t|dksJttt}t|dksJdS)Nzjax.numpy.logaddexp(a, b)zjax.numpy.logaddexp2(a, b))rrrr%r1r)laelae2r5r5r6test_jax_logaddexp5s  r:c ststdttttttf}ttttf|d}d\}}tj dddtj |||t fddt ||dDs@Jttttttf}ttttf|d}d\}}tj dddtj |||t fd dt ||dDszJdS) NJAX not installedr,r r-r=c3s|]}|VqdSNr5.0i_x_r5r6 Eztest_jax_sum..c3s|]}|VqdSr?r5r@rCr5r6rELrF) r,r)rr r rrr numpylinspaceallclosesumrangesfa_b_r5rCr6 test_jax_sum<s6:rQc ststdtttttttfttt f}t tttt tf|d}d\}}d\tj dddtj |||tfdd t||dDsPJdS) Nr;r,r<) r>r-r=c3s0|]}tdD] }||Vq qdS)r-N)rK)rArBj_c_d_rDr5r6rEZs.z)test_jax_multiple_sums..)r,r)rr rr rrrrr rGrHrIrJrKrLr5rUr6test_jax_multiple_sumsOs " rXcCststdtddd}tddd}t||}t||f|d}tjddgddgg}tjdd gd dgg}|||tj||ksHJdS) Nr;MNr,r-r>) r,r)r r(r rGarraymatmulall)rYr[cgrNmambr5r5r6test_jax_codegen_einsum]s   $rec Cststdtddd}tddd}tddd}tddd}tjddgdd gg}tjdd gd dgg}tjdd gddgg}tjdd gd d gg}t||}t||f|d} | ||tj|d dg|ddgksqJt ||}t||f|d} | ||||ksJt |||}t|||f|d} | ||||||ksJt ||||}t||||f|d} | ||||||||ksJt |dd g}t|f|d} | ||j ksJt t||gd}t||f|d} | ||tj tj|d dg|ddgdksJt t||d}t||f|d} | ||tjjtj|d dg|ddgdddksCJdS)Nr;rYrZr[PQr-r\r]r^r>rr,)r-rZr\rr-rZ)axis1axis2)r,r)r rGr_r!r einsumrar"r#T transposer$diagonal) rYr[rfrgrcrdmcmdrbrNr5r5r6test_jax_codegen_extrals>     ,   &8@rrcCs&tstdttd}ttf|d}tjgd}tj||gds'Jttd}ttf|d}tjgd}tj||gdsHJtdk}ttf|d}tjgd}tj||gdshJtdk}ttf|d}tjgd}tj||gdsJtdk}ttf|d}tjgd}tj||gd sJtdk}ttf|d}tjgd}tj||gd sJtdktd k@}ttf|d}tjgd}tj||gdsJtdktd kB}ttf|d}tjgd}tj||gd sJdS) Nr;r-r,)rr-rZ)FTF)TFT)TFF)TTF)FFT)FTTrZ)TTT) r,r)rr r rGr_ array_equalr)erNrDr5r5r6test_jax_relationalsD   rucCststdttt}tttf|d}tjgd}d}tj|||gds+Jtjgd}tjgd}tj|||gdsIJtjgd}tjgd}tj|||gdsgJdS) Nr;r,)rr-rZr\rZ)rr-rr-)rZrZrZrZ)rZr\r])rrrr) r,r)rrrr rGr_rs)rtrNrOrPr5r5r6 test_jax_mods  rwcCs8tstdtdddd}tg|d}|dksJdS)Nr;rZr>F)evaluater,g?)r,r)rr )exprrNr5r5r6 test_jax_pows  rzcCsBtstdttfttd}t|ddddtksJdS)Nr;r,g|=g#B ;)r,r)r rrabsJAX_DEFAULT_EPSILONrNr5r5r6test_jax_expm1s$r~cCs>tstdttfttd}t|dddtksJdS)Nr;r,g>N}a+)r,r)r rrr{r|r}r5r5r6test_jax_log1ps rcCs<tstdttttftttddddtksJdS)Nr;r,r\r]rv)r,r)r{r rrrr|r5r5r5r6test_jax_hypots0rcC6tstdtttfttdddtksJdS)Nr;r,drZ)r,r)r{r rrr|r5r5r5r6test_jax_log10s*rcCr)Nr;r,rv )r,r)r{r rrr|r5r5r5r6 test_jax_exp2*rcCr)Nr;r,)r,r)r{r rrr|r5r5r5r6 test_jax_log2rrcCrNr;r,r]rZ)r,r)r{r rrr|r5r5r5r6 test_jax_Sqrt rrcCrr)r,r)r{r rrr|r5r5r5r6 test_jax_sqrtrrcCststdtddd}tddd}|d||}t|||}t||f|d}t||f|d}tjgdgd gd g}tjj|dksKJtjgd }tj ||||||sbJdS) Nr;rYr\r r-r>r,)r-rZr\)r\rZrv)rvrh)r\r]rv) r,r)r rr rGr_linalg matrix_rankrI)rYr ry matsolve_exprrN f_matsolvem0x0r5r5r6test_jax_matsolves  "rcCsztstdtddd}tddd}tddd}tddd}t||g||gg}|jdks/Jt}||d ks;JdS) Nr;a_1r=r\a_2a_3a_4)rz)jax.numpy.block([[a_1, a_2], [a_3, a_4]]))r,r)r rshaper%r1)rrrrAr3r5r5r6 test_16857.s    rcststdtddd}t||tdd}tjddgddgg}tjddgddgg}|||ks8Jd d lm }|d d d tdt t fdddS)Nr;rYrZr,r-r\r]rvr)symbolsnT)integercsttdS)Nr,)r r r5r[rr5r6Lsz"test_issue_17006..) r,r)r r r rGr_rasympy.core.symbolrr*NotImplementedError)rYrNrcmrrr5rr6test_issue_17006=s    rcCs4ttddks JttddksJdS)N)ri)r\rvz!jax.numpy.array([[1, 2], [3, 5]])rizjax.numpy.array((1, 2)))r%r1r r5r5r5r6test_jax_arrayOsrcCsDtddksJtddksJtddksJtddks JdS) NNaNz jax.numpy.nan EulerGammazjax.numpy.euler_gammaacoszjax.numpy.arccoslogz jax.numpy.log)r&r'r5r5r5r6test_jax_known_funcs_constsTsrcCs&t}t|ds Jt|dsJdS)N _print_acos _print_log)r%hasattr)prntrr5r5r6test_jax_print_methods\rcCs&t}t|ds J|jdksJdS)N printmethod_jaxcode)r%rr)r3r5r5r6test_jax_printmethodbrrcCs4Gdddt}t}||tddksJdS)Nc@seZdZddZdS)z+test_jax_custom_print_method..expm1cSs.|j\}d||d}||jd|S)Nzexpm1().)args_print_module_format_module)selfr3r functionr5r5r6rlsz4test_jax_custom_print_method..expm1._jaxcodeN)__name__ __module__ __qualname__rr5r5r5r6rjs rr zjax.numpy.expm1(x))rr%r1r)rr3r5r5r6test_jax_custom_print_methodhsrN)_sympy.concrete.summationsrsympy.core.modrsympy.core.relationalrr(sympy.functions.elementary.miscellaneousr$sympy.functions.elementary.piecewiser&sympy.matrices.expressions.blockmatrixr"sympy.matrices.expressions.matexprr "sympy.matrices.expressions.specialr sympy.utilities.lambdifyr sympy.abcr r rrrrr sympy.corerrrsympy.codegen.matrix_nodesrsympy.codegen.numpy_nodesrrsympy.codegen.cfunctionsrrrrrrrsympy.tensor.arrayr 0sympy.tensor.array.expressions.array_expressionsr!r"r#r$sympy.printing.numpyr%r&r'3sympy.tensor.array.expressions.from_matrix_to_arrayr(sympy.testing.pytestr)r*sympy.externalr+r,rGfinfor_dtypedeafult_float_infoepsr|r7r:rQrXrerrrurwrzr~rrrrrrrrrrrrrrrr5r5r5r6sb       $ $     *3