Tenney–Euclidean metrics: Difference between revisions

From Xenharmonic Wiki
Jump to navigation Jump to search
Wikispaces>genewardsmith
**Imported revision 209349086 - Original comment: **
TE norm: clarify with a link
 
(50 intermediate revisions by 12 users not shown)
Line 1: Line 1:
<h2>IMPORTED REVISION FROM WIKISPACES</h2>
The '''Tenney-Euclidean metrics''' are {{w|metric (mathematics)|metrics}} defined in Tenney-Euclidean space. These consist of the TE norm, which measures the [[complexity]] of an [[interval]] in [[just intonation]], the TE temperamental norm, which measures the complexity of an interval ''as mapped by a temperament'', and the octave-equivalent TE seminorms of both.  
This is an imported revision from Wikispaces. The revision metadata is included below for reference:<br>
: This revision was by author [[User:genewardsmith|genewardsmith]] and made on <tt>2011-03-10 14:58:49 UTC</tt>.<br>
: The original revision id was <tt>209349086</tt>.<br>
: The revision comment was: <tt></tt><br>
The revision contents are below, presented both in the original Wikispaces Wikitext format, and in HTML exactly as Wikispaces rendered it.<br>
<h4>Original Wikitext content:</h4>
<div style="width:100%; max-height:400pt; overflow:auto; background-color:#f8f9fa; border: 1px solid #eaecf0; padding:0em"><pre style="margin:0px;border:none;background:none;word-wrap:break-word;white-space: pre-wrap ! important" class="old-revision-html">[[toc|flat]]
=The weighting matrix=
Let us define the val weighting matrix W to be the diagonal matrix with values 1, 1/log2(3), 1/log2(5) ... 1/log2(p) along the diagonal. Given a val "a" expressed as a row vector, the corresponding vector in weighted coordinates is aW, with transpose Wa* where the * denotes the transpose. Then the dot product of weighted vals is aW^2a*, which makes the Euclidean metric on vals, a measure of complexity, to be || &lt;a2 a3 ... ap| || = sqrt(a2^2 + a3^2/log2(3)^2 + ... + ap^2/log2(p)^2). Similarly, if b is a monzo, then in weighted coordinates the monzo becomes bW^(-1), and the dot product is bW^(-2)b*, leading to sqrt(b2^2 + log2(3)^2b3^2 + ... + log2(p)^2bp^2) as the norm on monzos, a measure of complexity we may call the Tenney-Euclidean, or TE, complexity.


