\h%SrSSKJr SSKJr SSKJr SSKJr SSK J r SSK J r SSK JrJr SS KJr SS KJr SS KJr SS KJr SS KJrJr SSKJrJrJrJr SSK J!r! S/r"Sr#Sr$Sr%Sr&g)z}Logic for applying operators to states. Todo: * Sometimes the final result needs to be expanded, we should do this by hand. )Sum)Add) NumberKind)Mul)Pow)S)sympify_sympify)AntiCommutator) Commutator)Dagger) InnerProduct) OuterProductOperator)StateKetBaseBraBase Wavefunction) TensorProductqapplyc0UR[S5$)zETransform the inner products in an expression by calling ``.doit()``.c.[U6R5$N)rdoitargss T/var/www/auris/envauris/lib/python3.13/site-packages/sympy/physics/quantum/qapply.pyip_doit_func..#st1D1I1I1K)replaceres r ip_doit_funcr$!s 99\#K LLr c0UR[S5$)z;Transform the sums in an expression by calling ``.doit()``.c.[U6R5$r)rrrs rrsum_doit_func..(sT (9r )r!rr"s r sum_doit_funcr(&s 99S9 ::r c SSKJn URSS5nURSS5nURSS5n[U5nUR[ :XaU(a [ U5$U$URSSS9n[U[5(aU$[U[5(a4SnURHnU[U40UD6- nM UR5$[X5(a2URVV s/sHunn [U40UD6U 4PM n nn U"U 6$[U[5(a-[URV s/sHn [U 40UD6PM sn 6$[U[5(aB[[UR40UD6/UR Q76nU(a [#U5nU$UnU$[U[$5(a#[UR&40UD6UR(-$[U[*5(aUR-5up[+U 6n[+U 6nU (dUnO4[U[*5(aU[/U40UD6-nOU[U40UD6-nX`:Xa%U(a[1[/[1U540UD65nU(a [ U5OUnU(a [#U5nU$UnU$U$s sn nfs sn f) a]Apply operators to states in a quantum expression. Parameters ========== e : Expr The expression containing operators and states. This expression tree will be walked to find operators acting on states symbolically. options : dict A dict of key/value pairs that determine how the operator actions are carried out. The following options are valid: * ``dagger``: try to apply Dagger operators to the left (default: False). * ``ip_doit``: call ``.doit()`` in inner products when they are encountered (default: True). * ``sum_doit``: call ``.doit()`` on sums when they are encountered (default: False). This is helpful for collapsing sums over Kronecker delta's that are created when calling ``qapply``. Returns ======= e : Expr The original expression, but with the operators applied to states. Examples ======== >>> from sympy.physics.quantum import qapply, Ket, Bra >>> b = Bra('b') >>> k = Ket('k') >>> A = k * b >>> A |k>>> qapply(A * b.dual / (b * b.dual)) |k> >>> qapply(k.dual * A / (k.dual * k)) QXX>*2v& 9?  As  aff((!%%// As  **,V gF  $ $:f888F6&4G44F ;6Jvayqapply_Mul..s4-{rzknjxPSUX>Y.Z.f^aef^f.frz9;c3z# UH1n[U[[[[45=(d US:Hv M3 g7frIrJrKs rrMrNs73Awps:cHeUXZ]C^3_3kcfjkck3kwrOT)r/z4Duplicated dummy indices in separate sums in qapply._apply_operator_apply_from_right_to).listrrOnelenr4rpoprr is_commutativerr8 is_Integerappendr7rketbrar r rrrfuncrallranger3r:rset variables intersection ValueErrorr6r5getattrNotImplementedErrorrrrintcomplexfloatr ) r#r;rextrar<rhslhscommnr6_apply _apply_rights rr:r:sT g> >LL 0g> >#sVCLL,88F3::F!&&$-.:'::S+T 2F  C+7+F ~s$:DA  # %c5W5~ c7 # # 3(@(@!#+F&3/00  t9>Haffte|5AA#EeK K FL ) )j:'::5@@affdmF*6g6u<rs] & "0?7/;AMM=  M ; tne=r