-i_ xSr/rSSKrSSKJrJr "SS\RR5r SSK J r S Sjr g) zz Matrix square root for general matrices and for upper triangular matrices. This module exists to avoid cyclic imports. N)ztrsyldtrsylc\rSrSrSrg) SqrtmErrorN)__name__ __module__ __qualname____firstlineno____static_attributes__r O/var/www/html/venv/lib/python3.13/site-packages/scipy/linalg/_matfuncs_sqrtm.pyrrsrr)within_block_loopc [R"U5n[R"U5=(a [R"USS9S:nU(dH[R"U[R SS9n[R"U[R S9nOG[R"U[R SS9n[R"U[R S9n[R"[R"U55nURu n[XQ-S5n[XV5upxUS-n Xh- n X-X--U:wa [S5e/n Sn X4X44H0up[U 5HnU RXU-45 X- n M M2 [X@X5 [U5HnU Uunn[US- S S 5HnXunnUUU2UU24nUU- S:a&UUUU2UU24R%UUU2UU245- nUUU2UU24nUUU2UU24nU(a['UUU5unnnO[)UUU5unnnUU-UUU2UU24'M M U$![an[!UR"6UeS nAff=f) a Matrix square root of an upper triangular matrix. This is a helper function for `sqrtm` and `logm`. Parameters ---------- T : (N, N) array_like upper triangular Matrix whose square root to evaluate blocksize : int, optional If the blocksize is not degenerate with respect to the size of the input array, then use a blocked algorithm. (Default: 64) Returns ------- sqrtm : (N, N) ndarray Value of the sqrt function at `T` References ---------- .. [1] Edvin Deadman, Nicholas J. Higham, Rui Ralha (2013) "Blocked Schur Algorithms for Computing the Matrix Square Root, Lecture Notes in Computer Science, 7782. pp. 171-182. g)initialrC)dtypeorder)rrzinternal inconsistencyN)npdiag isrealobjminasarray complex128float64sqrtshapemaxdivmod Exceptionrangeappendr RuntimeErrorrargsdotrr)T blocksizeT_diag keep_it_realRnnblocksbsmallnlargeblargensmallstart_stop_pairsstartcountsizeiejjstartjstopistartistopSRiiRjjxscaleinfos r _sqrtm_triurEs4WWQZF<<?Frvvfb'AQ'FL  JJq S 9F"--8 JJq # 6F"**5  A 77DAq!.!$GA'NF aZF  F (A-011 E(6*:; uA  # #UDL$9 : ME< )! 0: 7^(+ qsB#A,/MFE&,u ,-A1uqy&,f 4599!E&LrNs: # && 4W r