o X۷i @s^ddlmZddlZddlZddlZddlmZddZgdZdZ dd Z ej d d Z dS) ) annotationsN)_utilcCs|t|t||jd|jdkrtd|dddtjddf|dtjddddf}|d|jddS)av Khatri-rao product A column-wise Kronecker product of two matrices Parameters ---------- a : (n, k) array_like Input array b : (m, k) array_like Input array Returns ------- c: (n*m, k) ndarray Khatri-rao product of `a` and `b`. See Also -------- .. seealso:: :func:`scipy.linalg.khatri_rao` z6The number of columns for both arrays should be equal..N))r _assert_2dshape ValueErrorcupynewaxisreshape)abcrR/home/ubuntu/vllm_env/lib/python3.10/site-packages/cupyx/scipy/linalg/_matfuncs.py khatri_rao s  4r)gD`lCgD`\Cg`=Hb;Cg eCgJXBg"5Bg/cBg\L8BgpķAgsyAgS-Ag@gf@g?g{Gz@cCs|jdkr tjd|jdS|jd}t|jtjr|jntj}t| |}|tj ||d|}tj j |dd }|tk}|rUtttt|td}nd}|d|}||}||} || } tj |jd|d} tjt|d} t| || | | \} }}}|| | |}| ||}tj | |||}|}t|D]}||}qt|jtjrt|nt|}||9}|S)aCompute the matrix exponential. Parameters ---------- a : ndarray, 2D Returns ------- matrix exponential of `a` Notes ----- Uses (a simplified) version of Algorithm 2.3 of [1]_: a [13 / 13] Pade approximant with scaling and squaring. Simplifications: * we always use a [13/13] approximate * no matrix balancing References ---------- .. [1] N. Higham, SIAM J. MATRIX ANAL. APPL. Vol. 26(4), p. 1179 (2005) https://doi.org/10.1137/04061101X r)rr)dtyper)ordr)sizer zerosrr issubdtypeinexactfloat64diagsumeyelinalgnormitemth13intmathceillog2floatasarrayr _expm_innersolverangecomplexfloatingcmathexp)r na_dtypemuAnrmAscalesA2A4A6Ebbu1u2v1v2uvr13x_emurrrexpm@sF   "    rCc Cs|d||d||d|}|d||d||d||d|}|d||d ||d |}|d ||d ||d ||d|}||||fS)N r rrr) r7r4r5r6rr9r:r;r<rrrr's $0$0 r') __future__rr"r+r cupy.linalgrrrr rCfuser'rrrrs  $N