In this post, I shall describe a fitness function that can be locally maximized without gradient computations. This fitness function is my own. I initially developed this fitness function in order to evaluate block ciphers for cryptocurrency technologies, but I later found that this fitness function may be used to solve other problems such as the clique problem (which is NP-complete) in the average case and some natural language processing tasks as well. After describing algorithms for locally maximizing this fitness function, I conclude that such a fitness function is inherently interpretable and mathematical which is what we need for AI safety.

Let $K$ denote either the field of real numbers, complex numbers, or the division ring of quaternions. Given $(A_1,\dots ,A_r)\in M_n(K{)}^r$ and $1\le d<n$, the goal is to find a tuple $(X_1,\dots ,X_r)\in M_d(K{)}^r$ most similar to $(A_1,\dots ,A_r)$. In other words, we want to approximate the $n\times n$-matrices $(A_1,\dots ,A_r)$ with $d\times d$-matrices.

Suppose that $A_1,\dots ,A_r\in M_n\left(K\right),X_1,\dots ,X_r\in M_d\left(K\right).$ Define a function $\Gamma (A_1,\dots ,A_r;X_1,\dots ,X_r):M_{n,d}\left(K\right)\to M_{n,d}\left(K\right)$ by setting

$\Gamma (A_1,\dots ,A_r;X_1,\dots ,X_r)\left(X\right)={\sum }_j{A}_{j}X{X}_j^{\ast }$, and define$\Phi (A_1,\dots ,A_r)=\Gamma (A_1,\dots ,A_r;A_1,\dots ,A_r)$. 

If $X$ is a matrix, then let $\rho \left(X\right)$ denote the spectral radius of $X$.

Define the $L_2$-spectral radius similarity between $(A_1,\dots ,A_r)$ and $(X_1,\dots ,X_r)$

by setting $\parallel (A_1,\dots ,A_r)\simeq (X_1,\dots ,X_r){\parallel }_2$

$=\frac{\rho \left(\Gamma \right(A_1,\dots ,A_r;X_1,\dots ,X_r\left)\right)}{\rho \left(\Phi \right(A_1,\dots ,A_r\left){)}^{1/2}\rho \right(\Phi (X_1,\dots ,X_r){)}^{1/2}}.$ 

The quantity $\parallel (A_1,\dots ,A_r)\simeq (X_1,\dots ,X_r){\parallel }_2$ is always a real number in the interval $[0,1]$ (the proof is a generalization of the Cauchy-Schwarz inequality).

If $1\le d<n$, and $A_1,\dots ,A_r\in M_n\left(K\right)$, then we say that $(X_1,\dots ,X_r)\in M_d(K{)}^r$ is an $L_{2,d}$-spectral radius dimensionality reduction (LSRDR) of $(A_1,\dots ,A_r)$ if the similarity $\parallel (A_1,\dots ,A_r)\simeq (X_1,\dots ,X_r){\parallel }_2$ is locally maximized.

One can produce an LSRDR of $(A_1,\dots ,A_r)$ simply by locally maximizing $(X_1,\dots ,X_r)$ using gradient ascent. But there is another way to obtain LSRDRs since they behave mathematically.

Let $Z\left(K\right)$ denote the center of the algebra $K$. In particular,$Z\left(R\right)=R,Z\left(C\right)=C,Z\left(H\right)=R$.

Empirical observation: If $(X_1,\dots ,X_r)\in M_d(K{)}^r$ is an LSRDR of $(A_1,\dots ,A_r)$, then typically exists $\lambda \in Z\left(K\right)$ along with matrices $R,S$ with $X_j=\lambda R{A}_{j}S$ for all $j$ and where if $P=SR$, then $P$ is a (typically non-orthogonal) projection matrix. The projection matrix $P$ is typically unique in the sense that if we train an LSRDR of $(A_1,\dots ,A_r)$ with rank $d$ multiple times, then we generally end up with the same projection matrix $P$. If the projection matrix $P$ is unique, then we shall call $P$ the canonical LSRDR projection matrix of rank $d$. Let $H$ denote the dominant eigenvector of $\Gamma (A_1,\dots ,A_r;P{A}_{1}P,\dots ,P{A}_{r}P)$ with $\text{Tr}\left(H\right)=1$ (in the quaternionic case, we should define the trace as the real part of the sum of diagonal entries in the matrix). Let $G$ denote the dominant eigenvector of $\Gamma (A_1,\dots ,A_r;P{A}_{1}P,\dots ,P{A}_{r}P{)}^{\ast }$ with $\text{Tr}\left(G\right)=1$. Then $G,H$ are positive semidefinite matrices with $\text{Im}\left(G\right)=\text{Im}\left(P^{\ast }\right)=\ker (P{)}^{\perp }$ and $\text{Im}\left(H\right)=\text{Im}\left(P\right)=\ker (P^{\ast }{)}^{\perp }$.

