o
    oh,                     @   s  d dl mZmZmZmZmZmZmZmZm	Z	m
Z
mZmZmZmZmZmZmZmZ d dlmZmZmZmZmZmZmZmZmZ d dlmZ d dlm Z  d dl!m"Z" d dl#m$Z$m%Z% d dl&m'Z'm(Z(m)Z) d dl*m+Z+ ed	\Z,Z-Z.d
d Z/dd Z0dd Z1dd Z2dd Z3dd Z4dd Z5dd Z6dd Z7dd Z8dd Z9d d! Z:d"d# Z;d$d% Z<d&d' Z=d(d) Z>d*d+ Z?d,d- Z@d.d/ ZAd0d1 ZBd2d3 ZCd4d5 ZDd6d7 ZEd8d9 ZFd:d; ZGd<S )=    )SpioosymbolsRationalIntegerGoldenRatio
EulerGammaCatalanLambdaDummyEqNeLeLtGtGeMod)		PiecewisesincosAbsexpceilingsqrtsignfloor)ITE)raises)implemented_function)IndexedBaseIdx)MatrixSymbolSparseMatrixMatrix	rust_codezx,y,zc                   C   s,   t tddks
J t tddksJ d S )N*   42iz-56)r&   r    r)   r)   r/var/www/html/construction_image-detection-poc/venv/lib/python3.10/site-packages/sympy/printing/tests/test_rust.pytest_Integer   s   r+   c                   C   s   t tttdksJ t tttdksJ t tttdks!J t tttdks,J t tttdks7J t tttdksBJ d S )Nzx == yzx != yzx <= yzx < yzx > yzx >= y)	r&   r   xyr   r   r   r   r   r)   r)   r)   r*   test_Relational   s   r.   c                   C   s   t tdddksJ t tdddksJ t tdddks!J t td	ddks,J t ttdd d
ks9J t tddt dksFJ d S )N      z	3_f64/7.0   	   2iz
-3_f64/7.0zx + 3_f64/7.0z(3_f64/7.0)*x)r&   r   r,   r)   r)   r)   r*   test_Rational   s   r5   c                   C   sf   t tt dks
J t tt dksJ t tt dksJ t tt dks(J t t dks1J d S )Nzx + yzx - yzx*yzx/yz-x)r&   r,   r-   r)   r)   r)   r*   test_basic_ops(   s
   r6   c                  C   sH   G dd dt } t| tdksJ tddd}t|d dks"J d S )	Nc                   @   s   e Zd Zdd ZdS )ztest_printmethod.<locals>.fabsc                 S   s   d| | jd  S )Nz	%s.fabs()r   )_printargs)selfprinterr)   r)   r*   
_rust_code2   s   z)test_printmethod.<locals>.fabs._rust_codeN)__name__
__module____qualname__r;   r)   r)   r)   r*   fabs1   s    r?   zx.fabs()a   r/   )r   r   za[0])r   r&   r,   r"   )r?   r@   r)   r)   r*   test_printmethod0   s   rB   c                   C   sr   t tttt dksJ t ttdksJ t ttdks"J t ttdks,J t ttddks7J d S )Nzx.sin().powf(x.cos())zx.abs()x.ceil()z	x.floor()r/   zx - 3*((1_f64/3.0)*x).floor())r&   r   r,   r   absr   r   r   r)   r)   r)   r*   test_Functions9   s
   rE   c                  C   sJ  t dt dks