=Temperamental complexity=
== TE norm ==
Suppose now A is a matrix whose rows are vals defining a p-limit regular temperament. Then the corresponding weighted matrix is V = AW. The [[Tenney-Euclidean tuning|TE]] tuning projection matrix is then V`V, where V` is the [[Tenney-Euclidean Tuning|pseudoinverse]]. If the rows of V (or equivalently, A) are linearly independent, then we have V` = V*(VV*)^(-1), where V* denotes the transpose. In terms of vals, the tuning projection matrix is P = V`V = V*(VV*)^(-1)V = WA*(AW^2A*)^(-1)AW. P is a [[http://en.wikipedia.org/wiki/Positive-definite_matrix|positive semidefinite matrix]], so it defines a [[http://en.wikipedia.org/wiki/Definite_bilinear_form|positive semidefinite bilinear form]]. In terms of weighted monzos m1 and m2, m1Pm2* defines the semidefinite form on weighted monzos, and hence b1W^(-1)PW^(-1)b2* defines a semidefinite form on unweighted monzos, in terms of the matrix W^(-1)WA*(AW^2A*)^(-1)AWW^(-1) = A*(AW^2A*)^(-1)A = **P**. From the semidefinite form we obtain an associated [[http://en.wikipedia.org/wiki/Definite_quadratic_form|semidefinite quadratic form]] b**P**b* and from this the [[http://en.wikipedia.org/wiki/Norm_%28mathematics%29|seminorm]] sqrt(b**P**b*).  
The '''Tenney–Euclidean norm''' ('''TE norm''') or '''Tenney–Euclidean complexity''' ('''TE complexity''') applies to [[val]]s as well as to [[monzo]]s.  


It may be noted that (VV*)^(-1) = (AW^2A*)^(-1) is the inverse of the [[http://en.wikipedia.org/wiki/Gramian_matrix|Gram matrix]] used to compute [[Tenney-Euclidean temperament measures|TE complexity]], and hence is the corresponding Gram matrix for the dual space. Hence **P** represents a change of basis defined by the mapping given by the vals combined with an [[http://en.wikipedia.org/wiki/Inner_product_space|inner product]] on the result. Given a monzo b, bA* represents the tempered interval corresponding to b in a basis defined by the mapping A, and //P// = (AW^2A*)^(-1) defines a positive-definite quadratic form, and hence a norm, on the tempered interval space with basis defined by A.
Let us define the val weighting matrix ''W'' to be the {{w|diagonal matrix}} with values 1, 1/log<sub>2</sub>3, 1/log<sub>2</sub>5 … 1/log<sub>2</sub>''p'' along the diagonal. For the [[harmonic limit|''p''-limit]] prime basis ''Q'' = {{val| 2 3 5 … ''p'' }},  


Denoting the temperament-defined, or temperamental, seminorm by T(x), the subspace of interval space such that T(x) = 0 contains a sublattice of the lattice of monzos consisting of the commas of the temperament. The [[http://en.wikipedia.org/wiki/Quotient_space_%28linear_algebra%29|quotient space]] of the full vector space by the commatic subspace such that T(x) = 0 is now a [[http://en.wikipedia.org/wiki/Normed_vector_space|normed vector space]] with norm given by T, in which the intervals of the regular temperament define a lattice. The norm T on these lattice points is the //temperamental norm// or //temperamental complexity// of the intervals of the regular temperament; in terms of the basis defined by A, it is sqrt(t//P//t*) where t is the image of a monzo b by t = bA*
$$ W = \operatorname {diag} (1/\log_2 (Q)) $$


=The OETES=
Right-multiplying a row vector by this matrix scales each entry by the corresponding entry of the diagonal matrix.  
Instead of starting from a matrix of vals, we may start from a matrix of monzos. If B is a matrix with rows of monzos spanning the commas of a regular temperament, then M = BW^(-1) is the corresponding weighted matrix. Q = M`M is a projection matrix dual to P = I-Q, where I is the identity matrix, and P is the same symmetric matrix as in the previous section. If the rows define a basis for the commas of the temperament, and are therefor linearly independent, then P = I - M*(MM*)^(-1)M = I - W^(-1)B*(BW^(-2)B*)^(-1)BW^(-1), and mPm* = bW^(-1)PW^(-1)b*, or b(W^(-2) - W^(-2)B*(BW^(-2)B*)^(-1)BW^(-2))b*, so that the terms inside the parenthesis define a formula for **P** in terms of the matrix of monzos B.


To define the OETES, or Tenney-Euclidean octave equivalent seminorm, we simply add a row |1 0 0 ... 0&gt; representing 2 to the matrix B. An alternative proceedure is to find the [[normal lists|normal val list]], and remove the first val from the list, corresponding to the octave or some fraction thereof, and proceed as in the previous section on temperamental complexity. This seminorm is a measure of the octave-equivalent complexity of a given p-limit rational interval in terms of thenp-limit regular temperament given by A.
Given a val ''V'' expressed as a row vector, the corresponding row vector in weighted coordinates is {{nowrap| ''V''<sub>''W''</sub> {{=}} ''VW'' }}, with transpose {{nowrap| {{subsup|''V''|''W''|T}} {{=}} ''WV''{{t}} }} where {{t}} denotes the transpose. The {{w|dot product}} of a weighted val with itself, or the sum of the squares of its entries, is the squared Euclidean metric of the val, {{nowrap| {{subsup|‖''V''<sub>''W''</sub>‖|2|2}} {{=}} ''V''<sub>''W''</sub>{{subsup|''V''|''W''|T}} {{=}} ''VW''<sup>2</sup>''V''{{t}} }}. Thus the Euclidean metric on the val, a measure of complexity, is {{nowrap| ‖''V''<sub>''W''</sub>‖<sub>2</sub> {{=}} sqrt(''V''<sub>''W''</sub>{{subsup|''V''|''W''|T}}) }} {{nowrap| {{=}} sqrt({{subsup|''v''|1|2}} + {{subsup|''v''|2|2}}/(log<sub>2</sub>3)<sup>2</sup> + … + {{subsup|''v''|''n''|2}}/(log<sub>2</sub>''p'')<sup>2</sup>) }}, where {{nowrap|''n'' {{=}} π(''p'')}} is the {{w|prime-counting function}} which records the number of primes to ''p''; dividing this by sqrt(''n'') gives the TE norm of a val.  


=Logflat TE badness=
Similarly, if '''m''' is a monzo, then in weighted coordinates the monzo becomes {{nowrap| '''m'''<sub>''W''</sub> {{=}} ''W''{{inv}}'''m''' }}, and the dot product is {{nowrap| {{subsup|'''m'''|''W''|T}}'''m'''<sub>''W''</sub> {{=}} '''m'''{{t}}''W''<sup>-2</sup>'''m''' }}, leading to {{nowrap| sqrt({{subsup|'''m'''|''W''|T}}'''m'''<sub>''W''</sub>) {{=}} sqrt({{subsup|''m''|1|2}} + (log<sub>2</sub>3)<sup>2</sup>{{subsup|''m''|2|2}} + … + (log<sub>2</sub>''p'')<sup>2</sup>{{subsup|''m''|''n''|2}}) }}; multiplying this by sqrt(''n'') gives the dual RMS norm on monzos which serves as a measure of complexity.
Given a matrix A whose rows are linearly independent vals defining a regular temperament, then the rank r of the temperament is the number of rows, which equals the number of linearly independent vals. The dimension of the temperament is the number of primes it covers; if p is the largest such prime, then the dimension d is pi(p), the number of primes to p. If we define S(A) to be the simple badness (relative error) of A, and C(A) to be the complexity of A, then //logflat badness// is defined by the formula S(A) C(A)^(r/(d-r)). If we set a cutoff margin for logflat badness, there are still infinite numbers of new temperaments appearing as complexity goes up, at a lower rate which is approximately logarithmic in terms of complexity.


==Examples==
== TE temperamental norm ==
Consider the temperament defined by the 5-limit [[Patent val|patent vals]] for 15 and 22 equal. From the vals, we may contruct a 2x3 matrix A = [&lt;15 24 35|, &lt;22 35 51|]. From this we may obtain the matrix **P** as A*(AW^2A*)^(-1)A, approximately
Suppose now ''V'' is a matrix whose rows are vals defining a ''p''-limit regular temperament. Then the corresponding weighted matrix is {{nowrap| ''V''<sub>''W''</sub> {{=}} ''VW'' }}. The [[Tenney–Euclidean tuning|TE tuning]] [[projection matrix]] is then {{nowrap| ''P''<sub>''W''</sub> {{=}} {{subsup|''V''|''W''|+}}''V''<sub>''W''</sub> }}, where {{+}} denotes the {{w|Moore–Penrose pseudoinverse}}. If the rows of ''V''<sub>''W''</sub> (or equivalently, ''V'') are linearly independent, then we have {{nowrap| {{subsup|''V''|''W''|+}} {{=}} {{subsup|''V''|''W''|T}}(''V''<sub>''W''</sub>{{subsup|''V''|''W''|T}}){{inv}} }}. In terms of vals, the tuning projection matrix is {{nowrap| {{subsup|''V''|''W''|+}}''V''<sub>''W''</sub> {{=}} {{subsup|''V''|''W''|T}}(''V''<sub>''W''</sub>{{subsup|''V''|''W''|T}}){{inv}}''V''<sub>''W''</sub> }} {{nowrap| {{=}} ''WV''{{t}}(''VW''<sup>2</sup>''V''{{t}}){{inv}}''VW'' }}. ''P''<sub>''W''</sub> is a {{w|positive-definite matrix|positive semidefinite matrix}}, so it defines a {{w|definite bilinear form|positive semidefinite bilinear form}}. In terms of weighted monzos ('''m'''<sub>''W''</sub>)<sub>1</sub> and ('''m'''<sub>''W''</sub>)<sub>2</sub>, {{subsup|('''m'''<sub>''W''</sub>)|1|T}}''P''<sub>''W''</sub>('''m'''<sub>''W''</sub>)<sub>2</sub> defines the semidefinite form on weighted monzos, and hence {{subsup|'''m'''|1|T}}''W''{{inv}}''P''<sub>''W''</sub>''W''{{inv}}'''m'''<sub>2</sub> defines a semidefinite form on unweighted monzos, in terms of the matrix {{nowrap| ''P'' {{=}} ''W''{{inv}}''P''<sub>''W''</sub>''W''{{inv}} }} {{nowrap| {{=}} ''V''{{t}}(''VW''<sup>2</sup>''V''{{t}}){{inv}}''V''}}. From the semidefinite form we obtain an associated {{w|definite quadratic form|semidefinite quadratic form}} '''m'''{{t}}''P'''''m''' and from this the {{w|norm (mathematics)|seminorm}} sqrt('''m'''{{t}}''P'''''m''').


[0.9911 0.1118 -0.1440]
It may be noted that {{nowrap|(''V''<sub>''W''</sub>{{subsup|''V''|''W''|T}}){{inv}} {{=}} (''VW''<sup>2</sup>''V''{{t}}){{inv}}}} is the inverse of the {{w|Gramian matrix}} used to compute [[TE complexity]], and hence is the corresponding Gram matrix for the dual space. Hence ''P'' represents a change of basis defined by the mapping given by the vals combined with an {{w|inner product space|inner product}} on the result. Given a monzo '''m''', ''V'''''m''' represents the tempered interval corresponding to '''m''' in a basis defined by the mapping ''V'', and {{nowrap|''P''<sub>''T''</sub> {{=}} (''VW''<sup>2</sup>''V''{{t}}){{inv}}}} defines a positive-definite quadratic form, and hence a norm, on the tempered interval space with basis defined by ''V''.
[0.1118 1.1075 1.8086]
[-0.1440 1.8086 3.0624]


If we want to find the temperamental seminorm T(250/243) of 250/243, we convert it into a monzo as |1 -5 3&gt;. Now we may multiply **P** by this on the left, obtaining the zero vector. Taking the dot product of the zero vector with |1 -5 3&gt; gives zero, and taking the square root of zero we get zero, the temperametal seminorm T(250/243) of 250/243. All of this is telling us that 250/243 is a comma of this temperament, which is 5-limit [[Porcupine family|porcupine]].
Denoting the temperament-defined, or temperamental, seminorm by ''T''(''x''), the subspace of interval space such that {{nowrap|''T''(''x'') {{=}} 0}} contains a lattice consisting of the commas of the temperament, which is a sublattice of the lattice of monzos. The {{w|quotient space (linear algebra)|quotient space}} of the full vector space by the commatic subspace such that {{nowrap|''T''(''x'') {{=}} 0}} is now a {{w|normed vector space}} with norm given by ''T'', in which the intervals of the regular temperament define a lattice. The norm ''T'' on these lattice points is the '''TE temperamental norm''' or '''TE temperamental complexity''' of the intervals of the regular temperament; in terms of the basis defined by ''V'', it is sqrt('''t'''{{t}}''P''<sub>''T''</sub>'''t''') where '''t''' is the image of a monzo '''m''' by {{nowrap| '''t''' {{=}} ''V'''''m''' }}.


Similarly, starting from the monzo |-1 1 0&gt; for 3/2, we may multiply this by **P**, obtaining &lt;-0.8793 0.9957 1.9526|, and taking the dot product of this with |-1 1 0&gt; gives 1.875 with square root 1.3693, which is T(3/2).
== Octave-equivalent TE seminorm ==
Instead of starting from a matrix of vals, we may start from a matrix of monzos. If ''M'' is a matrix with columns of monzos spanning the commas of a regular temperament, then {{nowrap| ''M''<sub>''W''</sub> {{=}} ''W''{{inv}}''M''}} is the corresponding weighted matrix. {{nowrap| ''Q''<sub>''W''</sub> {{=}} ''M''<sub>''W''</sub>{{subsup|''M''|''W''|+}} }} is a projection matrix dual to {{nowrap| ''P''<sub>''W''</sub> {{=}} ''I'' − ''Q''<sub>''W''</sub> }}, where ''I'' is the identity matrix, and ''P''<sub>''W''</sub> is the same symmetric matrix as in the previous section. If the rows define a basis for the commas of the temperament, and are therefore linearly independent, then {{nowrap| ''P''<sub>''W''</sub> {{=}} ''I'' − ''M''<sub>''W''</sub>({{subsup|''M''|''W''|T}}''M''<sub>''W''</sub>){{inv}}{{subsup|''M''|''W''|T}} }} {{nowrap| {{=}} ''I'' − ''W''{{inv}}''M''(''M''{{t}}''W''<sup>−2</sup>''M''){{inv}}''M''{{t}}''W''{{inv}} }}, and {{nowrap| {{subsup|'''m'''|''W''|T}}''P''<sub>''W''</sub>'''m'''<sub>''W''</sub> {{=}} '''m'''{{t}}''W''{{inv}}''P''<sub>''W''</sub>''W''{{inv}}'''m''' }}, or {{nowrap| '''m'''{{t}}(''W''{{inv|2}} − ''W''{{inv|2}}''M''(''M''{{t}}''W''{{inv|2}}''M''){{inv}}''M''{{t}}''W''{{inv|2}})'''m''' }}, so that the terms inside the parenthesis define a formula for ''P'' in terms of the matrix of monzos ''M''.


We can, however, map the monzos to elements of a rank r abelian group (where r is the rank of the temperament) which abstractly represents the elements of the temperament without regard to tuning, the [[abstract regular temperament]]. If b is a monzo, this mapping is given by bA*. Hence we have |1 -5 3&gt;A* maps to [0 0] for the interval associated to 250/243, and |-1 1 0&gt;A* maps to [9 13] for the interval assciated to 3/2. This is the number of steps needed to get to 3/2 in 15et and 22et respectively. We now may obtain a matrix defining the temperamental norm on this abstract temperament by //P// = (AW^2A*)^(-1), which is approximately
To define the '''octave-equivalent Tenney–Euclidean seminorm''', or '''OETES''', we simply add a column {{monzo| 1 0 0 … 0 }} representing 2 to the matrix ''M''. An alternative procedure is to find the [[normal lists #Normal val list|normal val list]], and remove the first val from the list, corresponding to the octave or some fraction thereof, and proceed as in the previous section on temperamental complexity. This seminorm is a measure of the octave-equivalent complexity of a given ''p''-limit rational interval in terms of the ''p''-limit regular temperament given by ''V''.


[175.3265 -120.0291]
== Examples ==
[-120.0291 82.1730]
Consider the temperament defined by the 5-limit [[patent val]]s for 15 and 22 equal. From the vals, we may construct a 2×3 matrix {{nowrap|''V'' {{=}} {{mapping| 15 24 35 | 22 35 51 }}}}. From this we may obtain the matrix ''P'' as ''V''{{t}}(''VW''<sup>2</sup>''V''{{t}}){{inv}}''V'', approximately


Using this, we find the temperamental norm of [9 13] to be sqrt([9 13]//P//[9 13]*) ~ sqrt(1.875) ~ 1.3693, identical to the temperamental seminorm of 3/2. Note however that while **P** does not depend on the choice of basis vals for the temperament, //P// does; if we choose [&lt;1 2 3|, &lt;0 -3 -5|] for our basis instead, then 3/2 is represented by [1 -3] and //P// changes coordinates to produce the same final result of temperamental complexity.
<math>
\left[\begin{matrix}
0.9911 & 0.1118 & -0.1440 \\
0.1118 & 1.1075 & 1.8086 \\
-0.1440 & 1.8086 & 3.0624 \\
\end{matrix}\right]
</math>


If instead we want the OETES, we may remove the first row of [&lt;1 2 3|, &lt;0 -3 -5|], leaving just [&lt;0 -3 -5|]. If we now call this 1x3 matrix A, then (AW^2A*)^(-1) is a 1x1 matrix; in effect a scalar, with value [&lt;0.1215588|]. Multiplying a monzo b times A* gives a 1x1 matrix bA* whose value is the number of generator steps of porcupine (of size a tempered 10/9) it takes to get to the octave class to which b belongs. Performing the multiplication and taking the square root, we conclude the OE complexity is simply proportional to this number of steps.
If we want to find the temperamental seminorm ''T''(250/243) of 250/243, we convert it into a monzo as {{monzo| 1 -5 3 }}. Now we may multiply ''P'' by this on the left, obtaining the zero vector. Taking the dot product of the zero vector {{monzo| 1 -5 3 }} gives zero, and taking the square root of zero we get zero, the temperametal seminorm ''T''(250/243) of 250/243. All of this is telling us that 250/243 is a comma of this temperament, which is 5-limit [[porcupine]].


For a more substantial example we need to consider at least a rank three temperament, so let us turn to marvel, the 7-limit temperament tempering out 225/224. The 2x4 matrix of monzos whose first row represents 2 and whose second row 225/224 is [|1 0 0 0&gt;, |-5 2 2 -1&gt;]. If we denote log2 of the odd primes by p3, p5, p7 etc, then the monzo weighting of this matrix is M = [|1 0 0 0&gt;, |-5 2p3 2p5 -p7&gt;], and P = I - M`M = [|1 0 0 0&gt;, |0 4p5^2+p7^2 -4p3p5 2p3p7&gt;/H,  
Similarly, starting from the monzo {{monzo| -1 1 0 }} for 3/2, we may multiply this by ''P'', obtaining {{val| -0.8793 0.9957 1.9526 }}, and taking the dot product of this with {{monzo| -1 1 0 }} gives 1.875 with square root 1.3693, which is ''T''(3/2).
|0 -4p3p5 4p3^2+p7^2 2p5p7&gt;/H, |0 2p3p7 2p5p7 4(p3^2+p5^2)&gt;/H], where H = 4p3^2+4p5^2+p7^2. On the other hand, we may start from the normal val list for the temperament, which is [&lt;1 0 0 -5|, &lt;0 1 0 2|, &lt;0 0 1 2|]. Removing the first row gives [&lt;0 1 0 2|, &lt;0 0 1 2], and val weighting this gives C = [&lt;0 1/p3 0 2/p7|, &lt;0 0 1/p5 2/p7|]. Then P = C`C is precisely the same matrix we obtained before.
 
We can, however, map the monzos to elements of a rank-''r'' abelian group (where ''r'' is the rank of the temperament) which abstractly represents the elements of the temperament without regard to tuning, the [[abstract regular temperament]]. If '''m''' is a monzo, this mapping is given by ''V'''''m'''. Hence we have ''V''{{monzo| 1 -5 3 }} maps to {{monzo| 0 0 }} for the interval associated to 250/243, and ''V''{{monzo| -1 1 0 }} maps to {{monzo| 9 13 }} for the interval assciated to 3/2. This is the number of steps needed to get to 3/2 in 15et and 22et respectively. We now may obtain a matrix defining the temperamental norm on this abstract temperament by ''P''<sub>''T''</sub> = (''VW''<sup>2</sup>''V''{{t}}){{inv}}, which is approximately
 
<math>
\left[\begin{matrix}
175.3265 & -120.0291 \\
-120.0291 & 82.1730 \\
\end{matrix}\right]
</math>
 
Using this, we find the temperamental norm of {{monzo| 9 13 }} to be {{nowrap| sqrt([9 13]''P''<sub>''T''</sub>[9 13]{{t}}) ~ sqrt(1.875) | ~ 1.3693 }}, identical to the temperamental seminorm of 3/2. Note however that while ''P'' does not depend on the choice of basis vals for the temperament, ''P''<sub>''T''</sub> does; if we choose {{mapping| 1 2 3 | 0 -3 -5 }} for our basis instead, then 3/2 is represented by {{monzo| 1 -3 }} and ''P''<sub>''T''</sub> changes coordinates to produce the same final result of temperamental complexity.
 
If instead we want the OETES, we may remove the first row of {{mapping| 1 2 3 | 0 -3 -5 }}, leaving just {{mapping| 0 -3 -5 }}. If we now call this 1×3 matrix ''V'', then {{nowrap| ''P''<sub>''T''</sub> {{=}} (''VW''<sup>2</sup>''V''{{t}}){{inv}} }} is a 1×1 matrix; in effect a scalar, with value {{mapping| 0.1215588 }}. Multiplying a monzo '''m''' by ''V'' on the left gives a 1×1 matrix ''V'''''m''' whose value is the number of generator steps of porcupine (of size a tempered 10/9) it takes to get to the octave class to which '''m''' belongs. Performing the multiplication and taking the square root, we conclude the OE complexity is simply proportional to this number of generator steps.
 
For a more substantial example we need to consider at least a rank-3 temperament, so let us turn to 7-limit marvel, the 7-limit temperament tempering out 225/224. The 2×4 matrix of monzos whose first row represents 2 and whose second row 225/224 is {{monzo list| 1 0 0 0 | -5 2 2 -1 }}. If we denote log<sub>2</sub> of the odd primes by p3, p5, p7, etc., then the monzo weighting of this matrix is {{nowrap|''M''<sub>''W''</sub> {{=}} {{monzo list| 1 0 0 0 | -5 2p3 2p5 -p7 }}}}, and {{nowrap|''P''<sub>''W''</sub> {{=}} ''I'' − ''M''<sub>''W''</sub>{{subsup|''M''|''W''|+}} }} = [{{monzo| 1 0 0 0 }}, {{monzo| 0 4(p5)<sup>2</sup> + (p7)<sup>2</sup> -4(p3)(p5) 2(p3)(p7) }}/''H'', {{monzo| 0 -4(p3)(p5) 4(p3)<sup>2</sup> + (p7)<sup>2</sup> 2(p5)(p7) }}/''H'', {{monzo| 0 2(p3)(p7) 2(p5)(p7) 4((p3)<sup>2</sup> + (p5)<sup>2</sup>) }}/''H''], where {{nowrap| ''H'' {{=}} 4(p3)<sup>2</sup> + 4(p5)<sup>2</sup> + (p7)<sup>2</sup> }}. On the other hand, we may start from the normal val list for the temperament, which is {{mapping| 1 0 0 -5 | 0 1 0 2 | 0 0 1 2 }}. Removing the first row gives {{mapping| 0 1 0 2 | 0 0 1 2 }}, and val weighting this gives {{nowrap| ''C''<sub>''W''</sub> {{=}} {{mapping| 0 1/p3 0 2/p7 | 0 0 1/p5 2/p7 }} }}. Then {{nowrap|''P''<sub>''W''</sub> {{=}} ''C''<sub>''W''</sub>{{+}}''C''<sub>''W''</sub>}} is precisely the same matrix we obtained before.


Octaves are now projected to the origin as well as commas. We can as before form the quotient space with respect to the seminorm, and obtain a normed space in which octave-equivalent interval classes of the intervals of the temperament are the lattice points. This seminorm applied to monzos gives the OE complexity.
Octaves are now projected to the origin as well as commas. We can as before form the quotient space with respect to the seminorm, and obtain a normed space in which octave-equivalent interval classes of the intervals of the temperament are the lattice points. This seminorm applied to monzos gives the OE complexity.


If we start from a normal val list and remove the first val, the remaining vals map to the octave classes of the notes of the temperament. If we call this reduced list of vals R, then the inner product on note classes in this basis is defined by the symmetric matrix S = (RW^2R*)^(-1). In the case of marvel, we obtain S = [[p3^2(4p5^2+p7^2) -4p3^2p5^2], [-4p3^2p5^2 p5^2(4p3^2+p7^2)]]/H. If k = [a b] is a note class of marvel in the coordinates defined by the truncated val list R, which in this case has a basis corresponding to tempered 3 and 5, then sqrt(kSk*) gives the OE complexity of the note class.
If we start from a normal val list and remove the first val, the remaining vals map to the octave classes of the notes of the temperament. If we call this reduced list of vals ''R'', then the inner product on note classes in this basis is defined by the symmetric matrix {{nowrap| ''S'' {{=}} (''RW''<sup>2</sup>''R''{{t}}){{inv}} }}. In the case of marvel, we obtain {{nowrap| ''S'' {{=}} {{lbrack}}[(p3)<sup>2</sup>(4(p5)<sup>2</sup> + (p7)<sup>2</sup>) -4(p3)<sup>2</sup>(p5)<sup>2</sup>]}}, {{nowrap| [-4(p3)<sup>2</sup>(p5)<sup>2</sup> (p5)<sup>2</sup>(4(p3)<sup>2</sup> + (p7)<sup>2</sup>)]]/''H'' }}. If {{nowrap| '''k''' {{=}} {{monzo| ''k''<sub>1</sub> ''k''<sub>2</sub> }} }} is a note class of marvel in the coordinates defined by the truncated val list ''R'', which in this case has a basis corresponding to tempered 3 and 5, then sqrt('''k'''{{t}}''S'''''k''') gives the OE complexity of the note class.
</pre></div>
 
<h4>Original HTML content:</h4>
[[Category:Math]]
<div style="width:100%; max-height:400pt; overflow:auto; background-color:#f8f9fa; border: 1px solid #eaecf0; padding:0em"><pre style="margin:0px;border:none;background:none;word-wrap:break-word;width:200%;white-space: pre-wrap ! important" class="old-revision-html">&lt;html&gt;&lt;head&gt;&lt;title&gt;Tenney-Euclidean metrics&lt;/title&gt;&lt;/head&gt;&lt;body&gt;&lt;!-- ws:start:WikiTextTocRule:10:&amp;lt;img id=&amp;quot;wikitext@@toc@@flat&amp;quot; class=&amp;quot;WikiMedia WikiMediaTocFlat&amp;quot; title=&amp;quot;Table of Contents&amp;quot; src=&amp;quot;/site/embedthumbnail/toc/flat?w=100&amp;amp;h=16&amp;quot;/&amp;gt; --&gt;&lt;!-- ws:end:WikiTextTocRule:10 --&gt;&lt;!-- ws:start:WikiTextTocRule:11: --&gt;&lt;a href="#The weighting matrix"&gt;The weighting matrix&lt;/a&gt;&lt;!-- ws:end:WikiTextTocRule:11 --&gt;&lt;!-- ws:start:WikiTextTocRule:12: --&gt; | &lt;a href="#Temperamental complexity"&gt;Temperamental complexity&lt;/a&gt;&lt;!-- ws:end:WikiTextTocRule:12 --&gt;&lt;!-- ws:start:WikiTextTocRule:13: --&gt; | &lt;a href="#The OETES"&gt;The OETES&lt;/a&gt;&lt;!-- ws:end:WikiTextTocRule:13 --&gt;&lt;!-- ws:start:WikiTextTocRule:14: --&gt; | &lt;a href="#Logflat TE badness"&gt;Logflat TE badness&lt;/a&gt;&lt;!-- ws:end:WikiTextTocRule:14 --&gt;&lt;!-- ws:start:WikiTextTocRule:15: --&gt;&lt;!-- ws:end:WikiTextTocRule:15 --&gt;&lt;!-- ws:start:WikiTextTocRule:16: --&gt;
[[Category:Interval space]]
&lt;!-- ws:end:WikiTextTocRule:16 --&gt;&lt;!-- ws:start:WikiTextHeadingRule:0:&amp;lt;h1&amp;gt; --&gt;&lt;h1 id="toc0"&gt;&lt;a name="The weighting matrix"&gt;&lt;/a&gt;&lt;!-- ws:end:WikiTextHeadingRule:0 --&gt;The weighting matrix&lt;/h1&gt;
[[Category:Interval complexity measures]]
Let us define the val weighting matrix W to be the diagonal matrix with values 1, 1/log2(3), 1/log2(5) ... 1/log2(p) along the diagonal. Given a val &amp;quot;a&amp;quot; expressed as a row vector, the corresponding vector in weighted coordinates is aW, with transpose Wa* where the * denotes the transpose. Then the dot product of weighted vals is aW^2a*, which makes the Euclidean metric on vals, a measure of complexity, to be || &amp;lt;a2 a3 ... ap| || = sqrt(a2^2 + a3^2/log2(3)^2 + ... + ap^2/log2(p)^2). Similarly, if b is a monzo, then in weighted coordinates the monzo becomes bW^(-1), and the dot product is bW^(-2)b*, leading to sqrt(b2^2 + log2(3)^2b3^2 + ... + log2(p)^2bp^2) as the norm on monzos, a measure of complexity we may call the Tenney-Euclidean, or TE, complexity.&lt;br /&gt;
[[Category:Tenney-weighted measures]]
&lt;br /&gt;
 
&lt;!-- ws:start:WikiTextHeadingRule:2:&amp;lt;h1&amp;gt; --&gt;&lt;h1 id="toc1"&gt;&lt;a name="Temperamental complexity"&gt;&lt;/a&gt;&lt;!-- ws:end:WikiTextHeadingRule:2 --&gt;Temperamental complexity&lt;/h1&gt;
{{Todo| reduce mathslang | improve readability }}
Suppose now A is a matrix whose rows are vals defining a p-limit regular temperament. Then the corresponding weighted matrix is V = AW. The &lt;a class="wiki_link" href="/Tenney-Euclidean%20tuning"&gt;TE&lt;/a&gt; tuning projection matrix is then V`V, where V` is the &lt;a class="wiki_link" href="/Tenney-Euclidean%20Tuning"&gt;pseudoinverse&lt;/a&gt;. If the rows of V (or equivalently, A) are linearly independent, then we have V` = V*(VV*)^(-1), where V* denotes the transpose. In terms of vals, the tuning projection matrix is P = V`V = V*(VV*)^(-1)V = WA*(AW^2A*)^(-1)AW. P is a &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Positive-definite_matrix" rel="nofollow"&gt;positive semidefinite matrix&lt;/a&gt;, so it defines a &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Definite_bilinear_form" rel="nofollow"&gt;positive semidefinite bilinear form&lt;/a&gt;. In terms of weighted monzos m1 and m2, m1Pm2* defines the semidefinite form on weighted monzos, and hence b1W^(-1)PW^(-1)b2* defines a semidefinite form on unweighted monzos, in terms of the matrix W^(-1)WA*(AW^2A*)^(-1)AWW^(-1) = A*(AW^2A*)^(-1)A = &lt;strong&gt;P&lt;/strong&gt;. From the semidefinite form we obtain an associated &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Definite_quadratic_form" rel="nofollow"&gt;semidefinite quadratic form&lt;/a&gt; b&lt;strong&gt;P&lt;/strong&gt;b* and from this the &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Norm_%28mathematics%29" rel="nofollow"&gt;seminorm&lt;/a&gt; sqrt(b&lt;strong&gt;P&lt;/strong&gt;b*). &lt;br /&gt;
&lt;br /&gt;
It may be noted that (VV*)^(-1) = (AW^2A*)^(-1) is the inverse of the &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Gramian_matrix" rel="nofollow"&gt;Gram matrix&lt;/a&gt; used to compute &lt;a class="wiki_link" href="/Tenney-Euclidean%20temperament%20measures"&gt;TE complexity&lt;/a&gt;, and hence is the corresponding Gram matrix for the dual space. Hence &lt;strong&gt;P&lt;/strong&gt; represents a change of basis defined by the mapping given by the vals combined with an &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Inner_product_space" rel="nofollow"&gt;inner product&lt;/a&gt; on the result. Given a monzo b, bA* represents the tempered interval corresponding to b in a basis defined by the mapping A, and &lt;em&gt;P&lt;/em&gt; = (AW^2A*)^(-1) defines a positive-definite quadratic form, and hence a norm, on the tempered interval space with basis defined by A.&lt;br /&gt;
&lt;br /&gt;
Denoting the temperament-defined, or temperamental, seminorm by T(x), the subspace of interval space such that T(x) = 0 contains a sublattice of the lattice of monzos consisting of the commas of the temperament. The &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Quotient_space_%28linear_algebra%29" rel="nofollow"&gt;quotient space&lt;/a&gt; of the full vector space by the commatic subspace such that T(x) = 0 is now a &lt;a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Normed_vector_space" rel="nofollow"&gt;normed vector space&lt;/a&gt; with norm given by T, in which the intervals of the regular temperament define a lattice. The norm T on these lattice points is the &lt;em&gt;temperamental norm&lt;/em&gt; or &lt;em&gt;temperamental complexity&lt;/em&gt; of the intervals of the regular temperament; in terms of the basis defined by A, it is sqrt(t&lt;em&gt;P&lt;/em&gt;t*) where t is the image of a monzo b by t = bA*&lt;br /&gt;
&lt;br /&gt;
&lt;!-- ws:start:WikiTextHeadingRule:4:&amp;lt;h1&amp;gt; --&gt;&lt;h1 id="toc2"&gt;&lt;a name="The OETES"&gt;&lt;/a&gt;&lt;!-- ws:end:WikiTextHeadingRule:4 --&gt;The OETES&lt;/h1&gt;
Instead of starting from a matrix of vals, we may start from a matrix of monzos. If B is a matrix with rows of monzos spanning the commas of a regular temperament, then M = BW^(-1) is the corresponding weighted matrix. Q = M`M is a projection matrix dual to P = I-Q, where I is the identity matrix, and P is the same symmetric matrix as in the previous section. If the rows define a basis for the commas of the temperament, and are therefor linearly independent, then P = I - M*(MM*)^(-1)M = I - W^(-1)B*(BW^(-2)B*)^(-1)BW^(-1), and mPm* = bW^(-1)PW^(-1)b*, or b(W^(-2) - W^(-2)B*(BW^(-2)B*)^(-1)BW^(-2))b*, so that the terms inside the parenthesis define a formula for &lt;strong&gt;P&lt;/strong&gt; in terms of the matrix of monzos B.&lt;br /&gt;
&lt;br /&gt;
To define the OETES, or Tenney-Euclidean octave equivalent seminorm, we simply add a row |1 0 0 ... 0&amp;gt; representing 2 to the matrix B. An alternative proceedure is to find the &lt;a class="wiki_link" href="/normal%20lists"&gt;normal val list&lt;/a&gt;, and remove the first val from the list, corresponding to the octave or some fraction thereof, and proceed as in the previous section on temperamental complexity. This seminorm is a measure of the octave-equivalent complexity of a given p-limit rational interval in terms of thenp-limit regular temperament given by A.&lt;br /&gt;
&lt;br /&gt;
&lt;!-- ws:start:WikiTextHeadingRule:6:&amp;lt;h1&amp;gt; --&gt;&lt;h1 id="toc3"&gt;&lt;a name="Logflat TE badness"&gt;&lt;/a&gt;&lt;!-- ws:end:WikiTextHeadingRule:6 --&gt;Logflat TE badness&lt;/h1&gt;
Given a matrix A whose rows are linearly independent vals defining a regular temperament, then the rank r of the temperament is the number of rows, which equals the number of linearly independent vals. The dimension of the temperament is the number of primes it covers; if p is the largest such prime, then the dimension d is pi(p), the number of primes to p. If we define S(A) to be the simple badness (relative error) of A, and C(A) to be the complexity of A, then &lt;em&gt;logflat badness&lt;/em&gt; is defined by the formula S(A) C(A)^(r/(d-r)). If we set a cutoff margin for logflat badness, there are still infinite numbers of new temperaments appearing as complexity goes up, at a lower rate which is approximately logarithmic in terms of complexity.&lt;br /&gt;
&lt;br /&gt;
&lt;!-- ws:start:WikiTextHeadingRule:8:&amp;lt;h2&amp;gt; --&gt;&lt;h2 id="toc4"&gt;&lt;a name="Logflat TE badness-Examples"&gt;&lt;/a&gt;&lt;!-- ws:end:WikiTextHeadingRule:8 --&gt;Examples&lt;/h2&gt;
Consider the temperament defined by the 5-limit &lt;a class="wiki_link" href="/Patent%20val"&gt;patent vals&lt;/a&gt; for 15 and 22 equal. From the vals, we may contruct a 2x3 matrix A = [&amp;lt;15 24 35|, &amp;lt;22 35 51|]. From this we may obtain the matrix &lt;strong&gt;P&lt;/strong&gt; as A*(AW^2A*)^(-1)A, approximately &lt;br /&gt;
&lt;br /&gt;
[0.9911 0.1118 -0.1440]&lt;br /&gt;
[0.1118 1.1075 1.8086]&lt;br /&gt;
[-0.1440 1.8086 3.0624]&lt;br /&gt;
&lt;br /&gt;
If we want to find the temperamental seminorm T(250/243) of 250/243, we convert it into a monzo as |1 -5 3&amp;gt;. Now we may multiply &lt;strong&gt;P&lt;/strong&gt; by this on the left, obtaining the zero vector. Taking the dot product of the zero vector with |1 -5 3&amp;gt; gives zero, and taking the square root of zero we get zero, the temperametal seminorm T(250/243) of 250/243. All of this is telling us that 250/243 is a comma of this temperament, which is 5-limit &lt;a class="wiki_link" href="/Porcupine%20family"&gt;porcupine&lt;/a&gt;.&lt;br /&gt;
&lt;br /&gt;
Similarly, starting from the monzo |-1 1 0&amp;gt; for 3/2, we may multiply this by &lt;strong&gt;P&lt;/strong&gt;, obtaining &amp;lt;-0.8793 0.9957 1.9526|, and taking the dot product of this with |-1 1 0&amp;gt; gives 1.875 with square root 1.3693, which is T(3/2).&lt;br /&gt;
&lt;br /&gt;
We can, however, map the monzos to elements of a rank r abelian group (where r is the rank of the temperament) which abstractly represents the elements of the temperament without regard to tuning, the &lt;a class="wiki_link" href="/abstract%20regular%20temperament"&gt;abstract regular temperament&lt;/a&gt;. If b is a monzo, this mapping is given by bA*. Hence we have |1 -5 3&amp;gt;A* maps to [0 0] for the interval associated to 250/243, and |-1 1 0&amp;gt;A* maps to [9 13] for the interval assciated to 3/2. This is the number of steps needed to get to 3/2 in 15et and 22et respectively. We now may obtain a matrix defining the temperamental norm on this abstract temperament by &lt;em&gt;P&lt;/em&gt; = (AW^2A*)^(-1), which is approximately&lt;br /&gt;
&lt;br /&gt;
[175.3265 -120.0291]&lt;br /&gt;
[-120.0291 82.1730]&lt;br /&gt;
&lt;br /&gt;
Using this, we find the temperamental norm of [9 13] to be sqrt([9 13]&lt;em&gt;P&lt;/em&gt;[9 13]*) ~ sqrt(1.875) ~ 1.3693, identical to the temperamental seminorm of 3/2. Note however that while &lt;strong&gt;P&lt;/strong&gt; does not depend on the choice of basis vals for the temperament, &lt;em&gt;P&lt;/em&gt; does; if we choose [&amp;lt;1 2 3|, &amp;lt;0 -3 -5|] for our basis instead, then 3/2 is represented by [1 -3] and &lt;em&gt;P&lt;/em&gt; changes coordinates to produce the same final result of temperamental complexity.&lt;br /&gt;
&lt;br /&gt;
If instead we want the OETES, we may remove the first row of [&amp;lt;1 2 3|, &amp;lt;0 -3 -5|], leaving just [&amp;lt;0 -3 -5|]. If we now call this 1x3 matrix A, then (AW^2A*)^(-1) is a 1x1 matrix; in effect a scalar, with value [&amp;lt;0.1215588|]. Multiplying a monzo b times A* gives a 1x1 matrix bA* whose value is the number of generator steps of porcupine (of size a tempered 10/9) it takes to get to the octave class to which b belongs. Performing the multiplication and taking the square root, we conclude the OE complexity is simply proportional to this number of steps.&lt;br /&gt;
&lt;br /&gt;
For a more substantial example we need to consider at least a rank three temperament, so let us turn to marvel, the 7-limit temperament tempering out 225/224. The 2x4 matrix of monzos whose first row represents 2 and whose second row 225/224 is [|1 0 0 0&amp;gt;, |-5 2 2 -1&amp;gt;]. If we denote log2 of the odd primes by p3, p5, p7 etc, then the monzo weighting of this matrix is M = [|1 0 0 0&amp;gt;, |-5 2p3 2p5 -p7&amp;gt;], and P = I - M`M = [|1 0 0 0&amp;gt;, |0 4p5^2+p7^2 -4p3p5 2p3p7&amp;gt;/H, &lt;br /&gt;
|0 -4p3p5 4p3^2+p7^2 2p5p7&amp;gt;/H, |0 2p3p7 2p5p7 4(p3^2+p5^2)&amp;gt;/H], where H = 4p3^2+4p5^2+p7^2. On the other hand, we may start from the normal val list for the temperament, which is [&amp;lt;1 0 0 -5|, &amp;lt;0 1 0 2|, &amp;lt;0 0 1 2|]. Removing the first row gives [&amp;lt;0 1 0 2|, &amp;lt;0 0 1 2], and val weighting this gives C = [&amp;lt;0 1/p3 0 2/p7|, &amp;lt;0 0 1/p5 2/p7|]. Then P = C`C is precisely the same matrix we obtained before.&lt;br /&gt;
&lt;br /&gt;
Octaves are now projected to the origin as well as commas. We can as before form the quotient space with respect to the seminorm, and obtain a normed space in which octave-equivalent interval classes of the intervals of the temperament are the lattice points. This seminorm applied to monzos gives the OE complexity.&lt;br /&gt;
&lt;br /&gt;
If we start from a normal val list and remove the first val, the remaining vals map to the octave classes of the notes of the temperament. If we call this reduced list of vals R, then the inner product on note classes in this basis is defined by the symmetric matrix S = (RW^2R*)^(-1). In the case of marvel, we obtain S = [[p3^2(4p5^2+p7^2) -4p3^2p5^2], [-4p3^2p5^2 p5^2(4p3^2+p7^2)]]/H. If k = [a b] is a note class of marvel in the coordinates defined by the truncated val list R, which in this case has a basis corresponding to tempered 3 and 5, then sqrt(kSk*) gives the OE complexity of the note class.&lt;/body&gt;&lt;/html&gt;</pre></div>

Latest revision as of 09:54, 15 April 2025

The Tenney-Euclidean metrics are metrics defined in Tenney-Euclidean space. These consist of the TE norm, which measures the complexity of an interval in just intonation, the TE temperamental norm, which measures the complexity of an interval as mapped by a temperament, and the octave-equivalent TE seminorms of both.

TE norm

The Tenney–Euclidean norm (TE norm) or Tenney–Euclidean complexity (TE complexity) applies to vals as well as to monzos.

Let us define the val weighting matrix W to be the diagonal matrix with values 1, 1/log23, 1/log25 … 1/log2p along the diagonal. For the p-limit prime basis Q = 2 3 5 … p],

$$ W = \operatorname {diag} (1/\log_2 (Q)) $$

Right-multiplying a row vector by this matrix scales each entry by the corresponding entry of the diagonal matrix.

Given a val V expressed as a row vector, the corresponding row vector in weighted coordinates is VW = VW, with transpose V T
W
 
= WV T
where  T denotes the transpose. The dot product of a weighted val with itself, or the sum of the squares of its entries, is the squared Euclidean metric of the val, VW 2
2
 
= VWV T
W
 
= VW2V T
. Thus the Euclidean metric on the val, a measure of complexity, is VW2 = sqrt(VWV T
W
 
)
= sqrt(v 2
1
 
+ v 2
2
 
/(log23)2 + … + v 2
n
 
/(log2p)2)
, where n = π(p) is the prime-counting function which records the number of primes to p; dividing this by sqrt(n) gives the TE norm of a val.

Similarly, if m is a monzo, then in weighted coordinates the monzo becomes mW = W−1m, and the dot product is m T
W
 
mW = m T
W-2m, leading to sqrt(m T
W
 
mW) = sqrt(m 2
1
 
+ (log23)2m 2
2
 
+ … + (log2p)2m 2
n
 
)
; multiplying this by sqrt(n) gives the dual RMS norm on monzos which serves as a measure of complexity.

TE temperamental norm

Suppose now V is a matrix whose rows are vals defining a p-limit regular temperament. Then the corresponding weighted matrix is VW = VW. The TE tuning projection matrix is then PW = V +
W
 
VW
, where + denotes the Moore–Penrose pseudoinverse. If the rows of VW (or equivalently, V) are linearly independent, then we have V +
W
 
= V T
W
 
(VWV T
W
 
)−1
. In terms of vals, the tuning projection matrix is V +
W
 
VW = V T
W
 
(VWV T
W
 
)−1VW
= WV T(VW2V T)−1VW. PW is a positive semidefinite matrix, so it defines a positive semidefinite bilinear form. In terms of weighted monzos (mW)1 and (mW)2, (mW) T
1
 
PW(mW)2 defines the semidefinite form on weighted monzos, and hence m T
1
 
W−1PWW−1m2 defines a semidefinite form on unweighted monzos, in terms of the matrix P = W−1PWW−1 = V T(VW2V T)−1V. From the semidefinite form we obtain an associated semidefinite quadratic form m TPm and from this the seminorm sqrt(m TPm).

It may be noted that (VWV T
W
 
)−1 = (VW2V T
)−1 is the inverse of the Gramian matrix used to compute TE complexity, and hence is the corresponding Gram matrix for the dual space. Hence P represents a change of basis defined by the mapping given by the vals combined with an inner product on the result. Given a monzo m, Vm represents the tempered interval corresponding to m in a basis defined by the mapping V, and PT = (VW2V T)−1 defines a positive-definite quadratic form, and hence a norm, on the tempered interval space with basis defined by V.

Denoting the temperament-defined, or temperamental, seminorm by T(x), the subspace of interval space such that T(x) = 0 contains a lattice consisting of the commas of the temperament, which is a sublattice of the lattice of monzos. The quotient space of the full vector space by the commatic subspace such that T(x) = 0 is now a normed vector space with norm given by T, in which the intervals of the regular temperament define a lattice. The norm T on these lattice points is the TE temperamental norm or TE temperamental complexity of the intervals of the regular temperament; in terms of the basis defined by V, it is sqrt(t TPTt) where t is the image of a monzo m by t = Vm.

Octave-equivalent TE seminorm

Instead of starting from a matrix of vals, we may start from a matrix of monzos. If M is a matrix with columns of monzos spanning the commas of a regular temperament, then MW = W−1M is the corresponding weighted matrix. QW = MWM +
W
 
is a projection matrix dual to PW = IQW, where I is the identity matrix, and PW is the same symmetric matrix as in the previous section. If the rows define a basis for the commas of the temperament, and are therefore linearly independent, then PW = IMW(M T
W
 
MW)−1M T
W
 
= IW−1M(M TW−2M)−1M TW−1, and m T
W
 
PWmW = m T
W−1PWW−1m, or m T(W−2W−2M(M TW−2M)−1M TW−2)m, so that the terms inside the parenthesis define a formula for P in terms of the matrix of monzos M.

To define the octave-equivalent Tenney–Euclidean seminorm, or OETES, we simply add a column [1 0 0 … 0 representing 2 to the matrix M. An alternative procedure is to find the normal val list, and remove the first val from the list, corresponding to the octave or some fraction thereof, and proceed as in the previous section on temperamental complexity. This seminorm is a measure of the octave-equivalent complexity of a given p-limit rational interval in terms of the p-limit regular temperament given by V.

Examples

Consider the temperament defined by the 5-limit patent vals for 15 and 22 equal. From the vals, we may construct a 2×3 matrix V = [15 24 35], 22 35 51]]. From this we may obtain the matrix P as V T(VW2V T)−1V, approximately

[math]\displaystyle{ \left[\begin{matrix} 0.9911 & 0.1118 & -0.1440 \\ 0.1118 & 1.1075 & 1.8086 \\ -0.1440 & 1.8086 & 3.0624 \\ \end{matrix}\right] }[/math]

If we want to find the temperamental seminorm T(250/243) of 250/243, we convert it into a monzo as [1 -5 3. Now we may multiply P by this on the left, obtaining the zero vector. Taking the dot product of the zero vector [1 -5 3 gives zero, and taking the square root of zero we get zero, the temperametal seminorm T(250/243) of 250/243. All of this is telling us that 250/243 is a comma of this temperament, which is 5-limit porcupine.

Similarly, starting from the monzo [-1 1 0 for 3/2, we may multiply this by P, obtaining -0.8793 0.9957 1.9526], and taking the dot product of this with [-1 1 0 gives 1.875 with square root 1.3693, which is T(3/2).

We can, however, map the monzos to elements of a rank-r abelian group (where r is the rank of the temperament) which abstractly represents the elements of the temperament without regard to tuning, the abstract regular temperament. If m is a monzo, this mapping is given by Vm. Hence we have V[1 -5 3 maps to [0 0 for the interval associated to 250/243, and V[-1 1 0 maps to [9 13 for the interval assciated to 3/2. This is the number of steps needed to get to 3/2 in 15et and 22et respectively. We now may obtain a matrix defining the temperamental norm on this abstract temperament by PT = (VW2V T)−1, which is approximately

[math]\displaystyle{ \left[\begin{matrix} 175.3265 & -120.0291 \\ -120.0291 & 82.1730 \\ \end{matrix}\right] }[/math]

Using this, we find the temperamental norm of [9 13 to be sqrt([9 13]PT[9 13] T) ~ sqrt(1.875), ~ 1.3693, identical to the temperamental seminorm of 3/2. Note however that while P does not depend on the choice of basis vals for the temperament, PT does; if we choose [1 2 3], 0 -3 -5]] for our basis instead, then 3/2 is represented by [1 -3 and PT changes coordinates to produce the same final result of temperamental complexity.

If instead we want the OETES, we may remove the first row of [1 2 3], 0 -3 -5]], leaving just [0 -3 -5]]. If we now call this 1×3 matrix V, then PT = (VW2V T)−1 is a 1×1 matrix; in effect a scalar, with value [0.1215588]]. Multiplying a monzo m by V on the left gives a 1×1 matrix Vm whose value is the number of generator steps of porcupine (of size a tempered 10/9) it takes to get to the octave class to which m belongs. Performing the multiplication and taking the square root, we conclude the OE complexity is simply proportional to this number of generator steps.

For a more substantial example we need to consider at least a rank-3 temperament, so let us turn to 7-limit marvel, the 7-limit temperament tempering out 225/224. The 2×4 matrix of monzos whose first row represents 2 and whose second row 225/224 is [[1 0 0 0, [-5 2 2 -1]. If we denote log2 of the odd primes by p3, p5, p7, etc., then the monzo weighting of this matrix is MW = [[1 0 0 0, [-5 2p3 2p5 -p7], and PW = IMWM +
W
 
= [[1 0 0 0, [0 4(p5)2 + (p7)2 -4(p3)(p5) 2(p3)(p7)/H, [0 -4(p3)(p5) 4(p3)2 + (p7)2 2(p5)(p7)/H, [0 2(p3)(p7) 2(p5)(p7) 4((p3)2 + (p5)2)/H], where H = 4(p3)2 + 4(p5)2 + (p7)2. On the other hand, we may start from the normal val list for the temperament, which is [1 0 0 -5], 0 1 0 2], 0 0 1 2]]. Removing the first row gives [0 1 0 2], 0 0 1 2]], and val weighting this gives CW = [0 1/p3 0 2/p7], 0 0 1/p5 2/p7]]. Then PW = CW+CW is precisely the same matrix we obtained before.

Octaves are now projected to the origin as well as commas. We can as before form the quotient space with respect to the seminorm, and obtain a normed space in which octave-equivalent interval classes of the intervals of the temperament are the lattice points. This seminorm applied to monzos gives the OE complexity.

If we start from a normal val list and remove the first val, the remaining vals map to the octave classes of the notes of the temperament. If we call this reduced list of vals R, then the inner product on note classes in this basis is defined by the symmetric matrix S = (RW2R T)−1. In the case of marvel, we obtain S = [[(p3)2(4(p5)2 + (p7)2) -4(p3)2(p5)2], [-4(p3)2(p5)2 (p5)2(4(p3)2 + (p7)2)]]/H. If k = [k1 k2 is a note class of marvel in the coordinates defined by the truncated val list R, which in this case has a basis corresponding to tempered 3 and 5, then sqrt(k TSk) gives the OE complexity of the note class.