This post gives an example of some calculations that I did using my own machine learning algorithm. These calculations work out nicely which indicates that the machine learning algorithm I am using is interpretable (and it is much more interpretable than any neural network would be). These calculations show that one can begin with old mathematical structures and produce new mathematical structures, and it seems feasible to completely automate this process to continue to produce more mathematical structures. The machine learning models that I use are linear, but it seems like we can get highly non-trivial results simply by iterating the procedure of obtaining new structures from old using machine learning.
I made a similar post to this one about 7 months ago, but I decided to revisit this experiment with more general algorithms and I have obtained experimental results which I think look nice.
To illustrate how this works, we start off with the octonions. The octonions consists of an 8-dimensional inner product space V together with a bilinear operation ∗ and a unit 1∈V where 1∗v=v∗1=v for all v∈V and where ∥u∗v∥=∥u∥⋅∥v∥ for all u,v∈V. The octonions are uniquely determined up to isomorphism from these properties. The operation ∗ is non-associative, but the ∗ is closely related to the quaternions and complex numbers. If we take a single element in v∈V∖Span(1), then {1,v} generates a subalgebra of (V,∗) isomorphic to the field of complex numbers, and if u,v∈V and {1,u,v} are linearly independent, then {1,u,v,u∗v} spans a subalgebra of V isomorphic to the division ring of quaternions. For this reason, one commonly thinks of the octonions as the best way to extend the division ring of quaternions to a larger algebraic structure in the same way that the quaternions extend the field of complex numbers. But since the octonions are non-associative, they cannot be used to construct matrices, so they are not as well-known as the quaternions (and the construction of the octonions is more complicated too)
Suppose now that {e0,e1,…,e7} is an orthonormal basis for the division ring of octonions with e0=1. Then define matrices A0,…,A7:V→V by setting Ajv=ej∗v for all j. Our goal is to transform (A0,…,A7) into other tuples of matrices that satisfy similar properties.
If (A1,…,Ar),(B1,…,Br) are matrices, then define the L2
-spectral radius similarity between (A1,…,Ar) and (B1,…,Br) as
where ρ denotes the spectral radius, ⊗ is the tensor product, and ¯¯¯¯¯X is the complex conjugate of X applied elementwise.
Let d∈{1,…,8}, and let Fd,Gd,Hd denote the maximum value of the fitness level 8⋅∥(A0,…,A7)≃(X0,…,X7)∥2 such that each Xj is a complex d×d anti-symmetric matrix (X=−XT), a complex d×d symmetric matrix (X=XT), and a complex d×d-Hermitian matrix (X=X∗) respectively.
The following calculations were obtained through gradient ascent, so I have no mathematical proof that the values obtained are actually correct.
G1=2, H1=1
G2=3, H2=3
F3=1+√3 , G3=3.5, H3=1+2√2
F4=4 , G4=4, H4=1+3√2
F5=(5+√13)/2, G5=4.5, H5≈5.27155841
F6=5 , G6=5, H6=3+2√2
F7=6 , G7=2+2√3≈5.4641, H7=1+2√7
F8=7 , G8=6, H8=7
Observe that with at most one exception, all of these values Fd,Gd,Hd are algebraic half integers. This indicates that the fitness function that we maximize to produce Fd,Gd,Hd behaves mathematically and can be used to produce new tuples (X1,…,Xr) from old ones (A1,…,Ar). Furthermore, an AI can determine whether something notable is going on with the new tuple (X1,…,Xr) in several ways. For example, if ∥(A1,…,Ar)≃(X1,…,Xr)∥2 has low algebraic degree at the local maximum, then (X1,…,Xr) is likely notable and likely behaves mathematically (and is probably quite interpretable too).
The good behavior of Fd,Gd,Hd demonstrates that the octonions are compatible with the L2-spectral radius similarity. The operators (A0,…,A7) are all orthogonal, and one can take the tuple (A0,…,A7) as a mixed unitary quantum channel that is very similar to the completely depolarizing channel. The completely depolarizing channel completely mixes every quantum state while the mixture of orthogonal mappings (A0,…,A7) completely mixes every real state. The L2-spectral radius similarity works very well with the completely depolarizing channel, so one should expect for the L2-spectral radius similarity to also behave well with the octonions.
This post gives an example of some calculations that I did using my own machine learning algorithm. These calculations work out nicely which indicates that the machine learning algorithm I am using is interpretable (and it is much more interpretable than any neural network would be). These calculations show that one can begin with old mathematical structures and produce new mathematical structures, and it seems feasible to completely automate this process to continue to produce more mathematical structures. The machine learning models that I use are linear, but it seems like we can get highly non-trivial results simply by iterating the procedure of obtaining new structures from old using machine learning.
I made a similar post to this one about 7 months ago, but I decided to revisit this experiment with more general algorithms and I have obtained experimental results which I think look nice.
To illustrate how this works, we start off with the octonions. The octonions consists of an 8-dimensional inner product space V together with a bilinear operation ∗ and a unit 1∈V where 1∗v=v∗1=v for all v∈V and where ∥u∗v∥=∥u∥⋅∥v∥ for all u,v∈V. The octonions are uniquely determined up to isomorphism from these properties. The operation ∗ is non-associative, but the ∗ is closely related to the quaternions and complex numbers. If we take a single element in v∈V∖Span(1), then {1,v} generates a subalgebra of (V,∗) isomorphic to the field of complex numbers, and if u,v∈V and {1,u,v} are linearly independent, then {1,u,v,u∗v} spans a subalgebra of V isomorphic to the division ring of quaternions. For this reason, one commonly thinks of the octonions as the best way to extend the division ring of quaternions to a larger algebraic structure in the same way that the quaternions extend the field of complex numbers. But since the octonions are non-associative, they cannot be used to construct matrices, so they are not as well-known as the quaternions (and the construction of the octonions is more complicated too)
Suppose now that {e0,e1,…,e7} is an orthonormal basis for the division ring of octonions with e0=1. Then define matrices A0,…,A7:V→V by setting Ajv=ej∗v for all j. Our goal is to transform (A0,…,A7) into other tuples of matrices that satisfy similar properties.
If (A1,…,Ar),(B1,…,Br) are matrices, then define the L2
-spectral radius similarity between (A1,…,Ar) and (B1,…,Br) as
∥(A1,…,Ar)≃(B1,…,Br)∥2=
ρ(A1⊗¯¯¯¯¯¯B1+⋯+Ar⊗¯¯¯¯¯¯Br)ρ(A1⊗¯¯¯¯¯¯A1+⋯+Ar⊗¯¯¯¯¯Ar)1/2⋅ρ(B1⊗¯¯¯¯¯¯B1+⋯+Br⊗¯¯¯¯¯¯Br)1/2
where ρ denotes the spectral radius, ⊗ is the tensor product, and ¯¯¯¯¯X is the complex conjugate of X applied elementwise.
Let d∈{1,…,8}, and let Fd,Gd,Hd denote the maximum value of the fitness level 8⋅∥(A0,…,A7)≃(X0,…,X7)∥2 such that each Xj is a complex d×d anti-symmetric matrix (X=−XT), a complex d×d symmetric matrix (X=XT), and a complex d×d-Hermitian matrix (X=X∗) respectively.
The following calculations were obtained through gradient ascent, so I have no mathematical proof that the values obtained are actually correct.
G1=2, H1=1
G2=3, H2=3
F3=1+√3 , G3=3.5, H3=1+2√2
F4=4 , G4=4, H4=1+3√2
F5=(5+√13)/2, G5=4.5, H5≈5.27155841
F6=5 , G6=5, H6=3+2√2
F7=6 , G7=2+2√3≈5.4641, H7=1+2√7
F8=7 , G8=6, H8=7
Observe that with at most one exception, all of these values Fd,Gd,Hd are algebraic half integers. This indicates that the fitness function that we maximize to produce Fd,Gd,Hd behaves mathematically and can be used to produce new tuples (X1,…,Xr) from old ones (A1,…,Ar). Furthermore, an AI can determine whether something notable is going on with the new tuple (X1,…,Xr) in several ways. For example, if ∥(A1,…,Ar)≃(X1,…,Xr)∥2 has low algebraic degree at the local maximum, then (X1,…,Xr) is likely notable and likely behaves mathematically (and is probably quite interpretable too).
The good behavior of Fd,Gd,Hd demonstrates that the octonions are compatible with the L2-spectral radius similarity. The operators (A0,…,A7) are all orthogonal, and one can take the tuple (A0,…,A7) as a mixed unitary quantum channel that is very similar to the completely depolarizing channel. The completely depolarizing channel completely mixes every quantum state while the mixture of orthogonal mappings (A0,…,A7) completely mixes every real state. The L2-spectral radius similarity works very well with the completely depolarizing channel, so one should expect for the L2-spectral radius similarity to also behave well with the octonions.