J t td t td   krdks J  J t ttdks*J t ttj t td   kr>dksAJ  J t dtt dksMJ t ttj  t td   krbdkseJ  J t dt d	ksoJ t td t td   krd	ksJ  J t td d
ksJ t ttdd dksJ t dt dksJ t ttdksJ t td dksJ t ttd  dksJ t ttdd dksJ t	dt
tdt } t d| td ttt    td t  dksJ dd ddfdd ddfg}t td d|iddksJ t td d|iddks#J d S ) NrA   z	x.recip()g      zx.sqrt()g      ?zx.sqrt().recip()g      z
PI.recip()zPI.sqrt().recip()r/   zx.cbrt()   zx.exp2()zx.exp()z	x.powi(3)zx.powf(y.powi(3))zx.powf(2_f64/3.0)gg      @z.(3.5*2*x).powf(-x + y.powf(x))/(x.powi(2) + y)c                 S   s   |j S N
is_integerbaser   r)   r)   r*   <lambda>Z       ztest_Pow.<locals>.<lambda>dpowic                 S   s   |j  S rI   rJ   rL   r)   r)   r*   rN   [       powPowuser_functionsz
x.dpowi(3)g	@z
x.pow(3.2))r&   r,   r   r   Halfr   r   r   r-   r   r   )rH   _cond_cfuncr)   r)   r*   test_PowC   s0   ,.0,*"rX   c                   C   s   t tdksJ t tdksJ t tjdksJ t t dks"J t tjdks+J t tjdks4J t tddks>J t tjdksGJ d S )NPIINFINITYNEG_INFINITYNANrA   E)	r&   r   r   r   InfinityNegativeInfinityNaNr   Exp1r)   r)   r)   r*   test_constants`   s   rb   c                   C   s^   t dt dtd ksJ t dt dtd ksJ t dt dtd ks-J d S )NrG   z*const GoldenRatio: f64 = %s;
2*GoldenRatio   z"const Catalan: f64 = %s;
2*Catalanz(const EulerGamma: f64 = %s;
2*EulerGamma)r&   r   evalfr
   r	   r)   r)   r)   r*   test_constants_otherk   s   "re   c                   C   s   t ddksJ t tjdksJ t ddksJ t tjdks"J t tt@ dks,J t ttB dks6J t t dks?J t tt@ t@ dksKJ t ttB tB d	ksWJ t tt@ tB d
kscJ t ttB t@ dksoJ d S )NTtrueFfalsezx && yzx || yz!xzx && y && zzx || y || zzz || x && yzz && (x || y))r&   r   rf   rg   r,   r-   zr)   r)   r)   r*   test_booleanr   s   ri   c                      s  t ttdk ftd df t dksJ t dddksJ t dddd	ks*J t ttdk ftd td
k ftd df t dddksHJ t ddddksSJ t dddks]J dt ttdk ftd td
k ftd df  t dddks}J dt ttdk ftd td
k ftd df d  t dddksJ t ttdk ftd tdkftttdkf tt fdd d S )NrA   rG   Tz'if (x < 1) {
    x
} else {
    x + 2
}r	assign_toz,r = if (x < 1) {
    x
} else {
    x + 2
};)rl   inlinez$r = if (x < 1) { x } else { x + 2 };   )rm   z9if (x < 1) { x } else if (x < 5) { x + 1 } else { x + 2 }z>r = if (x < 1) { x } else if (x < 5) { x + 1 } else { x + 2 };zJr = if (x < 1) {
    x
} else if (x < 5) {
    x + 1
} else {
    x + 2
};z;2*if (x < 1) { x } else if (x < 5) { x + 1 } else { x + 2 }r'   z@2*if (x < 1) { x } else if (x < 5) { x + 1 } else { x + 2 } - 42r   c                      s   t  S rI   r%   r)   exprr)   r*   rN      rQ   z test_Piecewise.<locals>.<lambda>)r   r,   r&   r   r   
ValueErrorr)   r)   ro   r*   test_Piecewise   s<   
(

,
0
,rr   c                  C   s.   t t tt t } t| tgddksJ d S )Ndereferencezx + y + (*z) + (*z).sin())r,   r-   r   rh   r&   ro   r)   r)   r*   test_dereference_printing   s   ru   c                  C   s   t tt } t| dksJ t| dddksJ t tt d } t| dks(J t| dddks2J t tt} t| dks@J d S )	Nzy*x.signum()rj   rk   zr = y*x.signum();r'   z(x + y).signum() + 42zr = (x + y).signum() + 42;zx.cos().signum())r   r,   r-   r&   r   ro   r)   r)   r*   	test_sign   s   rv   c                  C   s   t d\} }t|}t|dksJ t||gddksJ t|dddks'J tt t|dd	 W d    d S 1 s=w   Y  d S )
Nzx ifz	if_.sin()rs   z(*if_).sin()_unreserved)reserved_word_suffixzif_unreserved.sin()T)error_on_reserved)r   r   r&   r   rq   )r,   r-   rp   r)   r)   r*   test_reserved_words   s   
"rz   c                  C   s$   t tdk tt} t| dksJ d S )NrA   z#if (x < 1) {
    y
} else {
    z
})r   r,   r-   rh   r&   ro   r)   r)   r*   test_ITE   s   r{   c            	      C   s   t ddd\} }}td| td|td|}}}td| }t|dks(J td	||f }t|d
ks8J td|||f }t|dksIJ d S )Nzn m oTintegerijkr,   zx[j]Az
A[m*i + j]BzB[m*o*i + o*j + k])r   r!   r    r&   )	nmor~   r   r   r,   r   r   r)   r)   r*   test_Indexed   s   "r   c                  C   sL   t ddtd\} }td}td}t| |} t||  ||  ddks$J d S )Nzi mT)r}   clsr,   r-   rk   z"for i in 0..m {
    y[i] = x[i];
})r   r   r    r!   r&   )r~   r   r,   r-   r)   r)   r*   test_dummy_loops   s   
r   c                  C   s   t ddd\} }td}td}td}td}td| }td	|}t|||f ||  || d
dks6J t|||f ||  ||  ||  || d
dksRJ d S )Nzm nTr|   r   r,   r-   rh   r~   r   rk   zrfor i in 0..m {
    y[i] = 0;
}
for i in 0..m {
    for j in 0..n {
        y[i] = A[n*i + j]*x[j] + y[i];
    }
}z|for i in 0..m {
    y[i] = x[i] + z[i];
}
for i in 0..m {
    for j in 0..n {
        y[i] = A[n*i + j]*x[j] + y[i];
    }
}r   r    r!   r&   )r   r   r   r,   r-   rh   r~   r   r)   r)   r*   
test_loops   s   

.
r   c                  C   s   t ddd\} }}}td}td}td}td|}td| }td	|}	td
|}
t|||	|
f ||||	|
f  || dd||  | | || |  |	|  |
 || | |	|  |
 f d kseJ d S )Nzn m o pTr|   r@   br-   r~   r   r   lrk   zfor i in 0..m {
    y[i] = 0;
}
for i in 0..m {
    for j in 0..n {
        for k in 0..o {
            for l in 0..p {
                y[i] = a[%s]*b[%s] + y[i];
            }
        }
    }
}r   )r   r   r   pr@   r   r-   r~   r   r   r   r)   r)   r*    test_loops_multiple_contractions	  s   



(>r   c                  C   s  t ddd\} }}}td}td}td}td}td| }td	|}	td
|}
td|}t|||	|
|f |||	|
|f  ||	|
|f  || d}|d|| | | |	| |  |
|  | || | | |	| |  |
|  | |	| | |
|  | f d ksJ d S )Nzm n o pTr|   r@   r   cr-   r~   r   r   r   rk   zfor i in 0..m {
    y[i] = 0;
}
for i in 0..m {
    for j in 0..n {
        for k in 0..o {
            for l in 0..p {
                y[i] = (a[%s] + b[%s])*c[%s] + y[i];
r   r   )r   r   r   r   r@   r   r   r-   r~   r   r   r   coder)   r)   r*   test_loops_addfactor"  s"   



:dr   c                   C   s   t tdd  d S )Nc                   S   s   t ttddS )Ngarbage)method)r&   r   r,   r)   r)   r)   r*   rN   >  s    ztest_settings.<locals>.<lambda>)r   	TypeErrorr)   r)   r)   r*   test_settings=  s   r   c                  C   s   t d} tdt| d|  }t|| dksJ tdt| d|  t }t|| dtd ks3J td}tdt d	d
d}tdt| | d|   d|   }t||| || ddks`J d S )Nr,   rH   rG   z2*xz$const Catalan: f64 = %s;
2*x/Catalanrc   r   r~   r   Tr|   rA   rk   z8for i in 0..n {
    A[i] = (A[i] + 1)*(A[i] + 2)*A[i];
})r   r   r   r&   r
   rd   r    r!   )r,   rH   r   r~   r)   r)   r*   test_inline_functionA  s   
 r   c                  C   s   t ddd} t ddd}ddd d	d
