o i @spdZddlmZmZmZmZddlmZdgZ d dee de ded eeee dfd ee e ff d dZ dS) z0Determines if a contraction can use BLAS or not.)ListSequenceTupleUnion)ArrayIndexTypecan_blasNinputsresult idx_removedshapesreturnc Cst|dkrdS|\}}t||D]-}||||}}|dks-|dks-||dkr0dS||dt||vkr?dSq|dur^|D]}|d|||d||kr]dSqFt|dkrfdSdd|D} | d|} | d|} t|} |d|dkrd S| d| dkrd S|| d|d| krd S|d| || dkrd S|| d|| dkrd S|d| |d| krd St| dkst| dkrd Sd S)aChecks 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 FNrz OUTER/EINSUMcSsg|]}t|qS)set).0xrrE/mnt/sdb/aimis/docanh/lib/python3.10/site-packages/opt_einsum/blas.py Wszcan_blas..DOTz DOT/EINSUMGEMMz GEMV/EINSUMTDOT)lenrcountintfind) rr r r input_left input_rightcnlnrsets keep_left keep_rightrsrrrr sH /$   )N) __doc__typingrrrropt_einsum.typingr__all__strrboolrrrrrs