Example #1
0
/* Calculates the appropriate power of the next-even transformation matrix
 * and returns the coefficients in the linear equation it represents.
 */
func calcTMat(n int) (aa, ab, ba, bb int) {
	t1Mat := mat64.NewDense(2, 2, []float64{3, 2, 2, 1})
	tMat := new(mat64.Dense)
	tMat.Pow(t1Mat, n)

	return round(tMat.At(0, 0)),
		round(tMat.At(0, 1)),
		round(tMat.At(1, 0)),
		round(tMat.At(1, 1))
}