Said differently, $P,G,H$ along with the eigenvalue $\lambda$ is a solution to the following equations:

1. $P^2=P$.
2. $\Gamma (A_1,\dots ,A_r;P{A}_{1}P,\dots ,P{A}_{r}P)\left(H\right)=\lambda H$.
3. $\Gamma (A_1,\dots ,A_r;P{A}_{1}P,\dots ,P{A}_{r}P{)}^{\ast }(G)=\overline{\lambda }G$.
4. $\text{Im}\left(H\right)=\text{Im}\left(P\right)$.
5. $\text{Im}\left(G\right)=\text{Im}\left(P^{\ast }\right)$.
6. $\text{Tr}\left(G\right)=\text{Tr}\left(H\right)=1$.

One may often solve an equation or system of equations using iteration. For example, if $(X,d)$ is a complete metric space, $\alpha <1$ and $f:X\to X$ is a function with $d\left(f\right(x),f(y\left)\right)\le \alpha \cdot d(x,y)$ for all $x,y\in X$, and $y_0={lim}_{n\to \infty }{f}^n\left(x_0\right)$, then $f\left(y_0\right)=y_0$ by the contraction mapping theorem. We shall also apply this idea to produce the $P,G,H$ in an LSRDR since the $P,G,H$ satisfy a system of equations.

We need a function $f$ that can be easily applied to matrices but where each projection matrix $P$ is an attractive fixed point for $f$. Consider the function $f\left(x\right)=3x^2-2x^3$. We observe that $f\left(0\right)=0,f\left(1\right)=1,f^{\prime }\left(0\right)=0=f^{\prime }\left(1\right)$. Furthermore, if $x\in [0,1/2)$, then the sequence $f^n\left(x\right)$ converges to $0$ very quickly and if $x\in (1/2,1]$, then $f^n\left(x\right)$ converges to $1$ very quickly. Actually, there are neighborhoods $U$ of $0$ and $V$ of $1$ in the complex plane where if $w\in U$, then $f^n\left(w\right)$ converges to $0$ quickly and if $z\in V$, then $f^n\left(z\right)$ converges to $1$ quickly. As a consequence, if $P$ is a projection matrix, then there is a neighborhood $O$ of $P$ where if $Q\in O$, then $f^n\left(Q\right)$ converges to some projection matrix very quickly. Let Define a partial function $f^{\infty }$ by setting $f^{\infty }\left(Q\right)={lim}_{n\to \infty }{f}^n\left(Q\right)$.

Iterative algorithm for computing LSRDRs: Let $t>0$ but set $t$ to a sufficiently small value. Let $s\in N\cup \left\{\infty \right\}\setminus \left\{0\right\}$. Let $P_0$ be an random orthonormal projection of rank $d$. Let $G_0,H_0$ be random matrices in $M_n\left(K\right)$. Define $G_N,H_N,P_N,G_N^{♯},H_N^{♯}$ for all $N\ge 0$ recursively by setting

1. $G_N^{♯}=\Gamma (A_1,\dots ,A_r;P_N{A}_1{P}_N,\dots ,P_N{A}_r{P}_N{)}^{\ast }(G_N)$,
2. $H_N^{♯}=\Gamma (A_1,\dots ,A_r;P_N{A}_1{P}_N,\dots ,P_N{A}_r{P}_N)\left(H_N\right)$,
3. $G_{N+1}=G_N^{♯}/\text{Tr}\left(G_N^{♯}\right)$,
4. $H_{N+1}=H_N^{♯}/\text{Tr}\left(H_N^{♯}\right)$, and
5. $P_{N+1}=f^s(P_N+t\cdot (P_N{G}_N^{\ast }+H_N{P}_N\left)\right)$.

Then $(P_N{)}_N$ typically converges to a matrix $P$ of rank $d$ at an exponential rate. Furthermore, the matrix $P$ typically does not depend on the initialization $G_0,H_0,P_0$, but $P$ does depend on $t$ and $s$. Therefore set $P={\text{Pro}}_s(A_1,\dots ,A_r)$ whenever $P$ does not depend on the initialization (just assume that the limit converges). If $s=\infty$, then $P$ is typically the canonical rank $d$ LSRDR projection operator. 

