Vh(RddlmZddlmZmZddlZddlmZddl m Z ddl m Zddlm cmZddlmZddl mZmZmZmZmZmZmZddlmZmZmZgdZdeej<d ej<d efd Z ddl!Z!dd l"m#Z#dd l$m%Z%e#ejLjNjPjRee%d1dedede*fdZ+dddddede*de*de*def dZ,e#ejLjNjPjZed d2dddedee.e*fdeed e/deej<def d!Z0d"Z1d#Z2ed d3ddd$edee.efded e/deej<def d%Z3ed d4ddd$edeee.efdeed e/deej<def d&Z4ed'd(d)dd5d$ed*e/de5eeeffd+Z6edd$edefd,Z7eed-ejp.d1ded/ede*defd0Z9y)6)partial)OptionalUnionN)Tensor)check_fp_or_complexcheck_is_matrixDimDimsTypeELEMENTWISE_TYPE_PROMOTION_KINDIntLikeTensorLikeType)_maybe_convert_to_dtype"elementwise_type_promotion_wrapper out_wrapper)diagonal matrix_normnormsvdsvdvals vector_normvecdotcrossdtypex_dtypefn_namectjtjxstjfdtjtjtjk(fdtjtj k(fdyy)zG Checks related to the dtype kwarg in `linalg.*norm` functions NcdS)Nz1: dtype should be floating point or complex. Got )rrsK/home/dcms/DCMS/lib/python3.12/site-packages/torch/_refs/linalg/__init__.pyz#_check_norm_dtype...swiPQVPWXcfdjtjrdSdS)Nz:{fn_name}: dtype should be {d} for {d} inputs. Got {dtype}complexreal)rdr)formatutilsis_complex_dtype)rrrsrr z#_check_norm_dtype..2sBPWW$55g>)XDJXr!cddS)Nz: the dtype of the input (zJ) should be convertible without narrowing to the specified dtype ({dtype})r)rrsrr z#_check_norm_dtype..:swi9'CAAr!)torch_checkr'is_float_dtyper(get_higher_dtype)rrrs```r_check_norm_dtyper.'s    ' H5+A+A%+H X    " "5 )U-C-CG-L L     " "5' 2e ; A r!)register_decomposition)pw_cast_for_opmathabdimcdtjjjk(dtjjdk(xrjdk(fdtj\t j jtjdj}j|dzdzj|dzdzzj|dzdzj|dzdzzz S)Ncy)Nz=linalg.cross: inputs must have the same number of dimensions.rrr!rr zcross..Lr!cVddjdjS)Nzlinalg.cross: inputs dim z must have length 3, got z and )size)r2r3r4srr zcross..Ps4+C50I!&&QT+V[\]\b\bcf\g[hir!)device) r*r+ndimr:broadcast_tensorsr'canonicalize_dimaranger; index_select)r2r3r4idxs``` rrrFs LL !&&O LL s q-QVVC[A-i  " "1a (DAq   -C ,,q *C >>#a1} - cAg]1  sS1WM*Q^^C#'Q-OO PPr!offsetdim1dim2inputrFrGrHreturnc4tj||||S)NrE)r*r)rIrFrGrHs rrrZs >>%T EEr!T) exact_dtyperxordkeepdimc  ddlm}t|jdt |t r|g}||j dk(rudkstdk(rbtj|duxrt|dk7fd|j}|J|D]" tj| dk7 fd$t||jdtj|tjj |\}}t#t$|} dk(r-tj&tj(|d||| Stdk(r1| tj*tj,||| Std k(r1| tj.tj,||| St%||}t#tj&|| } t t0rd zdk(nd zdk(} | rtj2|jstj,|}| tj4| tj4|dz S)Nr)guard_size_obliviouszlinalg.vector_norminfcddS)N&linalg.vector_norm cannot compute the zH norm on an empty tensor because the operation does not have an identityr)rOsrr zvector_norm..ys>r!cdddS)NrVz norm on the dimension zL because this dimension is empty and the operation does not have an identityr)r%rOsrr zvector_norm..s!@FC 66r!rMr4rPr)r4rPz-infr=@?)%torch.fx.experimental.symbolic_shapesrRrr isinstancer numelfloatr*r+lenshaper.r'reduction_dtypesREDUCTION_OUTPUT_TYPE_KINDCOMPLEX_TO_FLOATrrsumneamaxabsaminr r,pow) rNrOr4rPrrRr`computation_dtype result_dtypeto_result_dtype reduce_sum is_ord_evenr%s ` @rrrds K!56#seAGGIN+scU5\>Q tO -CA  >  A LLaA 6  eQWW&:;&+&<&< 5 + + < cAglC#IQTDT  4 4QWW = ! AuyyEIIa4E)Fc RSSr!c~t|Dcgc]}||k7s ||k7s|}}|j||f|Scc}w)N)rangeextend)dim0rGr>irets r_backshift_permutationrusADk =Q$Y191 =C =JJd| J >s :::ctt|tjdDcgc]\}}| c}}Scc}}w)Nr<)key)sorted enumerateoperator itemgetter)permrsjs r_inverse_permutationr~s0 4h6I6I!6LM N$!QA NN Ns =Ac t|dtj|j|}t |t r|f}t jt|dk(dt j|d|dk7dt||jdt |trt j|dvdt|jd|d k7 |d k(rt|d||| S| t||}t|d|d|j}t j t#t%j&||d |}|r5t)|}t%j&t j*|d |}|St-|}t j|ddt/dfvdt|jd|dk7 t1|dkDrt j2nt j4|} |dk(r| t||}t|d|d|j}| t#t%j&||d }|r5t)|}t%j&t j*|d |}|S|\} } |t/dk(r| | } } |s | | kr| dz} | t|d| ||| S)Nzlinalg.matrix_normr=cy)Nz4linalg.matrix_norm: dim must be a 2-tuple. Got {dim}rrr!rr zmatrix_norm..r7r!rr<cy)NzDlinalg.matrix_norm: dims must be different. Got ({dim[0]}, {dim[1]})rrr!rr zmatrix_norm..r7r!)fronuccyNz.linalg.matrix_norm: Order {ord} not supported.rrr!rr zmatrix_norm..r7r!r)allow_low_precision_dtypesrrMr1rTcyrrrr!rr zmatrix_norm..r7r!rS)rPrYr4rZrX)rr'canonicalize_dimsr>r\r r*r+r_r.rstrrrrrurdrprims transposer~ unsqueezergr^rrfrh) rrOr4rPrr|resultinv_permabs_ordmax_minrrrGs rrrsA+,  ! !!&&# .C#sf LL CA U LL A#a&V eQWW&:;#s > ! D  GG)cUl  %<q!S'? ? +Au5)#a&#a&!&&ADYYwuq$'?@"gNF/5)DhOMc( 1eEl+ + D  GG)cQh c %**uzz7S c> +Au5)#a&#a&!&&ADWU__Q%=>BGF/5)DhOMJD$%,&!4dt  AsgUKT r!cR|6t|tr|f}tjt |dvdn%|#tj|j dvd|5|t |dk(s|#|j dk(r|d}t |||||S|d}t|||||S)N)r<r=cy)NzHlinalg.norm: If dim is specified, it must be of length 1 or 2. Got {dim}rrr!rr znorm.. r7r!cy)NzZlinalg.norm: If dim is not specified but ord is, the input must be 1D or 2D. Got {A.ndim}Drrr!rr znorm..r7r!r=)rr<rMrY)r\r r*r+r_r>rr)rrOr4rPrs rrrs  c3 &C H  ^   FFf  p   SX] ! ;C1c3u== ;C1c3u==r!USVh full_matricesc0tj||S)Nr)rr)rrs rrr"s 99Qm 44r!c"t|ddS)NFrr<)r)rs rrr(s q &q ))r!)rNy)type_promoting_argstype_promotion_kindrctt|jd|j|zj|S)Nz linalg.vecdotr)rrconjrd)rNrr4s rrr.s0 1 FFHqL  #  &&r!)r1)r=NF)r)rDr1F)NNF)T): functoolsrtypingrrr* torch._prims_primsrtorch._prims_common _prims_commonr' torch._refs_refsrefstorch._refs.linalglinalgrrrr r r r r torch._prims_common.wrappersrrr__all__rrr.rz torch._decompr/torch._decomp.decompositionsr0_opsopsaten linalg_crossintrrlinalg_vector_normr^boolrrur~rrtuplerrDEFAULTrrr!rrs" ###  Xekk2 U[[ SV 01; ++889 PVPPSP:P( F F F  F  F  F ++>>? " 6T $( 6T6T ucz 6T ( 6T 6T EKK 6T 6T@6TrO  # M $( MM ucz M M M EKK MMMb (," > $( >> %s # $> ( > > EKK >>>D S#t.5>5$5%PV@V:W5/5  *~*&**  #"7??'f''c'6'  'r!