Math::PlanePath::SierpinskiArrowheadCentres -- self-similar triangular path traversal
use Math::PlanePath::SierpinskiArrowheadCentres; my $path = Math::PlanePath::SierpinskiArrowheadCentres->new; my ($x, $y) = $path->n_to_xy (123);
This is a version of the Sierpinski arrowhead path taking the centres of each triangle represented by the arrowhead segments. The effect is to traverse the Sierpinski triangle.
... ... / / . 30 . . . . . 65 . ... / \ / 28----29 . . . . . . 66 68 9 \ \ / 27 . . . . . . . 67 8 \ 26----25 19----18----17 15----14----13 7 / \ \ / / 24 . 20 . 16 . 12 6 \ / / 23 21 . . 10----11 5 \ / \ 22 . . . 9 4 / 4---- 5---- 6 8 3 \ \ / 3 . 7 2 \ 2---- 1 1 / 0 <- Y=0 -9 -8 -7 -6 -5 -4 -3 -2 -1 X=0 1 2 3 4 5 6 7
The base figure is the N=0 to N=2 shape. It's repeated up in mirror image as N=3 to N=6 then across rotated as N=6 to N=9. At the next level the same is done with the N=0 to N=8 shape, up mirrored as N=9 to N=17 and across rotated as N=18 to N=26, etc.
The X,Y coordinates are on a triangular lattice using every second integer X, per "Triangular Lattice" in Math::PlanePath.
The base pattern is a triangle like
.-------.-------. \ / \ / \ 2 / m \ 1 / \ / \ / .- - - -. \ / \ 0 / \ / .
Higher levels replicate this within the triangles 0,1,2 but the middle "m" is not traversed. The result is the familiar Sierpinski triangle by connected steps 2 across or 1 diagonal.
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
See the SierpinskiTriangle path to traverse by rows instead.
Counting the N=0,1,2 part as level 1, each replication level goes from
Nstart = 0 Nlevel = 3^level - 1 inclusive
For example level 2 from N=0 to N=3^2-1=9. Each level doubles in size,
0 <= Y <= 2^level - 1 - (2^level - 1) <= X <= 2^level - 1
The Nlevel position is alternately on the right or left,
Xlevel = / 2^level - 1 if level even \ - 2^level + 1 if level odd
The Y axis ie. X=0, is crossed just after N=1,5,17,etc which is is 2/3 through the level, which is after two replications of the previous level,
Ncross = 2/3 * 3^level - 1 = 2 * 3^(level-1) - 1
See "FUNCTIONS" in Math::PlanePath for the behaviour common to all path classes.
$path = Math::PlanePath::SierpinskiArrowheadCentres->new ()
Create and return a new path object.
($x,$y) = $path->n_to_xy ($n)
Return the X,Y coordinates of point number $n on the path. Points begin at 0 and if $n < 0 then the return is an empty list.
$n
$n < 0
If $n is not an integer then the return is on a straight line between the integer points.
An easy over-estimate of the range can be had from inverting the Nlevel formulas in "Level Ranges" above.
level = floor(log2(Ymax)) + 1 Nmax = 3^level - 1
For example Y=5, level=floor(log2(11))+1=3, so Nmax=3^3-1=26, which is the left end of the Y=7 row, ie. rounded up to the end of the Y=4 to Y=7 replication.
Math::PlanePath, Math::PlanePath::SierpinskiArrowhead, Math::PlanePath::SierpinskiTriangle
http://user42.tuxfamily.org/math-planepath/index.html
Copyright 2011 Kevin Ryde
Math-PlanePath is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3, or (at your option) any later version.
Math-PlanePath is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with Math-PlanePath. If not, see <http://www.gnu.org/licenses/>.
To install Math::PlanePath, copy and paste the appropriate command in to your terminal.
cpanm
cpanm Math::PlanePath
CPAN shell
perl -MCPAN -e shell install Math::PlanePath
For more information on module installation, please visit the detailed CPAN module installation guide.