0;ji>dZddlmZmZmZmZmZmZddlm Z ddl Z ddl Z ddl mZddlmZddlmZdd lmZed d 5Zed 5Ze jed Zdddn #1swxYwYdddn #1swxYwYeddddddededeeeefdededede jfdZedddddedede d de jfd Z eddd!dedede ddeeeeffd"Z edddddedededee jeeeefffd#Z edddd ddedededee jeeeefffd$Z d%Z!d&Z"eddddd'ededede d de jf d(Z#eddd!d'ededede ddeeeeff d)Z#eddddd'ededededee jeeeefff d*Z#edddd dd'ededededee jeeeefff d+Z#dS),z#Functions for interval construction) CollectionDictListUnionoverloadIterable)LiteralN) ArrayLike)cache) _FloatLike_co)_resource_filez librosa.corezintervals.msgpackrbF)strict_map_key )level T)bins_per_octavetuningsortn_binsfmin intervalsrrrreturnct|tr|dkr&d|tjd|tz|z z}n|dkrt ||}ns|dkrt dg|| }nY|d krt dd g|| }n>|d krt gd || }n#tj|}t|}tj ||z }tj dtj|z| d|}|rtj |}||zS)aConstruct a set of frequencies from an interval set Parameters ---------- n_bins : int The number of frequencies to generate fmin : float > 0 The minimum frequency intervals : str or array of floats in [1, 2) If `str`, must be one of the following: - `'equal'` - equal temperament - `'pythagorean'` - Pythagorean intervals - `'ji3'` - 3-limit just intonation - `'ji5'` - 5-limit just intonation - `'ji7'` - 7-limit just intonation Otherwise, an array of intervals in the range [1, 2) can be provided. bins_per_octave : int > 0 If `intervals` is a string specification, how many bins to generate per octave. If `intervals` is an array, then this parameter is ignored. tuning : float Deviation from A440 tuning in fractional bins. This is only used when `intervals == 'equal'` sort : bool Sort the intervals in ascending order. Returns ------- frequencies : array of float The frequencies Examples -------- Generate two octaves of Pythagorean intervals starting at 55Hz >>> librosa.interval_frequencies(24, fmin=55, intervals="pythagorean", bins_per_octave=12) array([ 55. , 58.733, 61.875, 66.075, 69.609, 74.334, 78.311, 82.5 , 88.099, 92.812, 99.112, 104.414, 110. , 117.466, 123.75 , 132.149, 139.219, 148.668, 156.621, 165. , 176.199, 185.625, 198.224, 208.828]) Generate two octaves of 5-limit intervals starting at 55Hz >>> librosa.interval_frequencies(24, fmin=55, intervals="ji5", bins_per_octave=12) array([ 55. , 58.667, 61.875, 66. , 68.75 , 73.333, 77.344, 82.5 , 88. , 91.667, 99. , 103.125, 110. , 117.333, 123.75 , 132. , 137.5 , 146.667, 154.687, 165. , 176. , 183.333, 198. , 206.25 ]) Generate three octaves using only three intervals >>> intervals = [1, 4/3, 3/2] >>> librosa.interval_frequencies(9, fmin=55, intervals=intervals) array([ 55. , 73.333, 82.5 , 110. , 146.667, 165. , 220. , 293.333, 330. ]) equalg@rdtype pythagoreanrrji3)primesrrji5ji7)r#r&N) isinstancestrnparangefloatpythagorean_intervalsplimit_intervalsarraylenceilmultiplyouterflattenr) rrrrrrratios n_octaves all_ratioss P/root/voice-cloning/.venv/lib/python3.11/site-packages/librosa/core/intervals.pyinterval_frequenciesr:sP)S!!&   ")AeDDDDWFF- ' '*?QUVVVFF %  %sO$FF%  %1vTFF%  % yy/F)$$f++011I""3")I*>*>#>GGOOQQJ )WZ((  .rrreturn_factorsr=cdSNr<s r9r.r.s Cr;r!cdSr?r@r<s r9r.r. Cr;cdSr?r@r<s r9r.r.rBr;ctj|tjtjdz\}|dk}||xxdz cc<|xxdz cc<t |rtj|}||}nt|}|rtfd|DStj d|S)aPythagorean intervals Intervals are constructed by stacking ratios of 3/2 (i.e., just perfect fifths) and folding down to a single octave:: 1, 3/2, 9/8, 27/16, 81/64, ... Note that this differs from 3-limit just intonation intervals in that Pythagorean intervals only use positive powers of 3 (ascending fifths) while 3-limit intervals use both positive and negative powers (descending fifths). Parameters ---------- bins_per_octave : int The number of intervals to generate sort : bool If `True` then intervals are returned in ascending order. If `False`, then intervals are returned in circle-of-fifths order. return_factors : bool If `True` then return a list of dictionaries encoding the prime factorization of each interval as `{2: p2, 3: p3}` (meaning `3**p3 * 2**p2`). If `False` (default), return intervals as an array of floating point numbers. Returns ------- intervals : np.ndarray or list of dictionaries The constructed interval set. All intervals are mapped to the range [1, 2). See Also -------- plimit_intervals Examples -------- Generate the first 12 intervals >>> librosa.pythagorean_intervals(bins_per_octave=12) array([1. , 1.067871, 1.125 , 1.201355, 1.265625, 1.351524, 1.423828, 1.5 , 1.601807, 1.6875 , 1.802032, 1.898437]) >>> # Compare to the 12-tone equal temperament intervals: >>> 2**(np.arange(12)/12) array([1. , 1.059463, 1.122462, 1.189207, 1.259921, 1.33484 , 1.414214, 1.498307, 1.587401, 1.681793, 1.781797, 1.887749]) Or the first 7, in circle-of-fifths order >>> librosa.pythagorean_intervals(bins_per_octave=7, sort=False) array([1. , 1.5 , 1.125 , 1.6875 , 1.265625, 1.898437, 1.423828]) Generate the first 7, in circle-of-fifths other and factored form >>> librosa.pythagorean_intervals(bins_per_octave=7, sort=False, return_factors=True) [ {2: 0, 3: 0}, {2: -1, 3: 1}, {2: -3, 3: 2}, {2: -4, 3: 3}, {2: -6, 3: 4}, {2: -7, 3: 5}, {2: -9, 3: 6} ] r#rc3<K|]}| |dVdS))r r#Nr@).0ipow2pow3s r9 z(pythagorean_intervals..s5;;!aT!W--;;;;;;r;r ) r+r,modflog2astypeintargsortrangelistpower)rrr= log_ratios too_smallidxrIrJs @@r9r.r.sL 9_ % %Dwtbgajj011JQIyQOOOqOOO ;;s  D %j$$_ O$$<;;;;;s;;;;;; 8Az " ""r;crtj|}tj|}tj|d}||z }tj|d}||z }tj||tj||z }tj|||zd|zz dS)uCompute the harmonic distance between ratios a and b. Harmonic distance is defined as `log2(a * b) - 2*log2(gcd(a, b))` [#]_. Here we are expressing a and b as prime factorization exponents, and the prime basis are provided in their log2 form. .. [#] Tenney, James. "On ‘Crystal Growth’ in harmonic space (1993–1998)." Contemporary Music Review 27.1 (2008): 47-56. rr )r+r0maximumminimumarounddot)logsaba_numa_denb_numb_dengcds r9__harmonic_distanceres  A  A Jq!  E AIE Jq!  E AIE *UE " "RZu%=%= =C 9TXXa!ea#go.. 2 22r;c|tj||tj|kS)z-Given two tuples of prime powers, break ties.)r\r+abs)r^r_r]s r9_crystal_tie_breakrh"s5 88BF1II  "&))!4!4 44r;r$cdSr?r@r$rrr=s r9r/r/' Cr;cdSr?r@rjs r9r/r/2rkr;cdSr?r@rjs r9r/r/=rkr;c Xtj|}tj|tj}g}t t |D]c}dgt |z}d||<|t|d||<|t|d|}t} t} tdgt |z} | | t | |kr.tj } d} t|D]z\}}d}| D]:}||f| vr%t|||| ||f<| ||f| ||f<|| ||fz };|| ks,tj|| rt||| |r|} |} {|| }| ||D]U}ttj|tj|z}|| vr||vr||Vt | |k.tjt| t$}tj||\}}|dk}||xxdz cc<||xxdzcc<|t,}|rtj|}||}nt |}|rwg}|D]p}t}||dkr || |d<|dt3|||D||q|Stjd|S)u Construct p-limit intervals for a given set of prime factors. This function is based on the "harmonic crystal growth" algorithm of [#1]_ [#2]_. .. [#1] Tenney, James. "On ‘Crystal Growth’ in harmonic space (1993–1998)." Contemporary Music Review 27.1 (2008): 47-56. .. [#2] Sabat, Marc, and James Tenney. "Three crystal growth algorithms in 23-limit constrained harmonic space." Contemporary Music Review 27, no. 1 (2008): 57-78. Parameters ---------- primes : array of odd primes Which prime factors are to be used bins_per_octave : int The number of intervals to construct sort : bool If `True` then intervals are returned in ascending order. If `False`, then intervals are returned in crystal growth order. return_factors : bool If `True` then return a list of dictionaries encoding the prime factorization of each interval as `{2: p2, 3: p3, ...}` (meaning `3**p3 * 2**p2`). If `False` (default), return intervals as an array of floating point numbers. Returns ------- intervals : np.ndarray or list of dictionaries The constructed interval set. All intervals are mapped to the range [1, 2). See Also -------- pythagorean_intervals Examples -------- Compare 3-limit tuning to Pythagorean tuning and 12-TET >>> librosa.plimit_intervals(primes=[3], bins_per_octave=12) array([1. , 1.05349794, 1.125 , 1.18518519, 1.265625 , 1.33333333, 1.40466392, 1.5 , 1.58024691, 1.6875 , 1.77777778, 1.8984375 ]) >>> # Pythagorean intervals: >>> librosa.pythagorean_intervals(bins_per_octave=12) array([1. , 1.06787109, 1.125 , 1.20135498, 1.265625 , 1.35152435, 1.42382812, 1.5 , 1.60180664, 1.6875 , 1.80203247, 1.8984375 ]) >>> # 12-TET intervals: >>> 2**(np.arange(12)/12) array([1. , 1.05946309, 1.12246205, 1.18920712, 1.25992105, 1.33483985, 1.41421356, 1.49830708, 1.58740105, 1.68179283, 1.78179744, 1.88774863]) Create a 7-bin, 5-limit interval set >>> librosa.plimit_intervals(primes=[3, 5], bins_per_octave=7) array([1. , 1.125 , 1.25 , 1.33333333, 1.5 , 1.66666667, 1.875 ]) The same example, but now in factored form >>> librosa.plimit_intervals(primes=[3, 5], bins_per_octave=7, ... return_factors=True) [ {}, {2: -3, 3: 2}, {2: -2, 5: 1}, {2: 2, 3: -1}, {2: -1, 3: 1}, {3: -1, 5: 1}, {2: -3, 3: 1, 5: 1} ] rrrErr c@i|]\}}|dk |t|S)r)rO)rGprSs r9 z$plimit_intervals..s*WWW5EUVJJaUJJJr;)r+ atleast_1drMfloat64rQr1appendtuplecopydictrRinf enumeratereiscloserhpopr0r-rLr\rNrOrPupdateziprS)r$rrr=r]seedsrHseedfrontier distancesrrootscorebest_ffpointHDs new_point_new_seedpowsrTrIrUrVfactorsvs r9r/r/Hsh]6 " "F 76 , , ,D E 3v;;  ""sS[[ Q U4[[!!!Q U4[[!!!!zz||HII !s6{{" # #D T i..? * *!(++  HAuB * *u:Y..*=dE1*M*MIah'*3AuH*=IeQh'i5))Ezz 2u%%&uhv.>EELL(( ### * *ARXi0028A;;>??Hy((XX-E-E))); i..? * *> 8DOO5 1 1 1Dwtxx~~..JQIyQOOOqOOO ;;s  D %j$$_ O$$   AAAw!||Qx! HHWWCQ4H4HWWW X X X NN1     8Az " ""r;)$__doc__typingrrrrrrtyping_extensionsr msgpacknumpyr+ numpy.typingr _cacher _typingr util.filesrimsgpackopen_fdescload INTERVALSrOr*r-boolndarrayr:r.rerhr/r@r;r9rs*)DDDDDDDDDDDDDDDD%%%%%%""""""######''''''^N$788?H t  ? GL>>> ?????????????????????????????? R iii i iS*U++, i  i  i iZiiiiX %(  EN  Z   !)-ELT] $sCx.   !C)-EI 2:tDcN+ +,  R tEd#d#d#d#(,d#EId# 2:tDcN+ +,d#d#d#d#N333@555  %(     EN  Z        DM   $sCx.           2:tDcN+ +,  R u#u#u# u#u#  u#  u#  2:tDcN+ +, u#u#u#u#u#u#s6BA:. B:A> >BA> BBB