o 0 i0@sLddlZddlZddlmZddlmZ ddZ ddZ dd Z d d Z dS) N)cublas)devicec Cs0ddlm}|jdkrtd|j|jdvr td|j|jd|jdkr.td|jd|jdkrBtd |j|j|j|jkrRtd |j|j|j}|d kr\d }n|d krcd }n|dkrjd}n|dkrqd}ntd|jt||d}t||d}t||d}|jd}|jdkr|jdnd} |j rt j } n |j rt j } ntd|j stdt} tj|tjd} tjdtjd} || |||jj|}tj||jd}|| |||jj||jj| jj| jjtjj|| || | || |jj|| jj|jj|| jj tjj|| dS)aSolve a linear matrix equation using cusolverDngetr[fs](). Computes the solution to a system of linear equation ``ax = b``. Args: a (cupy.ndarray): The matrix with dimension ``(M, M)``. b (cupy.ndarray): The matrix with dimension ``(M)`` or ``(M, K)``. Returns: cupy.ndarray: The matrix with dimension ``(M)`` or ``(M, K)``. Note: ``a`` and ``b`` will be overwritten. rcusolvera.ndim must be 2 (actual: {}))r"b.ndim must be 1 or 2 (actual: {})rza must be a square matrix.shape mismatch (a: {}, b: {}).zdtype mismatch (a: {}, b: {})fsdFcDzzunsupported dtype (actual:{})Zgetrf_bufferSizegetrfgetrsz'a must be F-contiguous or C-contiguous.zb must be F-contiguous.dtypeN)cupy_backends.cuda.libsrndim ValueErrorformatshaper TypeErrorgetattrZ _f_contiguous_cublas CUBLAS_OP_NZ _c_contiguous CUBLAS_OP_T_deviceget_cusolver_handle_cupyempty_numpyint32dataptrlinalg_util3_check_cusolver_dev_info_if_synchronization_allowed)ab _cusolverrthelperrrnnrhstranshandleZdipivZdinfoZlworkZdworkr4X/home/app/PaddleOCR-VL-test/.venv_paddleocr/lib/python3.10/site-packages/cupyx/lapack.pygesvsl       r6cCs\ddlm}|jdkrtd|j|jdkrd}n|jdkr&|jd}ntd|j|jd|jdkrBtd|j|j|j|jkrRtd|j|j|j}|d kr\d }n|d krcd }n|d krjd }n|dkrqd}ntd|t||d}t||d}tt|d}|dvrt||d} t||d} nt||d} t||d} tj } |j dvrtj } ntj } |j\} }t | |}tjdtjd}tj||d}t}t}tjd|d}| |kr|jd d}|jd d}||| ||jj| }tj||d}||| ||jj| |jj|jj||jj tjj||| |tj| | |||jj| |jj|jj| }tj||d}| |tj| | |||jj| |jj|jj| |jj||jjtjj| |||tjtj| tj|||jj|jj| |jj| |d|S| j!jd d}|}|jdkr|fn||f}tj"||d d}||d| <|||| |jj|}tj||d}|||| |jj||jj|jj||jj tjj||||tjtj| tj| ||jj|jj||jj| | |tj| ||||jj||jj|jj| }tj||d}| |tj| ||||jj||jj|jj||jj||jjtjj| ||S)aSolves over/well/under-determined linear systems. Computes least-square solution to equation ``ax = b` by QR factorization using cusolverDngeqrf(). Args: a (cupy.ndarray): The matrix with dimension ``(M, N)``. b (cupy.ndarray): The matrix with dimension ``(M)`` or ``(M, K)``. Returns: cupy.ndarray: The matrix with dimension ``(N)`` or ``(N, K)``. rrrrrr r zdtype mismatch (a: {}, b: {}).r r r rrrrzunsupported dtype (actual: {})Zgeqrf_bufferSizegeqrftrsmsdZormqr_bufferSizeormqrZunmqr_bufferSizeZunmqrfdrg?)orderN)rr<)#rrrrrrrrrrcharrZ CUBLAS_OP_Cminr"r#r$r%r r!Zget_cublas_handlearraycopyr&r'r(r)r*ZCUBLAS_SIDE_LEFTZCUBLAS_FILL_MODE_UPPERZCUBLAS_DIAG_NON_UNITctypesconjTZzeros)r+r,r-r1rr.Z geqrf_helperr7r8Z ormqr_helperr:Zno_transr2mr0Zmn_mindev_infotauZcusolver_handleZ cublas_handleoneZws_size workspaceZbbZ out_shaper4r4r5gelsPs                 rIc Csddlm}ddl}|jdstdt|j|j}t|d}|dkr-|j}|j }n+|dkr8|j }|j }n |dkrC|j }|j }n|dkrN|j}|j}n d |j}t||j|d d d }tj|}|jd d\} } tt|jdd } t} tj} tj| tjd}|| | | |jj | |jj | tj!j"#|||j}|$%| | dj|d d d }tj|}|jd d\}}tjdtjd}|| | | ||jj | |jj ||jj | tj!j"#|||$%|S)Nrr potrsBatchedzpotrsBatched is not availabler r rrDdtype must be float32, float64, complex64 or complex128 (actual: {})CTr<r@rr)&rrZcupyx.cusolverZcheck_availability RuntimeErrorr$ promote_typesrZ spotrfBatchedZ spotrsBatchedZ dpotrfBatchedZ dpotrsBatchedZ cpotrfBatchedZ cpotrsBatchedZ zpotrfBatchedZ zpotrsBatchedrrastyper"Z_coreZ _mat_ptrsrintprodr r!rCUBLAS_FILL_MODE_LOWERr#r%r&r'r(r)r*rBreshape)r+r,r-ZcupyxrZ potrfBatchedrJmsgZapldar0Z batch_sizer3uplorEb_shapebpldbr1r4r4r5 _batched_posvs^     r]c Csddlm}tjj}||||||||jdkr$t ||St |j |j }t |d}|dkr@|j }|j}|j}n4|dkrN|j}|j}|j}n&|dkr\|j}|j}|j}n|dkrj|j}|j}|j}n d|j }t||j|dd d }|j\} } t} tj } tj!d t j"d } || | | |j#j$| }tj!||d }|| | | |j#j$| |j#j$|| j#j$tjj%|| |j}|&| d j|dd d }|j\}}|| | | ||j#j$| |j#j$|| j#j$ tjj%|| t'|&|S)aPSolve the linear equations A x = b via Cholesky factorization of A, where A is a real symmetric or complex Hermitian positive-definite matrix. If matrix ``A`` is not positive definite, Cholesky factorization fails and it raises an error. Note: For batch input, NRHS > 1 is not currently supported. Args: a (cupy.ndarray): Array of real symmetric or complex hermitian matrices with dimension (..., N, N). b (cupy.ndarray): right-hand side (..., N) or (..., N, NRHS). Returns: x (cupy.ndarray): The solution (shape matches b). rrrr r rrrKTrMrrrO)(rrr"r(r)Z_assert_cupy_arrayZ_assert_stacked_2dZ_assert_stacked_squarerr]r$rQrZspotrfZspotrf_bufferSizeZspotrsZdpotrfZdpotrf_bufferSizeZdpotrsZcpotrfZcpotrf_bufferSizeZcpotrsZzpotrfZzpotrf_bufferSizeZzpotrsrrrRrr r!rrUr#r%r&r'r*rVZascontiguousarray)r+r,r-r)rZpotrfZpotrf_bufferSizeZpotrsrWrXr0r3rYrEZworksizerHrZr\r1r4r4r5posvsh         r^)numpyr$Zcupyr"rrrZ cupy.cudarr r6rIr]r^r4r4r4r5s  H ;