The iterative algorithm for computing LSRDRs is a proper extension of the notion of an LSRDR in some cases. For example, if we simply count the number of parameters, the matrices $R,S$ have $2\cdot n\cdot d$ parameters while the matrices $(X_1,\dots ,X_r)$ have $d^2\cdot r$ parameters. Therefore, if $d\cdot r\ll 2\cdot n$, then the matrices $R,S$ (and the projection matrix $P$) have more parameters than $(X_1,\dots ,X_r)$. This means that we cannot obtain a unique canonical LSRDR projection of dimension $d$ when $d\cdot r\ll 2\cdot n$. On the other hand, even when $d\cdot r\ll 2\cdot n$, the matrix ${\text{Pro}}_{\infty }(A_1,\dots ,A_r)$ typically exists, and if we set $P={\text{Pro}}_{\infty }(A_1,\dots ,A_r)$, there are $R,S$ where $P=SR$ and where $(R{A}_{1}S,\dots ,R{A}_{r}S)$ is an LSRDR of $(A_1,\dots ,A_r)$. This means that the iterative algorithm for computing LSRDRs gives more information than simply an LSRDR. The iterative algorithm for computing LSRDRs gives a projection operator.

Interpretability: Since there are several ways of computing LSRDRs, LSRDRs behave mathematically. A machine learning algorithm that typically produces the same output is more interpretable than one that does not for a few reasons including the conclusion that when there is only one output of a machine learning algorithm, that one output only depends on the input and it does not have any other source of random information contributing to it. Since we typically (but not always) attain the same local maximum when training LSRDRs multiple times, LSRDRs are both interpretable and mathematical. This sort of mathematical behavior is what we need to make sense of the inner workings of LSRDRs and other machine learning algorithms. There are several ways to generalize the notion of an LSRDR, but these generalized machine learning algorithms still tend to behave mathematically; they still tend to produce the exact same trained model after training multiple times.

Capabilities: Trained LSRDRs can solve NP-complete problems such as the clique problem. I have also trained LSRDRs to produce word embeddings for natural language processing and to analyze the octonions. I can generalize LSRDRs so that they behave more like deep neural networks, but we still have a way to go until these generalized LSRDRs perform as well as our modern AI algorithms, and I do not know how far we can push the performance of generalized LSRDRs while retaining their inherent interpretability and mathematical behavior. 

If mathematical generalized LSRDRs can compete in performance with neural networks, then we are closer to solving the problem of general AI interpretability. But if not, then mathematical generalized LSRDRs could possibly be used as a narrow AI tool or even as a component of general AI; in this case, generalized LSRDRs could still improve general AI interpretability a little bit. An increased usage of narrow AI will (slightly) decrease the need for general AI.

Added 5/28/2025

When $d\ll n$, the iterative algorithm for computing LSRDRs produces a low rank projection matrix, but a matrix $P$ of rank $d$ can be factored as $SR$ where $S\in M_{n,d}\left(K\right),R\in M_{d,n}\left(K\right)$. To save computational resources, we may just work with $S,R$ during the training.

Suppose that $P=SR$ with $S\in M_{n,d}\left(K\right),R\in M_{d,n}\left(K\right)$. Then there are several ways to easily factor $f\left(P\right)$ as the product of an $n\times d$-matrix with a $d\times n$-matrix including the following factorizations:

$f\left(P\right)=3(SR{)}^2-2(SR{)}^3=$$3SRSR-2SRSRSR=S\cdot (3RSR-2RSRSR)$

$=(3SRS-2SRSRS)\cdot R=$$SRS\cdot (3R-2RSR)=$$(3S-2SRS)\cdot RSR$.

Therefore, define 

$f_1(S,R)=(S,3RSR-2RSRSR),f_2(S,R)=(3SRS-2SRSRS,R)$

$f_3(S,R)=(SRS,3R-2RSR),f_4(S,R)=(3S-2SRS,RSR)$.

In particular, if $f_i(S,R)=(S_1,R_1)$, then $S_1{R}_1=f\left(SR\right)$.

Let $i_N\in \{1,2,3,4\}$ for all $N$, and define $F(S,R)={lim}_{N\to \infty }{f}_{i_N}\dots f_{i_1}(S,R)$.

