}Kgz>RdZddlmZddlmZmZmZmZmZm Z ddl m Z ddl Z ddl ZddlmZddlmZdd lmZej*d d 5Zej/d 5Ze j2ed Zddddddeddddddededeeeefdedededej>fdZ edddddedede d dej>fd Z!eddd!dedede ddeeeeffd"Z!edddddedededeej>eeeefffd#Z!edddd ddedededeej>eeeefffd$Z!d%Z"d&Z#eddddd'ededede d dej>f d(Z$eddd!d'ededede ddeeeeff d)Z$eddddd'ededededeej>eeeefff d*Z$edddd dd'ededededeej>eeeefff d+Z$y#1swYxYw#1swYxYw),z#Functions for interval construction) resources) CollectionDictListUnionoverloadIterable)LiteralN) ArrayLike)cache) _FloatLike_coz librosa.corezintervals.msgpackrbF)strict_map_key )level T)bins_per_octavetuningsortn_binsfmin intervalsrrrreturncVt|tr|dk(r&d|tjd|tz|z z}nt|dk(rt ||}na|dk(rt dg|| }nL|d k(rt dd g|| }n6|d k(r1t gd || }n tj|}t|}tj||z }tjjdtj|zjd|}|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 Z/home/alanp/www/video.onchill/myenv/lib/python3.12/site-packages/librosa/core/intervals.pyinterval_frequenciesr:s*P)S!  "))AeDDWF- '*?QUVF % %sO$F% %1vTF% % /F)$f+01I""3"))I*>#>GOOQJ WWZ(  .rrreturn_factorsr=cyNr<s r9r.r.sr;r!cyr?r@r<s r9r.r.r;cyr?r@r<s r9r.r.rBr;ctj|tjtjdz\}|dk}||xxdz cc<|xxdz cc<j t |rtj |}||}n t|}|rtfd|DStjd|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#rc36K|]}| |dyw))r r#Nr@).0ipow2pow3s r9 z(pythagorean_intervals..s!;s!aT!W-ssr ) r+r,modflog2astypeintargsortrangelistpower)rrr= log_ratios too_smallidxrIrJs @@r9r.r.sL 99_ %Dwwtbggaj01JQIyQOqO ;;s D jj$_ O$;s;;; 88Az ""r;cxtj|}tj|}tj|d}||z }tj|d}||z }tj||tj||z }tj|j ||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 JJq! E AIE JJq! E AIE **UE "RZZu%= =C 99TXXa!ea#go. 22r;c|jtj||jtj|kS)z-Given two tuples of prime powers, break ties.)r\r+abs)r^r_r]s r9_crystal_tie_breakrh"s/ 88BFF1I "&&)!4 44r;r$cyr?r@r$rrr=s r9r/r/'r;cyr?r@rjs r9r/r/2rkr;cyr?r@rjs r9r/r/=rkr;c tj|}tj|tj}g}t t |D]O}dgt |z}d||<|j t|d||<|j t|Q|j}t} t} tdgt |z} | j | t | |krtj} d} t|D]m\}}d}| D]0}||f| vrt|||| ||f<| ||f| ||f<|| ||fz }2|| ks(tj|| sYt||| |sj|} |} o|j!| }| j ||D]Q}ttj"|tj"|z}|| vs<||vsA|j |St | |krtj"t| t$}tj&|j)|\}}|dk}||xxdz cc<||xxdzcc<|j+t,}|rtj.|}||}n t |}|rxg}|D]o}t}||dk7r || |d<|j1t3|||Dcic]\}}|dk7s |t-|c}}|j |q|Stj4d|Scc}}w)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 )r+ atleast_1drMfloat64rQr1appendtuplecopydictrRinf enumeratereiscloserhpopr0r-rLr\rNrOrPupdateziprS)r$rrr=r]seedsrHseedfrontier distancesrrootscorebest_ffpointHDs new_point_new_seedpowsrTrIrUrVfactorsvprSs r9r/r/Hs!h]]6 "F 776 ,D E 3v; sS[ Q U4[!Q U4[! zz|HII !s6{" #D T i.? *!(+HAuBu:Y.*=dE1*MIah'*3AuH*=IeQh'i5)) Ez 2u%&uhv.>E!,$LL( #ARXXi0288A;>?Hy(XX-E)5 i.? *> 88DO5 1Dwwtxx~.JQIyQOqO ;;s D jj$_ O$AAAw!|Qx! HHCQ4HW4H5EUVJaUm4HW X NN1  88Az "" Xs < M M )%__doc__ importlibrtypingrrrrrr typing_extensionsr msgpacknumpyr+ numpy.typingr _cacher _typingrpathimsgpackopen_fdescload INTERVALSrOr*r-boolndarrayr:r.rerhr/r@r;r9rs*DD%"#Y^^N$78H t  GLL>  9 R i i iS*U++, i  i  i iZZiiX %(  EN  ZZ   !)-ELT] $sCx.   !C)-EI 2::tDcN+ +,  R tEd#d#(,d#EId# 2::tDcN+ +,d#d#N3@5  %(     EN  ZZ        DM   $sCx.           2::tDcN+ +,  R u# u#u#  u#  u#  2::tDcN+ +, u#u#s  98s$ HH3HH HH&