fdd dd
fgd}tt| |ddks)J tt| |ddks5J tt||ddksAJ d S )Nr,   Fr|   r   Tceilc                 S   s   | j  S rI   rJ   r,   r)   r)   r*   rN   X  rQ   z%test_user_functions.<locals>.<lambda>r?      c                 S   s   | j S rI   rJ   r   r)   r)   r*   rN   X  rO   rD   )r   r   rT   rC   zfabs(x)zabs(n))r   r&   r   r   )r,   r   custom_functionsr)   r)   r*   test_user_functionsS  s   r   c                   C   sX   t tg ddksJ tt t tg dg W d    d S 1 s%w   Y  d S )NrA   rG   r/   z	[1, 2, 3])r&   r$   r   rq   r)   r)   r)   r*   test_matrix_  s   
"r   c                   C   s@   t t ttg dg W d    d S 1 sw   Y  d S )Nr   )r   NotImplementedErrorr&   r#   r)   r)   r)   r*   test_sparse_matrixe  s   
"r   N)H
sympy.corer   r   r   r   r   r   r   r	   r
   r   r   r   r   r   r   r   r   r   sympy.functionsr   r   r   r   r   r   r   r   r   sympy.logicr   sympy.testing.pytestr   sympy.utilities.lambdifyr   sympy.tensorr    r!   sympy.matricesr"   r#   r$   sympy.printing.rustr&   r,   r-   rh   r+   r.   r5   r6   rB   rE   rX   rb   re   ri   rr   ru   rv   rz   r{   r   r   r   r   r   r   r   r   r   r   r)   r)   r)   r*   <module>   sD   P ,			
(
