o
    h                     @   s  d Z ddlm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mZ ddlmZ dd	lmZ dd
lmZ ddlmZ ddlmZ ddl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$ 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/m0Z0 ddl1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7 e2dZ8e3dZ9e4dZ:e4dZ;e4dZ<e4dZ=e5dZ>e5dZ?e6dZ@e7dZAed ZBed!ZCed"ZDed#d$d%ZEed&d$d$d'ZFe*d(ZGe*d)d*d+gid,ZHe&e8eAd"d-d.ZIe&e9eAd"d-d.ZJe&eAeAd"d-d.ZKe&e@eAd"d-d.ZLe&e=eAd"d-d.ZMe&e?eAd"d-d.ZNd/d0 ZOd1d2 ZPd3d4 ZQd5d6 ZRd7d8 ZSd9d: ZTd;d< ZUd=d> ZVd?S )@zTests for sho1d.py    )Sum)oo)IInteger)S)Symbolsymbols)	factorial)exp)sqrt)Abs)KroneckerDelta)Dagger)hbar)
Commutator)qapply)InnerProduct)XPx)ComplexSpace)	represent)simplify)import_module)IndexedBaseIdx)skipraises)	RaisingOp
LoweringOpSHOKetSHOBraHamiltonianNumberOpak      bHNomegam   pTintegerq)nonnegativer/   numpyscipyfromlistsparse)import_kwargssympybasisndimformatc                  C   sD  t ttksJ ttt tdksJ ttt tdt ks$J ttt t	tj
d ttj
d   ks<J ttt t	tj
d ttj
d   ksTJ ttt t	tj
d ttj
d   kslJ td tdt	tdt t t  tdt t tt t   ksJ tjttjksJ ttd D ]} t| d | f t	| d ksJ qtstd tttddd}ttd D ]} || d | f t t	| d ksJ qtstd t!std	 tttdd
dd}ttd D ]} || d | f t t	| d ksJ q|j"dksJ |j"dks J d S )Nr%   xp   znumpy not installed.r,   r2   r8   zscipy not installed.zscipy.sparselil)r9   r:   r;   spmatrixfloat64)#r   adr#   r   doitr   r)   r   r$   r   nr   expandkzkfrewriter   r+   r*   r   r   r   hilbert_spacer   r   Infinityranger:   ad_rep_sympynpr   r   floatr3   dtype)iad_rep_numpyad_rep_scipy rS   \/var/www/html/ai/venv/lib/python3.10/site-packages/sympy/physics/quantum/tests/test_sho1d.pytest_RaisingOp9   s4   000:"&(rU   c                  C   s4  t ttksJ ttt tdksJ ttt tks J ttt t	tj
ttj
td   ks8J ttt tdksDJ ttt t	tj
ttj
td   ks\J td tdt	tdt t t  tt tt t   ksJ ttd D ]} t| | d f t	| d ksJ qd S )Nr%   r   r=   r>   )r   r#   rB   r   rC   r   r)   r   r$   r   rD   r   rE   rF   rG   rH   r   r+   r*   r   r   r   rK   r:   a_reprP   rS   rS   rT   test_LoweringOpY   s   002"rX   c                  C   s@  t tt tksJ t tt tdt ksJ t tt tdks'J ttt tj	t 
 ks6J td tt ksCJ td tdtdt t t  td tt t d   tdtd  ksnJ td ttt  tdtd  ksJ ttD ]} t| | f | ksJ qttt ksJ d S )Nr<   r   r#   r=   r%   r>   r(   )r   r)   rB   rC   r#   r   r(   r   r$   rD   rE   rH   r+   r   r*   r   r   rK   r:   N_reprL   rV   rW   rS   rS   rT   test_NumberOpe   s   &
.rZ   c                  C   s,  t tt tdksJ ttt tt tj	tdtd   t 
 ks(J td tt tt tdtd   ksAJ td tdtdt  td tt t d   ks`J td tt ttdtd   kswJ ttD ]} t| | f tt | tdtd   ksJ q{d S )Nr   r%   r>   r#   r=   r)   )r   r(   r)   rC   r   r   r$   r   r*   rD   rE   rH   rB   r#   r+   r   r   rK   r:   H_reprW   rS   rS   rT   test_Hamiltonianr   s   62*..r\   c                   C   s   t d tks
J td t ksJ ttt ttjtjks$J tj	t
tjks.J ttjdf tdks;J tdtjf tdksHJ d S )Nr$   r'   r   r%   )r   
dual_classr    r   r'   r$   rC   r   rD   rI   r   r   rJ   k3_repk3r   b3_repb3rS   rS   rS   rT   test_SHOKet|   s    rb   c                  C   s  t tttt tddf} |  tdtd tdtd  ks$J tt tttt tdtftt ddttks@J td}t	dt
d}tttt || t| |dtf dd}||t kshJ tttt tttt tdtf t || t| |dtf dd}||t ksJ tt- tttt tttt tdtf t |t tt tdtf dd}W d    d S 1 sw   Y  d S )Nr   r%   Tsum_doitr#   rD   )cls)r   r   r-   r    rC   r   r   r0   r   r   r   r   
ValueError)e1r#   rD   resultrS   rS   rT   test_sho_sums   s.   ,",>
>"ri   c                  C   s   t ddd} tt| d  td t| t ttt tt tdt	f }t
tt| ddtt| d  td | t  ttt ksJJ tt
ttt | ddtt
tt|  | ddkshJ d S )NalphaT)
is_complexr>   r   rc   )r   r
   r   r   r   r-   r   r	   r   r   r   r    r0   r   r#   )rj   cstaterS   rS   rT   test_sho_coherant_state   s   BF@rm   c                  C   s   t dddd} t ddd}t dddd}t dddd}t| d	|   | t| t| |d
tf}tt|| t| dd}t|| | d	  | t|| ksRJ d S )NnbarT)realr1   rD   r.   rP   )r/   r1   jr%   r   rc   )r   r   r   r    r   r   r   r   )rn   rD   rP   rp   rhorh   rS   rS   rT   test_issue_26495   s   ,*rr   N)W__doc__sympy.concreter   
sympy.corer   sympy.core.numbersr   r   sympy.core.singletonr   sympy.core.symbolr   r   (sympy.functions.combinatorial.factorialsr	   &sympy.functions.elementary.exponentialr
   (sympy.functions.elementary.miscellaneousr   $sympy.functions.elementary.complexesr   (sympy.functions.special.tensor_functionsr   sympy.physics.quantumr   sympy.physics.quantum.constantsr   r   sympy.physics.quantum.qapplyr   "sympy.physics.quantum.innerproductr   sympy.physics.quantum.cartesianr   r   sympy.physics.quantum.hilbertr   sympy.physics.quantum.representr   sympy.simplifyr   sympy.externalr   sympy.tensorr   r   sympy.testing.pytestr   r   sympy.physics.quantum.sho1dr   r   r   r    r!   r"   rB   r#   r$   rF   rG   r_   r'   ra   r(   r)   r*   r+   r:   r-   r0   rM   r3   rL   rV   rY   r[   r^   r`   rU   rX   rZ   r\   rb   ri   rm   rr   rS   rS   rS   rT   <module>   sn      
