The computation of the spiral center of triangles and quadrilaterals
Tis Veugen
tis.veugen@gmail.com
27 March 2021
Introduction
Spiral tilings of triangles have been studied by Waldman [1] and Fathauer [2]. This short memo deals with the computation of the spiral center given a triangle suitable for tessellation. The benefit of knowing the spiral center in advance is simple computation of gnomon corners based on the spiral’s scale factor and rotation angle. Also, the drawing can then be centered at any wanted position.
Furthermore, it turns out that the spiral center of quadrilaterals as designed by Fathauer [3] can be computed in the same way as for triangles.
Triangles
The seed triangle
has corners A, B, C; the angles at these corners are :
,
,
.
The opposite sides have length a, b, and 1.
Without loss of generality the side from A to C is horizontal. This eases to express the corners with complex variables. From the context should be clear if a capital is a corner or its corresponding complex number. Now, B and C can be expressed relative to A.
(1)
(2)
The spiral is
constructed by adding a smaller triangle A’B’C’,
being the first gnomon, next to ABC. The scale factor is called s;
see [2] for its computation. The rotation angle
is
in clockwise direction. The number of gnomons until alignment with
the seed triangle equals n = 5 for the above figure. This
number is not needed for the way of computation of the spiral center.
Side A’C’ is part of AB, such that corner C’ is
equal to B. The corners of A’B’C’ can be expressed
in terms of ABC.
(3)
(4)
(5)
The expression for
B’ characterizes mostly the spiral construction, since it
includes the scale factor s and the (negative) rotation angle
.
Therefore, the characteristic spiral transformation factor is
introduced:
(6)
The 3 corners are now combined in a tuple T of 3 elements
(7)
The transformed tuple T’ can then be expressed in matrix notation:
(8)
with
(9)
With this transformation matrix M any subsequent iteration of the seed tuple T can be easily computed, yielding the corners of the successive gnomons. Moreover, the center of the spiral can be derived by multiplying M infinite times. To perform this computation the eigendecomposition of M is calculated. The diagonal matrix L of eigenvalues for M is:
(10)
The corresponding matrix V of eigenvectors for M is:
(11)
So, M can now be written as:
(12)
where
(13)
The matrix
after
infinite multiplications becomes:
(14)
or
(15)
The 3 rows of
are
the same because the 3 corners converge to the same spiral center.
For the computation of the center the tuple T is expressed relative
to A:
(16)
Now the spiral
center S can be computed using the first row
of
:
(17)
So, spiral center
relative to corner A equals:
(18)
If the spiral center
should lie at the origin, then take A such that
.
Changing A is in fact a translation of the whole spiral figure. The
advantage of having the spiral center at the origin is that the next
iterated corner can be easily computed by multiplication with the
characteristic spiral transformation factor
.
This can also be verified by comparing the elements of
defined
as
(19)
with the
corresponding 3 elements of
;
or in other words:
(20)
Hence, with induction holds for iteration k of the 3 corners:
(21)
Remark 1: Iterating
in the “other” direction with growing triangles can be
realized by using the inverse of M, or the reciprocal of
.
Remark 2: A dual
tessellation can be formed based on the centers of the triangles.
Also, any convex combination f , being a row vector of size 3,
of the 3 corners can be used to generate a tessellation. Formula (21)
can be used to show that the k-iterated convex point equals the
initial convex point that has been transformed k times with
.
Since:
(22)
Remark 3: Formula (18) can be written in an alternative way to focus on the geometry for computing the spiral center.
(23)
The numerator in
(23) equals piece A’A, or in fact complex variable (A’ –
A). The denominator is part of a geometric series with
as
ratio. So, the spiral center can be constructed by repeatedly scaling
and rotating A’A and adding it to the partial sum.
As Fathauer [2] mentions, Figure 1 can be interpreted as n (=5) counter-clockwise spiral arms. One spiral arm is highlighted in red in Figure 2.
The scale factor for
the successive triangles is now
and
the rotation angle is
in
counter-clockwise direction. Focussing on side AC, this piece is
repeatedly scaled and rotated to construct a spiral. So relatively to
C the spiral parts are:
(24)
where:
(25)
The angles
and
are
related as
,
see also [2]. So
(26)
Relatively to A the spiral center for the arm is:
(27)
From the figure is
clear that
.
We can also show this in an algebraic way using the properties of the
seed triangle (e.g. considering B as origin and rotating BC
horizontally):
(28)
.
Now:
(29)
So
(30)
Quadrilaterals
The nomenclature for
the seed quadrilateral is analoguous to the seed triangle in the
above Triangles section. For constructing the spiral the triangle ABC
within the quadrilateral is iterated in exactly the same way as
triangle ABC in the Triangles section, using the same variable names
b, s,
,
,
etc. For this reason the same formulas in the Triangles section can
be applied for computing the spiral center of the iterated
quadrilateral.
References
[1] Cye H. Waldman, “Gnomon is an Island Whorled Polygons”, http://old.nationalcurvebank.org/gnomon/Gnomon%20is%20an%20Island.pdf, 3 July 2016
[2] Robert W. Fathauer, “Logarithmic Spiral Tilings of Triangles”, https://www.mathartfun.com/FathauerBridges2021v1.pdf, Submitted to the 2021 Bridges Conference
[3] Robert W. Fathauer, “Spiral tilings of quadrilaterals”, https://twitter.com/RobFathauerArt/status/1360248744021618692, 12 Februari 2021