o °2úh¬ ã @spdZddlmZmZmZmZddlmZdgZ d dee de ded eeee dfd ee e ff d d„Z dS) z0Determines if a contraction can use BLAS or not.é)ÚListÚSequenceÚTupleÚUnion)ÚArrayIndexTypeÚcan_blasNÚinputsÚresultÚ idx_removedÚshapesÚreturnc Csºt|ƒdkrdS|\}}t||ƒD]-}| |¡| |¡}}|dks-|dks-||dkr0dS||dt||vƒkr?dSq|dur^|D]}|d| |¡|d| |¡kr]dSqFt|ƒdkrfdSdd„|Dƒ} | d|} | d|} t|ƒ} |d|dkr‡d S| d| dkr‘d S|| d…|d| …kr d S|d| …|| d…kr¯d S|| d…|| d…kr¿d S|d| …|d| …krÍd St| ƒdksÙt| ƒdkrÛd Sd S)aˆChecks if we can use a BLAS call. Parameters ---------- inputs : list of str Specifies the subscripts for summation. result : str Resulting summation. idx_removed : set Indices that are removed in the summation shapes : sequence of tuple[int], optional If given, check also that none of the indices are broadcast dimensions. Returns: ------- type : str or bool The type of BLAS call to be used or False if none. Notes: ----- We assume several operations are not efficient such as a transposed DDOT, therefore 'ijk,jki->' should prefer einsum. These return the blas type appended with "/EINSUM" to differentiate when they can still be done with tensordot if required, e.g. when a backend has no einsum. Examples: -------- >>> can_blas(['ij', 'jk'], 'ik', set('j')) 'GEMM' >>> can_blas(['ijj', 'jk'], 'ik', set('j')) False >>> can_blas(['ab', 'cd'], 'abcd', set()) 'OUTER/EINSUM' >>> # looks like GEMM but actually 'j' is broadcast: >>> can_blas(['ij', 'jk'], 'ik', set('j'), shapes=[(4, 1), (5, 6)]) False éFéNrz OUTER/EINSUMcSsg|]}t|ƒ‘qS©)Úset)Ú.0ÚxrrúJ/var/www/html/chatgem/venv/lib/python3.10/site-packages/opt_einsum/blas.pyÚ Wszcan_blas..ÚDOTz DOT/EINSUMÚGEMMz GEMV/EINSUMÚTDOT)ÚlenrÚcountÚintÚfind) rr r r Ú input_leftÚ input_rightÚcÚnlÚnrÚsetsÚ keep_leftÚ keep_rightÚrsrrrr sH /ÿ$ÿ   )N) Ú__doc__ÚtypingrrrrÚopt_einsum.typingrÚ__all__ÚstrrÚboolrrrrrÚs  üÿþýü û