In the following algorithm, we will need to sometimes replace a pair $S,R$ with a new pair $S_1,R_1$ such that $S_1{R}_1=SR$ but where $(S_1,R_1)$ has norm smaller than $(S,R)$ because if we do not do this, after much training, the norm of $(S,R)$ will grow very large while $SR$ is just a projection matrix. To do this, we decrease the norm of $(S,R)$ using gradient descent. In particular, we observe that $SR=S(I+X)(I+X{)}^{-1}R$. We are taking the gradient when $X=0$, so we may have $SR\approx S(I+X)(I-X)R$ when $X$ is small. We want to move in the direction that minimizes the sum $\parallel S(I+X){\parallel }_2^2+\parallel (I-X)R{\parallel }_2^2$, but 

${\nabla }_X(\parallel S(I+X){\parallel }_2^2+\parallel (I-X)R{\parallel }_2^2){|}_{X=0}=$$2(S^{\ast }S-R{R}^{\ast }).$ 

Iterative algorithm for computing LSRDRs with low rank factorization:

Let $t>0$ but $t$ needs to be sufficiently small. Let $r>0$. Set $S_0\in M_{n,d}\left(K\right),R_0\in M_{d,n}\left(K\right)$. Let $G_0,H_0$ be random rank $d$ matrices. Define $G_N,H_N,G_N^{♯},H_N^{♯},S_N,R_N,S_N^{♯},R_N^{♯},T_N$ recursively for all $N\ge 0$ by setting

1. $G_N^{♯}=\Gamma (A_1,\dots ,A_r;S_N{R}_N{A}_1{S}_N{R}_N,\dots ,S_N{R}_N{A}_r{S}_N{R}_N{)}^{\ast }(G_N)$,
2. $H_N^{♯}=\Gamma (A_1,\dots ,A_r;S_N{R}_N{A}_1{S}_N{R}_N,\dots ,S_N{R}_N{A}_r{S}_N{R}_N)\left(H_N\right)$,
3. $G_{N+1}=G_N^{♯}/\text{Tr}\left(G_N^{♯}\right)$,
4. $H_{N+1}=H_N^{♯}/\text{Tr}\left(H_N^{♯}\right)$,
5. $(S_N^{♯},R_N^{♯})=F(S_N+t\cdot H_N\cdot S_N,R_N+t\cdot R_N\cdot G_N^{\ast })$,
6. $T_N=(S_N^{♯}{)}^{\ast }{S}_N^{♯}-R_N^{♯}\cdot (R_N^{♯}{)}^{\ast }$,
7. $R_{N+1}=(I+r\cdot t\cdot T_N)\cdot R_N^{♯}$, and
8. $S_{N+1}=S_N^{♯}\cdot (I+r\cdot t\cdot T_N{)}^{-1}$.

Then if everything goes right, $S_N{R}_N$ would converge to ${\text{Pro}}_{\infty }(A_1,\dots ,A_r)$ at an exponential rate.

Added 5/31/2025

The iterative algorithm for computing LSRDRs can be written in terms of completely positive operators. We define a completely positive superoperator as an operator of the form $\Phi (A_1,\dots ,A_r)$ where $A_1,\dots ,A_r$ are square matrices over $K$.  Completely positive operators are typically defined when $K=C$ for quantum information theory, but we are using a more general context here. If $E=\Phi (A_1,\dots ,A_r)$, then 

$\Gamma (A_1,\dots ,A_r;P_N{A}_1{P}_N,\dots ,P_N{A}_r{P}_N{)}^{\ast }(G_N)=E^{\ast }(G_N\cdot P)\cdot P$ and

$\Gamma (A_1,\dots ,A_r;P_N{A}_1{P}_N,\dots ,P_N{A}_r{P}_N)\left(H_N\right)=E(H_N\cdot P^{\ast })\cdot P^{\ast }$.

Added 6/7/2025

Iterative algorithm for computing LSRDRs apply not just to completely positive superoperators, but it also applies to positive operators that are not completely positive as long as the superoperators do not stray too far away from complete positivity. A linear superoperator $E:M_m\left(C\right)\to M_n\left(C\right)$ is said to be positive if $E\left(P\right)$ is positive semidefinite whenever $P$ is positive semidefinite. Clearly, every completely positive superoperator is positive, but not every positive operator is completely positive. For example, the transpose map $T$ defined by $T\left(P\right)=P^{\top }$ is always positive but not completely positive whenever $m=n>1$. The difference of two completely positive superoperators from $M_m\left(C\right)$ to $M_n\left(C\right)$ is known as a Hermitian preserving map. Every positive map is Hermitian preserving. If $E:M_m\left(C\right)\to M_n\left(C\right)$ is linear but not too far from positivity, then use the iterative algorithm for computing LSRDRs, we typically get the same results every time we train, and if $E$ is Hermitian preserving, then the operators $G_N,H_N$ will be positive semidefinite after training.