o װ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:{})getrf_bufferSizegetrfgetrsz'a must be F-contiguous or C-contiguous.zb must be F-contiguous.dtypeN)cupy_backends.cuda.libsrndim ValueErrorformatshaper TypeErrorgetattr _f_contiguous_cublas CUBLAS_OP_N _c_contiguous CUBLAS_OP_T_deviceget_cusolver_handle_cupyempty_numpyint32dataptrlinalg_util3_check_cusolver_dev_info_if_synchronization_allowed)ab _cusolverrthelperrrnnrhstranshandledipivdinfolworkdworkr;@/home/ubuntu/.local/lib/python3.10/site-packages/cupyx/lapack.pygesvsl       r=cCs\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: {})geqrf_bufferSizegeqrftrsmsdormqr_bufferSizeormqrunmqr_bufferSizeunmqrfdrg?)orderN)rrG)#rrrrrrrrrr charr" CUBLAS_OP_Cminr%r&r'r(r#r$get_cublas_handlearraycopyr)r*r+r,r-CUBLAS_SIDE_LEFTCUBLAS_FILL_MODE_UPPERCUBLAS_DIAG_NON_UNITctypesconjTzeros)r.r/r0r4rr1 geqrf_helperr?r@ ormqr_helperrCno_transr5mr3mn_mindev_infotaucusolver_handle cublas_handleonews_size workspacebb out_shaper;r;r<gelsPs                 rcc 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: {})CTrGrMrr)&rrcupyx.cusolvercheck_availability RuntimeErrorr' promote_typesr spotrfBatched spotrsBatched dpotrfBatched dpotrsBatched cpotrfBatched cpotrsBatched zpotrfBatched zpotrsBatchedrrastyper%_core _mat_ptrsrintprodr#r$rCUBLAS_FILL_MODE_LOWERr&r(r)r*r+r,r-rRreshape)r.r/r0cupyxr potrfBatchedrdmsgapldar3 batch_sizer6uplorZb_shapebpldbr4r;r;r< _batched_posvs^     rc 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 rrreTrgrrri)(rrr%r+r,_assert_cupy_array_assert_stacked_2d_assert_stacked_squarerrr'rmrspotrfspotrf_bufferSizespotrsdpotrfdpotrf_bufferSizedpotrscpotrfcpotrf_bufferSizecpotrszpotrfzpotrf_bufferSizezpotrsrrrvrr#r$rr{r&r(r)r*r-r|ascontiguousarray)r.r/r0r,rpotrfpotrf_bufferSizepotrsrrr3r6rrZworksizer`rrr4r;r;r<posvsh         r)numpyr'cupyr%rrr cupy.cudarr#r=rcrrr;r;r;r<s  H ;