md`# Rotor Estimation and AD in GA. In this notebook we will use [ganja.js](https://github.com/enkimute/ganja.js) to estimate the transformation between two point clouds. To perform this task we will use a Clifford Algebra over the Dual Vectors. Ganja.js normally only supports Clifford Algebras over the reals, but with a little trickery, we can easily create this exotic algebra. To do automatic differentiation, specifically to calculate partial derivatives with just one function evaluation, one needs dual vectors, consisting out of a scalar and ${tex`n`} dual coefficients. For ${tex`n=3`} the Cayley table of such an algebra is :

1	e₀₁	e₀₂	e₀₃
e₀₁	0	0	0
e₀₂	0	0	0
e₀₃	0	0	0

Using the Custom Cayley table feature of [ganja.js](https://github.com/enkimute/ganja.js), we can easily create such an Algebra. Additionally we also create the Clifford Algebra ${tex`\mathbb R_{3,0}`}. We then use the ganja.js translator of the Dual Vectors class to translate the operators of ${tex`\mathbb R_{3,0}`}. While the original functions of ${tex`\mathbb R_{3,0}`} support only scalars, in these new versions each of the coefficients of a multivector can itself be a Dual Vector. For example, for ${tex`\mathbb R_{2,0}`} the body of the function that implements the geometric product normally looks like this : \`\`\`javascript res[0]=b[0]*a[0]+b[1]*a[1]+b[2]*a[2]-b[3]*a[3]; res[1]=b[1]*a[0]+b[0]*a[1]-b[3]*a[2]+b[2]*a[3]; res[2]=b[2]*a[0]+b[3]*a[1]+b[0]*a[2]-b[1]*a[3]; res[3]=b[3]*a[0]+b[2]*a[1]-b[1]*a[2]+b[0]*a[3]; \`\`\` Where **res**, **a** and **b** are javascript Arrays. After translating it using the **inline** function of the Dual Vector algebra, it will look like this: \`\`\`javascript res[0]=Dual.Sub(Dual.Add(Dual.Add(Dual.Mul(b[0],a[0]),Dual.Mul(b[1],a[1])),Dual.Mul(b[2],a[2])),Dual.Mul(b[3],a[3])); res[1]=Dual.Add(Dual.Sub(Dual.Add(Dual.Mul(b[1],a[0]),Dual.Mul(b[0],a[1])),Dual.Mul(b[3],a[2])),Dual.Mul(b[2],a[3])); res[2]=Dual.Sub(Dual.Add(Dual.Add(Dual.Mul(b[2],a[0]),Dual.Mul(b[3],a[1])),Dual.Mul(b[0],a[2])),Dual.Mul(b[1],a[3])); res[3]=Dual.Add(Dual.Sub(Dual.Add(Dual.Mul(b[3],a[0]),Dual.Mul(b[2],a[1])),Dual.Mul(b[1],a[2])),Dual.Mul(b[0],a[3])); \`\`\` It still operates on an array, but all coefficients can now be Dual Vectors. These functions can then be used to implement the automatic differentiation for rotation and cost functions needed to do the estimation. We use a simple Gradient Descent - much more optimal algorithms exist (and can be implemented without Linear Algebra). In the example below, we parametrize the rotor manifold using a pure bivector ${tex`R`}, and minimize the simple cost function ${tex`\frac 1 n \sum\limits_{i=1}\limits^n | e^{- \frac R 2}\mathbf xe^{\frac R 2} - y |^2`}. We will evaluate this function using our Algebra over Dual Vectors which will result in automatically producing the needed derivatives. `